hdu 1724 Ellipse——辛普森积分

题目：http://acm.hdu.edu.cn/showproblem.php?pid=1724

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define db double
using namespace std;
const db eps=1e-5;
int T;db a,b,l,r;
db f(db x){return sqrt((1-x*x/a)*b);}
db cal(db l,db r){return (f(l)+4*f((l+r)/2)+f(r))*(r-l)/6;}
db simp(db l,db r,db eps,db ret)
{db mid=(l+r)/2,vl=cal(l,mid),vr=cal(mid,r);if(fabs(vl+vr-ret)<=15*eps)return vl+vr+(vl+vr-ret)/15;//not ret+()/15return simp(l,mid,eps/2,vl)+simp(mid,r,eps/2,vr);
}
int main()
{scanf("%d",&T);while(T--){scanf("%lf%lf%lf%lf",&a,&b,&l,&r);a*=a;b*=b;printf("%.3lf\n",simp(l,r,eps,cal(l,r))*2);}return 0;
}

转载于:https://www.cnblogs.com/Narh/p/10141896.html

hdu 1724 Ellipse——辛普森积分相关推荐

HDU 1724 Ellipse ——Simpson积分
[题目分析] 一看题目,直接把椭圆积分起来就可以了嘛. 然后发现椭圆比较难积分,还是算了吧. 用Simpson积分硬上. 大概就是用二次函数去拟合面积. [代码] #include <cstdi ...
夜深人静写算法（四十八）- 辛普森积分
文章目录一.引例 1.简单函数积分 2.复杂函数积分二.梯形法则 1.梯形法则原理 2.梯形法则实现三.辛普森积分 1.线段拟合 2.抛物线拟合 3.抛物线的定积分 4.辛普森积分实现四.自适 ...
bzoj 2178 圆的面积并 —— 辛普森积分
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2178 先看到这篇博客:https://www.cnblogs.com/heisenberg- ...
[BZOJ1502]月下柠檬树(自适应辛普森积分)
1502: [NOI2005]月下柠檬树 Time Limit: 5 Sec Memory Limit: 64 MB Submit: 1387 Solved: 739 [Submit][Statu ...
数值分析——自适应辛普森积分
数值分析--自适应辛普森积分自适应辛普森积分法是利用二次函数对被积函数进行插值的一种模拟积分法,常用于估计积分的数值,利用递归自动适应所需精度. 公式推导基本思想公式: ∫abf(x)dx≈∫ab ...
定积分的计算与辛普森积分及龙贝格积分
定积分的计算与辛普森积分及龙贝格积分定积分的计算辛普森积分公式与自适应辛普森积分龙贝格积分(Romberg积分) 面积或者体积的计算 hdu1724 [NEERC 2014 D题 Damage ...
[NOI2005] 月下柠檬树（自适应辛普森积分）
题目原题链接:点这里总体思路–>问题转化先将原本的柠檬树分解成为多个圆台,再单独看圆台的投影一个圆在地面的投影,是等比例的,而一条竖线的投影长度d=htan⁡αd=\cfrac{h}{\ ...
[BZOJ1502][NOI2005]月下柠檬树（辛普森积分+解析几何）
题目: 我是超链接题解: 首先我们理解一下投影的性质,也就是投影出来的图形一定跟原图形全等. 考虑一下圆形投影下来是什么呢?和原来一样的圆形啊怎么转化竖着的树呢? 也就是树高*cot(α) 那么我 ...
[LOJ2267][SDOI2017]龙与地下城-FFT-自适应辛普森积分
龙与地下城 Description 小Q同学是一个热爱学习的人,但是他最近沉迷于各种游戏,龙与地下城就是其中之一. 在这个游戏中,很多场合需要通过掷骰子来产生随机数,并由此决定角色未来的命运,因此骰子 ...

hdu 1724 Ellipse——辛普森积分

hdu 1724 Ellipse——辛普森积分相关推荐

最新文章

热门文章