bzoj 1059: [ZJOI2007]矩阵游戏 最大流

Description

　　小Q是一个非常聪明的孩子，除了国际象棋，他还很喜欢玩一个电脑益智游戏——矩阵游戏。矩阵游戏在一个N

*N黑白方阵进行（如同国际象棋一般，只是颜色是随意的）。每次可以对该矩阵进行两种操作：行交换操作：选择

矩阵的任意两行，交换这两行（即交换对应格子的颜色）列交换操作：选择矩阵的任意行列，交换这两列（即交换

对应格子的颜色）游戏的目标，即通过若干次操作，使得方阵的主对角线(左上角到右下角的连线)上的格子均为黑

色。对于某些关卡，小Q百思不得其解，以致他开始怀疑这些关卡是不是根本就是无解的！！于是小Q决定写一个程

序来判断这些关卡是否有解。

Sample Input

2
2
0 0
0 1
3
0 0 1
0 1 0
1 0 0

【数据规模】
对于100%的数据，N ≤ 200

解：无论怎么操作，同行的永远同行，同列的永远同列--->所以同行同列的1，操作之后也是同行同列

题目转化为找到n个不同行的且均不同列的元素 or n个不同列且两两不同行的元素

 

考虑网络流  n行放左边  n列放右边

每行的流量是1，每列的流量是1

i行j列有点   add(i,j,1)  

最大流是n有解

(刚开始想到高斯消元，是不可行的 )

1.刚开始想求rank==n（出现两个全1向量就错了）

2.再想怎么修改消元过程（消元过程当中，没有调整的余地--->第一列有多个1选择的话，选了一个，那么这一列其他行的1直接舍弃了-->对后面的列进行选择时，发现不匹配无法向前调整了）-->又变成匹配问题了---->最大流就好了、

 

 

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#define mem(a,b) memset(a,b,sizeof a)
#define en '\n'
#define maxn 100005
using namespace std;
typedef long long ll;
template<class T>void rd(T &x)
{
    x=0;int f=0;char ch=getchar();
    while(ch<'0'||ch>'9')  {f|=(ch=='-');ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}
    x=f?-x:x;
    return;
}

const int inf =1<<29;
int head[maxn];
int n,m,tot,maxflow,s,t;
queue<int>q;
int d[maxn];
struct node{
int v,nxt,w;
}edge[2*maxn];

void add(int u,int v,int w){
    edge[++tot].nxt=head[u];head[u]=tot;edge[tot].v=v;edge[tot].w=w;
    edge[++tot].nxt=head[v];head[v]=tot;edge[tot].v=u;edge[tot].w=0;
}
bool bfs(){
mem(d,0);
    while(q.size())q.pop();
    q.push(s);
    d[s]=1;
    while(q.size()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=edge[i].nxt){
                if(edge[i].w and !d[edge[i].v]){
                    q.push(edge[i].v);
                    d[edge[i].v]=d[x]+1;
                    if(edge[i].v==t){
                        return 1;
                    }
                }
        }
    }
    return 0;
}
int dinic(int x,int flow){
    if(x==t)return flow;
    int rest=flow,k;
    for(int i=head[x];i and rest ;i=edge[i].nxt){
        if(edge[i].w and d[edge[i].v]==d[x]+1){
            k=dinic(edge[i].v,min(rest,edge[i].w));
            if(!k)d[edge[i].v]=0;
            edge[i].w-=k;
            edge[i^1].w+=k;
            rest-=k;
        }
    }
    return flow -rest;
}
signed main()
{
    #ifdef local
    freopen("input2.txt","r",stdin);
    #endif
    int T;
    cin>>T;
    while(T--){
        cin>>n;
        mem(head,0);
        tot=1;
        s=0;
        t=2*n+1;
        for(int i=1;i<=n;i++){
                add(s,i,1);
                add(i+n,t,1);
            for(int j=1;j<=n;j++){
                int tem;
                scanf("%d",&tem);
                if(tem)
                add(i,n+j,1);
            }
        }

    int flow=0;
    maxflow=0;
    while(bfs()){
        while(flow=dinic(s,inf))maxflow+=flow;
    }
    if(maxflow==n){
        puts("Yes");
    }else cout<<"No"<<en;
   }
    return 0;
}