思路
单调队列板子。
设 \(b_{i,j}\) 表示第 \(i\) 行,区间为 \(j\) 到 \(j+y-1\) 的最小值,不难发现这个可以用单调队列 \(O(nm)\) 预处理出来。
接下来我们的问题就转化成了有一个矩阵,求所有宽为 \(1\),长为 \(x\) 的子矩阵的最小值之和。
这个问题也十分好做,直接枚举列,然后在每一行单调队列即可。
代码
/*
A tree without skin will surely die.
A man without face is invincible.
*/
#include<bits/stdc++.h>
using namespace std;
int const N=3e3+10;
int ans,n,m,a[N][N],b[N][N],g[N*N];
struct node{int w,id;}q[N];
inline void go(int x,int y){
for (int i=1;i<=n;++i){
int l=1,r=0;
for (int j=1;j<=y;++j){
while (l<=r && q[r].w>=a[i][j]) --r;
q[++r].w=a[i][j];q[r].id=j;
}
b[i][1]=q[l].w;
for (int j=y+1;j<=m;++j){
while (q[l].id<=j-y && l<=r) ++l;
while (l<=r && q[r].w>=a[i][j]) --r;
q[++r].w=a[i][j];q[r].id=j;
b[i][j-y+1]=q[l].w;
}
}
for (int j=1;j<=m-y+1;++j){
int l=1,r=0;
for (int i=1;i<=x;++i){
while (l<=r && q[r].w>=b[i][j]) --r;
q[++r].w=b[i][j];q[r].id=i;
}
ans+=q[l].w;
for (int i=x+1;i<=n;++i){
while (q[l].id<=i-x && l<=r) ++l;
while (l<=r && q[r].w>=b[i][j]) --r;
q[++r].w=b[i][j];q[r].id=i;ans+=q[l].w;
}
}
}
signed main(){
//读入
for (int i=1;i<=n*m;++i) g[i]=(g[i-1]*xx%zz+yy)%zz;
for (int i=1;i<=n;++i)
for (int j=1;j<=m;++j) a[i][j]=g[(i-1)*m+j-1];
go(x,y);//输出
return 0;
}
原文地址:http://www.cnblogs.com/tx-lcy/p/16791321.html
1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长!
2. 分享目的仅供大家学习和交流,请务用于商业用途!
3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入!
4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解!
5. 如有链接无法下载、失效或广告,请联系管理员处理!
6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需!
7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员!
8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载
声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性