Online Judge | Problem Set | Authors | Online Contests | User | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Web Board Home Page F.A.Qs Statistical Charts | Current Contest Past Contests Scheduled Contests Award Contest |
跪拜求挑错,谢……#include <stdio.h> #include <math.h> #include <algorithm> using namespace std; struct point { int x,y; }; int dist(point a,point b) { return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y); } struct poly{ point p[70000]; int n; poly(){} void init(int nn){n=nn; for(int i=0;i<n;++i) scanf("%d %d",&p[i].x,&p[i].y); p[n]=p[0];return; }; }; int mmax(int a,int b) { if(a>b)return a; return b; } int cross(point a, point b, point c) { return (a.x - c.x)*(b.y - c.y)-(a.y - c.y)*(b.x - c.x); } int rotating_calipers(poly a) { int q=1,ans=0; for(int p=0;p<a.n;p++) { while(cross(a.p[p+1],a.p[q+1],a.p[p])>cross(a.p[p+1],a.p[q],a.p[p])) q=(q+1)%a.n; ans=max(ans,max(dist(a.p[p],a.p[q]),dist(a.p[p+1],a.p[q+1]))); } return ans; } bool operator<(const point &a,const point &b){ return a.y<b.y||a.y==b.y&&a.x<b.x; } void convex(poly a,poly &b){ b.n=0; sort(a.p,a.p+a.n); for(int i=0;i<a.n;++i){ while(b.n>1&&cross(a.p[i],b.p[b.n-2],b.p[b.n-1])<0)--b.n; b.p[b.n++]=a.p[i]; } int k=b.n; for(int i=a.n-2;i>=0;--i){ while(b.n>k&&cross(a.p[i],b.p[b.n-2],b.p[b.n-1])<0)--b.n; b.p[b.n++]=a.p[i]; } if(b.n>1)--b.n; b.p[b.n]=b.p[0]; } poly a; int main() { int n,i,j; while(scanf("%d",&n)!=EOF) { a.init(n); convex(a,a); int ans=rotating_calipers(a); printf("%d\n",ans); } return 0; } Followed by:
Post your reply here: |
All Rights Reserved 2003-2013 Ying Fuchen,Xu Pengcheng,Xie Di
Any problem, Please Contact Administrator