地狱飞龙

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 245    Accepted Submission(s): 83


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

Input
第一行为数据组数T(1<=T<=1000)
每组数据一行,包含4个实数,分别为 v1,v2,x,k(1v1,v2,x,k10)
 

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

Sample Input
1 1 1 1 1
 

Sample Output
2.36
 

实质为对k/(x^2-2*v2*x*t+(v1^2+v2^2)t^2)进行不定积分


#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
typedef long long ll;
double k,v1,x,v2;
double getans(ll t){
    double ans = k/(v1*x)*atan((v1*v1+v2*v2)/(v1*x)*(t-(v2*x)/(v1*v1+v2*v2)));//这里是求出来的关于t的积分公式,真是太强了。
    return ans;
}
int main(void)
{
    ll MAX = 999999999;
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%lf %lf %lf %lf",&v1,&v2,&x,&k);
        double ans = getans(MAX) - getans(0);//由于是不定积分,这里的无穷我们用MAX就能达到精度要求了。
        printf("%.2lf\n",ans);
    }
    return 0;
}