bzoj 3996 最小割

公式推出来后想了半天没思路，居然A是01矩阵。。。。。

如果一个问题是求最值，并那么尝试先将所有可能收益加起来，然后矛盾部分能否用最小割表达（本题有两个矛盾，第一个是选还是不选，第二个是i,j有一个不选，就不能获得bij的收益）。

 #include <cstdio>
 #include <cstring>
 #include <vector>
 #define N 510
 #define S N+N*N
 #define E S*10
 #define oo 0x3f3f3f3f
 #define fill(a,l,r,v) memset(a+l,v,sizeof(a[0])*(r-l+1))
 using namespace std;

 struct Edge {
     int u, v, f;
     Edge(){}
     Edge( int u, int v, int f ):u(u),v(v),f(f){}
 };
 struct Dinic {
     int n, src, dst;
     int head[S], next[E], etot;
     Edge edge[E];
     int dep[S], cur[S], qu[S], bg, ed;

     void init( int n, int src, int dst ) {
         this->n = n;
         this->src = src;
         this->dst = dst;
         fill( head, , n, 0x3f );
     }
     void adde( int u, int v, int f ) {
         next[etot]=head[u], head[u]=etot;
         edge[etot++] = Edge(u,v,f);
         next[etot]=head[v], head[v]=etot;
         edge[etot++] = Edge(v,u,);
     }
     bool bfs() {
         fill(dep,,n,);
         qu[bg=ed=] = src;
         dep[src] = ;
         while( bg<=ed ) {
             int u=qu[bg++];
             for( int t=head[u]; t!=oo; t=next[t] ) {
                 Edge &e = edge[t];
                 if( e.f && !dep[e.v] ) {
                     qu[++ed] = e.v;
                     dep[e.v] = dep[e.u]+;
                 }
             }
         }
         return dep[dst];
     }
     int dfs( int u, int a ) {
         if( u==dst || a== ) return a;
         int remain=a, past=, na;
         if( cur[u]==- ) cur[u]=head[u];
         for( int &t=cur[u]; t!=oo; t=next[t] ) {
             Edge &e = edge[t];
             Edge &ve = edge[t^];
             if( e.f && dep[e.v]==dep[e.u]+ && (na=dfs(e.v,min(e.f,remain))) ) {
                 remain -= na;
                 past += na;
                 e.f -= na;
                 ve.f += na;
                 if( !remain ) break;
             }
         }
         return past;
     }
     int flow() {
         int rt = ;
         while( bfs() ) {
             fill( cur, , n, - );
             rt += dfs(src,oo);
         }
         return rt;
     }
 }D;

 int n, src, dst, idc;
 int main() {
     scanf( "%d", &n );
     src = ;
     dst = n*n+n+;
     idc = ;
     D.init( dst, src, dst );
     int tot=;
     for( int i=; i<=n; i++ )
         for( int j=; j<=n; j++ ) {
             int bij;
             scanf( "%d", &bij );
             idc++;
             D.adde( src, idc, bij );
             D.adde( idc, i, oo );
             D.adde( idc, j, oo );
             tot += bij;
         }
     for( int i=; i<=n; i++ ) {
         int ci;
         scanf( "%d", &ci );
         D.adde( src, i,  );
         D.adde( i, dst, ci );
     }
     printf( "%d\n", tot-D.flow() );
 }

——————

我的最小割好慢。。。