CF #574 Div. 2 E. OpenStreetMap // 单调双端队列

CF #574 Div. 2  E. OpenStreetMap // 单调双端队列

https://codeforc.es/contest/1195/problem/E

题意：给m*n矩阵，求Σ(所有a*b的子矩阵中的最小元素)

思路：

先横着用单调双端队列扫一遍，得到1*b的所有矩阵的最小值。

然后纵着用单调双端队列扫一遍上述的最小值，就可以得到每一个小矩阵的最小元素。

单调双端队列和单调栈一样：元素单调，存的是下标。

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define FI first
#define SE second
#define MP make_pair
#define PII pair<int,int>
const LL mod = 1e9+7;
const int LG = 50;
const int MX = 3e3+10;
int g[MX][MX];
int G[MX][MX];
int main(){
    LL n,m,a,b,x,y,z,now;
    cin>>n>>m>>a>>b>>now>>x>>y>>z;
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)g[i][j]=now,now=(now*x+y)%z;
    for(int i=1;i<=n;i++){
        deque<int>q;
        for(int j=1;j<b;j++) {while(!q.empty()&&g[i][q.back()]>=g[i][j])q.pop_back();q.push_back(j);}
        for(int j=b;j<=m;j++){
            while(!q.empty()&&q.front()<=j-b) q.pop_front();
            while(!q.empty()&&g[i][q.back()]>=g[i][j])q.pop_back();q.push_back(j);
            G[i][j]=g[i][q.front()];
        }
    }
    /*for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)cout<<g[i][j]<<' ';
        puts("");
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++)cout<<G[i][j]<<' ';
        puts("");
    }*/
    LL ans=0;
    for(int i=b;i<=m;i++){
        deque<int>q;
        for(int j=1;j<a;j++) {while(!q.empty()&&G[q.back()][i]>=G[j][i]) q.pop_back();q.push_back(j);}
        for(int j=a;j<=n;j++){
            while(!q.empty()&&q.front()<=j-a)q.pop_front();
            while(!q.empty()&&G[q.back()][i]>=G[j][i]) q.pop_back();q.push_back(j);
            ans=ans+(LL)G[q.front()][i];
        }
    }
    cout<<ans<<'\n';
    return 0;
}