bzoj1059: [ZJOI2007]矩阵游戏

链接

　　http://www.lydsy.com/JudgeOnline/problem.php?id=1059

题解

　　最大流手生了….1A连击中止，连击数：2
　　这是道网络流神题，题目问你能否通过移动使矩阵的对角线上全是1。显然这个过程是可逆的，因此可以倒着考虑，考虑对角线上全是1的矩阵通过变换变出的矩阵应该满足什么条件。
　　显然，处于同一行的1怎么变换都还是在1行里；处于同一列的1怎么变换也都是在1列里。考虑目标状态中对角线上的1，显然那些1两两之间都互不同行、互不同列，而且怎么变换都是如此。
　　那问题就成了，看看能不能从矩阵里找出n个元素，两两之间不同行不同列。
　　问题变成：能否每一行选一个1，选出来1都不同列
　　这样就可以直接构建网络流模型了，一个坐标为(x,y)的1，对应左侧的第x个点到右侧的第y个点的一条边，在用S连向所有的左侧点，容量为1；右侧所有点连向T，容量为1。跑最大流看看等与不等于N。

代码

//网络流
#include <cstdio>
#include <algorithm>
#define maxn 1000
#define inf 0x3f3f3f3f
using namespace std;
int N, d[maxn], num[maxn], c[maxn*maxn], Exit, head[maxn], tot,
    nex[maxn*maxn], to[maxn*maxn], S, T, last[maxn];
void adde(int a, int b, int cc){to[++tot]=b;c[tot]=cc;nex[tot]=head[a];head[a]=tot;}
void adde2(int a, int b, int cc){adde(a,b,cc);adde(b,a,0);}
int isap(int pos, int in)
{
    if(pos==T)return in;
    int t, flow=0;
    for(int &p=last[pos];p;p=nex[p])
        if(c[p] and d[to[p]]+1==d[pos])
        {
            flow+= t=isap(to[p],min(in-flow,c[p]));
            c[p]-=t, c[p xor 1]+=t;
            if(Exit or flow==in)return flow;
        }
    num[d[pos]]--;
    Exit=num[d[pos]]==0;
    num[++d[pos]]++;
    last[pos]=head[pos];
    return flow;
}
bool work() 
{
    int i, j, x, flow=0;
    scanf("%d",&N);
    for(i=1;i<=T;i++)head[i]=last[i]=0;
    for(i=1;i<=tot;i++)nex[i]=0;
    tot=1;
    for(i=1;i<=N;i++)for(j=1;j<=N;j++)
    {
        scanf("%d",&x);
        if(x)adde2(i,N+j,1);
    }
    S=N+N+10, T=S+1;
    for(i=1;i<=N;i++)adde2(S,i,1),adde2(N+i,T,1);
    for(i=1;i<=T;i++)num[i]=d[i]=0;
    Exit=0;
    while(!Exit)flow+=isap(S,inf);
    return flow==N;
}
int main()
{
    int T;
    for(scanf("%d",&T);T;T--)
    {
        if(work())printf("Yes");else printf("No");
        putchar('\n');
    }
    return 0;
}