★ 最长递增子序列问题 （最多不相交路径）(分层思想) 网络流最大流

题目描述 Description

给定正整数序列x1,..... , xn  。
（1）计算其最长递增子序列的长度s。
（2）计算从给定的序列中最多可取出多少个长度为s的递增子序列。
（3）如果允许在取出的序列中多次使用x1和xn，则从给定序列中最多可取出多少个长
度为s的递增子序列。

输入描述 Input Description

第1 行有1个正整数n，表示给定序列的长度。接
下来的1 行有n个正整数x1.....xn 。

输出描述 Output Description

第1 行是最长递增子序列的长度s。第2行是可取出的长度为s 的递增子序列个数。第3行是允许在取出的序列中多次使用x1和xn时可取出的长度为s 的递增子序列个数。

样例输入 Sample Input

4
3 6 2 5

样例输出 Sample Output

2
2
3

数据范围及提示 Data Size & Hint

题目可以在这里提交：http://wikioi.com/homework/23/

注意：这题的后台数据第4,5,6,9,10组数据有误。

分析：首先还是感叹一下网络流的强大！ORZ。。。

第一问直接LIS求解即可。

第二问，对于最长递增子序列问题，这里怎么用网络流来写呢，首先的问题在于如何控制长度，这里用到了分层的思想：我们这len[i]为以a[i]开头的最长递增子序列的长度。预处理这个之后我们可以建边了。将每个点进行拆点，分成i和i+n两点，首先是从源点S出发的边，我们发现，只有len[i]==s的位置开始才存在长度为s的递增子序列，建边add(S,i,1)，然后是流入汇点T的边，当len[i]==1时，才可能是某个序列的终点，建边add(i+n,T,1)，对于上面所讲的控制长度的问题，当j>i时，如果a[j]>a[i]，并且len[i]==len[j]+1，则建边add(i+n,j,1)，到此第二问完成。

第三问，a[0]和a[n-1]可以被多次使用，则我们只需改变四条边的流量即可，分别将<0,0+n>,<n-1,n-1+n>,<S,0>,<n-1+n,T>的流量改为无穷大即可。

代码：

//Isap算法，复杂度O(n^2m)
#pragma comment(linker,"/STACK:102400000,102400000")
#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <vector>
#include <string>
#include <math.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
using namespace std;
typedef long long ll;   //记得必要的时候改成无符号
const int maxn=1005;
const int maxm=1000005;
const int INF=1000000000;
struct EdgeNode
{
    int from;
    int to;
    int cost;
    int next;
}edge[maxm];
int head[maxn],cnt;
void add(int x,int y,int z)
{
    edge[cnt].from=x;edge[cnt].to=y;edge[cnt].cost=z;edge[cnt].next=head[x];head[x]=cnt++;
    edge[cnt].from=y;edge[cnt].to=x;edge[cnt].cost=0;edge[cnt].next=head[y];head[y]=cnt++;
    //printf("%d %d %d\n",x,y,z);
}

void init()
{
    cnt=0;
    memset(head,-1,sizeof(head));
}

int S,T,n,m;
int d[maxn],gap[maxn],curedge[maxn],pre[maxn];
//curedge[]为当前弧数组，pre为前驱数组

int sap(int S,int T,int n)  //n为点数
{
    int cur_flow,flow_ans=0,u,tmp,neck,i;
    memset(d,0,sizeof(d));
    memset(gap,0,sizeof(gap));
    memset(pre,-1,sizeof(pre));
    for(i=0;i<=n;i++)curedge[i]=head[i]; //初始化当前弧为第一条邻接表
    gap[0]=n;
    u=S;
    while(d[S]<n)             //当d[S]>=n时，网络中肯定出现了断层
    {
        if(u==T)
        {
            cur_flow=INF;
            for(i=S;i!=T;i=edge[curedge[i]].to)
            {                           //增广成功，寻找瓶颈边
                if(cur_flow>edge[curedge[i]].cost)
                {
                    neck=i;
                    cur_flow=edge[curedge[i]].cost;
                }
            }
            for(i=S;i!=T;i=edge[curedge[i]].to)
            {                             //修改路径上的边容量
                tmp=curedge[i];
                edge[tmp].cost-=cur_flow;
                edge[tmp^1].cost+=cur_flow;
            }
            flow_ans+=cur_flow;
            u=neck;                     //下次增广从瓶颈边开始
        }
        for(i=curedge[u];i!=-1;i=edge[i].next)
            if(edge[i].cost&&d[u]==d[edge[i].to]+1)
               break;
        if(i!=-1)
        {
            curedge[u]=i;
            pre[edge[i].to]=u;
            u=edge[i].to;
        }
        else
        {
            if(0==--gap[d[u]])break;    //gap优化
            curedge[u]=head[u];
            for(tmp=n,i=head[u];i!=-1;i=edge[i].next)
                if(edge[i].cost)
                   tmp=min(tmp,d[edge[i].to]);
            d[u]=tmp+1;
            ++gap[d[u]];
            if(u!=S)u=pre[u];           //重标号并且从当前点前驱重新增广
        }
    }
    return flow_ans;
}

int a[maxn],dp[maxn],len[maxn];
int lis()
{
    int i,ans;
    fill(dp,dp+n,INF);
    for(i=0;i<n;i++){
        ans=lower_bound(dp,dp+n,a[i])-dp;
        dp[ans]=a[i];
    }
    for(i=n-1;i>=0;i--){
        if(dp[i]!=INF)return i+1;
    }
}
int solve(int st)
{
    int i,ans;
    fill(dp,dp+n+1,INF);
    for(i=st;i<n;i++){
        if(a[i]<a[st])continue;
        ans=lower_bound(dp+st,dp+n,a[i])-dp;
        dp[ans]=a[i];
    }
    for(i=n-1;i>=0;i--){
        if(dp[i]!=INF)return i-st+1;
    }
}

int main()
{
    int i,j,Length,c1,c2,c3;
    while(~scanf("%d",&n))
    {
        for(i=0;i<n;i++)scanf("%d",&a[i]);
        Length=c1=lis();
        for(i=0;i<n;i++)len[i]=solve(i);
        //for(i=0;i<n;i++)printf("%d ",len[i]);
        S=2*n; T=S+1; init();
        for(i=0;i<n;i++){
            add(i,i+n,1);
            if(len[i]==Length)add(S,i,1);
            if(len[i]==1)add(i+n,T,1);
            for(j=i+1;j<n;j++){
                if(a[j]>a[i]&&len[i]==len[j]+1)
                    add(i+n,j,1);
            }
        }
        c2=sap(S,T,T+1);

        //printf("\n\n\n");

        S=2*n; T=S+1; init();
        for(i=0;i<n;i++){
            if(i==0||i==n-1)add(i,i+n,INF);
            else add(i,i+n,1);
            if(len[i]==Length){
                if(i==0)add(S,i,INF);
                else add(S,i,1);
            }
            if(len[i]==1){
                if(i==n-1)add(i+n,T,INF);
                else add(i+n,T,1);
            }
            for(j=i+1;j<n;j++){
                if(a[j]>a[i]&&len[i]==len[j]+1)
                    add(i+n,j,1);
            }
        }
        c3=sap(S,T,T+1);
        printf("%d\n%d\n%d\n",c1,c2,c3);
    }
    return 0;
}