Codeforces Round #574 (Div. 2) E. OpenStreetMap （单调队列）

最新推荐文章于 2022-11-22 23:15:14 发布

_kikyou-

最新推荐文章于 2022-11-22 23:15:14 发布

阅读量82

点赞数

分类专栏：思维

本文链接：https://blog.csdn.net/qq_46030630/article/details/116125187

版权

思维专栏收录该内容

42 篇文章 1 订阅

订阅专栏

题目链接

题意：给定了一个最大为3000 * 3000的矩形，问里面每个a * b小矩形内的最小值之和。 (1<=a,b<=n,m)

分析：先对每行用单调队列维护长度为b的区间的最小值,存进row[i][j]里面,
row[i][j]表示第i行中,[j-b+1,j]区间内的最小值。
再对row数组的第b,b+1…m列每列都用单调队列维护长度为a的区间的最小值。

这样得到的就是一个a*b的小矩形的最小值。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
const int maxn = 3e3+10;
const int mx = 40;
const int mod = 1e9+7;
const ll inf = 34359738370;
const int INF = 1e9+7;
const double e = 2.71828;
int n,m,a,b;
ll h[maxn][maxn];
ll g[maxn*maxn];
int x,y,z;
int row[maxn][maxn];//每一行中 长度为b的最小值
ll ans=0;
void solve1(int x)//h矩阵 第x行的区间最小值
{
	deque<int> q;
	//先维护[1,b-1] ，这样从b位置开始 每次出队就是一个答案
	for(int i=1;i<b;i++)
	{
		while(q.size() && h[x][q.back()] > h[x][i]) q.pop_back();
		q.push_back(i);
	}
	for(int i=b;i<=m;i++)
	{
		while(q.size() && i-q.front() >= b) q.pop_front();
		while(q.size() && h[x][i] < h[x][q.back()]) q.pop_back();
		q.push_back(i);
		row[x][i]=h[x][q.front()];
	}
}
void solve2(int y)//列 对row的第y列进行维护
{
	deque<int> q;
	for(int i=1;i<a;i++)
	{
		while(q.size() && row[i][y] < row[q.back()][y]) q.pop_back();
		q.push_back(i);
	}
	for(int i=a;i<=n;i++)
	{
		while(q.size() && i-q.front()>= a) q.pop_front();
		while(q.size() && row[i][y] < row[q.back()][y]) q.pop_back();
		q.push_back(i);
		ans+=row[q.front()][y];
	}
}
int main()
{
    int t;t=1;
	while(t--)
	{
		cin>>n>>m>>a>>b;
		scanf("%I64d %d %d %d",&g[0],&x,&y,&z);
		for(int i=1;i<=2*m*n;i++) g[i]=(g[i-1]*x+y)%z;
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=m;j++)
			{
				h[i][j]=g[(i-1)*m+j-1];
			}
		}
		for(int i=1;i<=n;i++)
		{
			solve1(i);
		}
		for(int i=b;i<=m;i++) solve2(i);
		cout<<ans<<'\n';
	}
    return 0;
}