斜率优化题目（持续更新）

bzoj1597
关键词：排序+单调队列预处理、斜率优化
题意：将n个矩形分成若干组，每一组花费是最长边*最宽边（长宽不能互换），求分组的最小花费。
解法:
1.将n个矩形按长度从小到达排序，再按宽度求出单调递减队列，剩余的矩形可以并入这些留下的组中，而使花费不变！
2.对于剩下的矩形，x递增，y递减。dp[i]:前i个矩形分组的最小花费。dp[i]=min(dp[j]+y[j+1]*x[i]),0<=j<=i-1，然后进行斜率优化。
3.设k< j，且j优于k，即dp[j]+y[j+1]*x[i]< dp[k]+y[k+1]*x[i],可推出dp[j]-dp[k]<(y[k+1]-y[j+1])*x[i].
令up(j,k)=dp[j]-dp[k],down(j,k)=y[k+1-y[j+1]
则当k< j< i,且up(i,j) * down * (j,k)< down(i,j) *up(j,k)时，j不能取得最优值。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<limits.h>
#define mem(a,b) memset(a,sizeof(a),b)
#define ll long long
#define MP make_pair
const ll INF = 0x0fffffffffffffff;
using namespace std;
typedef long long lld;


const int maxn = 50000+10;

int n;
struct node{
    ll x,y;//x是长度，y是宽度
}a[maxn];
int tot;
ll x[maxn],y[maxn];
int que2[maxn],tail2,head2;
ll dp[maxn];

bool cmp(node a,node b){
    return (a.x==b.x)?a.y<b.y:a.x<b.x;
}

ll getup(int i,int j){
    return dp[i]-dp[j];
}

ll getdown(int i,int j){
    return y[j+1]-y[i+1];
}

int main(){
    scanf("%d",&n);
    a[0].y=INF;
    for(int i=1;i<=n;i++){
        scanf("%lld%lld",&a[i].x,&a[i].y);
    }
    sort(a+1,a+n+1,cmp);
    tot=0;
    //得到x递增，y递减的序列
    for(int i=1;i<=n;i++){
        while(tot>0&&y[tot]<=a[i].y) tot--;
        x[++tot]=a[i].x,y[tot]=a[i].y;
    }
    head2=0,tail2=-1;
    que2[++tail2]=0,dp[0]=0;
    for(int i=1;i<=tot;i++){
        while(head2+1<=tail2&&getup(que2[head2+1],que2[head2])<x[i]*getdown(que2[head2+1],que2[head2])) head2++;
        int t=que2[head2];
        dp[i]=dp[t]+y[t+1]*x[i];
        while(head2+1<=tail2&&getup(i,que2[tail2])*getdown(que2[tail2],que2[tail2-1])<=getdown(i,que2[tail2])*getup(que2[tail2],que2[tail2-1])) tail2--;
        que2[++tail2]=i;
    }
    printf("%lld\n",dp[tot]);
}

bzoj3675 序列分割
关键词：滚动数组、斜率优化
将n长序列分割k次，每次分割得到两个序列，增益是分割成的序列和之积，求最大增益
dp[i][k]:前i个数分割k次后的最大增益
dp[i][k]=max{ dp[j][k-1]+sum[j]*(sum[i]-sum[j]),1<=j< i },j！=0且j!=i（分割序列一定非空！）
简化：
1.滚动数组。省略第二维数组,k->1,(k-1)->0
2.斜率优化。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<limits.h>
#define mem(a,b) memset(a,sizeof(a),b)
#define ll long long
#define MP make_pair
const ll INF = 0x0fffffffffffffff;
using namespace std;
typedef long long lld;

const int maxn = 100000+10;

int n,k;
ll a[maxn];
int head,tail,que[maxn];
ll sum[maxn];
ll now[maxn],old[maxn];

ll up(int j,int k){
    return old[k]-old[j]-(sum[k]*sum[k]-sum[j]*sum[j]);
}

ll down(int j,int k){
    return sum[j]-sum[k];
}

ll dp(int i,int j){
    return old[j]+sum[j]*(sum[i]-sum[j]);
}

void DP(int x){
    head=0,tail=-1;
    //分割为非空序列，不能让j为0
   // que[++tail]=0;***
    //old[0]=0;***
    for(int i=x;i<=n;i++){
        while(head+1<=tail&&up(que[head+1],que[head])<sum[i]*down(que[head+1],que[head])) head++;
        now[i]=dp(i,que[head]);
        while(head+1<=tail&&up(i,que[tail])*down(que[tail],que[tail-1])<down(i,que[tail])*up(que[tail],que[tail-1])) tail--;
        que[++tail]=i;
    }
    for(int i=1;i<=n;i++) old[i]=now[i];
}

int main(){
    //freopen("a.txt","r",stdin);
    scanf("%d%d",&n,&k);
    sum[0]=0;
    for(int i=1;i<=n;i++) { scanf("%lld",&a[i]); }
    int tot=0;
    for(int i=1;i<=n;i++){
        if(a[i]!=0){
            a[++tot]=a[i];
            sum[tot]=sum[tot-1]+a[tot];
        }
    }
    n=tot;
    for(int i=1;i<=k;i++) DP(i);
    printf("%lld\n",old[n]);
}