标签:sdoi rip source include nbsp 斜率优化 amp 状态 机器
这个题一开始想着去处理一个人等了多长时间
然后就会发现等的时间是和前面分了多少批有关的,这样的话我们需要用二维状态表示到这个点,前面分了多少批
这样暴力会很萎
这时候我们会回想起一个叫做修车的题目,他对于每个点的处理相当于是这个点让后面的人多等了多久
那我们可以通过同样的方式思考,每分了一批其实就是让后面的所有人多等了一个S的时间,其余的并不影响
那么我们可以推出一维的状态
f[i]=min(f[j]+s*(F[n]-F[j])+T[i]*(F[i]-F[j]));
对于前百分之60的数据,T为正数,满足决策单调性,可以二分栈,但我懒得打了,就在codevs上AC了一个弱化版的n^2;
斜率优化在下面
// MADE BY QT666 #include<cstdio> #include<algorithm> #include<cmath> #include<iostream> #include<queue> #include<set> #include<cstdlib> #include<cstring> #include<string> #include<ctime> #define lson num<<1 #define rson num<<1|1 #define int long long using namespace std; typedef long long ll; const int N=100050; int gi() { int x=0,flag=1; char ch=getchar(); while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘) flag=-1;ch=getchar();} while(ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-‘0‘,ch=getchar(); return x*flag; } int f[N],t[N],F[N],s,n; int cal(int j,int i){ return f[j]+s*(F[n]-F[j])+t[i]*(F[i]-F[j]); } main() { n=gi(),s=gi(); for(int i=1;i<=n;i++) t[i]=gi()+t[i-1],F[i]=gi()+F[i-1]; for(int i=1;i<=n;i++){ f[i]=cal(0,i); for(int j=1;j<i;j++){ f[i]=min(f[i],cal(j,i)); } } printf("%lld",f[n]); }
我们来推一推斜率方程:
f[j]+s*(F[n]-F[j])+T[i]*(F[i]-F[j)<=f[k]+s*(F[n]-F[k])+T[i]*(F[i]-F[k]);
f[j]+s*F[n]-s*F[j]+T[i]*F[i]-T[i]*F[j]<=f[k]+s*F[n]-s*F[k]+T[i]*F[i]-T[i]*F[k];
f[j]-s*F[j]-T[i]*F[j]<=f[k]-s*F[k]-T[i]*F[k];
令Y(j)=f[j]-s*F[j];
令X(j)=F[j];
则变为:
Y(j)-T[i]*X(j)<=Y(k)-T[i]*X(k);
T[i]*(X(k)-X(j))<=Y(k)-Y(j);
若X(k)>=X(j):T[i]<=(Y(k)-Y(j))/(X(k)-X(j));
若X(k)<=X(j);T[i]>=(Y(k)-Y(j))/(X(k)-X(j));
因为F值都是正的,所以X(k)>=X(j);横坐标单调;
所以T[i]<=(Y(k)-Y(j))/(X(k)-X(j));
然而斜率T[i]不是单调的
然后我就不会做了
标签:sdoi rip source include nbsp 斜率优化 amp 状态 机器
原文地址:http://www.cnblogs.com/qt666/p/6524321.html