V=∫lrπf2(x)dxV=\int\limits_{l}^{r} {\pi f^2(x)dx}
在f(x)是多项式情况下
∫lrf(x)dx=∑i=0naii+1(ri+1−li+1)\int\limits_{l}^{r}{f(x)dx}=\sum\limits^{n}_{i=0}\frac{a_i}{i+1}(r^{i+1}-l^{i+1})
推出以上然后就可以做了,第二部分可以二分做
(代码很丑,是打广工训练的时候“抢”时间写的……)

#include<stdio.h>
#define eps 0.000001double a[23],b[23],xlow,xhigh,inc,l,r,mid;
int T,n,cs;double pow(double a,int b)
{double rt=1;while (b--) rt*=a;return rt;
}inline double V(double x)
{double rt=0;for (int i=0;i<=2*n;i++) rt+=3.1415926535*b[i]/(i+1)*(pow(x,i+1)-pow(xlow,i+1));return rt;
}inline void prepare()
{for (int i=0;i<=20;i++) b[i]=0;
}int main()
{while(~scanf("%d",&n)){cs++;prepare();for (int i=0;i<=n;i++) scanf("%lf",a+i);for (int i=0;i<=n;i++) for (int j=0;j<=n;j++) b[i+j]+=a[i]*a[j];scanf("%lf%lf%lf",&xlow,&xhigh,&inc);printf("Case %d: ",cs);printf("%.2f\n",V(xhigh));double fxxk=xlow;for (int i=1;i<=8;i++){l=xlow;r=xhigh;for (mid=(l+r)/2;l+eps<r;mid=(l+r)/2) if (V(mid)<inc) l=mid;else r=mid;if (l+eps<xhigh || inc<eps+V(l)) printf("%.2f ",(xlow=l)-fxxk);else{if (i==1) printf("insufficient volume");break;}}puts("");}
}

【CodeForces】gym-101205B Curvy Little Bottles (2012 ACM-ICPC World Finals B)相关推荐

  1. 【Codeforces】Gym 101173B Bipartite Blanket 霍尔定理+状压DP

    题意 给一张$n\times m$二分图,带点权,问有多少完美匹配子集满足权值和大于等于$t$ 这里有一个结论:对于二分图$\mathbb{A}$和$\mathbb{B}$集合,如果子集$A \in ...

  2. 【Codeforces】Round #488 (Div. 2) 总结

    [Codeforces]Round #488 (Div. 2) 总结 比较僵硬的一场,还是手速不够,但是作为正式成为竞赛生的第一场比赛还是比较圆满的,起码没有FST,A掉ABCD,总排82,怒涨rat ...

  3. 【Codeforces】925A Stairs and Elevators【贪心】

    [Codeforces]925A Stairs and Elevators [题目大意] 在一个n*m的矩阵里,有clcl个楼梯和cece个电梯,电梯和楼梯可以到任意一层,给出clcl个楼梯的位置和c ...

  4. 【二】gym初次入门一学就会---代码详细解析简明教程----平衡杆案例

    相关文章: [一]gym环境安装以及安装遇到的错误解决 [二]gym初次入门一学就会-简明教程 [三]gym简单画图 [四]gym搭建自己的环境,全网最详细版本,3分钟你就学会了! [五]gym搭建自 ...

  5. 【CodeForces】[546A]Soldier and Bananas

    直接运算 并没有什么弯 注意int就足以储存数据 不需要借钱时(n>res) 输出 0 #include<stdio.h> int main() {int k,n,w;while(s ...

  6. 【四】gym搭建自己的环境,全网最详细版本,3分钟你就学会了!

    相关文章: [一]gym环境安装以及安装遇到的错误解决 [二]gym初次入门一学就会-简明教程 [三]gym简单画图 [四]gym搭建自己的环境,全网最详细版本,3分钟你就学会了! [五]gym搭建自 ...

  7. 【CodeForces】445B A Lot of Games 字典树博弈

    传送门:[CodeForces]445B  A Lot of Games 题目大意:两人一起构造一个串,每人每次向串的末尾放一个字母,必须保证放了这个字母后能够成所给的N个串的前缀,如果某个人不能放时 ...

  8. 【Codeforces】【161Div2】

    [题目来源]http://www.codeforces.com/contest/263 [A. Beautiful Matrix] [解析]模拟即可.按照题目的意思,找到1所在的位置(x, y),然后 ...

  9. 【CodeForces】961 F. k-substrings 字符串哈希+二分

    [题目]F. k-substrings [题意]给定长度为n的串S,对于S的每个k-子串$s_ks_{k+1}...s_{n-k+1},k\in[1,\left \lceil \frac{n}{2} ...

  10. 【codeforces】【比赛题解】#940 CF Round #466 (Div. 2)

    人生的大起大落莫过如此,下一场我一定要回紫. [A]Points on the line 题意: 一个直线上有\(n\)个点,要求去掉最少的点,使得最远两点距离不超过\(d\). 题解: 暴力两重fo ...

最新文章

  1. ML之xgboost:利用xgboost算法(结合sklearn)训练mushroom蘑菇数据集(22+1,6513+1611)来预测蘑菇是否毒性(二分类预测)
  2. ESP32 OTA 接口简略说明
  3. 上传图片时,图片大小不一,设置宽高,图片拉伸,如何处理呢,我来告诉你...
  4. MySQL设置从库只读模式
  5. C# 参数签名字符串按 ASCII码排序,注意其中的坑
  6. mysql shell无法启动服务_[shell脚本] mysql服务启动脚本
  7. eclipse+mysql+tomcat配置JNDI
  8. AI驱动智能化日志分析 : 通过决策树给日志做聚类分析
  9. Java解码aac,Android MediaCodec解码AAC,AudioTrack播放PCM音频
  10. MTK modemlog分析
  11. CMMI 2.0 和 1.3
  12. 分享Echarts饼状图显示信息,内容,值,百分比都显示的代码
  13. mybatis处理xml大于小于号报异常
  14. 阿里怎么用DDD来拆分微服务?
  15. emmc/ufs, sd
  16. 11g中hanganalyze的格式
  17. push_back、emplace_back、std::move
  18. 基于TextField实现通用搜索组件
  19. BottomBarLayout——方便快捷地实现底部导航栏
  20. 如何将ppt转换为jpg图片

热门文章

  1. 软件优先级优先于硬件优先级_一切都是必须有优先级的
  2. 区块链毕业设计中期检查表模板
  3. 第5章第25节:如何寻找和下载视频、音频素材 [PowerPoint精美幻灯片实战教程]
  4. FileZilla Server目录乱码问题的解决
  5. RapidMiner遇到的问题
  6. mysql 纯真ip 导出dat_纯真IP数据库(qqwry.dat)转换成最新的IP数据库格式(ipwry.dat)
  7. vue仿微博评论回复_js模拟回帖/微博评论功能案例
  8. html5在线拍照,js html5 调用摄像头拍照
  9. 使用ExtendSim进行医疗仿真
  10. Ubuntu18.04 下载与安装(阿里云官方镜像站)