自适应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;
}
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;
}