題目描述

最近房地產(chǎn)商GDOI(Group of Dumbbells Or Idiots)從NOI(Nuts Old Idiots)手中得到了一塊開發(fā)土地。據(jù)了解，這塊土地是一塊矩形的區(qū)域，可以縱橫劃分為N×M塊小區(qū)域。GDOI要求將這些區(qū)域分為商業(yè)區(qū)和工業(yè)區(qū)來開發(fā)。根據(jù)不同的地形環(huán)境，每塊小區(qū)域建造商業(yè)區(qū)和工業(yè)區(qū)能取得不同的經(jīng)濟(jì)價(jià)值。更具體點(diǎn)，對于第i行第j列的區(qū)域，建造商業(yè)區(qū)將得到Aij收益，建造工業(yè)區(qū)將得到Bij收益。另外不同的區(qū)域連在一起可以得到額外的收益，即如果區(qū)域(I,j)相鄰（相鄰是指兩個(gè)格子有公共邊）有K塊（顯然K不超過4）類型不同于(I,j)的區(qū)域，則這塊區(qū)域能增加k×Cij收益。經(jīng)過Tiger.S教授的勘察，收益矩陣A,B,C都已經(jīng)知道了。你能幫GDOI求出一個(gè)收益最大的方案么？

輸入

輸入第一行為兩個(gè)整數(shù)，分別為正整數(shù)N和M，分別表示區(qū)域的行數(shù)和列數(shù)；第2到N+1列，每行M個(gè)整數(shù)，表示商業(yè)區(qū)收益矩陣A；第N+2到2N+1列，每行M個(gè)整數(shù)，表示工業(yè)區(qū)收益矩陣B；第2N+2到3N+1行，每行M個(gè)整數(shù)，表示相鄰額外收益矩陣C。第一行，兩個(gè)整數(shù)，分別是n和m（1≤n，m≤100）；

任何數(shù)字不超過1000”的限制

輸出

輸出只有一行，包含一個(gè)整數(shù)，為最大收益值。

樣例輸入

3 3
1 2 3
4 5 6
7 8 9
9 8 7
6 5 4
3 2 1
1 1 1
1 3 1
1 1 1

樣例輸出

81

---

題解

網(wǎng)絡(luò)流最小割

只考慮相鄰的兩個(gè)，問題轉(zhuǎn)化為：$i$和$j$各有兩種選法：選擇A可以獲得$a_i$或$a_j$的收益；選擇B可以獲得$b_i$或$b_j$的收益；如果選擇不同，則會獲得$c_i+c_j$的收益。問最大收益。

這是一個(gè)經(jīng)典的最小割模型，建圖方法：S連向i，容量為$a_i$，i連向T，容量為b_i；S連向j，容量為$b_j$，i連向T，容量為$a_j$（這兩步是反轉(zhuǎn)源匯的過程）。i和j之間連容量為$c_i+c_j$的雙向邊。

?

因此總的建圖為：黑白染色，黑點(diǎn)正常連，白點(diǎn)反轉(zhuǎn)源匯，然后相鄰的點(diǎn)之間連邊。答案為$\sum\limits a_i+\sum\limits b_i+\sum\limits(c_i+c_j)-mincut$。

#include <queue> #include <cstdio> #include <cstring> #define N 10010 #define M 1000010 #define pos(i , j) (i - 1) * m + j using namespace std; typedef long long ll; const int inf = 1 << 30; queue<int> q; int n , m , head[N] , to[M] , next[M] , cnt = 1 , s , t , dis[N]; ll a[110][110] , b[110][110] , c[110][110] , val[M]; void add(int x , int y , ll z) {to[++cnt] = y , val[cnt] = z , next[cnt] = head[x] , head[x] = cnt;to[++cnt] = x , val[cnt] = 0 , next[cnt] = head[y] , head[y] = cnt; } ll link(int x1 , int y1 , int x2 , int y2) {add(pos(x1 , y1) , pos(x2 , y2) , c[x1][y1] + c[x2][y2]);add(pos(x2 , y2) , pos(x1 , y1) , c[x1][y1] + c[x2][y2]);return c[x1][y1] + c[x2][y2]; } bool bfs() {int x , i;memset(dis , 0 , sizeof(dis));while(!q.empty()) q.pop();dis[s] = 1 , q.push(s);while(!q.empty()){x = q.front() , q.pop();for(i = head[x] ; i ; i = next[i]){if(val[i] && !dis[to[i]]){dis[to[i]] = dis[x] + 1;if(to[i] == t) return 1;q.push(to[i]);}}}return 0; } ll dinic(int x , ll low) {if(x == t) return low;ll temp = low , k;int i;for(i = head[x] ; i ; i = next[i]){if(val[i] && dis[to[i]] == dis[x] + 1){k = dinic(to[i] , min(temp , val[i]));if(!k) dis[to[i]] = 0;val[i] -= k , val[i ^ 1] += k;if(!(temp -= k)) break;}}return low - temp; } int main() {int i , j;ll ans = 0;scanf("%d%d" , &n , &m) , s = 0 , t = n * m + 1;for(i = 1 ; i <= n ; i ++ ) for(j = 1 ; j <= m ; j ++ ) scanf("%lld" , &a[i][j]);for(i = 1 ; i <= n ; i ++ ) for(j = 1 ; j <= m ; j ++ ) scanf("%lld" , &b[i][j]);for(i = 1 ; i <= n ; i ++ ) for(j = 1 ; j <= m ; j ++ ) scanf("%lld" , &c[i][j]);for(i = 1 ; i <= n ; i ++ ){for(j = 1 ; j <= m ; j ++ ){ans += a[i][j] + b[i][j];if((i & 1) ^ (j & 1)) add(s , pos(i , j) , b[i][j]) , add(pos(i , j) , t , a[i][j]);else{add(s , pos(i , j) , a[i][j]) , add(pos(i , j) , t , b[i][j]);if(i > 1) ans += link(i , j , i - 1 , j);if(i < n) ans += link(i , j , i + 1 , j);if(j > 1) ans += link(i , j , i , j - 1);if(j < m) ans += link(i , j , i , j + 1);}}}while(bfs()) ans -= dinic(s , inf);printf("%lld\n" , ans);return 0; }

?

?

轉(zhuǎn)載于:https://www.cnblogs.com/GXZlegend/p/7434706.html

總結(jié)

以上是生活随笔為你收集整理的【bzoj2132】圈地计划 网络流最小割的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。