2016寒假训练——二分（2）

最新推荐文章于 2022-08-24 17:19:06 发布

原创最新推荐文章于 2022-08-24 17:19:06 发布 · 338 阅读

0 ·

CC 4.0 BY-SA版权

ACM算法专栏收录该内容

179 篇文章

订阅专栏

本文介绍了解决POJ3122分饼问题的方法，通过估算饼的面积并模拟分配过程，确保每个朋友都能从同一饼中分得一块。文章提供了AC代码示例，并分享了一些关键技巧，例如使用 acos(-1.0) 计算 pi 值以提高精度。

来源：http://poj.org/problem?id=3122

POJ3122

就是说现在有N个饼，F个朋友来分，要求（F+1）个能分的一块饼上的一块，不能来自不同的饼，这个题目就明显多了，就是估计大概的饼的面积，然后分了试试看，模拟进行和人数进行比较。

技巧：

1.这里POJ的DISCUSS中提到，这里pi取3.1415926不够，所以这里可以用pi=acos(-1.0)这样，可以规避这个错误。

2.double类型不要用l<=r，而是l-r>1e-6来判断是否应该结束。当然，这里的精度要求没有这么高。

3.最后要用.f提交，这是题目要求。。。可怕的要求，不过好多人WA在这了

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
using namespace std;
const int Maxn=10010;
const double pi=acos(-1.0);
double a[Maxn];
int judge(int n,int f,double s){
    int cnt=0,i;
    for(i=0;i<n;i++){
        cnt+=(int)floor(a[i]/s);
        if(cnt>=f) return 1;
    }
    return 0;
}
bool cmp(double a,double b){
    return a>b;
}
int main(){
    int c;
    int n,f;
    double s,mp;
    double ans;
    scanf("%d",&c);
    while(c--){
        ans=s=0;
        memset(a,0,sizeof(a));
        scanf("%d%d",&n,&f);
        f++;
        for(int i=0;i<n;i++){
            scanf("%lf",&a[i]);
            a[i]=a[i]*a[i]*pi;
            s+=a[i];
        }
        double r,l;
        sort(a,a+n,cmp);
        //cout<<a[0]<<" "<<a[n-1]<<endl;
        r=s/(double)f;
        l=a[0]/(double)f;
        while(r>1e-5+l){
            double mid=(r-l)/2+l;
            if(judge(n,f,mid)){
                l=mid;
            }
            else{
                r=mid;
            }
        }
        printf("%.4f\n",l);
    }
    return 0;
}

AC代码（VJ提交，非POJ）：