标签：OpenStreetMap head CF1195E ll while tail && include

题目链接

题意

有一个\(n\times m\)的矩阵，询问其中所有大小为\(a \times b\)的子矩阵的最小值之和。
\(1\le n,m \le 3000\)

思路

因为是子矩阵的大小是固定的。所以想到先将其中一维的最小值求出来，然后在此基础上再去求另外一维的最小值。

看数据范围不能带\(log\)。在每一维上单独做的时候就是一个滑动窗口。所以直接单调队列做刚好。

代码

/*
* @Author: wxyww
* @Date:   2019-07-18 14:37:59
* @Last Modified time: 2019-07-18 14:55:00
*/
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<ctime>
using namespace std;
typedef long long ll;
const int N = 3010;
ll read() {
    ll x=0,f=1;char c=getchar();
    while(c<'0'||c>'9') {
        if(c=='-') f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9') {
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
}
ll ans;
int t[N * N],a[N][N],b[N][N],q[N];
int main() {
    int n = read(),m = read(),nn = read(),mm = read();
    t[0] = read();int x = read(),y = read(),mod = read();
    for(int i = 1;i <= n * m;++i) t[i] = (1ll * t[i - 1] * x % mod+ y) % mod;

    for(int i = 1;i <= n;++i) {
        int head = 1,tail = 0;
        for(int j = 1;j <= m;++j) {
            a[i][j] = t[(i - 1) * m + j - 1];
            while(tail >= head && a[i][q[tail]] >= a[i][j]) --tail;
            q[++tail] = j;
            while(tail >= head && q[head] <= j - mm) ++head;
            if(j >= mm) b[i][j - mm + 1] = a[i][q[head]];
        }
    }



    for(int j = 1;j <= m - mm + 1;++j) {
        int head = 1,tail = 0;
        for(int i = 1;i <= n;++i) {
            while(tail >= head && b[q[tail]][j] >= b[i][j]) --tail;
                q[++tail] = i;
            while(tail >= head && q[head] <= i - nn) ++head;
            if(i >= nn) ans += b[q[head]][j];
        }
    }

    cout<<ans<<endl;

    return 0;
}

标签：OpenStreetMap,head,CF1195E,ll,while,tail,&&,include
来源： https://www.cnblogs.com/wxyww/p/CF1195E.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。