HDU 2819(SWAP)
链接:http://acm.hdu.edu.cn/showproblem.php?pid=2819
题意:给出一个N*N的0,1矩阵问是否能通过行列交换使得矩阵的对角线全部唯一,可以输出在最后的矩阵中是那两行(列)进行了交换。
思路:二分图的最大匹配。
如果通过交换某些行没有办法的到解的话,那么只交换列 或者 既交换行又交换列 那也没办法得到解。
这个可以用矩阵的秩来解释,所有的对角线都是1,所以也就是矩阵的秩就是N,所以秩小于N就无解。另外,根据矩阵的性质,任意交换矩阵的两行 或者 两列,矩阵的秩不变,也就保证了如果通过 只交换行 或 只交换列 无法得到解的话,那么其他交换形式也必然无解。
明白这些后,开始建立二分图:
左边为横坐标x,右边为纵坐标y,
G[i][j] == 1 表示就从X的i向Y的j引一条边,那么这条边的含义就可以解释为可以将Y的第j列(因为Y表示的是列的集合)移到第i列,使得a[i][i]变成1,这样就相当于是第i行第i列就变成了1,也就是说对角线多了一个1。
然后求这个二分图的最大匹配(目的是为了让每一列只与X中的某一行匹配)。
code:
#include<bits/stdc++.h>
using namespace std;
const int maxn = 150;
int link[maxn];
vector<vector<int> >G(maxn);
bool vis[maxn];
bool dfs(int u){
for(int i = 0; i < G[u].size(); ++i){
int v = G[u][i];
if(!vis[v]){
vis[v] = 1;
if(link[v] == -1 || dfs(link[v])){
link[v] = u;
return true;
}
}
}
return false;
}
int hungary(int n){
memset(link,-1,sizeof(link));
int res = 0;
for(int i = 1; i <= n; ++i){
memset(vis,0,sizeof(vis));
if(dfs(i)) ++res;
}
return res;
}
int L[maxn],R[maxn];
int main(){
int n;
while(cin >> n){
int cnt = 0,x;
for(int i = 0; i <= n; ++i) G[i].clear();
for(int i = 1; i <= n; ++i){
for(int j = 1; j <= n; ++j){
cin >> x;
if(x) G[i].push_back(j);
}
}
if(hungary(n) < n){
puts("-1");
continue ;
}
else{
for(int i = 1; i <= n; ++i){
if(i != link[i]){
for(int j = 1; j <= n; ++j){
if(i == link[j]){
L[cnt] = i, R[cnt] = j,++cnt;
swap(link[i],link[j]);
}
}
}
}
printf("%d\n",cnt);
for(int i = 0; i < cnt; ++i){
printf("C %d %d\n",L[i],R[i]);
}
}
}
return 0;
}