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

原题见HDU 5531

依次输入n个点的坐标,分别为圆心。保证相邻圆心的距离是个正整数。第n个圆和第1个圆相邻。要求相邻两个圆要相切,求全部圆面积和的最小值,以及此时半径的取法。
先根据坐标求出相邻圆心距离依次为 a0,a1,...,an1
设圆半径分别为 x0,x1,...,xn1
化为方程,即:

x0+x1=a0x1+x2=a1x2+x3=a2...xn1+x0=an1
x0xi ,如下:
x1=a0x0x2=a1a0+x0x3=a2a1+a0x0...xn1=an1an2+...+(1)(n1)x0
xi=bi+cix0 ,则
b0=0,c0=1
bi=ai1bi1, ci=ci1

xn1=bn1+cn1x0 代入到①式,则
(cn1+1)x0=an1bn1

cn1=1x0xi
cn1=1an1bn10
考虑多解的情况。
x2i=(b2i+2bicix0+c2ix20)=nx20+2(bici)xi+b2i
化为二次方程的情形,则只需考虑 m=bicinx0 的区间范围。
x0 的范围由 xi 不小于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;
}
  • 0
    点赞
  • 0
    收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:编程工作室 设计师:CSDN官方博客 返回首页
评论

打赏作者

Danliwoo

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值