[ 题外话 : = =看了一眼题目就知道是求最值
然后就不会做了╮(╯▽╰)╭ 所以,数学渣就去学了一发拉格朗日乘数法 - -]
那么上正文TAT 由于公式太多,我就直接截图哈qwq
附代码:
#include "stdio.h"
#include "iostream"
#define rep(f,a,b) for(f=a;f<=b;f++)
using namespace std;
const double eps=1e-12;
const double inf=1e5;
const int N=10005;
void read(int &v){
char ch=getchar(); while (ch<'0'||ch>'9') ch=getchar();
while (ch<='9'&&ch>='0') { v=v*10+ch-'0'; ch=getchar();}
}
int n; double eu;
double s[N],k[N],vi[N];
double v[N],atime,lagra;
double speed(double y,int t,double &x){
double bot=max(0.0,y),top=inf;
while (x=(bot+top)/2,top-bot>eps){
if (2*lagra*k[t]*x*x*(x-y)<=1) bot=x; else top=x;
}
}
bool check(){
int i; double esum=0.0;
rep(i,1,n) {
speed(vi[i],i,v[i]);
esum+=k[i]*(v[i]-vi[i])*(v[i]-vi[i])*s[i];
}
return esum>=eu;
}
int main(){
//freopen("bicycling.in","r",stdin);
//freopen("bicycling.out","w",stdout);
read(n); int i;
scanf("%lf",&eu);
rep(i,1,n) {
scanf("%lf%lf%lf",s+i,k+i,vi+i);
} double bot=0.0,top=inf;
while (lagra=(bot+top)/2,top-bot>eps){
if (check()) bot=lagra;
else top=lagra;
}
rep(i,1,n) atime+=s[i]/v[i];
printf("%.10lf\n",atime);
return 0;
}