731. [网络流24题] 最长递增子序列
★★★☆ 输入文件:alis.in
输出文件:
alis.out
简单对比
时间限制:1 s 内存限制:128 MB
«问题描述:
给定正整数序列x1,..., xn。
(1)计算其最长递增子序列的长度s。
(2)计算从给定的序列中最多可取出多少个长度为s的递增子序列。
(3)如果允许在取出的序列中多次使用x1和xn,则从给定序列中最多可取出多少个长
度为s的递增子序列。
注意:这里的最长递增子序列即最长不下降子序列!!!
«编程任务:
设计有效算法完成(1)(2)(3)提出的计算任务。
«数据输入:
由文件alis.in提供输入数据。文件第1 行有1个正整数n(n<=500),表示给定序列的长度。接
下来的1 行有n个正整数x1,..., xn。
«结果输出:
程序运行结束时,将任务(1)(2)(3)的解答输出到文件alis.out中。第1 行是最长
递增子序列的长度s。第2行是可取出的长度为s 的递增子序列个数。第3行是允许在取出
的序列中多次使用x1和xn时可取出的长度为s 的递增子序列个数。
输入文件示例 输出文件示例
alis.in
4
3 6 2 5
alis.out
22
3
【问题分析】
第一问是LIS,动态规划求解,第二问和第三问用网络最大流解决。
【建模方法】
首先动态规划求出F[i],表示以第i位为开头的最长上升序列的长度,求出最长上升序列长度K。
1、把序列每位i拆成两个点<i.a>和<i.b>,从<i.a>到<i.b>连接一条容量为1的有向边。
2、建立附加源S和汇T,如果序列第i位有F[i]=K,从S到<i.a>连接一条容量为1的有向边。
3、如果F[i]=1,从<i.b>到T连接一条容量为1的有向边。
4、如果j>i且A[i] < A[j]且F[j]+1=F[i],从<i.b>到<j.a>连接一条容量为1的有向边。
求网络最大流,就是第二问的结果。把边(<1.a>,<1.b>)(<N.a>,<N.b>)(S,<1.a>)(<N.b>,T)这四条边的容量修改为无穷大,再求一次网络最大流,就是第三问结果。
【建模分析】
上述建模方法是应用了一种分层图的思想,把图每个顶点i按照F[i]的不同分为了若干层,这样图中从S出发到T的任何一条路径都是一个满足条件的最长上升子序列。由于序列中每个点要不可重复地取出,需要把每个点拆分成两个点。单位网络的最大流就是增广路的条数,所以最大流量就是第二问结果。第三问特殊地要求x1和xn可以重复使用,只需取消这两个点相关边的流量限制,求网络最大流即可。
转载自hzwer
ps:任务3若能取出无限多的序列,则输出任务2的答案
- //数据略坑
- #include<cstdio>
- #include<cstring>
- #include<algorithm>
- #define setfile(name) freopen(#name".in","r",stdin);freopen(#name".out","w",stdout);
- using namespace std;
- const int N=1005;
- const int inf=1e9;
- struct edge{int v,cap,next;}e[N*N*2];int tot=1,head[N];
- int n,S,T,K,ans,a[N],f[N],dis[N],q[N*N*2];
- inline int read(){
- int x=0,f=1;char ch=getchar();
- 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;
- }
- inline void add(int x,int y,int z){
- e[++tot].v=y;e[tot].cap=z;e[tot].next=head[x];head[x]=tot;
- e[++tot].v=x;e[tot].cap=0;e[tot].next=head[y];head[y]=tot;
- }
- inline bool bfs(){
- memset(dis,-1,sizeof dis);
- unsigned short h=0,t=1;q[t]=S;dis[S]=0;
- while(h!=t){
- int x=q[++h];
- for(int i=head[x];i;i=e[i].next){
- if(e[i].cap&&dis[e[i].v]==-1){
- dis[e[i].v]=dis[x]+1;
- if(e[i].v==T) return 1;
- q[++t]=e[i].v;
- }
- }
- }
- return 0;
- }
- int dfs(int x,int f){
- if(x==T) return f;
- int used=0,t;
- for(int i=head[x];i;i=e[i].next){
- if(e[i].cap&&dis[e[i].v]==dis[x]+1){
- t=dfs(e[i].v,min(e[i].cap,f));
- e[i].cap-=t;e[i^1].cap+=t;
- used+=t;f-=t;
- if(!f) return used;
- }
- }
- if(!used) dis[x]=-1;
- return used;
- }
- inline int dinic(){
- int res=0;
- while(bfs()) res+=dfs(S,inf);
- return res;
- }
- inline void DP(){
- for(int i=1;i<=n;i++){
- f[i]=1;
- for(int j=1;j<i;j++){
- if(a[i]>=a[j]){
- f[i]=max(f[i],f[j]+1);
- }
- }
- }
- K=*max_element(f+1,f+n+1);
- }
- void init(){
- n=read();
- for(int i=1;i<=n;i++) a[i]=read();
- }
- void ord_mapping(){
- tot=1;memset(head,0,sizeof head);
- for(int i=1;i<=n;i++){
- for(int j=1;j<i;j++){
- if(a[j]<=a[i]&&f[j]+1==f[i]){
- add(j+n,i,1);
- }
- }
- }
- }
- inline void work(){
- DP();printf("%d\n",K);
- S=0,T=n<<1|1;
- ord_mapping();
- for(int i=1;i<=n;i++){
- if(f[i]==1) add(S,i,1);
- if(f[i]==K) add(i+n,T,1);
- add(i,i+n,1);
- }
- ans=dinic();
- printf("%d\n",ans);
- ord_mapping();
- for(int i=1;i<=n;i++){
- int w=1;
- if(i==1||i==n) w=inf;
- if(f[i]==1) add(S,i,w);
- if(f[i]==K) add(i+n,T,w);
- add(i,i+n,w);
- }
- int lastans=ans;
- ans=dinic();
- if(ans>inf) printf("%d\n",lastans);
- else printf("%d\n",ans);
- }
- int main(){
- setfile(alis);
- init();
- work();
- return 0;
- }