计算几何:极角排序(poj 2007 Scrambled Polygon)与简单凸包(poj 1113 Wall)

本文主要是介绍计算几何:极角排序(poj 2007 Scrambled Polygon)与简单凸包(poj 1113 Wall),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

ps:好久没来写博客了..准备重新开始了、两道简单题

poj 2007:http://poj.org/problem?id=2007 

按照(0,0)逆时针排序,由于在-180 ~ 180之内,直接叉积极角排序即可

/*将p[1]到p[m-1]的点根据p[0]按逆时针方向输出排序*/
#include <iostream>
#include <algorithm>
#include <cstdio>
#define MAXN 60
using namespace std;struct point{int x,y;
}p[MAXN];
int cross(int x1,int y1,int x2,int y2){return x2*y1 - x1*y2;
}
bool cmp(const point &a,const point &b){int x1 = a.x - p[0].x,y1 = a.y - p[0].y;int x2 = b.x - p[0].x,y2 = b.y - p[0].y;return cross(x1,y1,x2,y2) < 0;
}
int main()
{int m =0;while(~scanf("%d%d",&p[m].x,&p[m].y)) m++;sort(p+1,p+m,cmp);for(int i = 0;i<m;i++)printf("(%d,%d)\n",p[i].x,p[i].y);return 0;
}

接下来一道凸包:

poj  1113  http://poj.org/problem?id=1113

求一个多边形距离为L的周长,一开始以为是求凸包然后等比例放大,后来一想每个顶点用一点园周即可,整个图形合起来就是凸包的周长加一个半径为L的圆的周长


#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cmath>
#include <stack>
#define PI acos(-1)
#define MAXN 1000+10
using namespace std;int N,L;struct point{int x,y;
}p[MAXN];
int cross(int x1,int y1,int x2,int y2){return x2*y1 - x1*y2;
}
double dis(const point &a,const point &b){return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}
double dis(int x1,int y1,int x2,int y2){return sqrt((x1 - x2)*(x1 - x2) + (y1 - y2) * (y1 - y2));
}
bool cmp(const point &a,const point &b){int x1 = a.x - p[0].x,y1 = a.y - p[0].y;int x2 = b.x - p[0].x,y2 = b.y - p[0].y;if(cross(x1,y1,x2,y2) != 0)return cross(x1,y1,x2,y2) < 0;elsereturn dis(a,p[0]) < dis(b,p[0]);
}
bool ok(point p1,point p2,point p3){int x1 = p2.x - p1.x;int y1 = p2.y - p1.y;int x2 = p3.x - p2.x;int y2 = p3.y - p2.y;int c = cross(x1,y1,x2,y2);if(c <= 0) return true;else return false;
}
//凸包扫描算法
void GrahamScan(){stack<point> s;s.push(p[0]);s.push(p[1]);int i = 2;while(i < N ){point p2 = s.top();s.pop();point p1 = s.top();point p3 = p[i];if(ok(p1,p2,p3)){s.push(p2);s.push(p[i]);i++;}}s.pop();double C = 0.0;point now = p[N-1];while(!s.empty()){//printf("pre (s.top):x= %d,y = %d\n",s.top().x,s.top().y);C += dis(now,s.top());now = s.top();s.pop();}C += dis(p[0],p[N-1]);long long ans = (long long )((C + (2.0*L*PI))+0.5);printf("%I64d\n",ans);
}int main()
{while(~scanf("%d%d",&N,&L)){int m = 1;scanf("%d%d",&p[0].x,&p[0].y);for(int i=1;i<N;i++){point pt;scanf("%d%d",&pt.x,&pt.y);if(pt.y < p[0].y || (pt.y == p[0].y && pt.x < p[0].x)){p[m].x = p[0].x;p[m++].y = p[0].y;p[0].x = pt.x;p[0].y = pt.y;}else{p[m].x = pt.x;p[m++].y = pt.y;}}sort(p+1,p+N,cmp);GrahamScan();}return 0;
}

此题推荐一个不错的排序方法(水平序): http://www.tuicool.com/articles/iyqiY3Z


这篇关于计算几何:极角排序(poj 2007 Scrambled Polygon)与简单凸包(poj 1113 Wall)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/855658

相关文章

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu2289(简单二分)

虽说是简单二分,但是我还是wa死了  题意:已知圆台的体积,求高度 首先要知道圆台体积怎么求:设上下底的半径分别为r1,r2,高为h,V = PI*(r1*r1+r1*r2+r2*r2)*h/3 然后以h进行二分 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#includ

【数据结构】——原来排序算法搞懂这些就行,轻松拿捏

前言:快速排序的实现最重要的是找基准值,下面让我们来了解如何实现找基准值 基准值的注释:在快排的过程中,每一次我们要取一个元素作为枢纽值,以这个数字来将序列划分为两部分。 在此我们采用三数取中法,也就是取左端、中间、右端三个数,然后进行排序,将中间数作为枢纽值。 快速排序实现主框架: //快速排序 void QuickSort(int* arr, int left, int rig

usaco 1.3 Prime Cryptarithm(简单哈希表暴搜剪枝)

思路: 1. 用一个 hash[ ] 数组存放输入的数字,令 hash[ tmp ]=1 。 2. 一个自定义函数 check( ) ,检查各位是否为输入的数字。 3. 暴搜。第一行数从 100到999,第二行数从 10到99。 4. 剪枝。 代码: /*ID: who jayLANG: C++TASK: crypt1*/#include<stdio.h>bool h

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n