标题 [ZJOI2007]棋盘制作 [NOIP1999]拦截导弹
[ZJOI2007]棋盘制作
ZJOI2007
悬线法求解最大矩型(DP)
1.有效竖线:除了两个端点以外,不覆盖任何一个障碍点的竖直线段1.有效竖线:除了两个端点以外,不覆盖任何一个障碍点的竖直线段
2.悬线:上端覆盖了一个障碍点或者到达整个矩形上边界的有效线段2.悬线:上端覆盖了一个障碍点或者到达整个矩形上边界的有效线段
如果把一条悬线向左右两个方向尽可能的移动,那么就得到了一个矩形。如果把一条悬线向左右两个方向尽可能的移动,那么就得到了一个矩形。
注意:悬线对应的矩型不一定是极大子矩阵,悬线定义中固定了悬线的下边界注意:悬线对应的矩型不一定是极大子矩阵,悬线定义中固定了悬线的下边界
因为取悬线要确定下标,所以悬线左右移动所得到的矩形无法向下扩展因为取悬线要确定下标,所以悬线左右移动所得到的矩形无法向下扩展
但是可以向上扩展,这句话现在可以不用理解。但是可以向上扩展,这句话现在可以不用理解。
有如下定义:有如下定义:
我觉得最妙的做法就是通过坐标奇偶性取反,这里要用到单调栈。
即通过取反将01问题转换成求最大的全1正方形和矩阵问题即通过取反将01问题转换成求最大的全1正方形和矩阵问题
正解应该就是悬线法了吧正解应该就是悬线法了吧
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e3+7;
typedef long long ll;
inline ll read() {
ll s = 0, w = 1;
char ch = getchar();
while (ch < 48 || ch > 57) {
if (ch == '-') w = -1;
ch = getchar();
}
while (ch >= 48 && ch <= 57) s = (s << 1) + (s << 3) + (ch ^ 48), ch = getchar();
return s * w;
}
int res[maxn][maxn],lt[maxn][maxn],rt[maxn][maxn],height[maxn][maxn],ans,ans2;
int main() {
int n=read(),m=read();
for(int i=1; i<=n; ++i) {
for(int j=1; j<=m; ++j) {
res[i][j]=read();
lt[i][j]=rt[i][j]=j;
height[i][j]=1;
}
}
for(int i=1;i<=n;++i)
for(int j=2;j<=m;++j)
if(res[i][j]!=res[i][j-1]) lt[i][j]=lt[i][j-1];
for(int i=1;i<=n;++i)
for(int j=m-1;j;--j)
if(res[i][j]!=res[i][j+1]) rt[i][j]=rt[i][j+1];
for(int i=1;i<=n;++i)
for(int j=1,a,b;j<=m;++j) {
if(i>1&&res[i][j] != res[i - 1][j]) {
lt[i][j]=max(lt[i][j],lt[i - 1][j]);
rt[i][j]=min(rt[i][j],rt[i - 1][j]);
height[i][j] = height[i - 1][j] + 1;
}
a=rt[i][j]-lt[i][j]+1;
b=min(a,height[i][j]);
ans=max(b*b,ans);
ans2=max(ans2,a*height[i][j]);
}
printf("%d\n%d",ans,ans2);
return 0;
}
题目二拦截导弹
[NOIP1999]拦截导弹
Dilworth定理:
定理1 令(X,≤)是一个有限偏序集,并令r是其最大链的大小。则X可以被划分成r个但不能再少的反链。定理1令(X,≤)是一个有限偏序集,并令r是其最大链的大小。则X可以被划分成r个但不能再少的反链。
其对偶定理称为Dilworth定理:其对偶定理称为Dilworth定理:
定理2 令(X,≤)是一个有限偏序集,并令m是反链的最大的大小。则X可以被划分成m个但不能再少的链。定理2令(X,≤)是一个有限偏序集,并令m是反链的最大的大小。则X可以被划分成m个但不能再少的链。
即:最长递增子序列的长度等于最少不递增序列的个数,最长不递增序列的长度等于最少递增子序列的个数即:最长递增子序列的长度等于最少不递增序列的个数,最长不递增序列的长度等于最少递增子序列的个数
最后结合一下二分和一个单调的数组最后结合一下二分和一个单调的数组
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+7;
typedef long long ll;
int a[maxn],que[maxn];
int main() {
int n=0,len,i;
while(~scanf("%d",&a[++n]));
--n;
for(i=n,len=0;i;--i) {
if(a[i]>=que[len]) que[++len]=a[i];
else {
int cnt=upper_bound(que+1,que+1+len,a[i])-que;
que[cnt]=a[i];
}
}
printf("%d\n",len);
for(i=1,len=0;i<=n;++i) {
if(a[i]>que[len]) que[++len]=a[i];
else {
int cnt=lower_bound(que+1,que+1+len,a[i])-que;
que[cnt]=a[i];
}
}
printf("%d\n",len);
}
二分 Dilworth定理 STL DP