【图论】【二分图匹配】[HDU2819]Swap
题目
分析:在对角线上全是1,逆向思考这个问题,将对角线上的1回到他们原来的地方后,发现这些1中没有任何两个在同一行或同一列,也就是这个图是一个置换矩形。
将每一行,每一列编号,分别作为x部和y部,将1所在的行列连边,求最大匹配。
至于输出方案,匹配完了之后,枚举每一行(或列,下同),将这一行的1置换到对角线上去即可。
代码:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 100
int n,c[MAXN+10],c2[MAXN+10],ans1[MAXN+10],ans2[MAXN+10],cnt;
bool vis[MAXN+10];
struct node{
int v;
node *next;
}edge[MAXN*MAXN+10],*adj[MAXN+10],*ecnt=edge;
void addedge(int u,int v){
node *p=++ecnt;
p->v=v;
p->next=adj[u];
adj[u]=p;
}
void Read(int &x){
char c;
while((c=getchar())&&c!=EOF)
if(c>='0'&&c<='9'){
x=c-'0';
while((c=getchar())&&c>='0'&&c<='9')
x=x*10+c-'0';
ungetc(c,stdin);
return;
}
exit(0);
}
void read(){
Read(n);
int i,j,t;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++){
Read(t);
if(t)
addedge(i,j);
}
}
bool dfs(int u){
int v;
for(node *p=adj[u];p;p=p->next){
v=p->v;
if(!vis[v]){
vis[v]=1;
if(!c[v]||dfs(c[v])){
c[v]=u;
c2[u]=v;
return 1;
}
}
}
return 0;
}
bool maxmatch(){
for(int i=1;i<=n;i++){
memset(vis,0,sizeof vis);
if(!dfs(i))
return 0;
}
return 1;
}
void print(){
int i;
for(i=1;i<=n;i++){
if(c[i]==i)
continue;
ans1[++cnt]=min(c[i],i),ans2[cnt]=max(c[i],i);
c2[c[i]]=c2[i];
c[c2[i]]=c[i];
c[i]=i;
c2[i]=i;
}
printf("%d\n",cnt);
for(i=1;i<=cnt;i++)
printf("R %d %d\n",ans1[i],ans2[i]);
}
int main()
{
while(1){
cnt=0;
memset(c,0,sizeof c);
memset(c2,0,sizeof c2);
memset(adj,0,sizeof adj);
ecnt=edge;
read();
if(maxmatch())
print();
else
puts("-1");
}
}