Looking for a challenge! 2 开坑…
网上好像只能找到 http://ruchiose.is-programmer.com/2014/6/18/main_amppz_2011_solution.48903.html 这个大佬写的一句话题解…
我太菜了,波兰题真的恐怖,我调这个东西调了三个半小时… 退役选手好像做这种稍微长一点的题就不行了
题解基本上是意译自 Looking for a challenge 2
遇到的非傻逼错误好像也就没滚动数组导致RE了… 我不知道我为啥写了这么久啊?
题意
给出一个n*m的子矩阵,求出其最大等差子矩阵的面积。
等差矩阵即每行每列都为等差序列之矩阵。
测试组数不超过10000,n,m <= 3000 元素值不超过1e9,每个测试文件大小不超过20MB
引导问题1:Air crashes
题意
给出一个数组a,对于每一个元素,找出其左侧最近的比它小的元素。
解法:following the arrows
在数组前添加一个哨兵节点,其值设为负无穷。用一个数组,称其为arrow,来维护每个元素左侧比它小的最近元素。
从左向右遍历数组a,若a_{i-1} < a_i 则a_{i-1} 即是a_i左侧最近的比它小的元素。否则,比较 arrow[a_{i-1}] 指向的值与 ai 的大小关系,不断循环,直到找到为止。
此解法的时间复杂度为线性,因为每一支箭都最多只被向前移动一次。
引导问题2:Plot
给出一个01矩阵A,找出其中面积最大的全1子矩阵之面积。
解法:O(nm)
建立一个辅助数组D,D[i,j] 表示(i,j)下方连续的1的个数,该数组可以方便地线性求得。
那么,对于每一个非零的D[i,j],我们在同一行里从(i,j)出发寻找左右两侧离其最近的D[i,k]使得D[i,k]<=D[i,j]。
计左侧点的标号为j’,右侧点的标号为j’’,那么 D[i,j]\times (j’’-j’-1) 即为一个符合条件的子矩阵面积。
对于矩阵中的每一个元素计算上式,最大的结果即为答案,时间复杂度为O(nm)
原问题解法
经过观察,我们可以将一个数列分解成若干个子序列,使得每一个子序列都是等差数列。显然,相邻的两个子序列有且仅有有一个公共元素。
于是我们可以在线性时间内处理出从宽度为1的最大等差子矩阵,相似地,我们也可以这样处理出长度为1,宽度为2和长度为2的最大等差子矩阵。
下面就只剩长度和宽度都大于3的情况待解决。对于每个3*3的子矩阵,我们判断它是否为等差矩阵,如果是,则标记该子矩阵的中心元素。经过进一步观察,我们发现,如果一个子矩阵是等差矩阵,当且仅当该矩阵的每一个3*3大小的子矩阵都为等差矩阵。
于是,利用Plot一题中的方法,求出标记过为3*3子矩阵中心的最大子矩阵,即能求出a矩阵的最大等差子矩阵。
总的时间复杂度为O(nm).
代码
#include <iostream>
#include <algorithm>
#include <climits>
const int maxn = 3000+5;
const int inf = INT_MIN;
int T;
int n, m;
int a[maxn][maxn], d[maxn][maxn];
int arrow_left[maxn], arrow_right[maxn];
int ans, now;
int diff, diff1, diff2;
int main(int argc, char *argv[]) {
std::ios::sync_with_stdio(false);
std::cin.tie(0);
std::cin >> T;
while (T--) {
std::cin >> n >> m;
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
std::cin >> a[i][j];
}
}
ans = 1;
if (m > 1) {
ans = std::max(ans, 2);
for (int i = 1; i <= n; i++) {
now = 2;
diff = a[i][2]-a[i][1];
for (int j = 3; j <= m; j++) {
if (a[i][j] - a[i][j-1] == diff) {
now++