hdu 1724 Ellipse —— 自适应辛普森积分

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef double db;
db const eps=1e-5;
db a,b,l,r;
db f(db x){return b*sqrt(1-(x*x)/(a*a));}
db simp(db l,db r){return (r-l)/6*(f(l)+4*f((l+r)/2)+f(r));}
db asr(db l,db r,db eps,db lst)
{
  db mid=(l+r)/2;
  db ls=simp(l,mid),rs=simp(mid,r);
  if(fabs(ls+rs-lst)<=15*eps)return (ls+rs+(ls+rs-lst)/15);//
  return asr(l,mid,eps/2,ls)+asr(mid,r,eps/2,rs);
}
db asme(db l,db r,db eps){return asr(l,r,eps,simp(l,r));}
int main()
{
  int T; scanf("%d",&T);
  while(T--)
    {
      scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
      printf("%.3f\n",asme(l,r,eps)*2);
    }
  return 0;
}