[HDU1724]Ellipse（辛普森积分）

题目描述

传送门

题解

辛普森积分裸题
答案要乘2吖，因为椭圆是两半

代码

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

const double eps=1e-10;
int dcmp(double x)
{
    if (x<=eps&&x>=-eps) return 0;
    return (x>0)?1:-1;
}
int T;
double a,b,l,r,ans;

double F(double x)
{
    return b*sqrt(1-((x*x)/(a*a)));
}
double calc(double l,double r)
{
    return (F(l)+F(r)+F((l+r)/2)*4)*(r-l)/6;
}
double simpson(double l,double r,double now)
{
    double mid=(l+r)/2;
    double p=calc(l,mid),q=calc(mid,r);
    if (dcmp(now-p-q)==0) return now;
    else return simpson(l,mid,p)+simpson(mid,r,q);
}
int main()
{
    scanf("%d",&T);
    while (T--)
    {
        scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
        ans=simpson(l,r,calc(l,r))*2;
        printf("%.3lf\n",ans);
    }
}