「网络流24题」最长递增子序列问题
题意有些奇怪
任务(2)是取出。。。且题中递增是非严格递增
我的代码任务3若能取出无限多的序列,则输出-1
输入
4
1 3 2 4
输出
3 1 2
样例2
输入
4
3 6 2 5
输出
2 2 -1
搬运byvoid神犇题解
「问题分析」
第一问是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可以重复使用,只需取消这两个点相关边的流量限制,求网络最大流即可。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 |
#include<cstdio> #include<cmath> #include<ctime> #include<cstring> #include<iostream> #include<algorithm> #include<queue> #include<set> #define ll long long #define mod 1000000007 #define inf 1000000000 using namespace std; int read() { int f=1,x=0;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; } int n,K,cnt=1,T; int a[1005],f[1005],last[1005]; int h[1005],q[1005]; struct edge{ int next,v,to; }e[100005]; void insert(int u,int v,int w) { e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;e[cnt].v=w; e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt;e[cnt].v=0; } bool bfs() { int head=0,tail=1; memset(h,-1,sizeof(h)); q[0]=0;h[0]=0; while(head!=tail) { int now=q[head];head++; for(int i=last[now];i;i=e[i].next) if(e[i].v&&h[e[i].to]==-1) { h[e[i].to]=h[now]+1; q[tail++]=e[i].to; } } return h[T]!=-1; } int dfs(int x,int f) { if(x==T)return f; int w,used=0; for(int i=last[x];i;i=e[i].next) if(h[e[i].to]==h[x]+1) { w=dfs(e[i].to,min(f-used,e[i].v)); e[i].v-=w;e[i^1].v+=w; used+=w;if(used==f)return f; } if(!used)h[x]=-1; return used; } int dinic() { int tmp=0; while(bfs())tmp+=dfs(0,inf); return tmp; } void solve1() { T=2*n+1; for(int i=1;i<=n;i++) { if(f[i]==1)insert(0,i,1); if(f[i]==K)insert(i+n,T,1); insert(i,i+n,1); } for(int i=1;i<=n;i++) for(int j=i+1;j<=n;j++) if(a[j]>=a[i]&&f[j]==f[i]+1) insert(i+n,j,1); printf("%d\n",dinic()); } void solve2() { cnt=1;memset(last,0,sizeof(last)); for(int i=1;i<=n;i++) { int v=1; if(i==1||i==n)v=inf; if(f[i]==1)insert(0,i,v); if(f[i]==K)insert(i+n,T,v); insert(i,i+n,v); } for(int i=1;i<=n;i++) for(int j=i+1;j<=n;j++) if(a[j]>=a[i]&&f[j]==f[i]+1) insert(i+n,j,1); int tmp=dinic(); if(tmp>=inf)puts("-1"); else printf("%d\n",tmp); } int main() { n=read(); for(int i=1;i<=n;i++)a[i]=read(); for(int i=1;i<=n;i++) { f[i]=1; for(int j=1;j<i;j++) if(a[j]<=a[i])f[i]=max(f[i],f[j]+1); K=max(K,f[i]); } printf("%d\n",K); solve1(); solve2(); return 0; } |
黄学长,好像你题解里面说的和你连边的方式不太一样。
应将solve2()中第二个添边部分中将i+n和j所连得边的流量修改为1而不是inf
恩是我写错了
然而对于
4
3 6 2 5
这个数据并没有输出正解
正解:
2
2
3
程序输出:
2
2
-1
改了一句(不然80分)
if(tmp>=inf) printf(“%dn”, n); // puts(“-1”);