栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 3913 Bob wants to pour water

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 3913 Bob wants to pour water

#include<stdio.h>#include<string.h>#include<math.h>#define PI acos(-1.0)#define eps 0.0000000001struct node{  double z,a,b,c;}p[100005];struct node1{ double z,r;}q[100005];double w,len,v;int n,m;double cal(double x){    int i,j;    double ans=0;    for(i=1;i<=m;i++)    {        if(x>=p[i].z+p[i].c/2.0)ans=ans+p[i].a*p[i].b*p[i].c;        else if(x<=p[i].z-p[i].c/2.0)ans=ans+0;        else ans=ans+p[i].a*p[i].b*(x-(p[i].z-p[i].c/2.0));    }    for(i=1;i<=n;i++)    {        if(x>=q[i].z+q[i].r)ans=ans+4.0/3*PI*q[i].r*q[i].r*q[i].r;        else if(x<=q[i].z-q[i].r)ans=ans+0;        else if(x>q[i].z-q[i].r&&x<q[i].z)ans=ans+PI*(q[i].r+x-q[i].z)*(q[i].r+x-q[i].z)*(q[i].r-(q[i].r+x-q[i].z)/3.0);        else ans=ans+4.0/3*PI*q[i].r*q[i].r*q[i].r-PI*(q[i].r+q[i].z-x)*(q[i].r+q[i].z-x)*(q[i].r-(q[i].r+q[i].z-x)/3.0);    }    return ans;}double erfen(double l,double r){    double mid,mianji;    while(fabs(l-r)>eps)    {        mid=(l+r)/2.0;        double x=cal(mid);        mianji=w*len*mid-x;        if(fabs(mianji-v)<=eps)return mid;        else if(mianji<v)l=mid;        else if(mianji>v)r=mid;    }    return mid;}int main(){    int T,i,j,k;    scanf("%d",&T);    while(T--)    {        scanf("%lf%lf%lf%d%d",&w,&len,&v,&m,&n);        for(i=1;i<=m;i++) scanf("%lf%lf%lf%lf",&p[i].z,&p[i].a,&p[i].b,&p[i].c);        for(i=1;i<=n;i++) scanf("%lf%lf",&q[i].z,&q[i].r);        double ans=erfen(0,10000000000005);        printf("%.8fn",ans);    }    return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/373882.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号