loj6005「网络流 24 题」最长递增子序列（dp+最大流）

首先第一问就是dp求就好啦，写了nlogn的。

现在我们已经有了dp[i]，表示以第i个数结尾的lis是多长。考虑如何建图实现第二问的限制，把每个点拆成两点，建边，容量为1，这样就满足了每个点最多被经过一次，然后源点向所有dp[i]为1的点建边，容量为1，所有dp[i]==ans的点向T建边，容量为1，对于每一个i，向所有能接在他后面的数j建边，容量为1。（因为会有相等的，所以一定要严格满足dp[i]+1=dp[j]）。然后跑最大流即可。

对于第三问，$a_1和a_n$可以随便用，因此我们把这两点内部的边的容量改为inf即可。还有s向1，n向t的边，如果有，也改成inf。

注意如果n=1,第三问要特判答案应该为1。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define N 510
#define inf 0x3f3f3f3f
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,a[N],h[N<<1],num=1,T=1001,lev[N<<1],d[N],dp[N],ans=0;
struct edge{
    int to,next,val;
}data[N*N<<1];
inline void add(int x,int y,int val){
    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].val=val;
    data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].val=0;
}
inline bool bfs(){
    queue<int>q;memset(lev,0,sizeof(lev));
    lev[0]=1;q.push(0);
    while(!q.empty()){
        int x=q.front();q.pop();
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(lev[y]||!data[i].val) continue;
            lev[y]=lev[x]+1;q.push(y);
        }
    }return lev[T];
}
inline int dinic(int x,int low){
    if(x==T) return low;int tmp=low;
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(lev[y]!=lev[x]+1||!data[i].val) continue;
        int res=dinic(y,min(tmp,data[i].val));
        if(!res) lev[y]=0;tmp-=res;data[i].val-=res;data[i^1].val+=res;
        if(!tmp) return low;
    }return low-tmp;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();for(int i=1;i<=n;++i) a[i]=read();
    if (n==1){
        printf("1\n1\n1\n");
        return 0;
    }
    for(int i=1;i<=n;++i){
        if(a[i]>=d[ans]) d[++ans]=a[i],dp[i]=ans;
        else{
            int id=upper_bound(d+1,d+ans+1,a[i])-d;
            d[id]=a[i];dp[i]=id;
        }
    }printf("%d\n",ans);
    for(int i=1;i<=n;++i){
        add(i,n+i,1);if(dp[i]==1) add(0,i,1);if(dp[i]==ans) add(n+i,T,1);
        for(int j=i+1;j<=n;++j) if(a[j]>=a[i]&&dp[j]==dp[i]+1) add(n+i,j,1);
    }int res=0;while(bfs()) res+=dinic(0,inf);
    printf("%d\n",res);num=1;memset(h,0,sizeof(h));
    for(int i=1;i<=n;++i){
        add(i,n+i,i==1||i==n?inf:1);
        if(dp[i]==1) add(0,i,i==1?inf:1);
        if(dp[i]==ans) add(n+i,T,i==n?inf:1);
        for(int j=i+1;j<=n;++j) if(a[j]>=a[i]&&dp[j]==dp[i]+1) add(n+i,j,1);
    }res=0;while(bfs()) res+=dinic(0,inf);
    printf("%d\n",res);
}