Beta Round ＃9 (酱油杯noi考后欢乐赛)最大伤害

题目：http://www.contesthunter.org/contest/Beta%20Round%20%EF%BC%839%20%28%E9%85%B1%E6%B2%B9%E6%9D%AFnoi%E8%80%83%E5%90%8E%E6%AC%A2%E4%B9%90%E8%B5%9B%29/%E6%9C%80%E5%A4%A7%E4%BC%A4%E5%AE%B3

题解：发现我最小割还是不怎么会T_T

          而且此题题意有些不明。。。

          看到出题人说就是选出一大坨，然后能加的都加。卧槽，这不是最大权闭合子图吗？

         最大权闭合子图中原来有边的现在连边为inf，表示i与j必须同割，那我们连边为x，就代表i，j不同割需要扣除多少。

         然后求出最小割，用总权值减去最小割，好象是对的？然后就A掉了。。。

代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 2000+5

#define maxm 1000000

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define mod 1000000007

using namespace std;

inline int read()

{

    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=10*x+ch-'0';ch=getchar();}

    return x*f;

}
int  n,m,s,t,sum,maxflow,tot=1,a[maxn],b[maxn][maxn],head[maxn],cur[maxn],h[maxn],q[maxn];

struct edge{int go,next,v;}e[maxm];

void ins(int x,int y,int z){e[++tot].go=y;e[tot].v=z;e[tot].next=head[x];head[x]=tot;}

void insert(int x,int y,int z){ins(x,y,z);ins(y,x,0);}

bool bfs()

{

    for(int i=s;i<=t;i++)h[i]=-1;

    int l=0,r=1;q[1]=s;h[s]=0;

    while(l<r)

    {

        int x=q[++l];

        for(int i=head[x];i;i=e[i].next)

         if(e[i].v&&h[e[i].go]==-1)

         {

            h[e[i].go]=h[x]+1;q[++r]=e[i].go;

         }

    }

    return h[t]!=-1;

}

int dfs(int x,int f)

{

    if(x==t) return f;

    int tmp,used=0;

    for(int i=cur[x];i;i=e[i].next)

     if(e[i].v&&h[e[i].go]==h[x]+1)

    {

        tmp=dfs(e[i].go,min(e[i].v,f-used));

        e[i].v-=tmp;if(e[i].v)cur[x]=i;

        e[i^1].v+=tmp;used+=tmp;

        if(used==f)return f;       

    }

    if(!used) h[x]=-1;

    return used;

}

void dinic()

{

    maxflow=0;

    while(bfs())

    {

        for (int i=s;i<=t;i++)cur[i]=head[i];maxflow+=dfs(s,inf);

    }

}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);

    n=read();m=read();s=0;t=n+m+1;
    for1(i,n+m)a[i]=read()*2;
    for1(i,n)for1(j,n)b[i][j]=read();
    for1(i,n)for1(j,m)b[i][j+n]=b[j+n][i]=read();
    for1(i,m)for1(j,m)b[i+n][j+n]=read();
    for1(i,n+m)
    {
        int t1=0,t2=0;
        if(a[i]>0)t1+=a[i],sum+=a[i];else t2+=abs(a[i]);
        for1(j,n+m)if(b[i][j])
        {
            t1+=b[i][j];
            if(b[i][j]>0)insert(i,j,b[i][j]),sum+=b[i][j];
        }
        if(t1)insert(s,i,t1);
        if(t2)insert(i,t,t2);
    }
    dinic();
    printf("%d
",(sum-maxflow)/2);

    return 0;

}