网络流24题之最长递增子序列问题

这题有点意思了。

三问：

①dp解决答案mxl

②网络流：

建图：

a.将一个位置i拆成两个点i1和i2

b.自源点向所有f[i]=mxl的i1连容量为1的边

c.自所有f[i]=1的i2向汇点连容量为1的边

d.每个i1向i2连容量为1的边

e.对于i<j&&a[i]<=a[j]&&f[i]+1==f[j]

由j2向i1连容量为1的边

跑最大流，即为第二问答案

③    将②中的s→1和n 的边容量扩为inf，将t→1和n的边容量扩为inf。

当然原图中如果没这条边就不管了

将1和n分别拆成的两点间容量扩为inf。

跑最大流，即为第三问答案

代码：

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>
#include<climits>
#include<ctime>
#include<queue>
#include<stack>
#include<map>
#include<set>
#define N 4002
#define M N*N
#define inf 1<<26
using namespace std;
int read()
{
	int x=0,f=1;char ch;
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
	return x*f;
}
int n,ls[N],lis[N],mxl;
int s,t;
int head[N],pos=-1,cur[N];
struct edge{int to,next,c;}e[M];
void add(int a,int b,int c)
{pos++;e[pos].to=b,e[pos].c=c,e[pos].next=head[a],head[a]=pos;}
queue<int>Q;bool vis[N];int dis[N];
bool bfs()
{
	for(int i=s;i<=t;i++)dis[i]=-1,vis[i]=0;
	vis[s]=1,dis[s]=0;Q.push(s);
	while(!Q.empty())
	{
		int u=Q.front();Q.pop();
		for(int i=head[u];i!=-1;i=e[i].next)
		{
			int v=e[i].to;
			if(v==t)
				v=t;
			if(vis[v]||e[i].c<=0)continue;
			dis[v]=dis[u]+1;
			vis[v]=1;Q.push(v);
		}
	}return vis[t];
}
int dfs(int u,int a)
{
	if(u==t||!a)return a;
	int f,flow=0;
	for(int &i=cur[u];i!=-1;i=e[i].next)
	{
		int v=e[i].to;
		if(dis[v]==dis[u]+1&&(f=dfs(v,min(a,e[i].c)))>0)
		{
			flow+=f,a-=f;
			e[i].c-=f,e[i^1].c+=f;
			if(!a)break;
		}
	}return flow;
}
int dinic()
{
	int ret=0;
	while(bfs())
	{
		for(int i=s;i<=t;i++)cur[i]=head[i];
		ret+=dfs(s,inf);
	}return ret;
}
void dp()
{
	for(int i=1;i<=n;i++)lis[i]=1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<i;j++)
			if(ls[j]<=ls[i])lis[i]=max(lis[i],lis[j]+1);
	for(int i=1;i<=n;i++)mxl=max(mxl,lis[i]);
}
void build()
{
	for(int i=1;i<=n;i++)
	{
		if(lis[i]==mxl)
			{add(s,i,1);add(i,s,0);}
		if(lis[i]==1)
			{add(i+n,t,1);add(t,i+n,0);}
	}
	for(int i=1;i<=n;i++)
		for(int j=i+1;j<=n;j++)
			if(lis[j]==lis[i]+1&&ls[i]<=ls[j])
				{add(j+n,i,1);add(i,j+n,0);}
	for(int i=1;i<=n;i++){add(i,i+n,1);add(i+n,i,0);}
}
void rebuild()
{
	memset(head,-1,sizeof(head));pos=-1;
	for(int i=1;i<=n;i++)
	{
		if(lis[i]==mxl)
			{add(s,i,1);add(i,s,0);}
		if(lis[i]==1)
			{add(i+n,t,1);add(t,i+n,0);}
	}
	for(int i=1;i<=n;i++)
		for(int j=1;j<i;j++)
			if(lis[j]+1==lis[i]&&ls[j]<=ls[i])
				{add(i+n,j,1);add(j,i+n,0);}
	for(int i=1;i<=n;i++){add(i,i+n,1);add(i+n,i,0);}
	add(1,1+n,inf);add(1+n,1,0);
	add(n,n+n,inf);add(n+n,n,0);
	if(lis[1]==mxl){add(s,1,inf);add(1,s,0);}
	else if(lis[1]==1){add(1,t,inf);add(t,1,0);}
	if(lis[n]==mxl){add(s,n,inf);add(n,s,0);}
	else if(lis[n]==1){add(n,t,inf);add(t,n,0);}
}
void init(){memset(head,-1,sizeof(head));}
int main()
{
	freopen("alis.in","r",stdin);
	freopen("alis.out","w",stdout);
	n=read();s=0,t=n*2+1;
	for(int i=1;i<=n;i++)
		ls[i]=read();dp();
	printf("%d\n",mxl);init();
	build();printf("%d\n",dinic());
	rebuild();printf("%d\n",dinic());
}