LA 3485 Bridge

原创 2013年12月04日 18:30:27

          自适应辛普森公式模板。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
#include<map>
#include<set>
#define LL long long
#define CLR(a, b) memset(a, b, sizeof(a))
using namespace std;

double a;

double F(double x)///需要积分的式子。
{
    return sqrt(1 + 4 * a * a * x * x);
}

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.0;
    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_arc_length(double w, double h)
{
    a = 4.0 * h / (w * w);
    return asr(0, w / 2, 1e-5) * 2;
}

int main()
{
    int t, cas = 1;
    scanf("%d", &t);
    while(t --)
    {
        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_arc_length(D1, m) < L1) x = m;
            else y = m;
        }
        if(cas > 1) puts("");
        printf("Case %d:\n%.2lf\n", cas ++, H - x);
    }
    return 0;
}


版权声明:本文为博主原创文章,未经博主允许不得转载。转载请标注:blog.csdn.net/ok_again

相关文章推荐

La3485 利用自适应辛普森公式 求积分的和

#include #include #include #include #include #include #include #define maxn 10005 #define INF 0xffff...

UVALive-3485 Bridge(二分答案+自适应辛普森积分)

前言本题是我在《训练指南》上发现的一道有趣的题目,有两种做法,一种是直接运用数学方法求积分,对这题而言较为繁琐;另一种则是运用自适应辛普森积分算法,简洁地求解此题。《训练指南》一书中还有很多比较好的题...
  • Ab_Ever
  • Ab_Ever
  • 2017年08月03日 09:02
  • 142

[ACM] POJ 3485 Highway (区间选点问题)

Highway Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 871   Accepte...

uva 3485 Optimal Array Multiplication Sequence

题目: I - Optimal Array Multiplication Sequence Time Limit:3000MS     Memory Limit:0KB     64bit IO ...

TTL与RS-485电平转换芯片MAX485/MAX3485

常用的转换芯片是MAX485(5V),MAX3485(3.3V),它们都是半双工的。        A端和B端分别为接收和发送的差分信号端,当A引...

uvalive 3485 - 数值积分

/** 抛物线积分求曲线长度。 首先确定出桥墩的个数。 题目要求桥墩数量要最少,也就是看D是不是能被总长度B整除, 若可以则桥墩数量n = B/D, 否则 n = B/D + 1, 一个简化的写法就是...

SP3485 pdf

  • 2011年11月22日 21:08
  • 512KB
  • 下载

sp3485+pcb+sch

  • 2013年10月30日 12:27
  • 4KB
  • 下载

社会化统计工具分享之Google分析、cnzz、51la、百度统计和Js比较

流量统计作为互联网站点最基础的数据分析工具,几乎与每个站长、互联网从业者相关。在互联网这个没有硝烟的战场上,流量统计工具就像打仗用的武器。能否选好适合自己的武器、用好武器,很大程度上决定了我们的胜败。...

MAX485和MAX3485模块资料

  • 2017年05月20日 14:55
  • 1.19MB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:LA 3485 Bridge
举报原因:
原因补充:

(最多只允许输入30个字)