标签：event stream include bsp 路径 重复 可重复 das 最长上升子序列
题解来自网络流24题：
【问题分析】
第一问时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可以重复使用，只需取消这两个点相关边的流量限制，求网络最大流即可。
还有这个题题意有些问题，不是递增，是不递减。
还有我这个题没有拆点在洛谷和codevs上都过了，究竟需不需要拆点啊。。
——代码（木有拆点）
1 #include <queue> 2 #include <cstdio> 3 #include <cstring> 4 #include <iostream> 5 #define N 1010 6 #define M 3000001 7 #define max(x, y) ((x) > (y) ? (x) : (y)) 8 #define min(x, y) ((x) < (y) ? (x) : (y)) 9 10 int n, ans, cnt, s, t, sum; 11 int a[N], f[N]; 12 int head[N], to[M], val[M], next[M], dis[N], cur[N]; 13 14 inline int read() 15 { 16 int x = 0, f = 1; 17 char ch = getchar(); 18 for(; !isdigit(ch); ch = getchar()) if(ch == ‘-‘) f = -1; 19 for(; isdigit(ch); ch = getchar()) x = (x << 1) + (x << 3) + ch - ‘0‘; 20 return x * f; 21 } 22 23 inline void add(int x, int y, int z) 24 { 25 to[cnt] = y; 26 val[cnt] = z; 27 next[cnt] = head[x]; 28 head[x] = cnt++; 29 } 30 31 inline bool bfs() 32 { 33 int i, u, v; 34 std::queue <int> q; 35 memset(dis, -1, sizeof(dis)); 36 q.push(s); 37 dis[s] = 0; 38 while(!q.empty()) 39 { 40 u = q.front(), q.pop(); 41 for(i = head[u]; i ^ -1; i = next[i]) 42 { 43 v = to[i]; 44 if(val[i] && dis[v] == -1) 45 { 46 dis[v] = dis[u] + 1; 47 if(v == t) return 1; 48 q.push(v); 49 } 50 } 51 } 52 return 0; 53 } 54 55 inline int dfs(int u, int maxflow) 56 { 57 if(u == t) return maxflow; 58 int i, v, d, ret = 0; 59 for(i = cur[u]; i ^ -1; i = next[i]) 60 { 61 v = to[i]; 62 if(val[i] && dis[v] == dis[u] + 1) 63 { 64 d = dfs(v, min(val[i], maxflow)); 65 ret += d; 66 cur[u] = i; 67 val[i] -= d; 68 val[i ^ 1] += d; 69 if(ret == maxflow) return ret; 70 } 71 } 72 return ret; 73 } 74 75 inline void clear() 76 { 77 int i, j; 78 sum = cnt = 0; 79 memset(head, -1, sizeof(head)); 80 for(i = 1; i <= n; i++) 81 { 82 if(f[i] == 1) add(s, i, 1), add(i, s, 0); 83 if(f[i] == ans) add(i, t, 1), add(t, i, 0); 84 } 85 for(i = 1; i <= n; i++) 86 for(j = 1; j < i; j++) 87 if(a[j] <= a[i] && f[j] + 1 == f[i]) 88 add(j, i, 1), add(i, j, 0); 89 } 90 91 int main() 92 { 93 int i, j, x; 94 n = read(); 95 s = 0, t = n + 1; 96 for(i = 1; i <= n; i++) 97 { 98 a[i] = read(); 99 x = 0; 100 for(j = 1; j < i; j++) 101 if(a[j] <= a[i]) 102 x = max(x, f[j]); 103 f[i] = x + 1; 104 ans = max(ans, f[i]); 105 } 106 printf("%d\n", ans); 107 clear(); 108 while(bfs()) 109 { 110 for(i = s; i <= t; i++) cur[i] = head[i]; 111 sum += dfs(s, 1e9); 112 } 113 printf("%d\n", sum); 114 clear(); 115 add(s, 1, 1e9), add(1, s, 0); 116 if(f[n] == ans) add(n, t, 1e9), add(t, n, 0); 117 while(bfs()) 118 { 119 for(i = s; i <= t; i++) cur[i] = head[i]; 120 sum += dfs(s, 1e9); 121 } 122 printf("%d\n", sum); 123 return 0; 124 }
标签：event stream include bsp 路径 重复 可重复 das 最长上升子序列
原文地址：http://www.cnblogs.com/zhenghaotian/p/6933397.html
