博客园 首页 私信博主 显示目录 隐藏目录 管理 动画

洛谷.4525.[模板]自适应辛普森法1(Simpson积分)

题目链接

Simpson积分公式:$$\int_a^bf(x)dx\approx\frac{b-a}{6}\left[f(a)+f(b)+4f(\frac{a+b}{2})\right]$$

推导过程大概就是,令\(f(x)=Ax^2+Bx+C\),代进去求一下积分就好了?
自适应是指根据区间大小控制精度。满足精度要求时直接返回。

那个,有人知道asr是指什么吗。。(什么的缩写?)

//0ms	1.7MB
#include <cmath>
#include <cstdio>
#include <algorithm>
#define eps (1e-6)

double a,b,c,d,L,R;

inline double f(double x){
	return (c*x+d)/(a*x+b);
}
inline double Simpson(double l,double r){
	return (r-l)*(f(r)+f(l)+4*f((l+r)*0.5))/6.0;
}
double asr(double l,double r,double Eps,double ans)
{
	double mid=(l+r)*0.5, lans=Simpson(l,mid), rans=Simpson(mid,r);
	if(fabs(lans+rans-ans)<Eps) return lans+rans;
	return asr(l,mid,Eps*0.5,lans)+asr(mid,r,Eps*0.5,rans);
}

int main()
{
	scanf("%lf%lf%lf%lf%lf%lf",&a,&b,&c,&d,&L,&R);
	printf("%.6lf",asr(L,R,eps,Simpson(L,R)));

	return 0;
}
posted @ 2018-06-11 08:49  SovietPower  阅读(261)  评论(0编辑  收藏  举报