3996: [TJOI2015]线性代数 思路题 最小割（最大权闭合图）

感觉是一道很巧妙的题。
开始没看到A是01矩阵哭瞎qwq。
首先先化简一下这个式子，手推一下还是比较好推的，由于我不会用公式编辑器所以就从网上贴一个啦，感谢TonyFang。

由于a是一个01矩阵，我们把ai看做“选”或“不选”，那么bi,j就可以看做同时选i和j两个物品的收益，ci可以看做选第i个物品的代价。
也就是说选bi,j必须选ci,cj，那就是最大权闭合图啦。

#include<iostream>
#include<cstdio>
#include<cstring>
#define inf 1000000007
#define N 251000
using namespace std;
int n,id,T,cnt=1,ans;
int head[N],q[N],dis[N],cur[N];
int list[N<<3],next[N<<3],key[N<<3];
inline int read()
{
    int a=0,f=1; char c=getchar();
    while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
    while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
    return a*f;
}
inline void insert(int x,int y,int z)
{
    next[++cnt]=head[x];
    head[x]=cnt;
    list[cnt]=y;
    key[cnt]=z;
}
inline bool BFS()
{
    memset(dis,-1,sizeof(dis));
    dis[0]=1; q[1]=0;
    int t=0,w=1,x;
    while (t<w)
    {
        int x=q[++t];
        for (int i=head[x];i;i=next[i])
            if (key[i]&&dis[list[i]]==-1)
                dis[list[i]]=dis[x]+1,q[++w]=list[i];
    }
    return dis[T]!=-1;
}
int find(int x,int flow)
{
    if (x==T) return flow;
    int w,used=0;
    for (int i=cur[x];i;i=next[i])
        if (key[i]&&dis[list[i]]==dis[x]+1)
        {
            w=find(list[i],min(key[i],flow-used));
            key[i]-=w; key[i^1]+=w; used+=w;
            if (key[i]) cur[x]=i;
            if (flow==used) return flow;
        }
    if (!used) dis[x]=-1;
    return used;
}
inline void dinic()
{
    while (BFS())
    {
        for (int i=0;i<=T;i++) cur[i]=head[i];
        ans-=find(0,inf);
    }
}
int main()
{
    n=read(); T=n*n+n+1;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
        {
            int x=read(); ++id; ans+=x;
            insert(0,id,x); insert(id,0,0);
            insert(id,i+n*n,inf); insert(i+n*n,id,0);
            insert(id,j+n*n,inf); insert(j+n*n,id,0);
        }
    for (int i=1;i<=n;i++)
    {
        int x=read();
        insert(i+n*n,T,x); insert(T,i+n*n,0);
    }
    dinic();
    cout << ans << endl;
    return 0;
}