洛谷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 ;

 }