HAUTOJ 1261: 地狱飞龙 自适应辛普森算法

1261: 地狱飞龙

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

题目描述

最近clover迷上了皇室战争，他抽到了一种地狱飞龙，很开心。假设地域飞龙会对距离为d的敌人每秒造成k/d^2伤害。假设地域飞龙位于坐标轴原点，以每秒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

本题其实就是求地狱飞龙对敌人造成的伤害总量，由于伤害是连续的，因此结果就是求k/d^2对时间t积分，即 

k是已知的，，求出这个积分就可以求出结果了。

但是这个积分式并不好求，因此考虑采用自适应辛普森算法求积分。

自适应辛普森算法原理可以参考辛普森算法

#include<bits/stdc++.h>
using namespace std;
double v1,v2,k,x;
double f(double t)//自定义函数（积分函数）
{
	double y;
	y=k/((x-v2*t)*(x-v2*t)+v1*v1*t*t);
	return y;
}

//自适应辛普森法
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));
}

int main()
{
	std::ios::sync_with_stdio(false);
	int T;
	double ans;
	cin>>T;
	while(T--)
	{
		cin>>v1>>v2>>x>>k;
		ans=asr(0,10000000,0.0001);
		printf("%.2f\n",ans);
	}
	return 0;
}