题目描述

如果玩过俄罗斯方块,应该见过如下图形:

我们称它为一个 \(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;
}

最新文章

  1. JSP页面元素构成
  2. c# 如何中List&lt;object&gt;中去掉object对象中的重复列数据?
  3. 情人节那点事,Power BI告诉你
  4. ORB
  5. 【网络】ssl与ssh
  6. 瞬态电压抑制二极管(TVS)选用原则
  7. oracle_有关表分区_查询
  8. 罪恶的SEO优化
  9. h5 meta学习
  10. 如何将外部的obj模型导入OpenGL
  11. Python:Day44 Javascript
  12. 使用sort函数进行排序
  13. Boost filessystem...
  14. InfluxDB、grafana、collect部署监控(centos6.8)
  15. spark RDD 常见操作
  16. 关于在ASP.NET中使用JavaScript的建议
  17. bean的实例化
  18. 输入N组父子对,求父子对所组成的二叉树的高度----17年某公司的笔试题
  19. day 44 django 学习入门
  20. Android系统file_contexts二进制与文本转换工具

热门文章

  1. CF1430 D. String Deletion(div 2)
  2. 图的全部实现(邻接矩阵 邻接表 BFS DFS 最小生成树 最短路径等)
  3. windows下nginx中配置端口转发 ----本文来自转载
  4. java程序练习:输入数字转换成中文输出(金额相关)
  5. Scala-1-字符处理
  6. 如何对List集合中的对象进行按某个属性排序
  7. monkey及其的日志管理和分析
  8. C语言100题集合005-删除一维数组中所有相同的数,使之只剩一个
  9. C++ 基础 6:模板
  10. vi&amp;vim 基本使用方法