1. 程式人生 > >hdu 5513 三分+找規律

hdu 5513 三分+找規律

和精度鬥智鬥勇的一晚上,以為保留兩位小數精度要求不會很大,結果卡了一晚上,這道題是說給你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;
}