http://www.lydsy.com/JudgeOnline/problem.php?id=1096
最基础的斜率DP
我的推导&草稿:
dp[i]:i建设工厂,那么之前所有费用的最小值
朴素方程:dp[i]=min{C[i]+dp[j]+P[j+1]*dis[j+1][i]+P[j+2]*dis[j+2][i]+...+P[i]*dis[i][i]}
求dp[n]
令d[i]=dis[1][i]
dp[i]=min{C[i]+dp[j]+P[j+1]*(d[i]-d[j+1])+P[j+2]*(d[i]-d[j+2])+...+P[i]*(d[i]-d[i])}
dp[i]=min{C[i]+dp[j]+(P[j+1]+P[j+2]+...+P[i])*d[i]-P[j+1]*d[j+1]-P[j+2]*d[j+2]-...-P[i]*d[i]}
令s[i]=P[1]+...+P[i] , f[i]=P[i]*d[i]
dp[i]=min{C[i]+dp[j]+(s[i]-s[j])*d[i]-(f[j+1]+f[j+2]+...+f[i])}
令t[i]=f[1]+...+f[i]
dp[i]=min{C[i]+dp[j]+(s[i]-s[j])*d[i]-(t[i]-t[j])}
dp[i]=min{dp[j]-s[j]*d[i]+t[j]}+C[i]+s[i]*d[i]-t[i]
讨论dp[j]-s[j]*d[i]+t[j]
设j < k,j 比 k 优
dp[j]-s[j]*d[i]+t[j] < dp[k]-s[k]*d[i]+t[k]
d[i] < ((dp[k]+t[k]) - (dp[j]+t[j])) / (s[k] - s[j])
Y[k]=dp[k]+t[k]
X[k]=s[k]
d[i] < (Y[k]-Y[j]) / (X[k]-X[j])
//#define _TEST _TEST #include <cstdio> #include <cstring> #include <cstdlib> #include <iostream> #include <cmath> #include <algorithm> using namespace std; /************************************************ Code By willinglive Blog:http://willinglive.cf ************************************************/ #define rep(i,l,r) for(int i=l,___t=(r);i<=___t;i++) #define per(i,r,l) for(int i=r,___t=(l);i>=___t;i--) #define MS(arr,x) memset(arr,x,sizeof(arr)) #define LL long long #define INE(i,u,e) for(int i=head[u];~i;i=e[i].next) inline const int getint() { int r=0,k=1;char c=getchar(); for(;c<'0'||c>'9';c=getchar())if(c=='-')k=-1; for(;c>='0'&&c<='9';c=getchar())r=r*10+c-'0'; return k*r; } ///////////////////////////////////////////////// int n; LL p[1000010],c[1000010]; LL d[1000010],s[1000010],t[1000010]; LL dp[1000010]; int l,r,q[1000010]; ///////////////////////////////////////////////// inline LL X(int k){return s[k];} inline LL Y(int k){return dp[k]+t[k];} double slope(int k,int j){return (double)(Y(k)-Y(j))/(X(k)-X(j));} ///////////////////////////////////////////////// void input() { scanf("%d",&n); rep(i,1,n) { d[i]=getint(), p[i]=getint(), c[i]=getint(); s[i]=s[i-1]+p[i]; t[i]=t[i-1]+p[i]*d[i]; } } void solve() { rep(i,1,n) { while(l<r&&d[i]>=slope(q[l],q[l+1]))l++; int j=q[l]; dp[i]=dp[j]-s[j]*d[i]+t[j]+c[i]+s[i]*d[i]-t[i]; while(l<r&&slope(q[r-1],q[r])>=slope(q[r],i))r--; q[++r]=i; } cout<<dp[n]<<endl; } ///////////////////////////////////////////////// int main() { #ifndef _TEST freopen("std.in","r",stdin); freopen("wl.out","w",stdout); #endif input(), solve(); return 0; }