HDU 1348 和 POJ 1113 wall

首先，这次换了个算法来写，我用了Melkman

问题就来了，不知道怎么的，POJ1113的话，代码上去倒是可以AC。。。

但是在HDU的时候自己就o(╯□╰)o。。。

到底是怎么疑惑ing。。。。

View Code 
 1     #include <iostream>
 2     #include <algorithm>
 3     #include <stdio.h>
 4     #include <math.h>
 5     #include <string.h>
 6     using namespace std;
 7     const int maxn = 10005;
 8     const double pi = 3.1415926;
 9     int rear,pre,ans[2*maxn];
10 
11     struct Point
12     {
13         float x,y;
14     }p[maxn];
15 
16     bool cmp(Point a,Point b)
17     {
18         if(a.x!=b.x)
19             return a.x<b.x;
20         else
21             return a.y<b.y;
22     }
23 
24     double xmult(Point p1,Point p2,Point p0)
25     {
26         return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
27     }
28 
29     double dis(Point p,Point q)
30     {
31         return (p.x-q.x)*(p.x-q.x)+(p.y-q.y)*(p.y-q.y);
32     }
33 
34     void Melkman(int n)
35     {
36         int i,t;
37         rear=n;
38         pre=n;
39         ans[pre]=0;    ans[++rear]=1;
40         for(i=2;i<n;i++)
41         {
42             if(xmult(p[ans[rear-1]],p[ans[rear]],p[i])!=0)
43                 break;
44             ans[rear]=i;
45         }
46         ans[++rear]=i;    ans[--pre]=i;
47         if(xmult(p[ans[rear-2]],p[ans[rear-1]],p[ans[rear]])<0)
48         {
49             t=ans[rear-2];
50             ans[rear-2]=ans[rear-1];
51             ans[rear-1]=t;
52         }
53         for(i++;i<n;i++)
54         {
55             if(xmult(p[ans[rear-1]],p[ans[rear]],p[i])>0&&xmult(p[ans[pre]],p[ans[pre+1]],p[i])>0)
56                 continue;
57             while(xmult(p[ans[rear-1]],p[ans[rear]],p[i])<=0)
58                 --rear;
59             ans[++rear]=i;
60             while(xmult(p[ans[pre]],p[ans[pre+1]],p[i])<=0)
61                 ++pre;
62             ans[--pre]=i;
63         }
64     }
65 
66     int main()
67     {
68         long n,m,i,t,j,len;
69         cin>>t;
70         while(t--)
71         {
72             memset(ans,0,sizeof(ans));
73             cin>>n>>m;
74             for(i=0;i<n;i++)
75                 cin>>p[i].x>>p[i].y;
76             sort(p,p+n,cmp);
77             Melkman(n);
78             double sum=0;
79             for(i=pre;i<rear;i++)
80             {
81                  sum+=sqrt(dis(p[ans[i]],p[ans[i+1]]));
82                  //cout<<p[ans[i]].x<<" "<<p[ans[i]].y<<endl;
83             }
84             sum+=2*pi*m;
85             printf("%.0lf\n",sum);
86             if(t)
87                 cout<<endl;
88         }
89         return 0;
90     }

烦...

posted @ 2012-10-01 21:06  YORU  阅读(338)  评论(0)  编辑  收藏  举报