题意:给一个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();
}
}
}
/*
4
1 0 1 0
0 1 0 0
0 1 0 0
0 0 0 1
4
0 0 1 0
1 0 0 0
0 1 0 0
0 0 0 1
3
0 1 0
1 0 0
0 0 1
*/