模拟赛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;
}
模拟赛38 B. T形覆盖 大模拟的更多相关文章
- noip模拟赛#38
我打开了#39的problem...想了半个小时多发现我一道题都不会写...于是我打开了#38的problem T1:循环数字的定义为能够将该数划分为若干相同长度的段并且都相同. n=2e18. =& ...
- [luogu#2019/03/10模拟赛][LnOI2019]长脖子鹿省选模拟赛赛后总结
t1-快速多项式变换(FPT) 题解 看到这个\(f(x)=a_0+a_1x+a_2x^2+a_3x^3+ \cdots + a_nx^n\)式子,我们会想到我们学习进制转换中学到的,那么我们就只需要 ...
- [模拟赛FJOI Easy Round #2][T1 sign] (模拟+求字符串重复字串)
[题目描述] 小Z在无意中发现了一个神奇的OJ,这个OJ有一个神奇的功能:每日签到,并且会通过某种玄学的算法计算出今日的运势.在多次试验之后,小Z发现自己的运势按照一定的周期循环,现在他找到了你,请通 ...
- 模拟赛小结：2017-2018 ACM-ICPC Nordic Collegiate Programming Contest (NCPC 2017)
比赛链接:传送门 本场我们队过的题感觉算法都挺简单的,不知道为啥做的时候感觉没有很顺利. 封榜后7题,罚时1015.第一次模拟赛金,虽然是北欧的区域赛,但还是有点开心的. Problem B Best ...
- 【2020.8.23NOIP模拟赛】失落
[ 2020.8.23 N O I P 模 拟 赛 ] 失 落 [2020.8.23NOIP模拟赛]失落 [2020.8.23NOIP模拟赛]失落 题目描述 出题人心情很失落,于是他直接告诉你让你求出 ...
- HDU 5920 Ugly Problem 高精度减法大模拟 ---2016CCPC长春区域现场赛
题目链接 题意:给定一个很大的数,把他们分为数个回文数的和,分的个数不超过50个,输出个数并输出每个数,special judge. 题解:现场赛的时候很快想出来了思路,把这个数从中间分为两部分,当位 ...
- 2016ACM-ICPC网络赛北京赛区 1001 （trie树牌大模拟）
[题目传送门] 1383 : The Book List 时间限制:1000ms 单点时限:1000ms 内存限制:256MB 描述 The history of Peking University ...
- 体育成绩统计——20180801模拟赛T3
体育成绩统计 / Score 题目描述 正所谓“无体育,不清华”.为了更好地督促同学们进行体育锻炼,更加科学地对同学们进行评价,五道口体校的老师们在体育成绩的考核上可谓是煞费苦心.然而每到学期期末时, ...
- [GRYZ]寒假模拟赛
写在前面 这是首次广饶一中的OIERS自编自导,自出自做(zuo)的模拟赛. 鉴于水平气压比较低,机(wei)智(suo)的WMY/XYD/HYXZC就上网FQ下海找了不少水(fei)题,经过他们优( ...
随机推荐
- maven 配置文件 settings.xml pom.xml
maven涉及的配置文件一般有三个:(作用范围依次减小,优先级依次升高) 全局配置(maven安装包中):${M2_HOME}/conf/settings.xml 用户配置(该文件是拷贝全局配置而来的 ...
- 论文解读《Plug-and-Play Priors for Model Based Reconstruction》
这篇论文主要概述了model-baesd的方法在解决图像恢复的逆问题的很好的效果,降噪问题其实就是前向模型的H是一个恒等算子,将state-of-the-art的降噪算法(先验模型)和相对应的逆问题的 ...
- 三分钟带你分清Mysql 和Oracle之间的误区
摘要:Mysql 和Oracle,别再傻傻分不清. mysql 和Oracle 在开发中的使用是随处可见的,那就简单去了解一下这俩款火的不行的数据库. 本质区别: Oracle数据库是一个对象关系数据 ...
- Linux 网络编程的5种IO模型：异步IO模型
Linux 网络编程的5种IO模型:异步IO模型 资料已经整理好,但是还有未竟之业:复习多路复用epoll 阅读例程, 异步IO 函数实现 背景 上一讲< Linux 网络编程的5种IO模型:信 ...
- NOIP2018 填坑记
Oct,22ed,2018 DAY -18 又是颓废的一天呢 我好菜啊,一个圆方树弄了一整天(点双怎么那么毒瘤).(铁人两项怎么那么多点) Oct,23rd,2018 DAY -17 又双叒叕颓了一天 ...
- 如何将别人Google云端硬盘中的数据进行保存
查了好久终于知道! 如何将别人Google云端硬盘中的数据进行copy,而不是右键发现只有添加快捷方式 只要shift+z就可以保存了! 之后等我弄清楚怎么将别人家的云盘中的数据集导到colab再来详 ...
- leetcode110:combination-sum-ii
题目描述 给出一组候选数C和一个目标数T,找出候选数中起来和等于T的所有组合. C中的每个数字在一个组合中只能使用一次. 注意: 题目中所有的数字(包括目标数T)都是正整数 组合中的数字 (a 1, ...
- leetcode25word-ladder
题目描述 给定两个单词(初始单词和目标单词)和一个单词字典,请找出所有的从初始单词到目标单词的最短转换序列的长度: 每一次转换只能改变一个单词 每一个中间词都必须存在单词字典当中 例如: 给定的初始单 ...
- CSS实现模拟百度分享侧边栏效果
在<JS模拟百度分享侧边栏效果>一文中对于Div区块的运动通过JS实现了鼠标移入滑出显示,鼠标移出滑入隐藏的效果.其实在CSS3中通过transition属性就可以较为轻松实现. < ...
- martini-md参数（mdp文件）
输入参数:一个典型的mdp文件 1 ; 2 ; STANDARD MD INPUT OPTIONS FOR MARTINI 2.x 3 ; Updated 02 feb 2013 by DdJ 4 ; ...