洛谷P4525 【模板】自适应辛普森法1(simpson积分…
2018-06-17 20:29:42来源:未知 阅读 ()
题目描述
计算积分
结果保留至小数点后6位。
数据保证计算过程中分母不为0且积分能够收敛。
输入输出格式
输入格式:
一行,包含6个实数a,b,c,d,L,R
输出格式:
一行,积分值,保留至小数点后6位。
输入输出样例
1 2 3 4 5 6
2.732937
说明
a,b,c,d∈[-10,10]
-100≤L<R≤100 且 R-L≥1
辛普森积分是用$y = Ax^2 +Bx +c$去拟合给定的函数
$$\int_a^bf(x)dx\approx\frac{(b-a)(f(a)+f(b)+4f(\frac{a+b}{2}))}{6}$$
自适应的含义是根据不同的区间大小选用不同的eps
// luogu-judger-enable-o2 #include<cstdio> #include<cmath> double a, b, c, d, L, R; double F(double x) { return (c * x + d) / (a * x + b); } double sim(double l, double r) { return (F(l) + F(r) + 4 * F((l + r) / 2)) * (r - l) / 6; } double asr(double L, double R, double eps, double ans) { double mid = (L + R) / 2; double LL = sim(L, mid), RR = sim(mid, R); if(fabs(LL + RR - ans) < eps) return LL + RR; else return asr(L, mid, eps / 2, sim(L, mid)) + asr(mid, R, eps / 2, sim(mid, R)); } main() { #ifdef WIN32 freopen("a.in", "r", stdin); #endif scanf("%lf %lf %lf %lf %lf %lf", &a, &b, &c, &d, &L, &R); printf("%lf", asr(L, R, 1e-6, sim(L, R))); }
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
- C++冒泡排序 (基于函数模板实现) 2020-05-31
- C++ 模板类vector 2020-05-31
- C++ 模板类array 2020-05-31
- C++ 模板类vector 2020-05-30
- 洛谷P1164->小A点菜 2020-05-18
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash