地狱飞龙 [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

#include <stdio.h> 
#include <math.h>
#include <algorithm>
#define eps 1e-8
using namespace std;
double v1,v2,x0,k;
double qwe(double x)
{
	return k/(v1*x*v1*x+(x0-v2*x)*(x0-v2*x));
}
double asd(double a,double b)
{
	return (b-a)*(3*qwe((a+2*b)/3)+3*qwe((b+2*a)/3)+qwe(a)+qwe(b))/8.0;//改用了 Simpson’s 3/8 rule
}
double sum(double a,double b)
{
	double k1=asd(a,b);
	double mid=(a+b)/2;
	double k2=asd(a,mid),k3=asd(mid,b);
	return fabs(k1-k2-k3)<eps?k1:(sum(a,mid)+sum(mid,b));
}
int main()
{
	int t;
	double l;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%lf%lf%lf%lf",&v1,&v2,&x0,&k);
		l=sum(0,1e9);
		printf("%.2f\n",l);
	}
}