博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
hdu 1724 Ellipse simpson积分
阅读量:5810 次
发布时间:2019-06-18

本文共 862 字,大约阅读时间需要 2 分钟。

/*hdu 1724 Ellipse simpson积分求椭圆的部分面积simpson积分法 http://zh.wikipedia.org/zh-tw/%E8%BE%9B%E6%99%AE%E6%A3%AE%E7%A7%AF%E5%88%86%E6%B3%95*/#include
#include
const double eps=1e-8;inline double jue(double a){ return a>0?a:-a;}int dbcmp(double a){ if(a>eps)return 1; if(a<(-eps)) return -1; return 0;}template
inline T sqr(T d){return d*d;}double a,b;inline double func(double x){ return 2*b*sqrt(1-sqr(x)/sqr(a));}double simpson(double l,double r){ return (func(l)+func(r)+4*func((l+r)/2))/6*(r-l);}double work(double l,double r){ double mid=(l+r)/2; double ans=simpson(l,r); if(dbcmp(ans-simpson(l,mid)-simpson(mid,r))==0) { return ans; }else { return work(l,mid)+work(mid,r); }}int main(){ int n; double l,r; scanf("%d",&n); while(n--) { scanf("%lf%lf%lf%lf",&a,&b,&l,&r); double ret=work(l,r); printf("%.3f\n",ret); } return 0;}

转载地址:http://mjjbx.baihongyu.com/

你可能感兴趣的文章
Android打开/data/目录以及导出文件
查看>>
黑苹果安装
查看>>
多叉树实现类目体系
查看>>
线程池的简单使用
查看>>
Thinkphp5.0关于自动生成模块
查看>>
Centos命令
查看>>
[WorldWind学习]23.TerrainAccessor
查看>>
关于SpringBoot的单元测试
查看>>
洛谷 P2015 二叉苹果树(codevs5565) 树形dp入门
查看>>
详解c++多线程(一)
查看>>
Windows Gdi & CDC和HDC的区别与转换
查看>>
EPEL源
查看>>
POJ - 3662 Telephone Lines
查看>>
如何在云服务器创建maven私有仓库
查看>>
委托的来由
查看>>
数据库开启对sys用户的审计
查看>>
sublimetext 自定义build
查看>>
郝小亮-读王坚《在线》:互联网能做的和没做的事还有很多
查看>>
(转)hasLayout与BFC的触发条件
查看>>
Leetcode Longest Substring Without Repeating Characters
查看>>