HDU5432-Pyramid Split[通俗易懂]

HDU5432-Pyramid Split[通俗易懂]最近以各种借口偷懒了,上次BC的题一直也没补,今天才补完第一题,非常好过,可当时就是想不通~#include#include#include#defineeps1e-5doublea[10005],b[10005];intn;intbSearch(doublebegin,doubleend,doublee){doublemid,left=begi

大家好,又见面了,我是你们的朋友全栈君。

最近以各种借口偷懒了,上次BC的题一直也没补,今天才补完第一题,非常好过,可当时就是想不通~

#include<stdio.h>
#include<string.h>
#include<math.h>
#define eps 1e-5
double a[10005],b[10005];
int n;
int bSearch(double begin,double end,double e)
{
    double mid, left = begin, right = end;
    while(right-left>eps)
    {
        mid = (left + right)/2;
        double sum=0;
        for(int i=0;i<n;i++){
            if(a[i]<mid) continue;
            double r=(a[i]-mid)/a[i]*b[i];
            double h=a[i]-mid;
            sum+=1.0/3*r*r*h;
        }
        if(fabs(sum-e)<eps) return mid;
        if(sum-e>eps) left = mid;
        else right = mid ;
    }
    return left;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {

        scanf("%d",&n);
        for(int i=0;i<n;i++)
            scanf("%lf",&a[i]);
        for(int i=0;i<n;i++)
            scanf("%lf",&b[i]);
        double sum=0;
        for(int i=0;i<n;i++){
            sum+=1.0/3*b[i]*b[i]*a[i];
        }
        sum/=2;
        double k=bSearch(0,5000,sum);
        printf("%d\n",(int)k);
    }
    return 0;
}

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/151015.html原文链接:https://javaforall.net

(0)
上一篇 2022年6月19日 上午6:00
下一篇 2022年6月19日 上午6:16


相关推荐

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号