NKOI 1941 最长递增子序列

【线性规划与网络流24题 6】最长递增子序列

Time Limit:10000MS  Memory Limit:65536K
Total Submit:19 Accepted:2 
Case Time Limit:1000MS

Description

给定正整数序列x1, .., xn。 
（1）计算其最长递增子序列的长度s。 
（2）计算从给定的序列中最多可取出多少个长度为s的递增子序列。 
（3）如果允许在取出的序列中多次使用x1和 xn，则从给定序列中最多可取出多少个长度为s的递增子序列。 

编程任务： 
设计有效算法完成（1）（2）（3）提出的计算任务。

Input

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

Output

程序运行结束时，将任务（1）（2）（3）的解答输出。 
第1 行是最长递增子序列的长度s。第2行是可取出的长度为s的递增子序列个数。第3行是允许在取出的序列中多次使用x1和xn时可取出的长度为s的递增子序列个数。

Sample Input

4
3 6 2 5

Sample Output

2
2
3

Source

感谢 Wo_ai_WangYuan 放上数据

对于第一个问题，用f[i]表示到第i个数能得到的递增子序列的个数，很容易dp求解,ans=max(ans,f[i])

对于第二个问题，就要用网络流解决了

源点连f[i]==1的i,汇点连f[j]==ans的j,这里注意拆点，两个点之间容量为1，保证对于每一个序列每个数最多用一次，求最大流就是答案

对于第三个问题，源点和1号点的容量和n号点和汇点的容量都变为inf,这样对于不同的序列就可以用多次1号点和n号点了

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=1005,inf=1e9;
inline void _read(int &x){
    char t=getchar();bool sign=true;
    while(t<'0'||t>'9')
    {if(t=='-')sign=false;t=getchar();}
    for(x=0;t>='0'&&t<='9';t=getchar())x=x*10+t-'0';
    if(!sign)x=-x;
}
struct edge{  
   int a,b,c,flow;  
   edge(){}  
   edge(int fr,int t,int c,int f)  
       :a(fr),b(t),c(c),flow(f){}  
};  
int s[maxn],f[maxn],ans=-inf,op,ed;
int dis[maxn],vd[maxn],g[maxn][maxn],n;
void dp(){
	int i,j;
	for(i=1;i<=n;i++){
		f[i]=1;
		for(j=1;j<i;j++)
		    if(s[i]>=s[j]&&f[i]<f[j]+1)
		        f[i]=f[j]+1;
		ans=max(ans,f[i]);
	}
}
int dfs(int u,int flow){    
    if(u==ed)return flow;    
    int v,temp,delta=0;    
    for(v=1;v<=ed;v++)    
        if(g[u][v]&&dis[u]==dis[v]+1){    
            temp=dfs(v,min(flow-delta,g[u][v]));    
            g[u][v]-=temp;    
            g[v][u]+=temp;    
            delta+=temp;    
            if(delta==flow||dis[op]>=ed)return delta;    
        }    
    if(dis[op]>=ed)return delta;    
    vd[dis[u]]--;    
    if(vd[dis[u]]==0)dis[op]=ed;    
    dis[u]++;    
    vd[dis[u]]++;    
    return delta;    
}
int main(){
	_read(n);
	int i,j;
	for(i=1;i<=n;i++)_read(s[i]);
	dp();
	cout<<ans<<endl;
	if(ans==1){
		cout<<n<<endl<<n<<endl;
		return 0;
	}
	op=2*n+1,ed=op+1;
	for(i=1;i<=n;i++){
		g[i][i+n]=1;
		if(f[i]==1)g[op][i]=1;
	    if(f[i]==ans)g[i+n][ed]=1;
	}
	for(i=2;i<=n;i++)
		for(j=1;j<i;j++)
		    if(s[i]>=s[j]&&f[i]==f[j]+1)
		        g[j+n][i]=1;
	int maxflow=0;
	while(dis[op]<ed)maxflow+=dfs(op,inf);
	cout<<maxflow<<endl;
	maxflow=0;
	memset(dis,0,sizeof(dis));
	memset(vd,0,sizeof(vd));
	memset(g,0,sizeof(g));
	for(i=1;i<=n;i++){
		int val=1;
		if(i==1||i==n)val=inf;
		g[i][i+n]=val;
		if(f[i]==1)g[op][i]=val;
	    if(f[i]==ans)g[i+n][ed]=val;
	}
	for(i=2;i<=n;i++)
		for(j=1;j<i;j++)
		    if(s[i]>=s[j]&&f[i]==f[j]+1)
		        g[j+n][i]=1;
	while(dis[op]<ed)maxflow+=dfs(op,inf);
	cout<<maxflow<<endl;
}