BZOJ3695: 滑行

题目大意：给一个平面，宽度为X，从上到下划分成了N个区域，每个区域有一个高度和速度，问从(0,0)到(X,Y)最快时间

这题我看了好长时间才看懂题面再说什么...

看一下发现和光的折射很像...

然后我们可以二分初始角度，然后根据折射定律可以直接算出终点

然后啊，这题十分非常特别的坑

二分的时候必须类似我这么写...

为什么呢？因为asin()函数里面可能不在定义域内，然后返回的值就是NaN，然后所有判断条件都是false的，我猜标程是这么写的，所以如果反过来就是错的....

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
double pi=acos(-1);
double ans;
double h[110],v[110];
int n,X;
double judge(double a)
{
	int i,j;
	ans=0;
	double nowx=0;
	for(i=1;i<=n;i++)
	{
		nowx+=h[i]*tan(a);
		ans+=(h[i]/cos(a))/v[i];
		if(i!=n)
		a=asin(v[i+1]*sin(a)/v[i]);
	}
	return nowx;
}
int main()
{
	scanf("%d%d",&n,&X);
	int i;
	for(i=1;i<=n;i++)
	scanf("%lf",&h[i]);
	for(i=1;i<=n;i++)
	scanf("%lf",&v[i]);
	double l=0,r=pi/2,mid;
	while(r-l>1e-13)
	{
		mid=(l+r)/2;
		if(judge(mid)>=X) r=mid;
		else l=mid;
	}
	printf("%.3lf",ans);
}