【LA3485】 Bridge

前言

哈哈哈，垃圾微积分哈哈哈
前置知识：自适应Simpson法与微积分初步，学会编程

Solution

考虑一下我们有的是什么：

一段桥梁的横向距离，悬线的长度，以及高度。

我们发现如果我们重新设一个元(h=H-x)，那么很显然就有(x=H-h)

这个时候感觉(h)好求一些

那么怎么做？

可以根据(h)和(d)算出悬线长度，然后直接二分答案就好了。（注意按照格式输出）

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<iostream>
using namespace std;
#define ll long long
#define re register
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
inline int gi(){
	int f=1,sum=0;char ch=getchar();
	while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
	return f*sum;
}
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 simpson(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)<=eps*15)return L+R+(L+R-A)/15;
    return simpson(a,c,eps/2,L)+simpson(c,b,eps/2,R);
}
double Ask(double a,double b,double eps){
    return simpson(a,b,eps,Simpson(a,b));
}
double check(double d,double h){
	a=4.0*h/(d*d);
	return Ask(0,d/2,1e-5)*2;
}
const double eps=1e-5;
int main(){
	int T=gi();
	for(int Case=1;Case<=T;Case++){
		printf("Case %d:
",Case);
		int D,H,B,L;
		scanf("%d%d%d%d",&D,&H,&B,&L);
		int n=(D+B-1)/D;
		double d=(double)B/n;
		double l=(double)L/n;
		double x=0,y=H;
		while(y-x>eps){
			double mid=x+(y-x)*0.5;
			if(check(d,mid)<l)x=mid;
			else y=mid;
		}
		printf("%.2lf
",H-x);
		if(Case!=T)puts("");
	}
	return 0;
}