单调队列优化dp

我的理解:

形如dp[i]=max/min(f[k])+g[i]     (k<i  且g[i]与k无关)时，可用单调队列优化。

优化的对象是f[k];

 

本来 i和k是要用两个for循环来套完，但是在此处，g[i]与k，无关，而f[k]又要滚动取最值，那么就可以在一个for循环中把f[k]得到最值，同时计算当前的g[i],(强调，for循环的遍历顺序要确定，要保证可能与g[i]产生结果的所有f[k]都已计算完毕)将两者的结果存入dp[i]中，就达到了降维的目的。

 

补充:上述中的k可以当成i来看，因为是与i一同遍历。

例题 CDOj 594

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>

using namespace std;
#define INF 0xfffffff
int dp[2][110]; //两个110的数组滚动优化，轮流代表上一个
int n,c;
int q[110];//存与上一个有关的最值
int head,tail,cur;


int main()
{
    int i,j,x,nowf;
    while(scanf("%d%d",&n,&c)==2)
    {
        scanf("%d",&x);
        cur=0;
        for(i=0;i<x;i++)
            dp[cur][i]=INF;
        for(i=x;i<=100;i++)
            dp[cur][i]=(x-i)*(x-i);
        for(i=1;i<n;i++){
            scanf("%d",&x);
            cur=1-cur;//切换数组

            //比上一个大
            head=tail=0;
            for(j=0;j<=100;j++){    //这个循环将两个互不干涉但对同一值起作用的两个循环套在一起，是单调队列优化的精髓（降维）
                nowf=dp[1-cur][j]-j*c;  //这个j表示的是遍历上一数组
                while(head<tail && q[tail-1]>nowf)
                    tail--;
                q[tail++]=nowf; //不断记录最值
                if(j<x){
                    dp[cur][j]=INF;
                }
                else {
                    dp[cur][j]=q[head]+j*c+(x-j)*(x-j);
                }
            }

            //比上一个小
            head=tail=0;
            for(j=100;j>=0;j--){
                nowf=dp[1-cur][j]+j*c;
                while(head<tail && q[tail-1]>nowf)
                    tail--;
                q[tail++]=nowf;
                if(j>=x){
                    dp[cur][j]=min(dp[cur][j],q[head]-j*c+(x-j)*(x-j));
                }
            }

        }
        int ans=INF;
        for(i=0;i<=100;i++){
            ans=min(ans,dp[cur][i]);
        }
        printf("%d\n",ans);
    }
    return 0;
}

前面是学的时候照着网上敲的，后来发现记录极值好像完全不用开始组。。。。。

此题可以直接只存极值，但有些题中对于k有额外的限制条件的时候们就不能简单的只存一个值，比如hdu3401，还要存之前的数量，就要使用这种队列优化

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>

using namespace std;
#define INF 0xfffffff
int dp[2][110];
int n,c;
int head,tail,cur;

int main()
{
    int i,j,x,nowf;
    int mi=INF;
    while(scanf("%d%d",&n,&c)==2)
    {
        scanf("%d",&x);
        cur=0;
        for(i=0;i<x;i++)
            dp[cur][i]=INF;
        for(i=x;i<=100;i++)
            dp[cur][i]=(x-i)*(x-i);
        for(i=1;i<n;i++){
            scanf("%d",&x);
            cur=1-cur;
            head=tail=0;
            mi=INF;
            for(j=0;j<=100;j++){
                nowf=dp[1-cur][j]-j*c;
                mi=min(mi,nowf);
                if(j<x){
                    dp[cur][j]=INF;
                }
                else {
                    dp[cur][j]=mi+j*c+(x-j)*(x-j);
                }
            }
            mi=INF;
            head=tail=0;
            for(j=100;j>=0;j--){
                nowf=dp[1-cur][j]+j*c;
                mi=min(mi,nowf);
                if(j>=x){
                    dp[cur][j]=min(dp[cur][j],mi-j*c+(x-j)*(x-j));
                }
            }
        }
        int ans=INF;
        for(i=0;i<=100;i++){
            ans=min(ans,dp[cur][i]);
        }
        printf("%d\n",ans);
    }
    return 0;
}