模拟赛38 B. T形覆盖 大模拟
题目描述
如果玩过俄罗斯方块,应该见过如下图形:
我们称它为一个 \(T\) 形四格拼板 。其中心被标记为\(×\)。
小苗画了一个 \(m\) 行 \(n\) 列的长方形网格。行从 \(0\) 至 \(m-1\) 编号,列从 \(0\) 至 \(n-1\) 编号。
她将网格中的一些格子标记为 特殊格子 。
然后,她想要她的朋友帮助她将 \(T\) 形四格拼板按找如下规则摆放:
\(1\)、 特殊格子的数量与 \(T\) 形四格拼板的数量相同,每个 \(T\) 形四格拼板的中心在网格上的位置必须是特殊格子。
\(2\)、 \(T\) 形四格拼板之间不能有重叠部分。
\(3\)、 所有拼板的部分均在网格内。
注意,\(T\) 形四格拼板有四种摆放方式:
如果方案不存在,输出 \(No\),否则请找出一种方案使得被拼板覆盖的数总和最大,求出这个最大值。
输入格式
从文件 \(t-covering.in\) 中读入数据。
第一行两个整数 \(m,n\) 分别表示行数和列数。
接下来 \(m\) 行,每行 \(n\) 个整数,第 \(i\) 行第 \(j\) 个数(从 \(0\) 开始编号)表示方格中第 \(i\) 行第 \(j\) 列的数 \(a_{i,j}\)
接下来一行,一个整数 \(k\) ,表示特殊格子的数量。
接下来 \(k\) 行,每行两个整数 \(r_i,c_i\),表示第 \(i\) 个被标记的特殊格子的位置。
输出格式
如果有方案,输出可能的被覆盖的格子内数总和的最大值,否则输出 \(No\)。
注意,本题 \(10\) 组测试数据均绑点, 只输出 \(1\) 个 \(No\) 不会得分.
样例
样例输入
5 6
7 3 8 1 0 9
4 6 2 5 8 3
1 9 7 3 9 5
2 6 8 4 5 7
3 8 2 7 3 6
3
1 1
2 2
3 4
样例输出
67
【样例 \(1\) 解释】
其中一种最优的方案如下:
\((1,1)\) 位置的摆放方式为 \(⊣\)。
\((2,2)\) 位置的摆放方式为\(⊢\)。
\((3,4)\) 位置的摆放方式为 \(⊥\)。
样例 \(2\) 输入
5 6
7 3 8 1 0 9
4 6 2 5 8 3
1 9 7 3 9 5
2 6 8 4 5 7
3 8 2 7 3 6
3
1 1
2 2
3 3
【样例 \(2\) 输出】
No
数据范围与提示
分析
把这道题打成了大模拟
思路是这样的:
\(1\)、把所有必须选的格子选上,这样的 \(T\) 有三种
第一种是在边界上,第二种是两个格子直接相邻,第三种是两个格子对角线相邻
在没有限制的情况下,每一个 \(T\) 可以在周围的四个格子中任意选三个
而给了限制之后就只能选固定的三个了
要注意在选择过程中我们选择了一个格子后会使其它 \(T\) 能选择的格子减少
也就是说只能选择固定格子的 \(T\) 增加了,这个要存一下
\(2\)、对于两个 \(T\) 相隔一个格子的情况用并查集维护
把所有相隔一个格子的 \(T\)合并,如果形成的联通块没有环,可以贪心去掉一个最小的
如果有一个环,必须全部选上
如果有两个环即以上,则不合法
\(3\)、剩下的贪心去掉最小的即可
代码
#include<cstdio>
#include<vector>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<iostream>
#define rg register
inline int read(){
rg int x=0,fh=1;
rg char ch=getchar();
while(ch<'0' || ch>'9'){
if(ch=='-') fh=-1;
ch=getchar();
}
while(ch>='0' && ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);
ch=getchar();
}
return x*fh;
}
const int maxn=1e6+5;
int n,m,k,cnt;
std::vector<int> a[maxn],id[maxn];
bool jud=0,vis[maxn],haha[maxn],inq[maxn],huan[maxn],s[maxn];
int tot[maxn],zxz[maxn],ans=0,js[maxn],sta[maxn],tp,jla[maxn],jlb[maxn];
struct asd{
int nx,ny;
asd(){}
asd(int aa,int bb){
nx=aa,ny=bb;
}
}b[maxn];
void xg(int nx,int ny){
if(nx<1 || ny<1 || nx>n || ny>m) return;
if(vis[id[nx][ny]]) return;
if(s[id[nx][ny]]==0) return;
js[id[nx][ny]]--;
if(!inq[id[nx][ny]]){
sta[++tp]=id[nx][ny];
inq[id[nx][ny]]=1;
}
}
int fa[maxn];
int zhao(int xx){
if(fa[xx]==xx) return xx;
return fa[xx]=zhao(fa[xx]);
}
int main(){
freopen("t_covering.in","r",stdin);
freopen("t_covering.out","w",stdout);
memset(zxz,0x3f,sizeof(zxz));
n=read(),m=read();
rg int aa,bb;
for(rg int i=1;i<=n;i++){
a[i].push_back(0);
id[i].push_back(0);
for(rg int j=1;j<=m;j++){
aa=read();
a[i].push_back(aa);
cnt++;
id[i].push_back(cnt);
}
}
k=read();
if(k*4>n*m){
printf("No\n");
for(rg int i=1;i<=k;i++){
aa=read(),bb=read();
}
return 0;
}
for(rg int i=1;i<=k;i++){
aa=read(),bb=read();
aa++,bb++;
if(s[id[aa][bb]]==1){
jud=1;
}
b[i].nx=aa,b[i].ny=bb;
jla[id[aa][bb]]=aa,jlb[id[aa][bb]]=bb;
s[id[aa][bb]]=1;
}
if(k){
if(n==1 || m==1) jud=1;
}
if(jud){
printf("No\n");
return 0;
}
for(rg int i=1;i<=k;i++){
js[id[b[i].nx][b[i].ny]]=4;
}
for(rg int i=1;i<=k;i++){
aa=b[i].nx,bb=b[i].ny;
if((aa==1 || aa==n) && (bb==1 || bb==m)){
printf("No\n");
return 0;
}
if(aa==1){
vis[id[aa][bb]]=1;
if(haha[id[aa][bb]] || haha[id[aa][bb-1]] || haha[id[aa][bb+1]] || haha[id[aa+1][bb]]){
printf("No\n");
return 0;
}
xg(aa,bb+2);
xg(aa,bb-2);
xg(aa+2,bb);
haha[id[aa][bb]]=haha[id[aa][bb-1]]=haha[id[aa][bb+1]]=haha[id[aa+1][bb]]=1;
ans+=a[aa][bb-1]+a[aa][bb+1]+a[aa+1][bb]+a[aa][bb];
} else if(aa==n){
vis[id[aa][bb]]=1;
if(haha[id[aa][bb]] || haha[id[aa][bb-1]] || haha[id[aa][bb+1]] || haha[id[aa-1][bb]]){
printf("No\n");
return 0;
}
xg(aa,bb+2);
xg(aa,bb-2);
xg(aa-2,bb);
haha[id[aa][bb]]=haha[id[aa][bb-1]]=haha[id[aa][bb+1]]=haha[id[aa-1][bb]]=1;
ans+=a[aa][bb-1]+a[aa][bb+1]+a[aa-1][bb]+a[aa][bb];
} else if(bb==1){
vis[id[aa][bb]]=1;
if(haha[id[aa][bb]] || haha[id[aa-1][bb]] || haha[id[aa+1][bb]] || haha[id[aa][bb+1]]){
printf("No\n");
return 0;
}
xg(aa-2,bb);
xg(aa+2,bb);
xg(aa,bb+2);
haha[id[aa][bb]]=haha[id[aa-1][bb]]=haha[id[aa+1][bb]]=haha[id[aa][bb+1]]=1;
ans+=a[aa-1][bb]+a[aa+1][bb]+a[aa][bb+1]+a[aa][bb];
} else if(bb==m){
vis[id[aa][bb]]=1;
if(haha[id[aa][bb]] || haha[id[aa-1][bb]] || haha[id[aa+1][bb]] || haha[id[aa][bb-1]]){
printf("No\n");
return 0;
}
xg(aa-2,bb);
xg(aa+2,bb);
xg(aa,bb-2);
haha[id[aa][bb]]=haha[id[aa-1][bb]]=haha[id[aa+1][bb]]=haha[id[aa][bb-1]]=1;
ans+=a[aa-1][bb]+a[aa+1][bb]+a[aa][bb-1]+a[aa][bb];
}
}
rg int mx,my;
for(rg int i=1;i<=k;i++){
aa=b[i].nx,bb=b[i].ny;
if(aa+1<=n){
mx=aa+1,my=bb;
if(s[id[mx][my]]){
if(!vis[id[mx][my]]){
vis[id[mx][my]]=1;
if(haha[id[mx][my+1]] || haha[id[mx][my-1]] || haha[id[mx+1][my]] || haha[id[mx][my]]){
printf("No\n");
return 0;
}
xg(mx,my+2);
xg(mx,my-2);
xg(mx+2,my);
haha[id[mx][my+1]]=haha[id[mx][my-1]]=haha[id[mx+1][my]]=haha[id[mx][my]]=1;
ans+=a[mx][my+1]+a[mx][my-1]+a[mx+1][my]+a[mx][my];
}
if(!vis[id[aa][bb]]){
vis[id[aa][bb]]=1;
if(haha[id[aa][bb-1]] || haha[id[aa][bb+1]] || haha[id[aa-1][bb]] || haha[id[aa][bb]]){
printf("No\n");
return 0;
}
xg(aa,bb-2);
xg(aa,bb+2);
xg(aa-2,bb);
haha[id[aa][bb-1]]=haha[id[aa][bb+1]]=haha[id[aa-1][bb]]=haha[id[aa][bb]]=1;
ans+=a[aa][bb-1]+a[aa][bb+1]+a[aa-1][bb]+a[aa][bb];
vis[id[aa][bb]]=1;
}
}
}
if(bb+1<=m){
mx=aa,my=bb+1;
if(s[id[mx][my]]){
if(!vis[id[mx][my]]){
vis[id[mx][my]]=1;
if(haha[id[mx+1][my]] || haha[id[mx-1][my]] || haha[id[mx][my+1]] || haha[id[mx][my]]){
printf("No\n");
return 0;
}
xg(mx-2,my);
xg(mx+2,my);
xg(mx,my+2);
haha[id[mx+1][my]]=haha[id[mx-1][my]]=haha[id[mx][my+1]]=haha[id[mx][my]]=1;
ans+=a[mx+1][my]+a[mx-1][my]+a[mx][my+1]+a[mx][my];
}
if(!vis[id[aa][bb]]){
vis[id[aa][bb]]=1;
if(haha[id[aa-1][bb]] || haha[id[aa+1][bb]] || haha[id[aa][bb-1]] || haha[id[aa][bb]]){
printf("No\n");
return 0;
}
xg(aa-2,bb);
xg(aa+2,bb);
xg(aa,bb-2);
haha[id[aa-1][bb]]=haha[id[aa+1][bb]]=haha[id[aa][bb-1]]=haha[id[aa][bb]]=1;
ans+=a[aa-1][bb]+a[aa+1][bb]+a[aa][bb-1]+a[aa][bb];
vis[id[aa][bb]]=1;
}
}
}
}
for(rg int i=1;i<=k;i++){
aa=b[i].nx,bb=b[i].ny;
if(aa+1<=n && bb+1<=m){
mx=aa+1,my=bb+1;
if(s[id[mx][my]]){
if(!vis[id[mx][my]]){
vis[id[mx][my]]=1;
if(haha[id[mx+1][my]] || haha[id[mx-1][my]] || haha[id[mx][my+1]] || haha[id[mx][my]]){
printf("No\n");
return 0;
}
xg(mx-2,my);
xg(mx+2,my);
xg(mx,my+2);
haha[id[mx+1][my]]=haha[id[mx-1][my]]=haha[id[mx][my+1]]=haha[id[mx][my]]=1;
ans+=a[mx+1][my]+a[mx-1][my]+a[mx][my+1]+a[mx][my];
}
if(!vis[id[aa][bb]]){
vis[id[aa][bb]]=1;
if(haha[id[aa+1][bb]] || haha[id[aa-1][bb]] || haha[id[aa][bb-1]] || haha[id[aa][bb]]){
printf("No\n");
return 0;
}
xg(aa-2,bb);
xg(aa+2,bb);
xg(aa,bb-2);
haha[id[aa+1][bb]]=haha[id[aa-1][bb]]=haha[id[aa][bb-1]]=haha[id[aa][bb]]=1;
ans+=a[aa+1][bb]+a[aa-1][bb]+a[aa][bb-1]+a[aa][bb];
vis[id[aa][bb]]=1;
}
}
}
if(bb-1>=1 && aa+1<=n){
mx=aa+1,my=bb-1;
if(s[id[mx][my]]){
if(!vis[id[mx][my]]){
vis[id[mx][my]]=1;
if(haha[id[mx+1][my]] || haha[id[mx-1][my]] || haha[id[mx][my-1]] || haha[id[mx][my]]){
printf("No\n");
return 0;
}
xg(mx-2,my);
xg(mx+2,my);
xg(mx,my-2);
haha[id[mx+1][my]]=haha[id[mx-1][my]]=haha[id[mx][my-1]]=haha[id[mx][my]]=1;
ans+=a[mx+1][my]+a[mx-1][my]+a[mx][my-1]+a[mx][my];
}
if(!vis[id[aa][bb]]){
vis[id[aa][bb]]=1;
if(haha[id[aa+1][bb]] || haha[id[aa-1][bb]] || haha[id[aa][bb+1]] || haha[id[aa][bb]]){
printf("No\n");
return 0;
}
xg(aa-2,bb);
xg(aa+2,bb);
xg(aa,bb+2);
haha[id[aa+1][bb]]=haha[id[aa-1][bb]]=haha[id[aa][bb+1]]=haha[id[aa][bb]]=1;
ans+=a[aa+1][bb]+a[aa-1][bb]+a[aa][bb+1]+a[aa][bb];
vis[id[aa][bb]]=1;
}
}
}
}
for(rg int i=1;i<=tp;i++){
if(vis[sta[i]]) continue;
if(js[sta[i]]<=2){
printf("No\n");
return 0;
} else {
mx=jla[sta[i]],my=jlb[sta[i]];
ans+=a[mx][my+1]+a[mx][my-1]+a[mx+1][my]+a[mx-1][my]+a[mx][my];
vis[sta[i]]=1;
if(haha[id[mx][my+1]]){
ans-=a[mx][my+1];
xg(mx-2,my);
xg(mx+2,my);
xg(mx,my-2);
} else if(haha[id[mx][my-1]]){
ans-=a[mx][my-1];
xg(mx-2,my);
xg(mx+2,my);
xg(mx,my+2);
} else if(haha[id[mx-1][my]]){
ans-=a[mx-1][my];
xg(mx,my-2);
xg(mx,my+2);
xg(mx+2,my);
} else {
ans-=a[mx+1][my];
xg(mx,my-2);
xg(mx,my+2);
xg(mx-2,my);
}
}
}
for(rg int i=1;i<=cnt;i++){
fa[i]=i;
}
rg int cc,dd;
for(rg int i=1;i<=k;i++){
if(vis[id[b[i].nx][b[i].ny]]) continue;
aa=b[i].nx,bb=b[i].ny;
cc=id[b[i].nx][b[i].ny];
tot[cc]+=a[aa][bb]+a[aa-1][bb]+a[aa+1][bb]+a[aa][bb-1]+a[aa][bb+1];
zxz[cc]=std::min(std::min(a[aa+1][bb],a[aa-1][bb]),std::min(a[aa][bb-1],a[aa][bb+1]));
}
for(rg int i=1;i<=k;i++){
if(vis[id[b[i].nx][b[i].ny]]) continue;
aa=b[i].nx,bb=b[i].ny;
mx=b[i].nx+2,my=b[i].ny;
if(mx<=n && s[id[mx][my]]){
cc=zhao(id[aa][bb]),dd=zhao(id[mx][my]);
if(cc==dd){
if(huan[cc]){
printf("No\n");
return 0;
}
huan[cc]=1;
tot[cc]-=a[aa+1][bb];
zxz[cc]=0;
} else {
tot[cc]+=tot[dd];
tot[cc]-=a[aa+1][bb];
zxz[cc]=std::min(zxz[cc],zxz[dd]);
fa[dd]=cc;
}
}
mx=b[i].nx,my=b[i].ny+2;
if(my<=m && s[id[mx][my]]){
cc=zhao(id[aa][bb]),dd=zhao(id[mx][my]);
if(cc==dd){
if(huan[cc]){
printf("No\n");
return 0;
}
huan[cc]=1;
tot[cc]-=a[aa][bb+1];
zxz[cc]=0;
} else {
tot[cc]+=tot[dd];
tot[cc]-=a[aa][bb+1];
zxz[cc]=std::min(zxz[cc],zxz[dd]);
fa[dd]=cc;
}
}
}
for(rg int i=1;i<=k;i++){
if(vis[id[b[i].nx][b[i].ny]]) continue;
mx=b[i].nx,my=b[i].ny;
if(fa[id[mx][my]]!=id[mx][my]) continue;
ans+=tot[id[mx][my]]-zxz[id[mx][my]];
}
printf("%d\n",ans);
return 0;
}