【二分】hdu 4033

#define LL long long
#define PI acos(-1.0)
#define N  105
#define MAX INT_MAX
#define MIN INT_MIN
#define eps 1e-8
double a[N];
int n;
double maxm;
double cal(double aa,double bb,double cc ){
    return (aa*aa+bb*bb-cc*cc)/(2*aa*bb);
}
int gao(double x){
    int i,j;
    double res=0.0;
    for(i=0;i<n;i++){
        if(a[i]+a[(i+1)%n]<=x)return 0;//big
        if(fabs( a[i]-a[(i+1)%n] )>=x)return -1;//small
        double tmp=cal(a[i],a[(i+1)%n],x);
        res+=acos(tmp);
    }
    if(fabs(res-2*PI)<eps)return 1;//equal
    else
    if(res<2*PI)return -1;//small
    else
    if(res>2*PI)
    return 0;//big
}
int main(){
    int t;
    scanf("%d",&t);
    int ca=1;
    while(t--){
        scanf("%d",&n);
        int i,j;
        maxm=0.0;
        for(i=0;i<n;i++){
            scanf("%lf",&a[i]);
            if(a[i]>maxm)maxm=a[i];
        }
        double l=0,r=2*maxm,mid;
        bool flag=0;
        while((r-l)>eps){
            mid=(l+r)/2.0;//cout<<mid<<endl;
            int tag=gao(mid);
            if(tag==1){
                flag=1;
                break;
            }
            if(tag==0){//big
                r=mid;
            }
            else{//small
                l=mid;
            }
        }
        printf("Case %d: ",ca++);
        if(!flag)printf("impossible\n");
        else{
            printf("%.3f\n",mid);
        }
    }
    return 0;
}




















  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值