「网络流 24 题」[6] 最长递增子序列

题意

对于一个序列a，解决下方三个问题

1.求出这个序列的最长递增子序列长度s

2.求出有多少个递增子序列长度为s

3.如果a[1]和a[n]可以多次使用，然后长度为s的递增子序列。

注：对于上述所有递增子序列均为不严格递增子序列。

分析

首先第一个问题用树状数组或者归并排序搞一搞就好啦。当然似乎由于数据范围比较小所以说n^2的求似乎也没有什么问题。

然后考虑第二个问题[虽然似乎也可以通过类似于树状数组一样的方法来解决啊xxx]

主要就是建图，考虑对于所有dp值为1的点，连向汇点t，所有dp值为s的点，都把源点s连向它，然后把所有点i+n连向点i，然后把所有满足(i<=j&&dp[i]+1==j)的点对，都将j连向n+i，然后跑一边最大流就好了....

真是神奇的建图。

建图真神奇。

code

#include<bits/stdc++.h>
#define inf 100000000
#define M 1020
using namespace std;
void read(int &x){
	x=0; char c=getchar();
	for (;c<48;c=getchar());
	for (;c>47;c=getchar())x=(x<<1)+(x<<3)+(c^48);
}
struct ed{
	int x,cap,nx;
}e[M*M];
int a[M],nx[M],ecnt,b[M];
void add(int x,int y,int cap){
	e[ecnt]=(ed){y,cap,nx[x]};
	nx[x]=ecnt++;
	e[ecnt]=(ed){x,0,nx[y]};
	nx[y]=ecnt++;
}
int ch(int l,int r,int x){
	for (;;){
		int mid=(l+r)>>1;
		if (b[mid]==x)return mid;
		if (b[mid]<x)l=mid+1;else r=mid-1;
	}
}
struct Dinic{
	int Q[M],level[M],nnx[M],s,t;
	bool bfs(int x){
		int l=0,r=0;
		memset(level,0,sizeof(level));
		level[Q[r++]=x]=1;
		for (;l<r;){
			x=Q[l++];
			for (int i=nx[x];~i;i=e[i].nx)if (e[i].cap>0&&!level[e[i].x]){
				level[Q[r++]=e[i].x]=level[x]+1;
				if (e[i].x==t)return 1;
			}
		}
		return level[t]>0;
	}
	int dfs(int x,int f){
		if (x==t)return f;
		int d,sum=0;
		for (int &i=nnx[x];~i;i=e[i].nx)if (e[i].cap>0&&level[x]+1==level[e[i].x]){
			d=dfs(e[i].x,min(e[i].cap,f-sum));
			e[i].cap-=d;
			e[i^1].cap+=d;
			sum+=d;
			if (sum==f)return f;
		}
		if (!sum)level[x]=0;
		return sum;
	}
	int solve(int S,int T){
		s=S; t=T;
		int res=0;
		for (;bfs(s);){
			memcpy(nnx,nx,sizeof(nx));
			res+=dfs(s,inf);
		}
		return res;
	}
}dinic;
void Max(int &x,int y){
	if (x<y)x=y;
}
int n;
struct Tree{
	#define lowbit(p) (p&(-p))
	int a[M];
	void add(int p,int x){
		for (;p<=n;p+=lowbit(p))Max(a[p],x);
	}
	int qu(int p){
		int res=0;
		for (;p;p-=lowbit(p))Max(res,a[p]);
		return res;
	}
}T;
int dp[M];
int main(){
//	freopen("1.txt","r",stdin);
	read(n);
	int i,res=0;
	for (i=1;i<=n;i++){read(a[i]); b[i]=a[i];}
	sort(b+1,b+n+1);
	int tot=unique(b+1,b+n+1)-b-1;
	for (i=1;i<=n;i++)a[i]=ch(1,tot,a[i]);//;,printf("%d ",a[i]); printf("\n");
	for (i=1;i<=n;i++){
		dp[i]=T.qu(a[i])+1;
		T.add(a[i],dp[i]);
		Max(res,dp[i]);
	}
	printf("%d\n",res);
	memset(nx,-1,sizeof(nx)); ecnt=0;
	int s=0,t=2*n+1,j;
	for (i=1;i<=n;i++){
		if (dp[i]==res)add(s,i,1);
		add(i+n,i,1);
		if (dp[i]==1)add(i,t,1);
		for (j=1;j<i;j++)if (a[j]<=a[i]&&dp[j]+1==dp[i])add(i,j+n,1);
	}
	int ans=dinic.solve(s,t);
	printf("%d\n",ans);
	add(1+n,1,inf);
	add(1,t,inf);
	if (dp[n]==res)add(s,n,inf);
	ans+=dinic.solve(s,t);
	printf("%d\n",ans);
	return 0;
}