HDU-#1724 Ellipse（自适应Simpson法+数值积分）

       题目大意：给出椭圆公式以及a,b值，再给定一个垂直x轴的两条直线l,r，求与椭圆围成的阴影面积。

      解题思路：这是一道数值积分的裸题，直接利用自适应Simpson公式就可求解。由于该题是模板题，code给出了详细分析，这里就不赘述了。具体的分析可以参见白书上的数学基础章节，如果需要学习理论知识【Simpson积分法】。因此利用该公式可以求出阴影部分的上半部分的阴影面积，乘以2即是结果。详见code。

      题目来源：http://acm.hdu.edu.cn/showproblem.php?pid=1724

      code：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

const double eps = 1e-6;
int a,b,l,r,n;
double ans;

double F(double x){ //数值积分函数 
    return sqrt(b*b*(1-x*x/(a*a)));
}

double simpson(double a,double b){ //三点simpson法
    double c=a+(b-a)/2;
    return (F(a)+4*F(c)+F(b))*(b-a)/6; //三点simpson公式
}

double asr(double a,double b,double eps,double A){// 自适应Simpson公式（递归过程）
    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; //eps为精度，用于算法自适应划分区间
    return asr(a,c,eps/2,L)+asr(c,b,eps/2,R);
}

double asr(double a,double b,double eps){ //自适应Simpson公式（主过程）
    return asr(a,b,eps,simpson(a,b)); //递归调用
}

int main(){
    scanf("%d",&n);
    while(n--){
        scanf("%d%d%d%d",&a,&b,&l,&r);
        ans=asr(l,r,eps); //调用自适应Simpson的主过程
        ans*=2; //该椭圆是上下两部分对称的，因此求出上部分需要乘以2
        printf("%.3lf\n",ans);
    }
    return 0;
}