Grid Repainting 3
题解
我们可以发现当我们对于一个点用它将它们这一行/列涂白时,我们肯定会尽可能先将在这一行/列的其它点用去涂它所在的列/行涂白,而这一列/行又会进行同样的操作,一直这样下去。
所以我们可以考虑将行/列看作一个点,将一个红点看作连接它所在的行与列的边,于是,我们可以得到一个图。
而对于这一个图的每一个连通块,我们都会有一个行与列不能被选择。
我们假设我们有
t
o
t
tot
tot个大于
1
1
1的连通块,总共包括了
a
a
a行与
b
b
b列。
由于选择不同的行与列涂白对于答案是没有区别的,我们可以假设我们选择了
x
(
x
∈
[
0
,
t
o
t
]
)
x(x\in[0,tot])
x(x∈[0,tot])个行,
t
o
t
−
x
tot-x
tot−x个列涂白,那么我们总涂白的格子数为,
n
m
−
(
n
−
a
+
x
)
(
m
−
b
+
t
o
t
−
x
)
nm-(n-a+x)(m-b+tot-x)
nm−(n−a+x)(m−b+tot−x)
这明显是一个下凹函数,所以最值只会在
0
0
0或
t
o
t
tot
tot处取得,我们可以枚举它们两个中哪个更大,确定
x
x
x的取值。
那么我们就选择
x
x
x个连通块,选择其中的一个节点开始遍历,开始节点的行不能选择,剩下
t
o
t
−
x
tot-x
tot−x个连通块开始遍历节点的列不能被选择。
然后就可以通过
d
f
s
dfs
dfs输出答案。
时间复杂度 O ( h w ) O\left(hw\right) O(hw)。
源码
#include<bits/stdc++.h>
using namespace std;
#define MAXN 2505
#define lowbit(x) (x&-x)
#define reg register
#define mkpr make_pair
#define fir first
#define sec second
typedef long long LL;
typedef unsigned long long uLL;
const int INF=0x7f7f7f7f;
const LL jzm=2333;
const double Pi=acos(-1.0);
typedef pair<int,int> pii;
const double PI=acos(-1.0);
template<typename _T>
_T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
void read(_T &x){
_T f=1;x=0;char s=getchar();
while(s>'9'||s<'0'){if(s=='-')f=-1;s=getchar();}
while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
x*=f;
}
template<typename _T>
void print(_T x){if(x<0){x=(~x)+1;putchar('-');}if(x>9)print(x/10);putchar(x%10+'0');}
int h,w,tot;bool flag,fb[MAXN<<1];
bool vis[MAXN<<1];
vector<int>P[MAXN<<1];
vector<pii>ans;
char maze[MAXN][MAXN];
void dosaka(int x,int y){
bool rx=vis[x],cy=vis[y+h];vis[x]=vis[y+h]=1;
if(!rx)for(int i=0;i<(int)P[x].size();i++)dosaka(x,P[x][i]);
if(!cy)for(int i=0;i<(int)P[y+h].size();i++)dosaka(P[y+h][i],y);
}
void print(int x,int y){
bool rx=vis[x],cy=vis[y+h];vis[x]=vis[y+h]=1;
if(!rx)for(int i=0;i<(int)P[x].size();i++)print(x,P[x][i]);
if(!cy)for(int i=0;i<(int)P[y+h].size();i++)print(P[y+h][i],y);
if(!rx&&!fb[x]&&(cy||fb[y+h]))printf("X %d %d\n",x,y);
if(!cy&&!fb[y+h]&&(rx||fb[x]))printf("Y %d %d\n",x,y);
}
signed main(){
read(h);read(w);
for(int i=1;i<=h;i++){
scanf("\n%s",maze[i]+1);
for(int j=1;j<=w;j++)if(maze[i][j]=='R')
P[i].push_back(j),P[j+h].push_back(i);
}
for(int i=1;i<=h;i++)
for(int j=1;j<=w;j++)
if(maze[i][j]=='R'&&!vis[i]&&!vis[j+h])
dosaka(i,j),ans.push_back(mkpr(i,j)),tot++;
int sumrow=0,sumcol=0,bst=0,maxx=0;
for(int i=1;i<=h;i++)sumrow+=vis[i];
for(int i=h+1;i<=h+w;i++)sumcol+=vis[i];
for(int i=0;i<=tot;i++)
if(h*w-(h-sumrow+i)*(w-sumcol+tot-i)>maxx)
bst=i,maxx=h*w-(h-sumrow+i)*(w-sumcol+tot-i);
for(int i=1;i<=h+w;i++)vis[i]=0;
printf("%d\n",sumcol+sumrow-tot);
for(int i=0;i<tot;i++){
if(i<bst)fb[ans[i].fir]=1;else fb[ans[i].sec+h]=1;
print(ans[i].fir,ans[i].sec);
}
return 0;
}