自适应辛普森(近似求积分)
2017-04-24 15:36
253 查看
其实就是写个板子,写一下这个公式的具体实现,可以做 hdu 1724 Ellipse、haut oj 1261地狱飞龙 练下手
具体代码:(代码中是求y=a*x+b在[L,R]上的积分)
具体代码:(代码中是求y=a*x+b在[L,R]上的积分)
#include <bits/stdc++.h> using namespace std; const double eps=1e-7; double a,b,l,r; int t; double F(double x) // F[x] { return a*x+b; } double simpson(double a,double b) { double c=a+(b-a)/2; return (F(a)+4*F(c)+F(b))*(b-a)/6; } double asr(double a,double b,double eps,double A) { double c=a+(b-a)/2; double L=simpson(a,c),R=simpson(c,b); if(fabs(L+R-A)<=15*eps) return L+R+(L+R-A)/15.0; return asr(a,c,eps/2,L)+asr(c,b,eps/2,R); } double asr(double a,double b,double eps) { return asr(a,b,eps,simpson(a,b)); } int main() { scanf("%d",&t); while(t--) { scanf("%lf%lf%lf%lf",&a,&b,&l,&r); printf("%.3f\n", asr(l,r,eps)); } return 0; }
相关文章推荐
- 自适应辛普森(近似求积分模板)
- 自适应辛普森积分
- 数值积分 (辛普森公式 辛普森自适应法则) UVA
- 自适应辛普森积分算法
- HDU 1724 Ellipse (自适应辛普森积分)
- HDU 1724 Ellipse 自适应辛普森积分
- UVALive - 3485 Bridge 二分+数学+自适应辛普森求积分
- hdu 1724 Ellipse (自适应辛普森求积分 模板题)
- 数值积分-自适应辛普森法
- 数值积分_自适应辛普森法
- ZOJ 2369 | SGU 217 | Acdream 1234 自适应辛普森积分
- 洛谷.4525.[模板]自适应辛普森法1(Simpson积分)
- BZOJ 1502 月下柠檬树 自适应辛普森积分
- 【自适应辛普森积分】hdu1724 Ellipse
- 自适应辛普森积分
- BZOJ 1502 NOI 2005 月下柠檬树 计算几何 自适应辛普森积分
- [LOJ2267][SDOI2017]龙与地下城-FFT-自适应辛普森积分
- SPOJ CIRU 圆的面积并 自适应辛普森积分
- LA 3485 (积分 辛普森自适应法) Bridge
- HDU 1724 Ellipse (自适应辛普森积分模板题)