国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發(fā)設計 > 正文

網絡流24題6. 最長遞增子序列問題

2019-11-06 08:19:40
字體:
來源:轉載
供稿:網友

最長遞增子序列問題

Description

給定正整數序列 x1,x2,x3,...,xn (1)計算其最長遞增子序列的長度 s。 (2)計算從給定的序列中最多可取出多少個長度為 s 的遞增子序列。 (3)如果允許在取出的序列中多次使用x1xn,則從給定序列中最多可取出多少個長度為 s 的遞增子序列。

Input

第 1 行有 1 個正整數 n,表示給定序列的長度。接下來的 1 行有 n 個正整數x1...xn

Output

第 1 行是最長遞增子序列的長度 s。第 2 行是可取出的長度為 s 的遞增子序列個數。第 3 行是允許在取出的序列中多次使用x1xn時可取出的長度為 s 的遞增子序列個數。

題解

題面冠冕堂皇的說著什么最長遞增子序列,數據竟給的是最長不上升子序列的,令人無語。 第一問可以用最簡單的O(n)的動態(tài)規(guī)劃求出,fi為以xi結束的最長遞增子序列長度,記最長遞增子序列長度為k。 下面考慮建圖,因為每個xi只能去一次,那么要限流,所以將每個xi拆成兩個點aibi,從aibi連一條容量為1的邊,就可以保證流過xi的流量最多為1。 然后從s點向每個fi=1ai連一條容量為1的邊,然后從每個fi=kbi向t連一條容量為1的邊。從滿足{i<j,xi<xj,fi+1=fj}biaj連一條容量為1的邊。這樣,一條從s到t的可行流就是一個最長遞增子序列,最大流即是答案。 對于第三問只需要將<a1,b1><an,bn><s,a1><bn,t>四條邊的容量賦為inf,再算一遍最大流即可。 【建模分析】 上述建模方法是應用了一種分層圖的思想,把圖每個頂點i按照fi的不同分為了若干層,這樣圖中從S出發(fā)到T的任何一條路徑都是一個滿足條件的最長上升子序列。由于序列中每個點要不可重復地取出,需要把每個點拆分成兩個點。單位網絡的最大流就是增廣路的條數,所以最大流量就是第二問結果。第三問特殊地要求x1xn可以重復使用,只需取消這兩個點相關邊的流量限制,求網絡最大流即可。

#include<cstdio>#include<cstring>#include<iostream>using namespace std;const int N = 1000 + 10, M = 2000000 + 10, inf = 0x3f3f3f3f;struct Edge{ int fr, to, cap, flow;}edg[M];int hd[N], nxt[M];int d[N], vis[N], q[N], dfn;int s, t;int n, ans, tot, k;int a[N], f[N];int e1, e2, e3, e4;void insert(int u, int v, int w){ edg[tot].fr = u, edg[tot].to = v, edg[tot].cap = w; nxt[tot] = hd[u], hd[u] = tot; if(u == 1 && v == n + 1) e1 = tot; if(u == n && v == n + n) e2 = tot; if(u == s && v == 1) e3 = tot; if(u == n + n && v == t) e4 = tot; tot++; edg[tot].fr = v, edg[tot].to = u; nxt[tot] = hd[v], hd[v] = tot; tot++;}bool bfs(){ int head = 1, tail = 1; q[1] = s; vis[s] = ++dfn; d[s] = 0; while(head <= tail){ int u = q[head++]; for(int i = hd[u]; i >= 0; i = nxt[i]){ Edge &e = edg[i]; if(vis[e.to] == dfn || e.cap <= e.flow) continue; vis[e.to] = dfn; d[e.to] = d[u] + 1; q[++tail] = e.to; } } return vis[t] == dfn;}int dfs(int x, int a){ if(x == t || a == 0) return a; int flow = 0, f; for(int i = hd[x]; i >= 0; i = nxt[i]){ Edge &e = edg[i]; if(d[e.to] == d[x] + 1 && (f = dfs(e.to, min(a, e.cap - e.flow))) > 0){ flow += f; e.flow += f; edg[i^1].flow -= f; a -= f; if(a == 0) break; } } return flow;}void dp(){ for(int i = 1; i <= n; i++){ for(int j = 1; j < i; j++) if(a[j] <= a[i] && f[j] > f[i]) f[i] = f[j]; f[i]++; k = max(k, f[i]); }}void build(){ s = 0, t = n * 2 + 1; int tmp[5]; for(int i = 1; i <= n; i++){ insert(i, i + n, 1); if(f[i] == 1) insert(s, i, 1); if(f[i] == k) insert(i + n, t, 1); } for(int i = 1; i <= n; i++) for(int j = i + 1; j <= n; j++) if(a[i] <= a[j] && f[i] + 1 == f[j]) insert(i + n, j, 1);}void init(){ memset(hd, -1, sizeof(hd)); scanf("%d", &n); for(int i = 1; i <= n; i++) scanf("%d", &a[i]);}void work(){ dp();
發(fā)表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發(fā)表
主站蜘蛛池模板: 酒泉市| 蓬安县| 内丘县| 宣威市| 柞水县| 和龙市| 察哈| 阿城市| 图木舒克市| 瓦房店市| 新巴尔虎左旗| 余庆县| 女性| 四川省| 遂昌县| 马边| 京山县| 沛县| 全椒县| 昔阳县| 和平县| 武穴市| 广宗县| 永修县| 水富县| 天台县| 金溪县| 红安县| 玉山县| 卓资县| 陈巴尔虎旗| 金堂县| 南雄市| 临夏县| 花垣县| 方正县| 金门县| 永仁县| 和林格尔县| 赤壁市| 蒙山县|