[NOI2012]骑行川藏（拉格朗日乘数法）

【题解】

显然最优情况下要用完所有能量，而且对于 Ei=ki*si*(xi-vi)^2，当abs(xi-vi)一定时，xi>vi比xi<vi用时更短 

转化为：寻找变量受一个条件所限制的多元函数的极值 的问题 
即：在 sigma(ki*si*(Xi-vi)^2) = Eu 的条件下，求 sigma(si/xi)的最小值 
用拉格朗日乘数法解决 

构造方程组 
发现：对于方程 -si/(Xi^2) + λ* 2*ki*si*(Xi-vi) =0，转化为 λ*2*ki * Xi^2 *(Xi-vi) = 1 (λ>0)后，Xi随λ的增大而减小 
因此二分λ解出所有Xi，带入 sigma(ki*si*(Xi-vi)^2) 与 Eu 比较大小即可 

【代码】

#include<stdio.h>
#include<stdlib.h>
#define INF 1e9
#define eps 1e-12
double s[10005]={0},k[10005]={0},v[10005]={0},x[10005]={0};
double Eu;
int n;
double jdz(double a)
{
	if(a<0) a=0-a;
	return a;
}
double max(double a,double b)
{
	if(a>b) return a;
	return b;
}
double pd(double lam)//λ
{
	double sum=0,left,right,mid;
	int i;
	for(i=1;i<=n;i++)
	{
		left=max(0,v[i]);
		right=INF;
		while(right-left>eps)//二分求x[i]：[ λ*2*ki*Xi^2*(Xi-vi) ]关于x[i]单调递增 
		{
			mid=(left+right)/2;
			if(jdz(lam*2*k[i]*mid*mid*(mid-v[i])-1)<=eps) break;
			if(lam*2*k[i]*mid*mid*(mid-v[i])>1) right=mid;
			else left=mid;
		}
		x[i]=mid;
		sum+=k[i]*s[i]*(x[i]-v[i])*(x[i]-v[i]);
	}
	return sum;
}
void erfen()
{
	double left=0,right=INF,mid;
	while(right-left>eps)
	{
		mid=(left+right)/2;
		if(pd(mid)>Eu) left=mid;
		else right=mid;
	}
}
int main()
{
	double t=0;
	int i;
	scanf("%d%lf",&n,&Eu);
	for(i=1;i<=n;i++)
		scanf("%lf%lf%lf",&s[i],&k[i],&v[i]);
	erfen();
	for(i=1;i<=n;i++)
		t+=s[i]/x[i];
	printf("%.8lf",t);
	return 0;
}