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)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 在 docker 之间导出导入镜像「建议收藏」

    在 docker 之间导出导入镜像「建议收藏」很喜欢玩docker,但最新遇到一个问题,公司给的新机器的dns有问题,导致pull不下来镜像。没办法了,没有镜像什么神马都干不了,又不能花很多时间去搭建私有的镜像库,只有另寻办法了。废话少说,经过探究,发现docker提供把镜像导出export(保存save)为文件的机制,这样就可以把镜像copy到任意地方了。1.导…

    2025年9月19日
    6
  • PyCharm安装库numpy失败的解决方法

    PyCharm安装库numpy失败的解决方法事情是这样的,博主初学python和机器学习,在跑一个代码的时候被提示出现以下错误:(能被提示出现这个错误,可见确实是初学了!)图1:跑代码时候的报错注:图1是安装好了numpy后出现的第二个错误,错误本质是一样的,都是缺少某个库百度查资料后得知在PyCharm中有一个安装库的方法是:Settings>>PythonInterpreter>>点击图2中红色圈起来的加号,出现图3图2:安装库的一个方法图3:点击图2的加号后出现的界面顺利的话,只

    2022年8月25日
    16
  • 为什么要进行分销?

    为什么要进行分销?

    2021年6月17日
    95
  • python中关于if-else使用性能的一点感悟

    python中关于if-else使用性能的一点感悟

    2021年11月19日
    48
  • datatable删除行

    先列出正确的写法,如果你只想马上改错就先复制吧,如果你有时间想学习一下就继续看下面列出可能出错的可能性吧。1.如果只是想删除datatable中的一行,可以用DataRow的delete,但是必须

    2021年12月24日
    36
  • java.math.BigDecimal 比较大小

    java.math.BigDecimal 比较大小BigDecimala=newBigDecimal(101);BigDecimalb=newBigDecimal(111);//使用compareTo方法比较//注意:a、b均不能为null,否则会报空指针if(a.compareTo(b)==-1){System.out.println("a小于b");}if(a.compareTo(b)==…

    2022年7月14日
    31

发表回复

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

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