【图论】【二分图匹配】[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");
    }
}
posted @ 2015-08-26 20:18  outer_form  阅读(111)  评论(0编辑  收藏  举报