地狱飞龙 【simpson 求积分】

题目描述

最近clover迷上了皇室战争，他抽到了一种地狱飞龙，很开心。假设地域飞龙会对距离为d的敌人每秒造成k/d2伤害。假设地域飞龙位于坐标轴原点，以每秒v1的速度向y轴正方向移动，敌人在（x,0）的位置，以每秒v2的速度向x轴负方向移动。问，敌人至少有多少血量永远才不会被地狱飞龙喷死。（伤害是连续造成的，不是一秒一秒间断的）

输入

第一行为数据组数T（1<=T<=1000）
每组数据一行，包含4个实数，分别为v1,v2,x,k(1≤v1,v2,x,k≤10)。

输出

每组数据输出一行，为敌人最小血量，结果保留2为有效数字.

样例输入

1
1 1 1 1
样例输出

2.36
以下是对 Simpson’s rule 的改进算法simpson

代码

#include <bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
const int maxN = 100000;
const int maxM = 1e6;
const double PI = acos(-1.0);
const double eps = 1e-8;
inline int read(){
    int x=0,f=1; char ch=getchar();
    while (ch<'0'||ch>'9') { if (ch=='-') f=-1; ch=getchar(); }
    while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
    return x*f;
}
double v1,v2,x0,k;
inline double f(double x)  // 被积函数，据题而定
{
    return k/(v1*v1*x*x+(x0-v2*x)*(x0-v2*x));
}

inline double getappr(double a,double b)
{
    return (b-a)*(f(a)+f(b)+3*f((2*a+b)/3)+3*f((a+2*b)/3))/8.0;
 } 
double simpson(double l,double r)  // 近似计算 积分
{
    double sum=getappr(l,r);
    double mid=(l+r)/2;
    double suml=getappr(l,mid);
    double sumr=getappr(mid,r);
    return fabs(sum-suml-sumr)<eps?sum:simpson(l,mid)+simpson(mid,r);
 } 
 // 以上都是小模版，可以用来求积分
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf%lf%lf%lf",&v1,&v2,&x0,&k);
        printf("%.2lf\n",simpson(0,1e8));
    }
    return 0;
}