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     }

烦...

转载于:https://www.cnblogs.com/yoru/archive/2012/10/01/2709919.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值