http://www.lydsy.com/JudgeOnline/problem.php?id=4873

最大權閉合子圖。。。

建圖：

1.d[i][j]:i->j區間的費用，d[i][j] > 0 ins(S,id(i,j),d[i][j]) 否則ins(id(i,j),T,-d[i][j]) 套路

2.對于壽司怎么搞，m=1，ins(種類,T,a[i]*a[i]),ins(壽司，種類,inf):必須割掉初始的費用,ins(壽司，T，a[i]),ins(區間,壽司, inf):每個區間割掉需要壽司的花費

3.ins(id(i,j),id(i+1,j)，inf),ins(id(i,j),id(i,j-1)，inf)：選了大的區間的必須選小的區間

但是上面有一步可以改進，就是2的最后。因為選了大的一定會選小的，那么我們只用將[i,i]這個區間向壽司連邊就行了。

記住a[i]有1000，并且匯點不要取太小。。。

#include<bits/stdc++.h> using namespace std; const int N = 100010, inf = 1 << 29; struct edge {int nxt, to, f; } e[N * 2]; int dis[N], used[N], head[N], q[N], iter[N], d[110][110], a[N]; int n, m, sum, T = 0, num = 0, cnt = 1; namespace maxflow {void link(int u, int v, int f){e[++cnt].nxt = head[u];head[u] = cnt;e[cnt].to = v;e[cnt].f = f;}void ins(int u, int v, int f){link(u, v, f); link(v, u, 0);}bool bfs(){int l = 1, r = 0; q[++r] = 0;memset(dis, 0, sizeof(dis)); dis[0] = 1;while(l <= r){int u = q[l++];for(int i = head[u]; i; i = e[i].nxt) if(!dis[e[i].to] && e[i].f){dis[e[i].to] = dis[u] + 1;q[++r] = e[i].to;}}return dis[T] > 0;}int dfs(int u, int delta){if(u == T) return delta;int ret = 0;for(int &i = iter[u]; i && delta; i = e[i].nxt) if(e[i].f && dis[e[i].to] == dis[u] + 1){int x = dfs(e[i].to, min(delta, e[i].f));e[i].f -= x; e[i ^ 1].f += x;ret += x; delta -= x;}return ret;}int id(int i, int j) { return (i - 1) * n + j; }void build(){//每個編號和T連邊，每個壽司和對應編號連邊 int D = n * n;T = N - 2;for(int i = 1; i <= n; ++i){ // i + D:壽司 a[i] + 2 * D: 種類 id(i, i): 區間 ins(i + D, T, a[i]); //每個壽司 if(d[i][i] < 0) ins(id(i, i), T, -d[i][i]);else ins(0, id(i, i), d[i][i]);ins(id(i, i), i + D, inf);if(!m) continue;if(!used[a[i]]){used[a[i]] = 1;ins(a[i] + 2 * D, T, a[i] * a[i]);}ins(i + D, a[i] + 2 * D, inf); } for(int i = 1; i <= n; ++i) for(int j = i + 1; j <= n; ++j){if(d[i][j] < 0) ins(id(i, j), T, -d[i][j]);else ins(0, id(i, j), d[i][j]); if(i < n) ins(id(i, j), id(i + 1, j), inf);if(j > 1) ins(id(i, j), id(i, j - 1), inf);} }int dinic(){int ret = 0;while(bfs()){for(int i = 0; i <= T; ++i) iter[i] = head[i];ret += dfs(0, inf);}return ret;} } using namespace maxflow; int main() {scanf("%d%d", &n, &m);for(int i = 1; i <= n; ++i) scanf("%d", &a[i]);for(int i = 1; i <= n; ++i)for(int j = i; j <= n; ++j) {scanf("%d", &d[i][j]);if(d[i][j] > 0) sum += d[i][j];}build();sum -= dinic();printf("%d\n", sum); return 0; } View Code

?

轉載于:https://www.cnblogs.com/19992147orz/p/6833105.html

總結

以上是生活随笔為你收集整理的bzoj4873的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。