ACM比赛用的简单模板.docVIP

邵伯仲的模板 凸包 以斜率为基准的gr法 #includeiostream #includealgorithm #includevector #includemath.h using namespace std; typedef struct { int x,y; double k; }Node; Node *a; bool cmp(const Node a,const Node b)// 根据斜率和x进行排序 { if (a.k==b.k) return a.xb.x; return a.kb.k; } bool mutiply(const Node one,const Node two,const Node three)//差集 { double p1x,p1y,p2x,p2y,result; p1x=two.x-one.x; p1y=two.y-one.y; p2x=three.x-one.x; p2y=three.y-one.y; result=p1x*p2y-p2x*p1y; if (result=0) return true; return false; } double Area(const Node one,const Node two)//计算面积 { double p1x,p1y,p2x,p2y,result; p1x=one.x-0; p1y=one.y-0; p2x=two.x-0; p2y=two.y-0; result=(p1x*p2y-p2x*p1y)/2; return result; } int main() { double sum; int last,count; vector Nodes; Node one,two,three; int n,i,min,t; cint; while(t--) { cinn; a=new Node[n]; sum=0; min=0; for(i=0;in;i++) { cina[i].xa[i].y; if (a[min].x=a[i].x) { if (a[min].xa[i].x||a[min].ya[i].y) min=i; } } if (n3) {cout0.0endl;continue;} for(i=0;in;i++)//计算斜率 { a[i].k=(double)(a[i].y-a[min].y)/(double)(a[i].x+1-a[min].x); } a[min].k=-2100000000; sort(a,a+n,cmp);//按斜率和x排序 s.push_back(a[0]);//开始形成凸包 s.push_back(a[1]); for(i=2;in;i++) { last=s.size(); two=s[last-1]; one=s[last-2]; while(mutiply(one,two,a[i])ins.size()=2) { s.pop_back(); if(s.size()==1) break; last=s.size(); two=s[last-1]; one=s[last-2]; } s.push_back(a[i]); }//形成了凸包 last=s.size();//计算面积 one=s[0];two=s[1]; for(i=2;ilast;i++) { sum+=Area(one,two); one=two;two=s[i]; }