hdu2819 Swap(二分图匹配+输出交换方案)

题目

N*N(N<=100)的01矩阵，每次交换两行或两列

问是否能将矩阵交换为主对角线全1的矩阵，

如果可以，输出交换方案，每一步都会影响后面

相当于一步一步把原矩阵换成目的矩阵

思路来源

https://www.cnblogs.com/gj-Acit/archive/2013/08/17/3265502.html

题解

根据满秩矩阵及三秩相等，

若交换列可行，则交换行列一定可行；

反之交换列不行，交换行列一定不行，

二分图左顶点代表行，右顶点代表列，

(x,y)为1的话就连一条边，匹配数小于n一定不行，否则一定可以

最后交换方案的时候，枚举行，link[j]=i代表(i,j)为1

交换i列和j列，使得(i,i)为1，每次使得link[i]=i，这样n次即可

心得

看很多人建图直接建1到n的，左顶点和右顶点共用一个点

这样的图可能有自环，但似乎好像没事……

代码

#include<iostream>
#include<cstdio> 
#include<cstring>
using namespace std;
const int N=105;
bool vis[N*2];
int head[N*2],cnt,n,ans;
int link[N*2]; 
int a[N],b[N],tot,v;
struct edge{int to,next,w;}e[N*N*2];
void init()
{
	memset(head,0,sizeof head);
	cnt=0;
	tot=0;
}
void add(int u,int v)
{
	e[++cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}
bool dfs(int u)
{
	for(int i=head[u];i;i=e[i].next)
	{
		int v=e[i].to;
		if(!vis[v])
		{
			vis[v]=1;
			if(link[v]==-1||dfs(link[v]))
			{
				link[v]=u;
				return 1;
			}
		}
	}
	return 0;
}
int hungary(int n)
{
	int res=0;
	memset(vis,0,sizeof vis);
	memset(link,-1,sizeof link);
	for(int u=1;u<=n;++u)
	{
		memset(vis,0,sizeof vis);
		res+=dfs(u);
	}
	return res;
}
int main()
{
	while(~scanf("%d",&n))
	{
		init();
		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=n;++j)
			{
				scanf("%d",&v);
				if(v)add(i,j);
			}
		}
		ans=hungary(n);
		if(ans<n){puts("-1");continue;}
		int j;
		for(int i=1;i<=n;++i)//枚举行 
		{
			for(j=1;j<=n;++j)
			if(link[j]==i)break;
			if(j!=i)//交换第i列和第j列 
			{
				a[++tot]=i;
				b[tot]=j;
				swap(link[i],link[j]);
			}
		}
		printf("%d\n",tot);
		for(int i=1;i<=tot;++i)
		printf("C %d %d\n",a[i],b[i]);
	} 
	return 0;
}