hdu 5144 NPY and shot(三分)

http://acm.hdu.edu.cn/showproblem.php?pid=5144
题目概要：一个人站在高度为H的地方以V的速度扔出球，以最好的角度扔，以便球能扔得够远，求扔出去的最远距离。
解析：物理题：已知高度H和速度V，最好的角度angle。
设小球竖直的速度为v1,水平速度为v2，v1=sin(angle)V,v2=cos(angle)V
设小球飞行的时间为t,t=t1+t2,t1为小球竖直速度向上的时间，t2为小球竖直速度向下的时间。
t1=v1/g;
t2=sqrt((v1t1-0.5gt1t1+h)2/g);
s=v2(t1+t2);
求最好的角度，速度与水平所成的角度为a(0<=a<=90(π/2))，在这个范围中a是单峰极值函数，可以使用三分查找。

#include<stdio.h>
#include<stack>
#include<string.h>
#include<algorithm>
#include<math.h>
const double pi=3.1415926;
const double g=9.8;
const double eps=1e-8;
double h,v;
using namespace std;
double calo(double t){
	double t1,t2,h1;
	double v1=sin(t)*v,v2=cos(t)*v;
	t1=v1/g;
	t2=sqrt((v1*t1-0.5*g*t1*t1+h)*2/g);
	double s=v2*(t1+t2);
	return s;
}
double san(double l,double r){
	double min,max,mid,mid_mid;
	while(l+eps<r){
		mid=(l+r)/2;
		mid_mid=(mid+r)/2;
		if(calo(mid)>=calo(mid_mid))
			r=mid_mid;
		else 
			l=mid;
	}
	return l;
}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%lf%lf",&h,&v);
		double angle=san(0,pi/2);
		printf("%.2lf\n",calo(angle));
	}
	
}