loj6005「网络流 24 题」最长递增子序列

http://www.elijahqi.win/2017/12/10/loj6005%e3%80%8c%e7%bd%91%e7%bb%9c%e6%b5%81-24-%e9%a2%98%e3%80%8d%e6%9c%80%e9%95%bf%e9%80%92%e5%a2%9e%e5%ad%90%e5%ba%8f%e5%88%97/
题目描述

给定正整数序列 x1∼xn x_1 \sim x_n x​1​​∼x​n​​，以下递增子序列均为非严格递增。

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

输入格式

文件第 1 1 1 行有 1 1 1 个正整数 n n n，表示给定序列的长度。接下来的 1 1 1 行有 n n n 个正整数 x1∼xn x_1 \sim x_n x​1​​∼x​n​​。
输出格式

第 1 1 1 行是最长递增子序列的长度 s s s。第 2 2 2 行是可取出的长度为 s s s 的递增子序列个数。第 3 3 3 行是允许在取出的序列中多次使用 x1 x_1 x​1​​ 和 xn x_n x​n​​ 时可取出的长度为 s s s 的递增子序列个数。
样例
样例输入

4
3 6 2 5

样例输出

2
2
3

数据范围与提示

1≤n≤500 1 \leq n \leq 500 1≤n≤500

觉得自己真心废柴啊 什么都想不出这个建图方法想了半天也没想到的x

然后n*log(n)的dp也写不对 好丧啊

首先这个n*log(n)的dp是 设b数组为长度为下标的最小值是多少

那么我每次二分找一下比我小的最小值 然后记下他的长度更新我当前的dp值

然后这条件二和条件三应该怎么满足

可以知道 dp[i]表示一定选i号这个 并且以i结尾的最长上升子序列长度最长是多少

那么所有dp[i]=1的地方一定是起点 所有dp[i]=ans1的地方一定是终点

然后把每个数拆成两个点 中间权值为1 表示限制着一个数只能用一次

然后循环的时候 因为最长上升子序列是递增的所以我每回需要向后扫描 观察什么时候我dp[j]=dp[i]+1 并且满足a[j]>=a[i]的时候表明我有可能能从这里转移过去我建一条权值为1 的点 跑最大流即可

最后 针对第三问只需要把和1或者n有关的点改成无穷的权值即可

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1100
#define inf 0x3f3f3f3f
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x;
}
int n,b[N],a[N],dp[N],num=1,h[N],level[N],T;
struct node{
    int y,z,next;
}data[N*N];
inline void insert1(int x,int y,int z){
    data[++num].y=y;data[num].next=h[x];data[num].z=z;h[x]=num;
    data[++num].y=x;data[num].next=h[y];data[num].z=0;h[y]=num;
}
inline bool bfs(){
    memset(level,0,sizeof(level));queue<int>q;q.push(0);level[0]=1;
    while(!q.empty()){
        int x=q.front();q.pop();
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z;
            if (level[y]||!z) continue;level[y]=level[x]+1;q.push(y);if (y==T) return 1;
        }
    }return 0;
}
inline int dfs(int x,int s){
    if (x==T) return s;int ss=s;
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (level[x]+1==level[y]&&z){
            int xx=dfs(y,min(z,s));s-=xx;data[i].z-=xx;data[i^1].z+=xx;if (!xx) level[y]=0;
            if (!s) return ss;
        }
    }return ss-s;
}
int main(){
    freopen("loj.in","r",stdin);
    n=read();int ans1=0;T=n*2+1;
    for (int i=1;i<=n;++i) a[i]=read();memset(b,0x3f,sizeof(b));b[1]=a[1];dp[1]=1;
    for (int i=2;i<=n;++i){
        int len=upper_bound(b+1,b+n+1,a[i])-b-1;dp[i]=len+1; 
        if(a[i]<b[len+1]) b[len+1]=a[i];ans1=max(ans1,len+1);
    }printf("%d\n",ans1);int ans=0;
    //for (int i=1;i<=n;++i) printf("%d ",dp[i]);printf("\n");
    for (int i=1;i<=n;++i) {
        if (dp[i]==1) insert1(0,i,1);insert1(i,i+n,1);
        if (dp[i]==ans1) insert1(i+n,T,1);
        for (int j=i+1;j<=n;++j) if (dp[j]==dp[i]+1&&a[j]>=a[i]) insert1(i+n,j,1);
    }
    while(bfs()) ans+=dfs(0,inf);printf("%d\n",ans);num=1;memset(h,0,sizeof(h));ans=0;
    for (int i=1;i<=n;++i) {
        if (dp[i]==1) insert1(0,i,i==1?inf:1);
        if(i==1||i==n) insert1(i,i+n,inf);else insert1(i,i+n,1);
        if (dp[i]==ans1) insert1(i+n,T,i==n?inf:1);
        for (int j=i+1;j<=n;++j) if (dp[j]==dp[i]+1&&a[j]>=a[i]) insert1(i+n,j,1);
    }while(bfs()) ans+=dfs(0,inf);printf("%d\n",ans);
    return 0;
}