这题看了题解一句话恍然大悟!!!
如果从前往后推的话dp方程是有后效性的,但是从后往前考虑第i个对最后结果的贡献,其为处理第i个所用的时间再乘上从第i个开始到最后的f的和。
所以设t[i]为从1到i的t的和,f[i]为1到i的f的和
方程为dp[i] = min(dp[i],dp[j] + (m+t[j-1]-t[i-1])*(f[n]-f[i-1]))。
用斜率优化,注意这个凹折线是从后往前维护的。
我的代码:
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const LL maxn = 10005;
const LL INF = 0x3f3f3f3f;
LL n,m,t[maxn],f[maxn];
LL dp[maxn],q[maxn],hd,tl;
LL getDP(LL i,LL j){
return dp[j] + (m+t[j-1]-t[i-1])*(f[n]-f[i-1]);
}
LL getY(LL j,LL k){
return (dp[k]+f[n]*t[k-1])-(dp[j]+f[n]*t[j-1]);
}
LL getX(LL j,LL k){
return t[k-1] - t[j-1];
}
void solve(){
hd = tl = 0;
q[tl++] = n+1;
for(LL i=n;i>=1;i--){
while(hd+1<tl && getY(q[hd+1],q[hd])
>= f[i-1]*getX(q[hd+1],q[hd])) hd++;
dp[i] = getDP(i,q[hd]);
while(hd+1<tl && getY(i,q[tl-1])*getX(q[tl-1],q[tl-2])
>= getY(q[tl-1],q[tl-2])*getX(i,q[tl-1])) tl--;
q[tl++] = i;
}
printf("%I64d\n",dp[1]);
}
int main(){
while(~scanf("%I64d%I64d",&n,&m)){
memset(dp,0,sizeof(dp));
for(LL i=1;i<=n;i++){
scanf("%I64d%I64d",&t[i],&f[i]);
t[i] += t[i-1];f[i] += f[i-1];
}
solve();
}
return 0;
}