P1886 P2216 单调队列模板

何为单调队列？

单调队列是一个队列(废话)

而且必须同时满足下标单调和值单调两个单调特性。

跟优先队列不同，优先队列直接使用堆(heap)来实现，如何删去特定下标元素？不明。

本人喜欢用单调队列存下标，这样比存值不知道高到哪里去了。

新来一个元素，进队。之后特判长度有没有超过。超过则把最前面的元素出队。

之后，如果不满足性质，就把前面的元素顶掉，直到满足性质为止。

然后才可以把队首元素拿来用。

有几个很坑的地方，具体看代码注释。

(为什么top是队尾而tail是队首?laughcry)

例题：洛谷P1886，高级模板题。还有一道剧毒的P2216。(排版很奇怪，都怪反人类的博客园☺)

AC代码：

#include <cstdio>
const int N = 1000010;

int pl[N], l_t, l_h = 1, ansl[N];
int ps[N], s_t, s_h = 1, anss[N];
int a[N];

int main() {
    int n, k;
    scanf("%d%d", &n, &k);
    for(int i = 1; i <= n; i++) {
        scanf("%d", &a[i]);
    }
    for(int i = 1; i <= n; i++) {
        pl[++l_t] = i;
        ps[++s_t] = i;
        if(pl[l_h] + k <= i) {
            l_h++;
        }
        if(ps[s_h] + k <= i) {
            s_h++;
        }
        while(l_t > l_h && a[pl[l_t - 1]] <= a[i]) {
            l_t--;
        }
        while(s_t > s_h && a[ps[s_t - 1]] >= a[i]) {
            s_t--;
        }
        pl[l_t] = i;
        ps[s_t] = i;
        ansl[i] = pl[l_h];
        anss[i] = ps[s_h];
    }
    for(int i = k; i <= n; i++) {
        printf("%d ", a[anss[i]]);
    }
    printf("
");
    for(int i = k; i <= n; i++) {
        printf("%d ", a[ansl[i]]);
    }
    return 0;
}

#include <cstdio>
#include <algorithm>
using namespace std;
const int N=1002;
struct ah{int large,small;}f[N][N];
int c[N][N];
int p1[N],p2[N],top1,tail1=1,top2,tail2=1;///存下标
int a,b,n;///p1求min,p2求max
long long int ans=999999999,s;

int main()
{
    scanf ("%d%d%d",&a,&b,&n);
    for(int i=1;i<=a;i++)
    {
        top1=top2=0;
        tail1=tail2=1;
        for(int j=1;j<=b;j++)
        {
            scanf ("%d",&c[i][j]);
            p1[++top1]=p2[++top2]=j;
            while(c[i][p1[top1]]<=c[i][p1[top1-1]] && top1>tail1)
            {///单调队列维护
                p1[top1-1]=p1[top1];
                top1--;
            }
            while(c[i][p2[top2]]>=c[i][p2[top2-1]] && top2>tail2)
            {
                p2[top2-1]=p2[top2];
                top2--;
            }
            if(p1[top1]-p1[tail1]>=n) tail1++;///区间长度n
            if(p2[top2]-p2[tail2]>=n) tail2++;
            if(j>=n)
            {///产生最大最小值
                f[i][j-n+1].large=c[i][p2[tail2]];
                f[i][j-n+1].small=c[i][p1[tail1]];
            }
        }
    }
    ///读入完毕，处理①完毕
/*
    for(int i=1;i<=a;i++)///调试
    {
        for(int j=1;j<=b-n+1;j++)
        {
            printf("%d-%d ",f[i][j].large,f[i][j].small);
        }
        printf("
");
    }*/


    for(int j=1;j<=b-n+1;j++)
    {
        top1=top2=0;
        tail1=tail2=1;
        for(int i=1;i<=a;i++)///p1 small  p2 large
        {
            p1[++top1]=p2[++top2]=i;
            while(f[p1[top1]][j].small<=f[p1[top1-1]][j].small && top1>tail1)
            {
                p1[top1-1]=p1[top1];
                top1--;
            }
            while(f[p2[top2]][j].large>=f[p2[top2-1]][j].large && top2>tail2)
            {
                p2[top2-1]=p2[top2];
                top2--;
            }
            if(p1[top1]-p1[tail1]>=n) tail1++;
            if(p2[top2]-p2[tail2]>=n) tail2++;
            if(i>=n)
            {///产生最值
                s=1ll*(f[p2[tail2]][j].large-f[p1[tail1]][j].small);
                //printf("%d ",s);
                ans=min(ans,s);
            }

        }
        //printf("
");
    }
    printf("%lld",ans);
    return 0;
}

博客园排版有剧毒。。。