HDU 5531 Rebuild ——— 2015ACM-ICPC亚洲区长春站

原题见HDU 5531

依次输入n个点的坐标，分别为圆心。保证相邻圆心的距离是个正整数。第n个圆和第1个圆相邻。要求相邻两个圆要相切，求全部圆面积和的最小值，以及此时半径的取法。
先根据坐标求出相邻圆心距离依次为$a_{0},a_{1},...,a_{n-1}$
设圆半径分别为$x_{0},x_{1},...,x_{n-1}$
化为方程，即：

$$\begin{cases}x_{0}+x_{1}=a_{0}\\ x_{1}+x_{2}=a_{1}\\ x_{2}+x_{3}=a_{2}\\ ...\\ x_{n-1}+x_{0}=a_{n-1} &\text{①}\end{cases}$$ $则x_{0}可以代替所有其他的x_{i}$,如下： $$\begin{cases}x_{1}=a_{0}-x_{0}\\ x_{2}=a_{1}-a_{0}+x_{0}\\ x_{3}=a_{2}-a_{1}+a_{0}-x_{0}\\ ...\\ x_{n-1}=a_{n-1}-a_{n-2}+...+(-1)^{(n-1)}x_{0} \end{cases}$$ $令x_{i}=b_{i}+c_{i}*x_{0}$,则 $$b_{0}=0,c_{0}=1$$ $$b_{i}=a_{i-1}-b_{i-1},~c_{i}=-c_{i-1}$$
将 $x_{n-1}=b_{n-1}+c_{n-1}*x_{0}$代入到①式，则 $$(c_{n-1}+1)x_{0}=a_{n-1}-b_{n-1}$$
$当c_{n-1}=1时，方程有唯一解x_{0}，只要在代入其他x_{i}过程中不产生负数，即是一组解。否则无解。$
$当c_{n-1}=-1时，若a_{n-1}-b_{n-1}不为0，则方程无解。否则有多解。$
考虑多解的情况。 $$\sum x_i^2=\sum (b_i^2+2b_ic_ix_0+c_i^2x_0^2)=nx_0^2+2(\sum b_ic_i)x_i+\sum b_i^2$$ 化为二次方程的情形，则只需考虑 $m=-{\sum b_ic_i\over n}可取的x_0$的区间范围。
而 $x_0$的范围由 $x_i$不小于0决定。

附code

/*--------------------------------------------
 * File Name: HDU 5531
 * Author: Danliwoo
 * Mail: Danliwoo@outlook.com
 * Created Time: 2015-11-01 22:38:00
--------------------------------------------*/

#include <bits/stdc++.h>
using namespace std;
#define LL long long
const double Pi = acos(-1);
int n;
struct point{
    LL x, y;
    void sc(){
        scanf("%I64d%I64d", &x, &y);
    }
    LL dis(point b){
        LL d = (x-b.x)*(x-b.x)+(y-b.y)*(y-b.y);
        LL c = sqrt(d);
        return c;
    }
}p[10010];
LL a[10010], b[10010], c[10010];
double x[10010];
void give(double y)
{
    double s = 0;
    if(y < 0){
        printf("IMPOSSIBLE\n");
        return;
    }
    x[0] = y;
    s += pow(x[0], 2);
    for(int i = 1;i < n;i++)
    {
        x[i] = b[i]+c[i]*x[0];
        if(x[i] < 0){
            printf("IMPOSSIBLE\n");
            return;
        }
        s += pow(x[i],2);
    }
    s *= Pi;
    printf("%.2f\n", s);
    for(int i = 0;i < n;i++)
        printf("%.2f\n", x[i]);
}
void solve()
{
    double L = 0, R = a[0];
    for(int i = 2;i < n;i++)
    {
        if(c[i] == 1)
            L = max(L, -1.0*b[i]);
        else if(c[i] == -1)
            R = min(R, 1.0*b[i]);
        if(L > R)
        {
            printf("IMPOSSIBLE\n");
            return;
        }
    }
    double B = 0;
    for(int i = 0;i < n;i++)
        B += 2*b[i]*c[i];
    double m = -B/n/2.0;
    if(R < m)
        give(R);
    else if(L > m)
        give(L);
    else
        give(m);
}
int main()
{
    int T;
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        for(int i = 0;i < n;i++)
            p[i].sc();
        for(int i = 0;i < n;i++)
            a[i] = p[i].dis(p[(i+1)%n]);
        b[0] = 0; c[0] = 1;
        for(int i = 1;i < n;i++)
        {
            c[i] = -c[i-1];
            b[i] = a[i-1]-b[i-1];
        }
        if(c[n-1] == 1)
        {
            give((a[n-1]-b[n-1])/2.0);
            continue;
        }
        if(b[n-1] != a[n-1])
        {
            printf("IMPOSSIBLE\n");
            continue;
        }
        solve();
    }
    return 0;
}