hdu 1724（Simpson积分）

传送门
模板题。。。
注意最后答案乘以2，因为阴影部分由上下两部分共同组成
公式：

图片来自：http://blog.csdn.net/linraise/article/details/18624685

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
#define eps 1e-10
#define dcmp(x) ((x>-eps)-(x<eps))
double a,b,l,r;
inline double f(double x) {//x^2/a^2+y^2/b^2=1
    return b*sqrt(1.0-x*x/(a*a));
}
inline double cal(double l,double r) {
    double mid=l+(r-l)/2.0;
    return (r-l)*(f(l)+f(mid)*4.0+f(r))/6.0;
}
double simpson(double l,double r) {
    double mid=(l+r)/2.0;
    double a=cal(l,mid),b=cal(mid,r),p=cal(l,r);
    if (dcmp(p-a-b)==0) return p;
    return simpson(l,mid)+simpson(mid,r);
}
int main() {
//  freopen("hdu 1724.in","r",stdin);
    int kase;
    scanf("%d",&kase);
    while (kase--) {
        scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
        printf("%.3lf\n",simpson(l,r)*2.0);
    }
    return 0;
}