UVA - 1356 Bridge （辛普森公式）

Bridge

Time Limit: 3000MS

64bit IO Format: %lld & %llu

Submit Status uDebug

Description

题意：

桥上等距离摆放着塔高为H的若干个塔，忽略宽度，相邻塔的距离不超过D，塔之间绳索形成全等的对称抛物线。桥长度为B，绳索总长为L，求绳索最下端离地的高度y。

分析：

蓝书上提供了两种解法，一种凭借高数功底解积分，另一种则是根据辛普森公式去逼近积分。第一种普适性不强，这里着重说一下辛普森公式（感觉在做数值分析课的实验）。

首先有一点必须明确的是可导函数f(x)在区间[a,b]的弧长为∫√(1+[f'(x)]^2 dx 上下限分别为a，b。通过基本概念可求出对于宽度为w，高度为h的抛物线f(x)=ax^2,a=4h/w^2。从而得到弧长公式。

之后二分高度即可，将高度纳入计算弧长的函数计算，精度足够时输出即可。

计算弧长的函数是通过递归调用自适应辛普森公式对精确值进行逼近。辛普森公式就不多说了。

#include <algorithm>
#include <iostream>
#include <sstream>
#include <cstring>
#include <cstdlib>
#include <string>
#include <vector>
#include <cstdio>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;
#define INF 0x3f3f3f3f
const long long N=10005;
const long long mod=1e9;
const double eps=1e-14;
const double PI=acos(-1.0);
typedef long long ll;

int n;
double a;

double f(double x) {
    double a2=a*a,x2=x*x;
    return sqrt(1+4*a2*x2);
}

double simpson(double a,double b) {
    double c=a+(b-a)/2;
    return (f(a)+4*f(c)+f(b))*(b-a)/6;
}

double asr(double a,double b,double eps,double A) {
    double c=a+(b-a)/2;
    double l=simpson(a, c),r=simpson(c, b);
    if (fabs(l+r-A)<=15*eps) {
        return l+r+(l+r-A)/15;
    }
    return asr(a, c, eps/2, l)+asr(c, b, eps/2, r);
}

double asr(double a,double b,double eps) {
    return asr(a, b, eps, simpson(a, b));
}

double parabola_arec_length(double w,double h) {
    a=4.0*h/(w*w);
    return asr(0, w/2, 1e-5)*2;
}

int main() {
    int T;
    cin>>T;
    for (int kase=1; kase<=T; kase++) {
        int d,h,b,l;
        scanf("%d%d%d%d",&d,&h,&b,&l);
        int n=(b+d-1)/d;
        double d1=(double)b/n;
        double l1=(double)l/n;
        double x=0,y=h;
        while (y-x>1e-5) {
            double m=x+(y-x)/2;
            if (parabola_arec_length(d1, m)<l1) {
                x=m;
            } else y=m;
        }
        if (kase>1) {
            printf("\n");
        }
        printf("Case %d:\n%.2lf\n",kase,h-y);
    }
    return 0;
}