codeforces 407D Largest Submatrix 3
codeforces 407D Largest Submatrix 3
题意
找出最大子矩阵,须满足矩阵内的元素互不相等。
题解
官方做法
http://codeforces.com/blog/entry/11333
- \(O(n^6)\)
- 枚举子矩阵,暴力check。
- \(O(n^4)\)
- 枚举上下边界,双指针。
- \(O(n^3log_2n)\)
- 假设当前上边界 \(up\), 下边界 \(down\),\(R_i\) 表示当 \(i\) 为左边界时,右边界最大是 \(R_i\)。
- 当 \(down->down+1\) 时,\(R_i\) 要么不变,要么减小。
- 减小的情况:\(a_{down+1, j}(i<=j<=R_i)\) 在矩阵 \((up, i, down, R_i)\) 中出现。
- 因此对于每个 \(a_{down+1, j}(i<=j<=R_i)\),我们需要寻找在矩阵 \((up, i, down, R_i)\) 中出现的,从左边最靠近 \(j\),和从右边最靠近 \(j\) 的两个位置。这个可以用set维护。
- \(O(n^3)\)
- 如何优化掉 \(O(n^3log_2n)\) 做法的 \(log_2n\) ?
- 从大到小枚举 \(up\)。当 \(up->up-1\) 时,\(O(n^2)\) 更新 \((up, 1, n, m)\) 中所有点的从左边最靠近,和从右边最靠近的两个位置。
区间dp做法
\(f_{i, l, r}\)表示下边界是 \(i\), 左右边界是 \(l, r\) 时,上边界的值。
\(f_{i, l, r}=max\{f_{i-1, l, r}, f_{i, l+1, r}, f_{i, l, r-1}, a_{i, l}与a_{1-i, r}的限制, a_{i, r}与a_{1-i, l}的限制\}\)
可以这样理解:\(f_{i, l, r}\) 需要维护A、B、C、D、E、F这六个块之间的关系,\(f_{i-1, l, r}\) 只维护了A、B、C之间的关系,\(f_{i, l+1, r}\) 只维护了B、C、E、F之间的关系,\(f_{i, l, r-1}\) 只维护了A、B、D、E之间的关系,剩下的还需维护A与F、C与D、D与F的关系。
// 代码实现可以省掉第一维。
代码
\(O(n^3log_2n)\)做法
T了 不知道有没有写错
#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define rep(i, a, b) for(int i=(a); i<(b); i++)
#define sz(x) (int)x.size()
#define de(x) cout<< #x<<" = "<<x<<endl
#define dd(x) cout<< #x<<" = "<<x<<" "
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
const int N=404;
int n,m;
int a[N][N], R[N], cnt[N*N];
set<int> pos[N*N];
int solve() {
int ans=1;
rep(up,1,n+1) {
rep(i,0,N*N) pos[i].clear();
rep(i,0,m+1) R[i]=m;
rep(down,up,n+1) {
int r=0;
rep(i,1,m+1) {
while(r<m&&!cnt[a[down][r+1]]) {
++r;
++cnt[a[down][r]];
}
R[i]=min(R[i], r);
--cnt[a[down][i]];
}
rep(j,1,m+1) {
int c=a[down][j];
auto it1=pos[c].upper_bound(j);
auto it2=pos[c].lower_bound(j);
if(it1!=pos[c].begin()) {
--it1;
if(*it1<=j) R[*it1]=min(R[*it1], j-1);
}
if(it2!=pos[c].end()) {
if(j<=*it2) R[j]=min(R[j], *it2-1);
}
}
for(int i=m-1;i;--i) R[i]=min(R[i], R[i+1]);
rep(i,1,m+1) if(i<=R[i]) ans=max(ans, (down-up+1)*(R[i]-i+1));
rep(i,1,m+1) pos[a[down][i]].insert(i);
}
}
return ans;
}
int main() {
while(~scanf("%d%d",&n,&m)) {
///read
rep(i,1,n+1) rep(j,1,m+1) scanf("%d",&a[i][j]);
///solve
printf("%d\n",solve());
}
return 0;
}
\(O(n^3)\) 做法
#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define rep(i, a, b) for(int i=(a); i<(b); i++)
#define sz(x) (int)x.size()
#define de(x) cout<< #x<<" = "<<x<<endl
#define dd(x) cout<< #x<<" = "<<x<<" "
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
const int N=404;
int n, m;
int a[N][N], p[N*N], l[N][N], r[N][N], R[N];
void Min(int &a, int b) {
if(a>b) a=b;
}
void Max(int &a, int b) {
if(a<b) a=b;
}
int solve() {
int ans=1;
rep(i,0,n+1) rep(j,0,m+1) l[i][j]=0, r[i][j]=m+1;
for(int u=n;u;--u) {
rep(i,1,m+1) p[a[u][i]]=0;
rep(i,1,m+1) l[u][i]=p[a[u][i]], p[a[u][i]]=i;
rep(i,1,m+1) p[a[u][i]]=m+1;
for(int i=m;i;--i) r[u][i]=p[a[u][i]], p[a[u][i]]=i;
memset(p,0,sizeof(p));
rep(i,1,m+1) {
p[a[u][i]]=i;
rep(v,u+1,n+1) Max(l[v][i], p[a[v][i]]);
}
rep(i,0,N*N) p[i]=m+1;
for(int i=m;i;--i) {
p[a[u][i]]=i;
rep(v,u+1,n+1) Min(r[v][i], p[a[v][i]]);
}
rep(i,1,m+1) R[i]=m+1;
rep(v,u,n+1) {
rep(i,1,m+1) {
Min(R[l[v][i]], i);
Min(R[i], r[v][i]);
}
for(int i=m-1;i;--i) Min(R[i], R[i+1]);
rep(i,1,m+1) {
Max(ans, (v-u+1)*(R[i]-i));
}
}
}
return ans;
}
int main() {
while(~scanf("%d%d",&n,&m)) {
///read
rep(i,1,n+1) rep(j,1,m+1) scanf("%d",&a[i][j]);
///solve
printf("%d\n",solve());
}
return 0;
}
区间dp做法
#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define rep(i, a, b) for(int i=(a); i<(b); i++)
#define sz(x) (int)x.size()
#define de(x) cout<< #x<<" = "<<x<<endl
#define dd(x) cout<< #x<<" = "<<x<<" "
typedef long long ll;
typedef pair<int, int> pii;
typedef vector<int> vi;
const int N=404;
int n,m;
int a[N][N], f[N][N], p[N][160004];
int main() {
while(~scanf("%d%d",&n,&m)) {
rep(i,1,n+1) rep(j,1,m+1) scanf("%d",&a[i][j]);
memset(f,0,sizeof(f));
memset(p,0,sizeof(p));
int ans=1;
rep(i,1,n+1) {
rep(len,1,m+1) {
for(int l=1, r=len;r<=m;++l, ++r) {
if(l==r) {
f[l][r]=max(f[l][r], p[l][a[i][l]]);
} else {
f[l][r]=max(f[l][r], max(f[l][r-1], f[l+1][r]));
f[l][r]=max(f[l][r], max(p[l][a[i][r]], p[r][a[i][l]]));
if(a[i][l]==a[i][r]) f[l][r]=i;
}
ans=max(ans, (i-f[l][r])*(r-l+1));
}
}
rep(j,1,m+1) p[j][a[i][j]]=i;
}
printf("%d\n",ans);
}
return 0;
}
codeforces 407D Largest Submatrix 3的更多相关文章
- Largest Submatrix(动态规划)
Largest Submatrix Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others ...
- POJ-3494 Largest Submatrix of All 1’s (单调栈)
Largest Submatrix of All 1’s Time Limit: 5000MS Memory Limit: 131072K Total Submissions: 8551 Ac ...
- hdu 2870 Largest Submatrix(平面直方图的最大面积 变形)
Problem Description Now here is a matrix with letter 'a','b','c','w','x','y','z' and you can change ...
- Largest Submatrix of All 1’s
Given a m-by-n (0,1)-matrix, of all its submatrices of all 1’s which is the largest? By largest we m ...
- Largest Submatrix of All 1’s(思维+单调栈)
Given a m-by-n (0,1)-matrix, of all its submatrices of all 1's which is the largest? By largest we m ...
- POJ 3494 Largest Submatrix of All 1’s 单调队列||单调栈
POJ 3494 Largest Submatrix of All 1’s Description Given a m-by-n (0,1)-matrix, of all its submatrice ...
- POJ - 3494 Largest Submatrix of All 1’s 单调栈求最大子矩阵
Largest Submatrix of All 1’s Given a m-by-n (0,1)-matrix, of all its submatrices of all 1’s which is ...
- HDU 2870 Largest Submatrix (单调栈)
http://acm.hdu.edu.cn/showproblem.php? pid=2870 Largest Submatrix Time Limit: 2000/1000 MS (Java/Oth ...
- MINSUB - Largest Submatrix
MINSUB - Largest Submatrix no tags You are given an matrix M (consisting of nonnegative integers) a ...
随机推荐
- sgu-203 Hyperhuffman(哈夫曼编码)
Hyperhuffman You might have heard about Huffman encoding - that is the coding system that minimizes ...
- redis的数据类型(二)string类型
下面讲解value,value包括String.List.Set.Sorted Set.Hash 一.String类型 1.string类型 String是最基本的类型,而且Stirng类型是二 ...
- Linq lambda 匿名方法
课程6 委托.匿名方法.Lambda表达式.Linq查询表达式 上课日志1 一.委托的基本认识 提问:能不能把方法作为参数传递??? 也即是能不能声明一个能存放方法的变量呢——委托. 委托是一种数据类 ...
- ios开发 学习积累20161101
20161101 XML的声明 <?XML version="1.0" encoding="UTF-8" ?> XML文档必须有根元素 XML 对大 ...
- 记一次tomcat自动退出问题
问题 环境: centos/tomcat8/jdk1.8 最近遇到部署在服务器的tomcat总是过一段时间就自动结束进程 ; 通过监控tomcat 日志文件(tail -f ./logs/catali ...
- golang rpc 简单范例
RPC(Remote Procedure Call Protocol)--远程过程调用协议,它是一种通过网络从远程计算机程序上请求服务,而不需要了解底层网络技术的协议. 它的工作流程如下图: go ...
- 洛谷P4165 [SCOI2007]组队(排序 堆)
题意 题目链接 Sol 跟我一起大喊:n方过百万,暴力踩标算! 一个很显然的思路是枚举\(H, S\)的最小值算,复杂度\(O(n^3)\) 我们可以把式子整理一下,变成 \[A H_i + B S_ ...
- css 超出部分以省略号的形式显示
想要实现文字超出部分以省略号的形式显示首先需要给此元素设置一个宽度,然后添加以下属性 overflow: hidden;/*内容超出后隐藏*/ text-overflow: ellipsis;/*超出 ...
- 使用nodeJs安装Vue-cli (win10 使用管理员身份)
TIP:win10下安装,使用管理员身份进行,否则会有权限限制. 1,安装完成node,node有自带的npm,可以直接在cmd中,找到nodeJs安装的路径下,进行命令行全局安装vue-cli.(n ...
- 50+ Useful Docker Tools
As containers take root, dozens of tools have sprung up to support them. Check out your options for ...