网络流24题--最长不下降子序列个数

 残余网络

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef int lint;
const lint inf = 0x3f3f3f3f;
const int maxn = 1005;
struct dinic{
    static const int N = 5005;
    static const int M = 5000005;
    static const lint inf = 0x3f3f3f3f;
    int he[N],ne[M],ver[M],d[N];
    lint edge[M];
    int s,t,tot;
    queue<int> que;
    void init( int n ){
        for( int i = 0;i <= n;i++ ){
            he[i] = 0;
        }
        tot = 1;
    }
    void add( int x,int y,lint z ){
        ver[++tot] = y; ne[ tot ] = he[x]; he[x] = tot; edge[ tot ] = z;
        ver[++tot] = x; ne[tot] = he[y]; he[y] = tot;edge[tot] = 0;
    }
    bool bfs(){
        memset( d,0,sizeof( d ) );
        while( que.size() ) que.pop();
        que.push( s );
        d[s] = 1;
        while( que.size() ){
            int x = que.front();
            que.pop();
            for( int cure = he[x];cure;cure = ne[cure] ){
                int y = ver[cure];
                if( edge[cure] && !d[y] ){
                    que.push( y );
                    d[ y ] = d[x] + 1;
                    if( y == t ) return 1;
                }
            }
        }
        return 0;
    }
    lint dfs( int x,lint flow ){
        if( x== t ) return flow;
        lint rest = flow,k;
        for( int cure = he[x];cure && rest;cure = ne[cure] ){
            int y = ver[cure];
            if( edge[cure] && d[ y ] == d[x] + 1 ){
                k = dfs( y,min( rest,edge[cure] ) );
                if( !k ) d[ y ] = 0;
                edge[cure] -= k;
                edge[ cure^1 ] += k;
                rest -= k;
            }
        }
        return flow - rest;
    }
    lint max_flow( int x,int y ){
        s = x; t = y;
        lint maxflow = 0;
        lint flow = 0;
        while( bfs() )
            while( flow = dfs( s,inf ) ) maxflow += flow;
        return maxflow;
    }
} g;
int a[maxn],dp[maxn];
int main(){
   int n;scanf("%d",&n);
   for( int i = 1;i <= n;i++ ) scanf("%d",&a[i]);
   int s = 0;
   for( int i = 1;i <= n;i++ ){
       dp[i] = 1;
       for( int j = 1;j < i;j++ ){
           if( a[j] <= a[i] ) dp[i] = max( dp[i],dp[j] +1 );
       }
       s = max( s,dp[i] );
   }
   printf("%d\n",s);
   int S = 0,T=  2*n+1;
   g.init(T);
   for( int i = 1;i <= n;i++ ){
       g.add( i,i+n,1 );
       if( dp[i] == 1 ){
           g.add( S,i,1 );
       }
       if( dp[i] == s ){
           g.add( i+n,T,1 );
       }
       for( int j = 1;j < i;j++ ){
           if( dp[j] == dp[i]-1 && a[j] <= a[i] ){
               g.add( j+n,i,1 );
           }
       }
   }
   int ans = g.max_flow(S,T);
   printf("%d\n",ans);
   g.add( 1,1+n,inf );
   g.add( n,n+n,inf );
   g.add( S,1,inf );
   if( dp[n] == s ){
       g.add( n+n,T,inf );
   }
   ans += g.max_flow(S,T);
   printf("%d\n",ans);
   return 0;
}