POJ 1819 Disks

最新推荐文章于 2019-08-07 19:02:15 发布

taozifish

最新推荐文章于 2019-08-07 19:02:15 发布

阅读量699

点赞数

CC 4.0 BY-SA版权

分类专栏：计算几何文章标签： ini

本文链接：https://blog.youkuaiyun.com/taozifish/article/details/6598991

计算几何专栏收录该内容

28 篇文章

订阅专栏

本文详细解析了计算几何领域的模拟题，包括技巧和难点，特别关注了最后一个圆的问题。通过实例，逐步指导读者理解并解决这类问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

//计算几何的模拟题，有trick，和恶心，注意最后的圆
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const double EP=1e-10;
int main(){
    //freopen("1.txt","r",stdin);
    int n,i,j,touch[1010],sum,tmax;
    bool flag[1010];
    double r[1010], x[1010],lmax;
    while(scanf("%d", &n)!=EOF){
        lmax=0;sum=0;tmax=n-1;
    memset(touch,-1,sizeof(touch));
    memset(flag,0,sizeof(flag));
    for(i=0; i<n; i++){
        scanf("%lf", &r[i]);
        x[i]=r[i];
    }
    for(i=1; i<n; i++){
        for(j=0; j<i; j++){
            if(x[j]+2*sqrt(r[i]*r[j])-x[i]>EP){
            x[i]=x[j]+2*sqrt(r[i]*r[j]);
            touch[i]=j;
            }
        }
        if(x[i]+r[i]-lmax>EP){
                lmax=x[i]+r[i];
                tmax=i;
            }
        }
    for(i=0; i<n; i++)
        for(j=touch[i]+1; j<i; j++)
        if(!flag[j]){
        flag[j]=1;
        sum++;
        }
    for(i=tmax+1;i<n;i++)
        if(!flag[i]){
        flag[i]=1;
        sum++;
        }
    printf("%d\n", sum);
    for(i=0; i<n; i++)
        if(flag[i])
        printf("%d\n",i+1);
    }
    return 0;
}