自适应辛普森积分

自适应Simpson积分

近几场训练赛碰到了有关积分的题，于是去看了看。

应用在于有的积分函数求原函数很难求，所以可以利用插值法将所求的整体分为若干个整体，最基本的化整为零的例子就是小学学过的黎曼法，分为若干小矩形，不过误差很大。而自适应就是根据你定义的精度来划分为更小的区间，递归求解，误差很小，我们有了积分函数和区间就可以利用这个simpson积分了。

来看两道例题吧：

Ellipse

 simpson裸题，求椭圆某一部分面积。

我们根据椭圆公式可以推出y的函数式，而且这是个高度对称的图形，我们求出上半部分的面积即可。

y=b*b-(x*x*b*b)/a*a。注意精度至少0.0001

double F(double x)      //这里自定义函数
{
    return sqrt(b-b*x*x/a);
}
double simpson(double a,double b)   //返回区间[a,b]的三点辛普森值
{
    double c = a + (b - a) / 2.0;
    return (F(a) + 4 * F(c) + F(b)) * (b - a) / 6.0;
}
double asr(double a,double b,double eps,double A)   //自适应辛普森递归过程
{
    double c = a + (b - a) / 2.0;                   //A为区间[a,b]的三点辛普森值
    double L = simpson(a,c), R = simpson(c,b);
    if(fabs(L + R - A) <= 15 * eps) return L + R + (L + R - A) / 15.0;
    return asr(a,c,eps/2.0,L) + asr(c,b,eps/2.0,R);
}
double asr(double a,double b,double eps)    //自适应辛普森主过程
{
    return asr(a,b,eps,simpson(a,b));
}

int main()
{
   int t;
   scanf("%d",&t);
   while(t--)
   {
       scanf("%lf%lf%lf%lf",&a,&b,&l,&r);
       a*=a,b*=b;
       printf("%.3f\n",2*asr(l,r,0.0001));
   }
   return 0;
}

                                                                            UVALive - 3485 

Bridge

 

 二分+数学+simpson积分。

题意：一个铁索桥，长度为B,桥墩之间的距离至少为D，桥墩的高度是H，铁索的总长度为L。题目保证所有数据合法。

求最少放置桥墩的情况下，铁索离地面的高度是多少。

二分铁索的深度h，可知h与铁索的长度是正相关的，然后每个桥墩之间的铁索长度都可以求出来，铁索自然态肯定是一个抛物线状的，即：F(x)=aX^2；

假设两个桥墩之间的距离W，那么h=a*(W/2)^2。在利用积分去验证长度关系。

怎么积分呢，不会，看题解的，还是感觉不严谨。吐槽一下UVA 没有PE，只有WA。

int d,h,b,L;
double l,tc;
double F(double x)
{
    return sqrt(1+4*tc*tc*x*x);
}
double simpon(double a,double b)
{
    double c=a+(b-a)/2.0;
    return (F(a)+4*F(c)+F(b))*(b-a)/6.0;
}
double asr(double a,double b,double eps,double A)
{
    double c=a+(b-a)/2.0;
    double li=simpon(a,c),R=simpon(c,b);
    if(fabs(li+R-A)<=15*eps) return li+R+(li+R-A)/15.0;
    return asr(a,c,eps/2.0,li)+asr(c,b,eps/2.0,R);
}
double asr(double a,double b,double eps)
{
    return 2*asr(a,b,eps,simpon(a,b));
}
int main()
{
    int t;
    scanf("%d",&t);
    int t1=t;
    while(t--){
      scanf("%d%d%d%d",&d,&h,&b,&L);
      int x=b/d;
      if(b%d) x++;
      l=L*1.0/x;
      double w=b*1.0/x;
      double left=0,rig=h,ans;
      while(left+eps<=rig)
      {
          double mid=(left+rig)/2.0;
          tc=4*mid/(w*w);
          if(asr(0,w/2.0,eps)>=l) ans=mid,rig=mid-eps;
          else left=mid+eps;
      }
      printf("Case %d:\n%.2f\n",t1-t,h-ans);
      if(t) puts("");
    }
    return 0;
}