【辛普森积分*模板】2017河工大校赛 B地狱飞龙

1261: 地狱飞龙

时间限制: 1 秒  内存限制: 64 MB
提交: 184  解决: 27
提交  状态 

题目描述

最近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为有效数字.

样例输入

11 1 1 1

样例输出

2.36

提示

来源

提交  状态 
终于知道了和队友的差距，打完比赛回到房间我是开始打王者，队友把校赛题目拿出来继续看题；队友过了一会儿说B题是辛普森积分，直到这个玩意儿但是没有模拟过，所以就没优先考虑这道题；

我是不懂什么是辛普森积分的，上百度了解了公式和代码； 辛普森积分公式就是对于不能够求原函数的积分近似求值的；

这道题被积函数也是比较容易列出来，代码完全就是套的模板； 代码wa了三次，是因为精度不够； 代码也是参考别人的，之后会作为我的模板；
代码：

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const double eps=0.000001;

double v1,v2,x,k;

double f(double t)
{
    double d=v1*t*v1*t+(x-v2*t)*(x-v2*t);  //被积公式，这个过程没必要分段；
    return k/d;
}

double get_sub(double l,double r)
{
    return (f(l)+4.0*f((l+r)/2)+f(r))*(r-l)/6.0;
}

double simpson(double l,double r)
{
    double sum=get_sub(l,r);
    double mid=(l+r)/2;
    double suml=get_sub(l,mid);
    double sumr=get_sub(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,&x,&k);
        double ans=simpson(0,20000);
        printf("%.2lf\n",ans);
    }
    return 0;
}