最长递增子序列

给定正整数序列 x1∼xn，以下递增子序列均为非严格递增。

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

开心~这题想了一下一遍很轻松的A了。第一个问题，最长不上升子序列 普及组难度的DP应该不用讲吧？
第二个问题，取出就不能用了。每个数只能用一次。我们想到航空路线问题的不相交路径，对于每个数我们都可以对他可以达到的递增的数连一条容量为1的边。然后问题来了。我们发现这样并没有办法做。。因为他们取数是有先后顺序的。
然后我灵光一现，想到了星际转移的分层图。然后就这样写辣。对于每个点拆成s个点。然后分层图。
在每一层中都连容量1的边。这就是一个星际转移+航空路线了。

对于第三问，只要把2条边容量限制取消就好了。

写完以后，我看了下题解。。发现和我的方法并不一样啊。。妙妙。
我先发我的方法，到时候研究一下题解。（Byvoid）

#include<bits/stdc++.h>
using namespace std;

const int MAXN=1e5+5;
const int INF=1e9;

struct edge{
    int to,next,w;
}e[MAXN<<1],e2[MAXN<<1];

int n,m,s,t;
int head[MAXN],cur[MAXN],a[1000],f[1000],cnt=1;
inline void add(int u,int v,int w){
    e[++cnt]=(edge){v,head[u],w},head[u]=cnt;
    e[++cnt]=(edge){u,head[v],0},head[v]=cnt;
}

int head2[MAXN],cur2[MAXN],cnt2=1;
inline void add2(int u,int v,int w){
    e2[++cnt2]=(edge){v,head2[u],w},head2[u]=cnt2;
    e2[++cnt2]=(edge){u,head2[v],0},head2[v]=cnt2;
}

queue<int>q;
int dep[MAXN];
bool bfs(int x){
    memset(dep,0,sizeof(dep));
    q.push(x);dep[x]=1;
    while(q.size()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].to,w=e[i].w;
            if(!dep[v]&&w){
                dep[v]=dep[u]+1;q.push(v);
            }
        }
    }
    if(!dep[t])return 0;
    return 1; 
}

int dfs(int u,int flow){
    if(u==t||flow==0)return flow;
    for(int &i=cur[u];i;i=e[i].next){
        int v=e[i].to,w=e[i].w;
        if(dep[v]==dep[u]+1&&w){
            int tem=dfs(v,min(w,flow));
            if(tem){
                e[i].w-=tem;
                e[i^1].w+=tem;
                return tem;
            }
        }
    }
    return 0;
}

int dinic(){
    int ans=0;
    while(bfs(s)){
        for(int i=s;i<=t;i++)cur[i]=head[i];
        while(int d=dfs(s,INF))ans+=d;
    }
    return ans;
}

bool bfs2(int x){
    memset(dep,0,sizeof(dep));
    q.push(x);dep[x]=1;
    while(q.size()){
        int u=q.front();q.pop();
        for(int i=head2[u];i;i=e2[i].next){
            int v=e2[i].to,w=e2[i].w;
            if(!dep[v]&&w){
                dep[v]=dep[u]+1;q.push(v);
            }
        }
    }
    if(!dep[t])return 0;
    return 1; 
}

int dfs2(int u,int flow){
    if(u==t||flow==0)return flow;
    for(int &i=cur2[u];i;i=e2[i].next){
        int v=e2[i].to,w=e2[i].w;
        if(dep[v]==dep[u]+1&&w){
            int tem=dfs2(v,min(w,flow));
            if(tem){
                e2[i].w-=tem;
                e2[i^1].w+=tem;
                return tem;
            }
        }
    }
    return 0;
}

int dinic2(){
    int ans=0;
    while(bfs2(s)){
        for(int i=s;i<=t;i++)cur2[i]=head2[i];
        while(int d=dfs2(s,INF))ans+=d;
    }
    return ans;
}

int main(){
    int ans=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    for(int i=n;i>=1;i--){
        f[i]=1;
        for(int j=n;j>i;j--){
            if(a[i]<=a[j])f[i]=max(f[i],f[j]+1);
        }
        ans=max(ans,f[i]);
    }
    s=0,t=n*ans+1;
    for(int i=1;i<n;i++){
        for(int j=i+1;j<=n;j++){
            if(a[i]<=a[j]){
                int tem=i,tem2=j;
                while(tem+n<=n*ans){
                    add(tem,tem2+n,1);
                    add2(tem,tem2+n,1);
                    tem+=n;tem2+=n;
                }
            }
        }
    }
    for(int i=1;i<=n;i++){
        add(s,i,1);
        if(i==1)add2(s,i,INF);
        else add2(s,i,1);   
    }
    for(int i=n*(ans-1)+1;i<=n*ans;i++){
        add(i,t,1); 
        if(i==n*ans)add2(i,t,INF);
        else add2(i,t,1);
    }
    printf("%d\n%d\n%d\n",ans,dinic(),dinic2());
    return 0;
}