luogu P2766 最长不下降子序列问题

关于构图：

首先，第一问dp（O(n^2) or O(nlogn)这里就不用二分优化了，详情请咨询我的blog）稳稳的。

很明显用网络流做，明明是看到标签。

很明显第三问是第二问的升级版，只需将超级源点和汇点到相应点的边权改为INF即可。

构图（p[i]表示原序列第i个的数，dp[i]表示以i为终点的最长不下降子序列的长度）：

1.超级源点和汇点：st=0,ed=n*2+1。

2.从i'向i建一条流量为1的边，这样它就可以流回去，进行下一次选择（另一个最长不下降子序列）。

3.判断，当dp[i]为1时，表示它是某些（个）最长不下降子序列的起点，由源点向它建一条流量为1的边。

4.判断，当dp[i]为maxx时，表示它是某些（个）最长不下降子序列的终点，由它向终点建一条流量为1的边。

5.判断，当p[i]<=p[j]&&(dp[j]==dp[i]+1||ans==1)时，表示它是某些（个）最长不下降子序列的中间点（不为起点和终点），由i向j'建一条流量为1的边。

注意：第二问和第三问分别求最大流。

关于代码：

#include
#include
#include
#include
#define INF 2147483647
using namespace std; 
queue f;
	struct node{int x,y,z,next;} a[100000];
	int last[10000],p[1000],dp[1000];
	int n,m,len,st,ed,ans=0;
void ins(int x,int y,int z)
{
	a[++len].x=x;a[len].y=y;a[len].z=z;a[len].next=last[x];last[x]=len;
}
void init(int x)
{
	len=-1;
	memset(last,-1,sizeof(last));
	for(int i=1;i<=n;i++)
	{
		ins(i+n,i,1),ins(i+n,i,0);
		if(dp[i]==1) ins(st,i,i==1?x:1),ins(i,st,0);
		if(dp[i]==ans) ins(i,ed,i==n?x:1),ins(ed,i,0);
		for(int j=i+1;j<=n;j++)
			if(p[i]<=p[j]&&(dp[j]==dp[i]+1||ans==1)) ins(i,n+j,1),ins(n+j,i,0);
	}
}
int h[100000];
bool bfs()
{
	memset(h,0,sizeof(h));
	h[st]=1;
	f.push(st);
	while(!f.empty())
	{
		int x=f.front();
		for(int i=last[x];i>=0;i=a[i].next)
		{
			int y=a[i].y;
			if(a[i].z>0&&h[y]==0)
			{
				h[y]=h[x]+1;
				f.push(y);
			}
		}
		f.pop();
	}
	if(h[ed]) return true; else return false;
}
int dfs(int x,int f)
{
	int s=0,t;
	if(x==ed) return f;
	for(int i=last[x];i>=0;i=a[i].next)
	{
		int y=a[i].y;
		if(a[i].z>0&&h[y]==h[x]+1&&f>s)
		{
			s+=(t=(dfs(y,min(f-s,a[i].z))));
			a[i].z-=t;
			a[i^1].z+=t;
		}
	}
	if(!s) h[x]=0;
	return s;
}
int dinic()
{
	int sum=0;
	while(bfs())
		sum+=dfs(st,INF);
	return sum;
}
int main()
{
	scanf("%d",&n);
	st=0,ed=2*n+1;
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&p[i]);
		dp[i]=1;
		for(int j=1;j