2180. 最长递增子序列问题 (网络流拆点)

题目
题意:
给定正整数序列 x1,⋯,xn。
1.计算其最长非递减子序列的长度 s。
2.计算从给定的序列中最多可取出多少个长度为 s 的非递减子序列。（给定序列中的每个元素最多只能被取出使用一次）
3.如果允许在取出的序列中多次使用 x1 和 xn，则从给定序列中最多可取出多少个长度为 s 的非递减子序列。
思路: n<=500.第一问可以用n^2的dp，这样便于下两问的求解。

第2问的话，首先特判一下s为1的特殊情况，就是递减子序列，很好判断。我们可以在两个能够dp转移的点之间建一条流量为1的边，源点向所有dp[i] == 1的点连一条流量为1的边，所有dp[i] == s的点向汇点连一条流量为1的边。除此之外记得要拆点，因为每个点只能用一次，形如1->3->4，2->3->5这样的连边也是有可能的，但是3这个点用了两次，就不符合题意了。所以每个点拆成入点和出点，入点和出点之间连一条流量为1的边，可以限制每个数只用一次。

第3问的话，就是把x1和xn的限制从1变成INF，其余不变，基于网络流的特性，在残留网络里跑仍然可以得到结果，不需要重新建图。具体操作就是入点和出点之间的边改成INF，然后源点到x1、xn的边改成INF，x1、xn到汇点的边改成INF.实际上，其中有些情况不会出现，但是都判断也没有问题。
代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 502*4;
const int M = 2*N;
const int INF = 1e9;
int f[502];
int a[502];
int cnt = 0;
int n,m,k,T; int st,ed;
int cur[N],d[N];
int h[N],e[M],ne[M],w[M],idx = 0;
void add(int a,int b,int c)
{
	e[idx] = b,w[idx] = c,ne[idx] = h[a],h[a] = idx++;
	e[idx] = a,w[idx] = 0,ne[idx] = h[b],h[b] = idx++;
}
int now = 0;
bool bfs()
{
	queue<int> q;
	memset(d,-1,sizeof(d));
	d[st] = 0,cur[st] = h[st],q.push(st);
	while(q.size())
	{
		int u = q.front(); q.pop();
		for(int i=h[u];~i;i=ne[i])
		{
			int j = e[i];
			if(d[j]==-1&&w[i])
			{
				d[j] = d[u]+1;
				cur[j] = h[j];
				if(j==ed) return 1;
				q.push(j);
			}
		}
	}
	return 0;
}
ll find(int u,ll limit)
{
	if(u==ed) return limit;
	ll flow = 0;
	for(int i=h[u];~i&&flow<limit;i=ne[i])
	{
		int j = e[i];
		cur[u] = i;
		if(d[j]==d[u]+1&&w[i])
		{
			ll t = find(j,min(1ll*w[i],limit-flow));
			if(!t) d[j] = -1;
			flow += t,w[i] -= t,w[i^1] += t;
		}
	}
	return flow;
}
ll dinic()
{
	ll ans = 0,flow;
	while(bfs()) while(flow = find(st,INF)) ans += flow;
	return ans;
}
void solve()
{
	memset(h,-1,sizeof(h)); idx = 0;
	cin>>n;
	// cout<<"??\n";
	for(int i=1;i<=n;++i) 
	{
		cin>>a[i];
	}
	int mx = 0;
	for(int i=1;i<=n;++i)
	{
		f[i] = 1;
		for(int j=1;j<i;++j)
		{
			if(a[j]<=a[i]) f[i] = max(f[i],f[j]+1);
		}
		mx = max(mx,f[i]);
	}
	cout<<mx<<"\n";
	if(mx==1)
	{
		cout<<n<<"\n"<<n;
		return ;
	}
	ll t1,t2;
	st = 0, ed = N-1;
	for(int i=1;i<=n;++i)
	{
		add(i,i+n,1); //拆点，入点向出点连一条流量为1的边. 
		if(f[i]==1)
		{
			add(st,i,1);	
			continue;
		}
		for(int j=1;j<i;++j)
		{
			if(a[j]<=a[i]&&f[i]==f[j]+1) 
			{
				add(j+n,i,1); //出点 
			}
		}
		if(f[i]==mx)
		{
			add(n+i,ed,1); //序列结尾,出点 
		}
	}
	t1 = dinic();
	for(int i=0;i<idx;i+=2)
	{
		int l = e[i^1],r = e[i];
		if((l==1||l==n)&&r==l+n) w[i] = INF;
		else if(l==st&&(r==1||r==n)) w[i] = INF;
		else if(r==ed&&(l==n+1||l==n+n)) w[i] = INF;	
	} 	
	t2 = t1 + dinic();
	cout<<t1<<"\n"<<t2<<"\n";
}
signed main(void)
{
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	solve();
	return 0;
}