【子矩阵——优先队列】

题目

代码

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int N = 1e3 + 10;
const int mod = 998244353;
int g[N][N], ymax[N][N], ymin[N][N];
deque<int> qxmax[N], qxmin[N], qymax[N], qymin[N];

int main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    
    int n, m, a, b;
    cin >> n >> m >> a >> b;

    // 初始化矩阵
    for (int i = 1; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            cin >> g[i][j];
        }
    }

    // 垂直方向的滑动窗口(按列计算)
    for (int j = 1; j <= m; j++) {
        for (int i = 1; i <= n; i++) {
            // 最大值队列
            while (!qymax[j].empty() && g[qymax[j].back()][j] <= g[i][j]) qymax[j].pop_back();
            qymax[j].push_back(i);

            // 最小值队列
            while (!qymin[j].empty() && g[qymin[j].back()][j] >= g[i][j]) qymin[j].pop_back();
            qymin[j].push_back(i);

            // 窗口内数据有效
            if (i >= a) {
                if (i - qymax[j].front() + 1 > a) qymax[j].pop_front();
                if (i - qymin[j].front() + 1 > a) qymin[j].pop_front();
                ymax[i][j] = g[qymax[j].front()][j]; // 存储具体值
                ymin[i][j] = g[qymin[j].front()][j]; // 存储具体值
            }
        }
    }

    // 水平方向的滑动窗口(按行计算)
    ll res = 0;
    for (int i = a; i <= n; i++) {
        for (int j = 1; j <= m; j++) {
            // 最大值队列
            while (!qxmax[i].empty() && ymax[i][qxmax[i].back()] <= ymax[i][j]) qxmax[i].pop_back();
            qxmax[i].push_back(j);

            // 最小值队列
            while (!qxmin[i].empty() && ymin[i][qxmin[i].back()] >= ymin[i][j]) qxmin[i].pop_back();
            qxmin[i].push_back(j);

            // 窗口内数据有效
            if (j >= b) {
                
                if (j - qxmax[i].front() + 1 > b) qxmax[i].pop_front();
                if (j - qxmin[i].front() + 1 > b) qxmin[i].pop_front();
                
                res = (res + (ll)ymax[i][qxmax[i].front()] * ymin[i][qxmin[i].front()]) % mod;
            }
        }
    }

    // 输出结果
    cout << res << endl;

    return 0;
}
相关推荐
老师我太想进步了20266 小时前
【数据结构-4】算法美学之顺序表
数据结构·算法
风影同志783246 小时前
编码器是激活函数吗?—— 一个看起来很像、其实不是的类比
算法
UIU1147 小时前
P1481 魔族密码
c++·算法·动态规划
加宝的起点雷达10 小时前
40 条候选里 38 条来自同一个源,主卡整屏重样
人工智能·算法·创业
秋天的怀念_David10 小时前
算法一千零一夜 · 第二十六夜 协议的母语
java·spring boot·算法
闻缺陷则喜何志丹13 小时前
【C++特性】四种类型转
c++·算法·类型转换
喵本喵叁肆13 小时前
用 Boosting 把弱隐私算法炼成强隐私算法(英文原题:Boosting and Differential Privacy (FOCS 2010))
算法·集成学习·boosting
青梅味猪大肠13 小时前
【操作系统-60】临近适应算法
java·linux·开发语言·嵌入式硬件·算法·操作系统