Codeforces Round #574 (Div. 2) E. OpenStreetMap 【单调队列】

一、题目

　　OpenStreetMap 

二、分析

　　对于二维空间找区间最小值，那么一维的很多好用的都无法用了，这里可以用单调队列进行查找。

　　先固定一个坐标，然后进行一维的单调队列操作，维护一个区间长度为$b$的最小值，可以确定一个数组$ans$。

　　对于新数组，进行另一个维度的单调队列操作，找到$a*b$的最小值。（由于英语太差，读题目时理解成了是点在区间左上角的意思）

　　写的时候一定注意第二次单调队列的对象数组是新数组。

三、AC代码

#include <bits/stdc++.h>

using namespace std;
#define ll long long
#define Min(a,b) ((a)>(b)?(b):(a))
#define Max(a,b) ((a)>(b)?(a):(b))
const int maxn = 3e3 + 13;
ll h[maxn][maxn];
ll ans[maxn][maxn];

int que[maxn];
int ql, qr;

int main()
{
    //freopen("input.txt", "r", stdin);
    int n, m, a, b;
    int x, y, z;
    while(scanf("%d%d%d%d", &n, &m, &a, &b) != EOF)
    {
        int i, j;
        ll g0;
        scanf("%I64d%d%d%d", &g0, &x, &y, &z);
        for(i = 1; i <= n; i++)
        {
            for(j = 1; j <= m; j++)
            {
                h[i][j] = g0;
                g0 = (1ll*g0*x + y)%z;
            }
        }        
        for(i = 1; i <= n; i++)
        {
            ql = 0, qr = 0;
            for(j = 1; j <= m; j++)
            {
                while(ql < qr && j - que[ql] >= b)
                    ql++;
                while(ql < qr && h[i][j] <= h[i][que[qr-1]])
                    qr--;
                que[qr++] = j;
                ans[i][j] = h[i][que[ql]];
            }
        }
        ll tmp = 0;
        for(j = b; j <= m; j++)
        {
            ql = 0, qr = 0;
            for(i = 1; i <= n; i++)
            {
                while(ql < qr && i - que[ql] >= a)
                    ql++;
                while(ql < qr && ans[i][j] <= ans[que[qr-1]][j])
                    qr--;
                que[qr++] = i;
                if(i >= a)
                    tmp += ans[que[ql]][j];
            }
        }
        printf("%I64d
", tmp);
    }
    return 0;
}