题意:
给你一个n阶矩阵,问你能不能通过变换行列式来达到主对角行全为1的情况。可以输出要变换的次数和具体的变换情况,不可以的话就输出-1。
题解:
可以用匈牙利算法算出能不能满足对角线全为1的情况,即是匹配数是否为n即可,然后难点就在于怎么变换了。现在再来说一下匈牙利匹配的时候得出的结果数组吧,这个数组在这里的含义相当于第几列与第几列进行交换的意思,那么我们就可以在匹配完之后找到满足cy[j]==i的时候就进行交换即可了,下面请看代码
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int MAXN=100+7;
int map[MAXN][MAXN];
bool vis[MAXN];
int cx[MAXN],cy[MAXN];
int a[MAXN],b[MAXN];
int n;
int dfs(int u)
{
for(int v=1;v<=n;v++)
{
if(map[u][v]&&!vis[v])
{
vis[v]=true;
if(cy[v]==-1||dfs(cy[v]))
{
cx[u]=v;
cy[v]=u;
return 1;
}
}
}
return 0;
}
int main()
{
while(~scanf("%d",&n))
{
memset(map,0,sizeof(map));
memset(cx,-1,sizeof(cx));
memset(cy,-1,sizeof(cy));
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
scanf("%d",&map[i][j]);
}
}
int sum=0;
for(int i=1;i<=n;i++)
{
memset(vis,false,sizeof(vis));
sum+=dfs(i);
}
if(sum<n)
{
printf("-1\n");
continue;
}
int cnt=0;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
for(int i=1;i<=n;i++)
{
int k=0;
for(int j=1;j<=n;j++)
if(cy[j]==i)
{
k=j;
break;
}
if(i!=k)
{
// printf("%d %d\n",i,k);
a[cnt]=i;
b[cnt]=k;
cnt++;
swap(cy[i],cy[k]);
}
}
printf("%d\n",cnt);
for(int i=0;i<cnt;i++)
printf("C %d %d\n",a[i],b[i]);
}
}