洛谷.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;

}