hdu 5531 Rebuild 三分

题目大意：

给你一个多边形，每个多边形的顶点都有一个圆，问是否存在一种情况，使得所有相邻的圆相切，切面积和最小。

比赛的时想到了三分，但是并没有发现其中的奥义，所以并不能AC！!

思路：

首先我们可以发现，因为相邻两个圆是紧贴着的，所以只要一个圆的大小改变了，其余所有圆都要改变。因为圆的面积是个二次函数，这样的话我们只需通过三分某一个边的长度来求解。

我们设dist[i]代表第i个点和第i+1个点的距离。P0所在圆的半径为r0，则P1为r1，以此类推。r0 + r1 = dist[0]，r1 + r2 = dist[1]。

当n为奇数的时候，可以根据这些式子直接算出r0。

当n为偶数的时候，可以三分r0的长度，进而求得总面积。至于左右区间，因为每个圆的半径都非负，即 r0 >= 0; r1 >= 0 -> dist[0] - r0 >= 0; 可以得到很多关于r0的不等式，这就是左右区间的范围。

最后有一个很难被注意到的细节，就是当n为偶数的时候，根据上面公式，可以发现奇数边长度和 == 偶数边长度和 如果忘记判断这一条 也会WA的。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cstdlib>
#include <cmath>
using namespace std;

const int maxn = 10000 + 10;
const double eps = 1e-8;
const double PI = acos(-1.0);

pair<int, int> point[maxn];
double dist[maxn];
double a[maxn], b[maxn];
int n;

int dcmp(double x){
    if(fabs(x) < eps) return 0;
    return x < 0 ? -1 : 1;
}

void solve1(){
    double ans = 0;
    for(int i = 0; i < n; i++){
        if(i&1) ans -= dist[i];
        else ans += dist[i];
    }
    ans /= 2;
    if(dcmp(ans) < 0){
        printf("IMPOSSIBLE\n");
        return;
    }
    int f = 0;
    double area = ans * ans, cur_r = ans;
    for(int i = 1; i < n; i++){
        area += (dist[i-1] - cur_r) * (dist[i-1] - cur_r);
        cur_r = dist[i-1] - cur_r;
        if(dcmp(cur_r) < 0){
            f = 1;
            break;
        }
    }
    if(f){
        printf("IMPOSSIBLE\n");
        return;
    }
    printf("%.2lf\n%.2lf\n", PI * area, ans);
    cur_r = ans;
    for(int i = 1; i < n; i++){
        printf("%.2lf\n", dist[i-1] - cur_r);
        cur_r = dist[i-1] - cur_r;
    }
}

double f(double r){
    double area = 0;
    for(int i = 0; i < n; i++) area += (a[i] * r + b[i]) * (a[i] * r + b[i]);
    return area;
}

void solve2(){
    a[0] = 1, b[0] = 0;
    double l = 0, r = 1e9;
    double re = 0;
    for(int i = 0; i < n; i++){
        if(i & 1) re += dist[i];
        else re -= dist[i];
    }
    if(dcmp(re) != 0){
        printf("IMPOSSIBLE\n");
        return;
    }
    for(int i = 1; i < n; i++){
        a[i] = -a[i-1];
        b[i] = dist[i-1] - b[i-1];
        if(a[i] > 0) l = max(l, -b[i]);
        else r = min(r, b[i]);
    }
    if(dcmp(l - r) > 0){
        printf("IMPOSSIBLE\n");
        return;
    }
    while(l + eps < r){
        double mid = (l + r) / 2;
        double midmid = (mid + r) / 2;
        double res1 = f(mid), res2 = f(midmid);
        if(res1 > res2) l = mid;
        else r = midmid;
    }
    printf("%.2lf\n", f(l) * PI);
    for(int i = 0; i < n; i++)
        printf("%.2lf\n", a[i] * l + b[i]);
}

void solve(){
    scanf("%d", &n);
    for(int i = 0; i < n; i++)
        scanf("%d%d", &point[i].first, &point[i].second);
    for(int i = 0; i < n; i++){
        int j = (i + 1) % n;
        dist[i] = sqrt((point[i].first - point[j].first) * (point[i].first - point[j].first) +
                       (point[i].second - point[j].second) * (point[i].second - point[j].second));
    }
    if(n & 1) solve1();
    else      solve2();
}

int main()
{
    int T;
    scanf("%d", &T);
    while(T--) solve();
    return 0;
}