bzoj 3996[TJOI2015]线性代数

把那个式子化简一下（这个矩阵不能用交换律2333）出来的大概就是∑∑ai*aj*bij-∑ai*ci貌似这就可以最小割了

从起点向不bij连边，流量bij，bij向ci，cj连边，流量inf，然后c[]向汇点连边，流量c[].

（连边写的太丑，勿喷233）

#include<bits/stdc++.h>
#define inf 0x7fffffff
#define LL long long
#define N 200005
using namespace std;
inline int ra()
{
    int x=0,f=1; char ch=getchar();
    while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
    while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
int S,T;
struct nodE{
    int to,next,v;
}e[N<<3];
int n,head[N<<1],cnt=1,h[N<<1];
int q[N<<2],ans;
void insert(int x, int y, int v)
{
    e[++cnt].to=y;
    e[cnt].next=head[x];
    e[cnt].v=v;
    head[x]=cnt;
}
bool bfs()
{
    int l=0,r=1; q[0]=S; 
    for (int i=0; i<=T; i++) h[i]=0; h[S]=1;
    while (l<r)
    {
        int x=q[l++];
        for (int i=head[x];i;i=e[i].next)
            if (e[i].v && h[e[i].to]==0)
            {
                h[e[i].to]=h[x]+1;
                q[r++]=e[i].to;
            }
    }
    return h[T];
}
int dfs(int x, int f)
{
    if (x==T) return f;
    int w,ww=0;
    for (int i=head[x];i;i=e[i].next)
        if (h[e[i].to]==h[x]+1 && e[i].v)
        {
            w=dfs(e[i].to,min(e[i].v,f-ww));
            e[i].v-=w; e[i^1].v+=w; ww+=w;
            if (ww==f) return f;
        }
    if (!ww) h[x]=0;
    return ww;
}
void dinic()
{
    while (bfs()) ans-=dfs(S,inf);
}
int main()
{
    n=ra(); S=0,T=n*n+n+2;
    for (int i=1; i<=n; i++)
        for (int j=1; j<=n; j++)
        {
            int x=ra();
            ans+=x;
            insert(S,(i-1)*n+j,x);
            insert((i-1)*n+j,S,0);
            insert((i-1)*n+j,n*n+i,inf);
            insert(n*n+i,(i-1)*n+j,0);
            insert((i-1)*n+j,n*n+j,inf);
            insert(n*n+j,(i-1)*n+j,0);
        }
    for (int i=1; i<=n; i++)
    {
        int x=ra();
        insert(n*n+i,T,x);
        insert(T,n*n+i,0);
    }
    dinic();
    cout<<ans;
    return 0;
}