CF1195E OpenStreetMap
2019-08-16 07:51:00来源:博客园 阅读 ()
CF1195E OpenStreetMap
题目链接
题意
有一个\(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;
}
原文链接:https://www.cnblogs.com/wxyww/p/CF1195E.html
如有疑问请与原作者联系
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash