网址:http://codeforces.com/problemset/problem/1195/E
给出一个矩阵,求以$(x,y)$为左下角,长度为$(a,b)$的子矩阵的元素的最大值。
见以下博客,同类型,可照搬:
https://www.cnblogs.com/Aya-Uchida/p/11332822.html
每行使用单调队列处理出行的最大值,然后得到的矩阵再对每列处理出列最大值,然后查询的时候直接输出即可。
AC代码:
#include <bits/stdc++.h> #pragma GCC optimize(2) using namespace std; int num[3005][3005]; int minn[3005][3005]; int ans[3005][3005]; int n,m,a,b,g,x,y,z; void data_gen() { num[0][0]=g; long long tmp=g; for(int i=0;i<n;++i) { int j=i?0:1; for(;j<m;++j) { tmp=(tmp*x+y)%z; num[i][j]=tmp; } } } int main() { cin>>n>>m>>a>>b>>g>>x>>y>>z; data_gen(); deque<pair<int,int> >que; for(int i=0;i<n;++i) { while(!que.empty()) que.pop_front(); for(int j=0;j<b;++j) { while(!que.empty()&&que.back().first>num[i][j]) que.pop_back(); que.push_back(make_pair(num[i][j],j)); } minn[i][0]=que.front().first; for(int j=b;j<m;++j) { while(!que.empty()&&que.back().first>num[i][j]) que.pop_back(); que.push_back(make_pair(num[i][j],j)); while(!que.empty()&&que.front().second<j-b+1) que.pop_front(); minn[i][j-b+1]=que.front().first; } } for(int i=0;i<m-b+1;++i) { while(!que.empty()) que.pop_front(); for(int j=0;j<a;++j) { while(!que.empty()&&que.back().first>minn[j][i]) que.pop_back(); que.push_back(make_pair(minn[j][i],j)); } ans[0][i]=que.front().first; for(int j=a;j<n;++j) { while(!que.empty()&&que.back().first>minn[j][i]) que.pop_back(); que.push_back(make_pair(minn[j][i],j)); while(!que.empty()&&que.front().second<j-a+1) que.pop_front(); ans[j-a+1][i]=que.front().first; } } long long res=0; for(int i=0;i<n-a+1;++i) for(int j=0;j<m-b+1;++j) res+=ans[i][j]; cout<<res<<endl; return 0; }
codeforces1195E OpenStreetMap 单调队列
原文:https://www.cnblogs.com/Aya-Uchida/p/11332856.html