凸包。。。我要刷光各大OJ的凸包。。。

從discuss學來的。。。剛開始沒看懂題意都。。想麻煩了。。

這個題是給你凸包上的點，給的無序，讓你排出來有序的。

這個方法遇到一種情況就壞了。。。

就是下圖：

但是，對於凸包來說，一定能找到一個點，過這個點平行於Y軸或者X軸的直線 能劃分成兩個區域。

不過具體怎麼求，還要想想。

這道題是直接告訴你了，按x最小的去求。

啊。。不對。。這題應該是資料弱了。。這麼做是不對滴。。。上面的圖就是反例啊。。。

我試試另一種做法~剛想起來~

下面這個程式碼是很弱的。。。看下個程式碼吧~^ ^

#include <queue> #include <stack> #include <math.h> #include <stdio.h> #include <stdlib.h> #include <iostream> #include <limits.h> #include <string.h> #include <string> #include <algorithm> using namespace std; const int MAX = 100010; struct point{ int x,y; }; point c[MAX]; int m; long long crossProduct(point a,point b,point c)//向量 ac 在 ab 的方向 { return (c.x - a.x*1ll)*(b.y - a.y) - (b.x - a.x*1ll)*(c.y - a.y); } long long disp2p(point a,point b) { return ( a.x - b.x*1ll ) * ( a.x - b.x ) + ( a.y - b.y*1ll ) * ( a.y - b.y ); } bool cmp1(point a,point b) // 排序 { long long len = crossProduct(c[0],a,b); if( len == 0 ) return disp2p(c[0],a) < disp2p(c[0],b); return len < 0; } bool cmp2(point a,point b) { long long len = crossProduct(c[0],a,b); if( len == 0 ) return disp2p(c[0],a) > disp2p(c[0],b); return len < 0; } void output(point *p,int n) { for(int i=0; i<n; i++) printf("%d %d/n",p[i].x,p[i].y); } void solve() { int tmp = 0; for(int i=1; i<m; i++) if( c[i].x < c[tmp].x || c[i].x == c[tmp].x && c[i].y < c[tmp].y ) tmp = i; swap(c[0],c[tmp]); point L[MAX],U[MAX]; int l = 0,u = 0; for(int i=1; i<m; i++) if( c[i].y <= c[0].y ) L[l++] = c[i]; else U[u++] = c[i]; sort(L,L+l,cmp1); sort(U,U+u,cmp2); printf("%d/n",m); printf("%d %d/n",c[0].x,c[0].y); output(L,l); output(U,u); } int main() { int ncases,n; int x,y; char ch; scanf("%d",&ncases); while( ncases-- ) { scanf("%d",&n); m = 0; for(int i=0; i<n; i++) { scanf("%d %d %c",&x,&y,&ch); if( ch == 'Y' ) { c[m].x = x; c[m++].y = y; } } solve(); } return 0; }

昨晚失眠了，3點半的時候床上躺了會，迷迷瞪瞪地，四點被姐的飛信吵醒了 = =。。遲鈍死了都。突然想起來怎麼劃分這題了。

直接找最遠點對，然後這個點對分成的凸包的兩個區域再按照上面的排序，就OK了。。。

然後下床實現~排除各種BUG。。A掉了，時間和上個差不多。都是N*LOGN的~~

#include <queue> #include <stack> #include <math.h> #include <stdio.h> #include <stdlib.h> #include <iostream> #include <limits.h> #include <string.h> #include <string> #include <algorithm> using namespace std; const int MAX = 100010; struct point{ int x,y; }; point c[MAX]; int stk[MAX],top; int m; long long crossProduct(point a,point b,point c)//向量 ac 在 ab 的方向 { return (c.x - a.x*1ll)*(b.y - a.y) - (b.x - a.x*1ll)*(c.y - a.y); } long long disp2p(point a,point b) { return ( a.x - b.x*1ll ) * ( a.x - b.x ) + ( a.y - b.y*1ll ) * ( a.y - b.y ); } bool cmp1(point a,point b) // 排序 { long long len = crossProduct(c[0],a,b); if( len == 0 ) return disp2p(c[0],a) < disp2p(c[0],b); return len < 0; } bool cmp2(point a,point b) { long long len = crossProduct(c[0],a,b); if( len == 0 ) return disp2p(c[0],a) > disp2p(c[0],b); return len < 0; } int Left,Right; void Rotating_calipers(int stk[],int n) { stk[n] = stk[0]; int q = 1; long long ans = 0ll; int a,b; for(int i=0; i<n; i++) { while( crossProduct(c[stk[i]],c[stk[i+1]],c[stk[q]]) < crossProduct(c[stk[i]],c[stk[i+1]],c[stk[q+1]]) ) q = (q+1)%n; if( disp2p(c[stk[i]],c[stk[q]]) > ans ) { ans = disp2p(c[stk[i]],c[stk[q]]); a = stk[i]; b = stk[q]; } } if( c[a].x < c[b].x || c[a].x == c[b].x && c[a].y < c[b].y ) { Left = a; Right = b; } else { Left = b; Right = a; } } void Graham(int n) { int tmp = 0; for(int i=1; i<n; i++) if( c[i].x < c[tmp].x || c[i].x == c[tmp].x && c[i].y < c[tmp].y ) tmp = i; swap(c[0],c[tmp]); sort(c+1,c+n,cmp1); stk[0] = 0; stk[1] = 1; top = 1; for(int i=2; i<n; i++) { while( crossProduct(c[stk[top]],c[stk[top-1]],c[i]) <= 0ll && top >= 1 ) top--; stk[++top] = i; } Rotating_calipers(stk,top+1); } vector<point> p; void solve() { Graham(m); point L[MAX],U[MAX]; int l = 0,u = 0; for(int i=0; i<m; i++) { if( i == Left || i == Right ) continue; if( crossProduct(c[Left],c[Right],c[i]) >= 0 ) L[l++] = c[i]; else U[u++] = c[i]; } sort(L,L+l,cmp1); sort(U,U+u,cmp2); printf("%d/n",m); p.clear(); p.push_back(c[Left]); for(int i=0; i<l; i++) p.push_back(L[i]); p.push_back(c[Right]); for(int i=0; i<u; i++) p.push_back(U[i]); int tmp = 0; for(int i=0; i<p.size(); i++) if( p[i].x < p[tmp].x || p[i].x == p[tmp].x && p[i].y < p[tmp].y ) tmp = i; int k = 0; while( k < m ) { printf("%d %d/n",p[tmp].x,p[tmp].y); tmp = (tmp+1)%m; k++; } } int main() { int ncases,n; int x,y; char ch; scanf("%d",&ncases); while( ncases-- ) { scanf("%d",&n); m = 0; for(int i=0; i<n; i++) { scanf("%d %d %c",&x,&y,&ch); if( ch == 'Y' ) { c[m].x = x; c[m++].y = y; } } solve(); } return 0; }

好啦，第三個程式碼~內點排序，用重心哈，凸包的重心一定是在內部的（也不可能在邊上）。

因為加了精度判斷，慢了不少。

寫比較函式很糾結啊。。。

#include <queue> #include <stack> #include <math.h> #include <stdio.h> #include <stdlib.h> #include <iostream> #include <limits.h> #include <string.h> #include <string> #include <algorithm> using namespace std; const int MAX = 100010; struct point{ double x,y; }; point c[MAX]; int m; const double eps = 1e-8; bool dy(double x,double y) { return x > y + eps;} // x > y bool xy(double x,double y) { return x < y - eps;} // x < y bool dyd(double x,double y) { return x > y - eps;} // x >= y bool xyd(double x,double y) { return x < y + eps;} // x <= y bool dd(double x,double y) { return fabs( x - y ) < eps;} // x == y // 叉積 (三點) double crossProduct(point a,point b,point c)//向量 ac 在 ab 的方向 { return (c.x - a.x)*(b.y - a.y) - (b.x - a.x)*(c.y - a.y); } point bary_center(point p[],int n) { point ans,t; double area = 0.0,t2; ans.x = 0.0; ans.y = 0.0; for(int i=1; i<n-1; i++) { t2 = crossProduct(p[i],p[0],p[i+1])/2.0; ans.x += (p[0].x + p[i].x + p[i+1].x)*t2; ans.y += (p[0].y + p[i].y + p[i+1].y)*t2; area += t2; } ans.x /= (3*area); ans.y /= (3*area); return ans; } point s; bool cmp(point a,point b) { double l1 = crossProduct(s,c[0],a); double l2 = crossProduct(s,c[0],b); if( dyd(l1,0.0) && dyd(l2,0.0) ) return xy(crossProduct(s,a,b),0.0); if( dyd(l1,0.0) && xyd(l2,0.0) ) return 0; if( xyd(l1,0.0) && xyd(l2,0.0) ) return xy(crossProduct(s,a,b),0.0); return 1; } void solve() { s = bary_center(c,m); int tmp = 0; for(int i=1; i<m; i++) if( xy(c[i].x,c[tmp].x) || dd(c[i].x,c[tmp].x) && xy(c[i].y,c[tmp].y) ) tmp = i; swap(c[tmp],c[0]); sort(c+1,c+m,cmp); printf("%d/n",m); for(int i=0; i<m; i++) printf("%.0lf %.0lf/n",c[i].x,c[i].y); } int main() { int ncases,n; double x,y; char ch; scanf("%d",&ncases); while( ncases-- ) { scanf("%d",&n); m = 0; for(int i=0; i<n; i++) { scanf("%lf %lf %c",&x,&y,&ch); if( ch == 'Y' ) { c[m].x = x; c[m++].y = y; } } solve(); } return 0; }