HDU 3998 Sequence 最长上升子序列+网络流 求不相交的最长上升子序列个数

题目大意：给你一个序列，求不相交的最长上升子序列的个数。

思路：我们找0为超级源点，2*n+1为超级汇点，我们先dp出LIS，在dp的过程中建图，把一个点拆成两个点i和i+n在它们之间连一条权值为1的边，如果a[i]>a[j]&&dp[i]<=dp[j]+1  那么我们在j+n和i之间建一条权值为1的边，最后我们求出最长上升子序列的长度后，把dp[i]=1的点和超级源点相连边的权值为1，代表这个点是上升子序列的起点，把dp[i]=len(len为最长上升子序列的长度)的点和超级汇点连一条权值为1的边，代表这个点是最长上升子序列的终点，之后跑一遍最大流，最多能流到终点的流也就是不相交的最长上升子序列的个数。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define inf 0x3fffffff
using namespace std;
const int maxn=1010;
int first[maxn],sign;
int s,t,d[maxn],dp[maxn],cur[maxn];
int n,a[maxn];
struct node
{
    int to,w,next;
}edge[maxn*20];
void creat()
{
    memset(first,-1,sizeof(first));
    sign=0;
}
void add_edge(int u,int v,int w)
{
    edge[sign].w=w;
    edge[sign].to=v;
    edge[sign].next=first[u];
    first[u]=sign++;
}
int LIS()
{
    int ans=1;
    for(int i=1;i<=n;i++)
    {
        dp[i]=1;
        for(int j=1;j<i;j++)
        {
            if(a[i]>a[j]&&dp[i]<=dp[j]+1)
            {
                dp[i]=dp[j]+1;
                add_edge(j+n,i,1);
                add_edge(i,j+n,0);
            }
        }
        ans=max(ans,dp[i]);
    }
    for(int i=1;i<=n;i++)
    {
        add_edge(i,i+n,1);
        add_edge(i+n,i,0);
        if(dp[i]==1)
        {
            add_edge(s,i,1);
            add_edge(i,s,0);
        }
        if(dp[i]==ans)
        {
            add_edge(i+n,t,1);
            add_edge(t,i+n,0);
        }
    }
    return ans;
}
int bfs()
{
    queue<int>q;
    q.push(s);
    memset(d,0,sizeof(d));
    d[s]=1;
    while(!q.empty())
    {
        int top=q.front();
        q.pop();
        for(int i=first[top];~i;i=edge[i].next)
        {
            int to=edge[i].to;
            if(edge[i].w>0&&d[to]==0)
            {
                d[to]=d[top]+1;
                if(to==t)return 1;
                q.push(to);
            }
        }
    }
    return d[t]!=0;
}
int dfs(int top,int flow)
{
    if(top==t)return flow;
    int ans=0,x=0;
    for(int i=cur[top];~i;i=edge[i].next)
    {
        int to=edge[i].to;
        if(edge[i].w>0&&d[to]==d[top]+1)
        {
            x=dfs(to,min(flow-ans,edge[i].w));
            ans+=x;
            edge[i].w-=x;
            edge[i^1].w+=x;
            if(edge[i].w)
                cur[top]=i;
            if(flow==ans)return flow;
        }
    }
    if(ans==0)
        d[top]=0;
    return ans;
}
int dinic()
{
    int ans=0;
    while(bfs())
    {
        for(int i=0;i<maxn;i++)
            cur[i]=first[i];
        ans+=dfs(s,inf);
    }
    return ans;
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        creat();
        memset(dp,0,sizeof(dp));
        s=0,t=2*n+1;
        for(int i=1;i<=n;i++)
            scanf("%d",&a[i]);
        int len=LIS();
        int ans=dinic();
        printf("%d\n",len);
        printf("%d\n",ans);
    }
    return 0;
}