hdu2819(二分构图)

来源:互联网 发布:java 文字转图片 编辑:程序博客网 时间:2024/06/05 05:55

       题意:给一个n*n的01矩阵,要求通过行变换或列变换使对角线上都是1,并给出变换步骤。若无解则输出-1。

       分析:就是一个简单的二分构图题,行和列匹配就行了,有解的话就是完备匹配。关键在于怎么给出步骤。如果把行数放在左边,列数放在右边,这样的匹配结果中match[i]就表示了第i行的1所在的位置,而为了使对角线上全是1,应有match[i]==i,所以用选择排序的方法对match[]进行排序,并记录下match的下标交换过程,然后输出就行了。 

                   还有一点就是,如果所给矩阵有解,则一定能只通过行变换或列变换求解,这想一想就能明白了。

代码如下:

#include<cstdio>#include<cstring>#include<iostream>#include<cstdlib>#include<cmath>#include<algorithm>#include<queue>#include<stack>#include<set>#include<map>using namespace std;const int M=1010;const int N=110;int g[N][N],vis[N],m[N][N];int n,ans,match[N];void build(){    memset(m,0,sizeof(m));    memset(g,0,sizeof(g));    for (int i=0;i<n;i++)        for (int j=0;j<n;j++)        {            scanf("%d",&g[i][j]);            if (g[i][j]){                m[i][j]=1;            }        }}bool dfs(int x){    int i;    for (i=0;i<n;i++) if (m[x][i] && !vis[i]) {        vis[i]=true;        if (match[i]==-1 || dfs(match[i])){            match[i]=x;            return true;        }    }    return false;}void print(){    int i,j,k=0,t=0;    pair<int,int> p[N];    for (i=0;i<n;i++) {        k=i;        for (j=i+1;j<n;j++) if (match[k]>match[j]){            k=j;        }        if (i!=k){            swap(match[i],match[k]);            p[t].first=match[i];            p[t].second=match[k];            t++;        }    }    printf("%d\n",t);    for (i=0;i<t;i++) printf("R %d %d\n",p[i].first+1,p[i].second+1);}int main(){    int i,j;    while(~scanf("%d",&n))    {        build();        ans=0;        memset(match,-1,sizeof(match));        for (i=0;i<n;i++) {            memset(vis,false,sizeof(vis));            if (dfs(i)) ans++;        }        if (ans<n){            printf("-1\n");        } else {            print();        }    }}/*41 0 1 00 1 0 00 1 0 00 0 0 140 0 1 01 0 0 00 1 0 00 0 0 130 1 01 0 00 0 1*/


0 0
原创粉丝点击