【題意】給出一個(gè)無(wú)向圖，每個(gè)點(diǎn)有一個(gè)標(biāo)號(hào)mark[i]，不同點(diǎn)可能有相同的標(biāo)號(hào)。對(duì)于一條邊(u, v)，它的權(quán)值定義為mark[u] xor mark[v]。現(xiàn)在一些點(diǎn)的標(biāo)號(hào)已定，請(qǐng)決定剩下點(diǎn)的標(biāo)號(hào)，使得總的邊權(quán)和最小。（0 < N <= 500, 0 <= M <= 3000, 0 <=?mark[i] <= 2^31-1） 胡伯濤神牛《最小割模型在信息學(xué)競(jìng)賽中的應(yīng)用》中的例題。非常好的一道題！非常推薦！ 【思路】 我們把問(wèn)題數(shù)學(xué)化就是： ?Minimum ?sigma(w_e) = sigma_{(u, v)∈E?}( mark(u) xor mark(v) ) 對(duì)于異或問(wèn)題，我們發(fā)現(xiàn)這樣的二進(jìn)制按位運(yùn)算各個(gè)二進(jìn)制位之間是互不影響的，所以我們可以一位一位的做這類題。 那么我們的式子又可以進(jìn)一步轉(zhuǎn)化為： Minimum ?sigma_{(u, v)∈E}?{ sigma_i=0~oo(2^i)?? sigma(mark(u, i) xor mark(v, i))?} 這樣我們就把mark的限制加強(qiáng)了：只可能是0或1。即這些點(diǎn)將分成兩類。 再觀察我們發(fā)現(xiàn)，xor運(yùn)算，只有當(dāng)u、v不同時(shí)結(jié)果才為1，即這樣的有效邊的兩端點(diǎn)一定屬于不同點(diǎn)集。這像什么？不就是割邊嘛！~而題目正好又是要求最小，這樣問(wèn)題便轉(zhuǎn)化為最小割了~ ? ?（要注意培養(yǎng)這種問(wèn)題轉(zhuǎn)化和模型發(fā)現(xiàn)的能力！） 那么具體的最小割網(wǎng)絡(luò)G_N模型：建一個(gè)源點(diǎn)，向每一個(gè)標(biāo)號(hào)為1的點(diǎn)連一條oo流量的邊（后面解釋為什么源點(diǎn)連標(biāo)號(hào)1的點(diǎn)）；建一個(gè)匯點(diǎn)，向每一個(gè)標(biāo)號(hào)為0的點(diǎn)連一條oo流量的邊；原圖中的邊容量設(shè)為1加入到G_N中。求出來(lái)的最小割便是該二進(jìn)制位下的標(biāo)號(hào)xor的和最小的情況。 然而題目還要求輸出所有點(diǎn)的標(biāo)號(hào)，并且需要標(biāo)號(hào)的和也最小。那么怎么保證標(biāo)號(hào)的和最小呢？無(wú)非就是盡可能的取0。那么又該怎么做？ 首先先看怎么給那些未標(biāo)號(hào)的點(diǎn)標(biāo)號(hào)：容易想到最小割把網(wǎng)絡(luò)分成了兩個(gè)點(diǎn)集，那么顯然每個(gè)點(diǎn)標(biāo)號(hào)應(yīng)該和它所在點(diǎn)集已標(biāo)號(hào)的點(diǎn)一致，所以當(dāng)然希望標(biāo)號(hào)為0的點(diǎn)集點(diǎn)更多一些。然后注意我們劃分點(diǎn)集是從源點(diǎn)開(kāi)始dfs，那么這樣劃出來(lái)的最小割邊集顯然更偏向源點(diǎn)，即這樣劃分出來(lái)的S集點(diǎn)是最少的。于是源點(diǎn)當(dāng)然連標(biāo)號(hào)為1的點(diǎn)吶~ 【代碼】 ? #include #include #include #include #include #include #define MID(x,y) ((x+y)/2) #define mem(a,b) memset(a,b,sizeof(a)) using namespace std; const int MAXV = 505; const int MAXE = 10005; const int oo = 0x3fffffff;/* Dinic-2.0-2013.07.21: adds template. double & int 轉(zhuǎn)換方便多了，也不易出錯(cuò) ~*/ template struct Dinic{struct node{int u, v;T flow;int opp;int next;}arc[2*MAXE];int vn, en, head[MAXV];int cur[MAXV];int q[MAXV];int path[2*MAXE], top;int dep[MAXV];void init(int n){vn = n;en = 0;mem(head, -1);}void insert_flow(int u, int v, T flow){arc[en].u = u;arc[en].v = v;arc[en].flow = flow;arc[en].next = head[u];head[u] = en ++;arc[en].u = v;arc[en].v = u;arc[en].flow = 0;arc[en].next = head[v];head[v] = en ++;}bool bfs(int s, int t){mem(dep, -1);int lq = 0, rq = 1;dep[s] = 0;q[lq] = s;while(lq < rq){int u = q[lq ++];if (u == t){return true;}for (int i = head[u]; i != -1; i = arc[i].next){int v = arc[i].v;if (dep[v] == -1 && arc[i].flow > 0){dep[v] = dep[u] + 1;q[rq ++] = v;}}}return false;}T solve(int s, int t){T maxflow = 0;while(bfs(s, t)){int i, j;for (i = 1; i <= vn; i ++) cur[i] = head[i];for (i = s, top = 0;;){if (i == t){int mink;T minflow = 0x3fffffff;for (int k = 0; k < top; k ++)if (minflow > arc[path[k]].flow){minflow = arc[path[k]].flow;mink = k;}for (int k = 0; k < top; k ++)arc[path[k]].flow -= minflow, arc[path[k]^1].flow += minflow;maxflow += minflow;top = mink;i = arc[path[top]].u;}for (j = cur[i]; j != -1; cur[i] = j = arc[j].next){int v = arc[j].v;if (arc[j].flow && dep[v] == dep[i] + 1)break;}if (j != -1){path[top ++] = j;i = arc[j].v;}else{if (top == 0) break;dep[i] = -1;i = arc[path[-- top]].u;}}}return maxflow;} }; Dinic dinic; int mark[MAXV]; bool if_mark[MAXV]; struct path{int u, v; }p[MAXE]; bool vis[MAXV]; int st[MAXV]; //ST集 void dfs(int u){vis[u] = 1;st[u] = 1;for (int i = dinic.head[u]; i != -1; i = dinic.arc[i].next){if (dinic.arc[i].flow <= 0) continue;int v = dinic.arc[i].v;if (!vis[v]){dfs(v);}}return ; } int main(){int t;scanf("%d", &t);while(t --){int n, m;scanf("%d %d", &n, &m);for (int i = 1; i <= m; i ++){scanf("%d %d", &p[i].u, &p[i].v);}int k;mem(mark, 0);mem(if_mark, false);scanf("%d", &k);int maxn = 0;for (int i = 0; i < k; i ++){int u;scanf("%d", &u);scanf("%d", &mark[u]);maxn = max(maxn, mark[u]);if_mark[u] = true;}int oi = ceil(log(maxn)/log(2));for (int k = 0; k < oi; k ++){dinic.init(n+2);for (int i = 1; i <= n; i ++){if (!if_mark[i])continue;if ((mark[i] & (1 << k))){dinic.insert_flow(n+1, i, oo);}else{dinic.insert_flow(i, n+2, oo);}}for (int i = 1; i <= m; i ++){dinic.insert_flow(p[i].u, p[i].v, 1);dinic.insert_flow(p[i].v, p[i].u, 1);}dinic.solve(n+1, n+2);mem(st, 0);mem(vis, 0);dfs(n+1); //殘留網(wǎng)絡(luò)中dfs確定點(diǎn)S、T集for (int i = 1; i <= n; i ++){if (st[i] == 1 && !if_mark[i]){mark[i] += (1 << k);}}}for (int i = 1; i <= n; i ++){printf("%d\n", mark[i]);}}return 0; }

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

總結(jié)

以上是生活随笔為你收集整理的SPOJ-OPTM Optimal Marks ★★(按位建图 最小割)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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