cf1195E. OpenStreetMap [单调队列]

传送门
题意：给出一个n * m的矩阵(n,m<=3000)，求其所有a*b的子矩阵内的最小值之和

题解：对于每一列倒序维护一个单调递增的单调队列可以得到以每个格子为起点的长度为a的子矩阵的最小值，然后对于每行倒序对这些得到的最小值维护一个单调递增的单调队列即可得到每个格子为左上角的大小为a*b的子矩阵的最小值，求和即可(维护单调队列的题目要注意正序还是倒序)

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
#define debug(x) cout<<#x<<" is "<<x<<endl;

const int maxn=3e3+5;

ll h[maxn][maxn],g[maxn*maxn],ac[maxn],pos[maxn],ak[maxn][maxn];

int main(){
    ll n,m,a,b,g0,x,y,z;
    scanf("%lld%lld%lld%lld%lld%lld%lld%lld",&n,&m,&a,&b,&g0,&x,&y,&z);
    g[0]=g0;
    for(int i=1;i<n*m;i++)g[i]=(g[i-1]*x+y)%z;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            h[i][j]=g[(i-1)*m+(j-1)];
        }
    }

    for(int i=1;i<=m;i++){
        int head,tail;
        head=1;
        tail=0;
        for(int j=n;j>=1;j--){
            while(head<=tail&&ac[tail]>=h[j][i])tail--;
            ac[++tail]=h[j][i];
            pos[tail]=j;
            if(j>n-a+1)continue;
            while(head<=tail&&pos[head]-j>=a)head++;
            ak[j][i]=ac[head];
        }
    }
    ll ans=0;
    for(int i=1;i<=n-a+1;i++){
        int head,tail;
        head=1;
        tail=0;
        for(int j=m;j>=1;j--){
            while(head<=tail&&ac[tail]>=ak[i][j])tail--;
            ac[++tail]=ak[i][j];
            pos[tail]=j;
            if(j>m-b+1)continue;
            while(head<=tail&&pos[head]-j>=b)head++;
            ans+=ac[head];
        }
    }
    printf("%lld\n",ans);
    return 0;
}