hdu 5513 三分+找規律
阿新 • • 發佈:2018-11-14
和精度鬥智鬥勇的一晚上,以為保留兩位小數精度要求不會很大,結果卡了一晚上,這道題是說給你N個點構成的一個凸多邊形,對於每個點畫一個圓,要求相鄰兩個圓必須相切,求所有圓的面積最小值為多少。這道題剛開始卡在fixed上,因為資料量比較大,所以每次用fixed會超時,改成scanf以後就好了,另一個坑點是求總面積,用∑(riripi)會wa,必須用∑(ri*ri)*pi,這又改了好長時間才發現。這道題的規律就是奇數個點的第一個r是一個區間,對這個區間三分。偶數個點的情況r是一個定值。
#include<bits/stdc++.h> using namespace std; const int maxn=1e5+5; const double pi =acos(-1.0); const double eps=1e-10; struct Node { double x,y; }node[maxn]; double cnt[maxn],res[maxn]; int n; double len(Node a,Node b) { return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y)); } double check(double r) { double ans=r*r; res[1]=r; for(int i=2;i<=n;i++) { res[i]=cnt[i-1]-res[i-1]; ans+=res[i]*res[i]; } return ans*pi; } int main() { int T;scanf("%d",&T); while(T--) { int flag=0;scanf("%d",&n); double ans=0; for(int i=1;i<=n;i++) scanf("%lf%lf",&node[i].x,&node[i].y); for(int i=1;i<n;i++) cnt[i]=len(node[i],node[i+1]); cnt[n]=len(node[1],node[n]); for(int i=1;i<=n;i++) { if(i%2==0) ans-=cnt[i]; else ans+=cnt[i]; } if(n%2==1) { res[1]=ans/2.0;ans=0; for(int i=2;i<=n;i++) res[i]=cnt[i-1]-res[i-1]; for(int i=1;i<=n;i++) ans+=res[i]*res[i]*pi; if(fabs(res[1]+res[n]-cnt[n])>eps) flag=1; } else if(n%2==0) { if(fabs(ans)>eps) flag=1; else { double l=0,r=min(cnt[1],cnt[n]),sum=0; for(int i=1;i<=n;i++) { sum=cnt[i]-sum; if(i%2==1) r=min(r,sum); else l=max(l,-sum); } for(int i=0;i<100;i++) { double lmid=(l+r)/2.0; double rmid=(lmid+r)/2.0; if(check(lmid)<check(rmid)) r=rmid; else l=lmid; } ans=check((l+r)/2.0); } } for(int i=1;i<=n;i++) if(res[i]<0.00) flag=1; if(flag==1) printf("IMPOSSIBLE\n"); else { printf("%.2f\n",ans); for(int i=1;i<=n;i++) printf("%.2f\n",res[i]); } } return 0; }