Libre 6005 「网络流 24 题」最长递增子序列 / Luogu 2766 最长递增子序列问题(网络流,最大流)

Libre 6005 「网络流 24 题」最长递增子序列 / Luogu 2766 最长递增子序列问题(网络流,最大流)

Description

问题描述:

给定正整数序列x1,...,xn 。

(1)计算其最长递增子序列的长度s。

(2)计算从给定的序列中最多可取出多少个长度为s的递增子序列。

(3)如果允许在取出的序列中多次使用x1和xn,则从给定序列中最多可取出多少个长度为s的递增子序列。

编程任务:

设计有效算法完成(1)(2)(3)提出的计算任务。

Input

第1 行有1个正整数n,表示给定序列的长度。接下来的1 行有n个正整数n:x1, ..., xn。

Output

第1 行是最长递增子序列的长度s。第2行是可取出的长度为s 的递增子序列个数。第3行是允许在取出的序列中多次使用x1和xn时可取出的长度为s 的递增子序列个数。

Sample Input

4
3 6 2 5

Sample Output

2
2
3

Http

Libre:https://loj.ac/problem/6005
Luogu:https://www.luogu.org/problem/show?pid=2766

Source

网络流,最大流

解决思路

看清题目,不是最长递增子序列是最长不下降子序列。
这道题目首先运用动态规划的方式求出最长不下降子序列,这也是第一问的内容。注意,本题不能使用单调队列的方式,因为要求出到每一个数的最长不下降子序列长度(后面记为F),这在后两问中要用。
那么如何求解第二问呢?
我们把每一个数拆成两个点入点和出点,在每一个数的入点和出点之间连容量为1的边,同时设置一个源点一个汇点。从前往后扫描每一个数,若发现第i个数的F[i]最长不下降子序列长度,则在源点与i的出点之间连一条容量为1的边。若F[i]1,则在其出点与汇点之间连一条容量为1的边。并且,对于任何数i,扫描其前面的每一个数j,若F[i]==F[j]+1且第j位的数<=第i位的数,则在i的出点与j的入点之间连一条容量为1的边。
这样建图,有点类似于分层图的思想,从最高层的F为最长不下降子序列长度,往下每一层长度减一,直到最底下一层长度为1。这样我们跑一边最短路就可以了。
至于第三问,我们只要在重新建图的时候把1到汇点,1的入点到出点,n的入点到出点,源点到n(如果存在的话)这几条边设置为无穷大即可。
但要注意,第三问要特判一下递减的情况,因为这样最长不下降子序列长度为1,跑最大流会出现无穷大的流的情况。
另:这里用Dinic求解最大流,具体请移步我的这篇文章

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

const int maxN=2000;
const int maxM=maxN*maxN*4;
const int inf=147483647;

class Edge
{
public:
	int u,v,flow;
};

int n;
int cnt=-1;
int F[maxN];
int Arr[maxN];
int Head[maxN];
int Next[maxM];
Edge E[maxM];
int depth[maxN];
int cur[maxN];
int Q[maxM];

void Add_Edge(int u,int v,int flow);
bool bfs();
int dfs(int u,int flow);

int main()
{
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
		scanf("%d",&Arr[i]);
	for (int i=1;i<=n;i++)//动态规划求出最长不下降子序列
	{
		F[i]=1;
		for (int j=1;j<i;j++)
		    if (Arr[j]<=Arr[i])
				F[i]=max(F[i],F[j]+1);
	}
	int maxlength=1;
	for (int i=1;i<=n;i++)
		maxlength=max(maxlength,F[i]);//得出第一问答案
	printf("%d\n",maxlength);
	memset(Head,-1,sizeof(Head));
	for (int i=1;i<=n;i++)//构造第二问的图
	{
		Add_Edge(i,i+n,1)//连接入点和出点
		if (F[i]==maxlength)//若与最长长度相同,则连接源点
			Add_Edge(0,i,1);
		if (F[i]==1)//若为最小长度1,则连接汇点
			Add_Edge(i+n,n*2+1,1);
		for (int j=1;j<i;j++)
			if ((F[j]==F[i]-1)&&(Arr[j]<=Arr[i]))//向前统计能连的
				Add_Edge(i+n,j,1);
	}
	int Ans=0;//求解最大流
	while (bfs())
	{
		for (int i=0;i<=2*n+1;i++)
			cur[i]=Head[i];
		while (int di=dfs(0,inf))
			Ans+=di;
	}
	cout<<Ans<<endl;
	memset(Head,-1,sizeof(Head));
	cnt=-1;
	for (int i=1;i<=n;i++)//构造第三问的图
	{
		int nowflow=1;
		if ((i==1)||(i==n))//1和n的流量为无穷大
			nowflow=inf;
        if (maxlength==1)//注意这里特判递减序列的情况
            Add_Edge(i,i+n,1);
        else
            Add_Edge(i,i+n,inf);
		if (F[i]==maxlength)
			Add_Edge(0,i,nowflow);
		if (F[i]==1)
			Add_Edge(i+n,n*2+1,nowflow);
		for (int j=1;j<i;j++)
			if ((F[j]==F[i]-1)&&(Arr[j]<=Arr[i]))
				Add_Edge(i+n,j,1);
	}
	Ans=0;//求解最大流
	while (bfs())
	{
		for (int i=0;i<=2*n+1;i++)
			cur[i]=Head[i];
		while (int di=dfs(0,inf))
			Ans+=di;
	}
	cout<<Ans<<endl;
	return 0;
}

void Add_Edge(int u,int v,int flow)
{
	cnt++;
	Next[cnt]=Head[u];
	Head[u]=cnt;
	E[cnt].u=u;
	E[cnt].v=v;
	E[cnt].flow=flow;

	cnt++;
	Next[cnt]=Head[v];
	Head[v]=cnt;
	E[cnt].v=u;
	E[cnt].u=v;
	E[cnt].flow=0;
}

bool bfs()
{
	memset(depth,-1,sizeof(depth));
	int h=1,t=0;
	Q[1]=0;
	depth[0]=1;
	do
	{
		t++;
		int u=Q[t];
		//cout<<u<<endl;
		for (int i=Head[u];i!=-1;i=Next[i])
		{
			int v=E[i].v;
			if ((depth[v]==-1)&&(E[i].flow>0))
			{
				depth[v]=depth[u]+1;
				h++;
				Q[h]=v;
			}
		}
	}
	while (t!=h);
	if (depth[n*2+1]==-1)
		return 0;
	return 1;
}

int dfs(int u,int flow)
{
	if (u==n*2+1)
		return flow;
	for (int &i=cur[u];i!=-1;i=Next[i])
	{
		int v=E[i].v;
		if ((depth[v]==depth[u]+1)&&(E[i].flow>0))
		{
			int di=dfs(v,min(flow,E[i].flow));
			if (di>0)
			{
				E[i].flow-=di;
				E[i^1].flow+=di;
				return di;
			}
		}
	}
	return 0;
}
posted @ 2017-08-03 17:22  SYCstudio  阅读(89)  评论(0编辑  收藏