hdu 1724 Ellipse simpson积分

时间:2024-01-16 22:52:08
/*
hdu 1724 Ellipse simpson积分 求椭圆的部分面积
simpson积分法 http://zh.wikipedia.org/zh-tw/%E8%BE%9B%E6%99%AE%E6%A3%AE%E7%A7%AF%E5%88%86%E6%B3%95
*/
#include<stdio.h>
#include<math.h>
const double eps=1e-8;
inline double jue(double a)
{
return a>0?a:-a;
}
int dbcmp(double a)
{
if(a>eps)return 1;
if(a<(-eps)) return -1;
return 0;
}
template<typename T>
inline T sqr(T d){return d*d;}
double a,b;
inline double func(double x)
{
return 2*b*sqrt(1-sqr(x)/sqr(a));
}
double simpson(double l,double r)
{
return (func(l)+func(r)+4*func((l+r)/2))/6*(r-l);
}
double work(double l,double r)
{
double mid=(l+r)/2;
double ans=simpson(l,r);
if(dbcmp(ans-simpson(l,mid)-simpson(mid,r))==0)
{
return ans;
}else
{
return work(l,mid)+work(mid,r);
}
}
int main()
{
int n;
double l,r;
scanf("%d",&n);
while(n--)
{
scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
double ret=work(l,r);
printf("%.3f\n",ret);
}
return 0;
}