[网络流24题] 洛谷P2766 最长不下降子序列 最大流

给出$n\le 5e2$个数，然后分别求最长不下降序列长度$s$，选出最多的长度为$s$的不下降子序列的个数，以及$1$和$n$可以选择无限次情形下最多的不下降子序列的个数。
这个题比较先进。对于第一个问题直接dp当前位置开始的最长的长度。
然后对于第二问，我们知道如果dp值为$s$，才有可能是开头，如果dp值为1，则可以成为末尾。建立一个源点向所有的dp值为$s$的点连流量$1$的边，所有dp值为$1$的点向汇点连流量为$1$的边，并且满足$i<j$且$a_i\le a_j$且$d{p}_i=d{p}_j+1$的点对$i$向$j$连一条流量为$1$的边。每个点只能用一次，所以要拆点，自己向自己连一条流量为$1$的边。
对于第三问，和第二问的区别在于$1$和$n$结点可以使用无限次。判断一下这两个点是否可以作为开头和结尾，如果是的话，源点连无限流量到$1$，$1$连无限流量到自身。结尾同理。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int N=5e2+7;
int a[N],dp[N]; 
namespace Dinic {
	const int N=1e4+3,M=2e5+3; 
	int front[N],head[N],dep[N],tot=-1,n,S,T;
	struct edge{ int v,nxt;ll w; }e[2*M];
	void init() {
		tot=-1;
		for(int i=1;i<=n;i++) front[i]=-1;
	}
	void add(int u,int v,ll w) {
		e[++tot].v=v,e[tot].w=w,e[tot].nxt=front[u],front[u]=tot;
		e[++tot].v=u,e[tot].w=0,e[tot].nxt=front[v],front[v]=tot;
	}
	bool bfs() {
		for(int i=1;i<=n;i++) dep[i]=0;
		queue<int> q;
		q.push(S);
		dep[S]=1;
		while(!q.empty()) {
			int u=q.front();
			q.pop();
			for(int i=front[u];i!=-1;i=e[i].nxt) {
				int v=e[i].v;
				if(dep[v]==0&&e[i].w) {
					dep[v]=dep[u]+1;
					q.push(v); 
				}
			}
		}
		return dep[T];
	}
	ll dfs(int cur,ll dis=1e18) {
		if(cur==T) return dis;
		ll flow=0,sum=0;
		for(int &i=head[cur];i!=-1;i=e[i].nxt) {
			int v=e[i].v;
			if(dep[v]==dep[cur]+1&&e[i].w) {
				flow=dfs(v,min(dis,e[i].w));
				dis-=flow,sum+=flow;
				e[i].w-=flow,e[i^1].w+=flow;
				if(!dis) break;
			}
		}
		return sum;
	}
	ll solve() {
		ll ans=0;
		while(bfs()) {
			for(int i=1;i<=n;i++) head[i]=front[i];
			ans+=dfs(S);
		}
		return ans;
	}
}
int n,mx=0;
void gao1() {
	int S=2*n+1;
	int T=2*n+2;
	Dinic::n=2*n+2;
	Dinic::S=S;
	Dinic::T=T;
	Dinic::init();
	for(int i=1;i<=n;i++) {
		if(dp[i]==mx) {
			Dinic::add(S,i,1);
		}
	} 
	for(int i=1;i<=n;i++) {
		if(dp[i]==1) {
			Dinic::add(i+n,T,1);
		}
	}
	for(int i=1;i<=n;i++) 
		Dinic::add(i,i+n,1);
	for(int i=1;i<=n;i++) {
		for(int j=i+1;j<=n;j++) {
			if(a[j]>=a[i]&&dp[i]==dp[j]+1) {
				Dinic::add(i+n,j,1);
			}
		}
	}
}
void gao2() {
	gao1();
	int S=2*n+1;
	int T=2*n+2; 
	if(dp[1]==mx) Dinic::add(1,1+n,1e9),Dinic::add(S,1,1e9);
	if(dp[n]==1) Dinic::add(n,n+n,1e9),Dinic::add(n+n,T,1e9);
}
int main() {
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int i=n;i>=1;i--) {
		dp[i]=1;
		for(int j=i+1;j<=n;j++) {
			if(a[j]>=a[i]) {
				dp[i]=max(dp[i],dp[j]+1);
			}		
		}
		mx=max(mx,dp[i]);
	}
	printf("%d\n",mx);
	gao1();
	printf("%lld\n",Dinic::solve());
	gao2();
	printf("%lld\n",Dinic::solve()); 
	return 0;
}