[ARC119D]Grid Repainting 3

Grid Repainting 3

题解

我们可以发现当我们对于一个点用它将它们这一行/列涂白时，我们肯定会尽可能先将在这一行/列的其它点用去涂它所在的列/行涂白，而这一列/行又会进行同样的操作，一直这样下去。
所以我们可以考虑将行/列看作一个点，将一个红点看作连接它所在的行与列的边，于是，我们可以得到一个图。
而对于这一个图的每一个连通块，我们都会有一个行与列不能被选择。

我们假设我们有$tot$个大于$1$的连通块，总共包括了$a$行与$b$列。
由于选择不同的行与列涂白对于答案是没有区别的，我们可以假设我们选择了$x(x\in [0,tot])$个行，$tot-x$个列涂白，那么我们总涂白的格子数为，
$$nm-(n-a+x)(m-b+tot-x)$$
这明显是一个下凹函数，所以最值只会在$0$或$tot$处取得，我们可以枚举它们两个中哪个更大，确定$x$的取值。

那么我们就选择$x$个连通块，选择其中的一个节点开始遍历，开始节点的行不能选择，剩下$tot-x$个连通块开始遍历节点的列不能被选择。
然后就可以通过$dfs$输出答案。

时间复杂度$O\left(hw\right)$。

源码

#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;
}

谢谢！！！