洛谷P2365 任务安排 [解法二 斜率优化]

解法一：http://www.cnblogs.com/SilverNebula/p/5926253.html

解法二：斜率优化

在解法一中有这样的方程：dp[i]=min(dp[i],dp[j]+(sumf[i]-sumf[j])*sumt[i]+s*(sumf[n]-sumf[j]) )

其中min的后半部分，也就是dp[j]+(sumf[i]-sumf[j])*sumt[i]+s*(sumf[n]-sumf[j]) 计算了将j~i分为一组的花费（以及提前计算的受影响花费）

设f(j)=dp[j]+(sumf[i]-sumf[j])*sumt[i]+s*(sumf[n]-sumf[j])，i不变时，若 f(j1)<f(j2) ，显然从j1到i分为一组比j2到i分为一组的答案更优，而如果j1<j2，显然j2可以被舍弃掉。由以上两个限制条件很容易联想到单调队列，进而想到斜率优化（并不）。

现在来考虑 j1<j2 ,f(j1)<f(j2) 的情况。把f()展开写再化简，可以得到(dp[j1]-dp[j2])/(sumf[j1]-sumf[j2])<=sumt[i]+s    (sumf和sumt分别是f、t的前缀和）

利用这个式子列斜率方程，维护一个下凸壳即可←然而并不能理解

我的想法：(dp[j1]-dp[j2])/(sumf[j1]-sumf[j2])显然是越小越好，我们可以据此维护斜率单调队列的队尾(具体看代码)，而上面那个式子用来维护队头，即可行：

斜率优化10ms，O（n^2）算法43ms

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cmath>
 #include<queue>
 #include<cstring>
 using namespace std;
 const int mxn=;
 int n;
 int s;
 int t[mxn],f[mxn];
 int sumt[mxn],sumf[mxn];
 int dp[mxn];
 int q[mxn];
 int gup(int j,int k){
     return (dp[j]-dp[k]);
 }
 int gdown(int j,int k){
     return sumf[j]-sumf[k];
 }
 int gdp(int i,int j){
     return dp[j]+(sumf[i]-sumf[j])*sumt[i]+s*(sumf[n]-sumf[j]);
 }
 int main(){
     scanf("%d%d",&n,&s);
     int i,j;
     for(i=;i<=n;i++){
         scanf("%d%d",&t[i],&f[i]);
         sumt[i]=sumt[i-]+t[i];
         sumf[i]=sumf[i-]+f[i];
     }
     memset(dp,0x3f,sizeof dp);
     dp[]=;
     int hd=,tl=;
     q[hd]=;
     for(i=;i<=n;i++){
         while(hd<tl && gup(q[hd],q[hd+])>=(sumt[i]+s)*gdown(q[hd],q[hd+]) )
             hd++;
         dp[i]=gdp(i,q[hd]);
         while(hd<tl && gup(i,q[tl])*gdown(q[tl],q[tl-])<=gup(q[tl],q[tl-])*gdown(i,q[tl]) )tl--;
         q[++tl]=i;
     }
     printf("%d",dp[n]);
     return ;
 }