hdu 5513 三分+找规律

版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/yqdjl6/article/details/83757396

和精度斗智斗勇的一晚上，以为保留两位小数精度要求不会很大，结果卡了一晚上，这道题是说给你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;
}