hdu 1724 Ellipse——辛普森积分

最新推荐文章于 2025-09-07 00:18:35 发布

weixin_34179762

最新推荐文章于 2025-09-07 00:18:35 发布

阅读量107

点赞数

CC 4.0 BY-SA版权

文章标签： php

原文链接：http://www.cnblogs.com/Narh/p/10141896.html

本文提供了一种使用C++实现的辛普森积分法解决方案，用于解决HDU 1724题目的数学积分问题。通过定义函数并采用递归方式优化积分精度，确保了计算结果的准确性。

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

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