1. 程式人生 > >HDU 1565 方格取數(1)(最大點權獨立集)

HDU 1565 方格取數(1)(最大點權獨立集)

題目大意:給你一個n*n的格子的棋盤,每個格子裡面有一個非負數。
從中取出若干個數,使得任意的兩個數所在的格子沒有公共邊,就是說所取的數所在的2個格子不能相鄰,並且取出的數的和最大。

解題思路:最大點權獨立集,關鍵是怎麼建圖了,我們可以採用染色的思想對這張圖進行染色,然後分成兩個點集
假設將第一個格子染成白色,然後將它相鄰的格子染成相反的顏色黑色,以此類推,這樣就可以將一張圖分成染成黑白兩種顏色的點集了
然後就是連邊了,連邊的話,我們只考慮白色格子的連向黑色格子的,因為兩者之間是相對的,所以只需要取一條就好了
這樣圖就建好了
最大權值就是:總權值-最小割了(最小割就是最小點權覆蓋了)具體的證明請看

胡伯濤:演算法合集之《最小割模型在資訊學競賽中å的應用》

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;

#define M 1000010
#define N 10010
#define INF 0x3f3f3f3f
struct Edge{
    int u, v, cap, flow, next;
    Edge() {}
    Edge(int u, int v, int cap, int flow, int
next): u(u), v(v), cap(cap), flow(flow), next(next) {} }E[M]; struct Dinic{ int head[N], d[N]; int tot, sink, source; void init() { memset(head, -1, sizeof(head)); tot = 0; } inline void AddEdge(int u, int v, int cap) { E[tot] = Edge(u, v, cap, 0, head[u]); head[u] = tot++; u = u ^ v; v = u ^ v; u = u ^ v; E[tot] = Edge(u, v, 0
, 0, head[u]); head[u] = tot++; } inline bool bfs(int s) { int u, v; memset(d, 0, sizeof(d)); queue<int> Q; Q.push(s); d[s] = 1; while (!Q.empty()) { u = Q.front(); Q.pop(); if (u == sink) return true; for (int i = head[u]; ~i; i = E[i].next) { v = E[i].v; if (!d[v] && E[i].cap - E[i].flow > 0) { d[v] = d[u] + 1; Q.push(v); } } } return false; } int dfs(int x, int a) { if (x == sink || a == 0) return a; int f, flow = 0; for (int i = head[x]; ~i; i = E[i].next) { int v = E[i].v; if (d[v] == d[x] + 1 && E[i].cap - E[i].flow > 0) { f = dfs(v, min(a, E[i].cap - E[i].flow)); E[i].flow += f; E[i^1].flow -= f; flow += f; a -= f; if (!a) break; } } if (flow == 0) d[x] = 0; return flow; } int Maxflow(int source, int sink) { int flow = 0; this->sink = sink; while (bfs(source)) flow += dfs(source, INF); return flow; } }; Dinic dinic; #define maxn 25 int n, source, sink, Sum; int num[maxn][maxn]; bool mark[maxn][maxn]; int dir[4][2] = {{-1, 0}, {1, 0}, {0, -1}, {0, 1}}; void init() { source = 0, sink = n * n + 1, Sum = 0; dinic.init(); for (int i = 0; i < n; i++) for (int j = 1; j <= n; j++) { scanf("%d", &num[i][j]); Sum += num[i][j]; } memset(mark, 0, sizeof(mark)); for (int i = 0; i < n; i++) for (int j = 1; j <= n; j++) if (!mark[i][j]) { dinic.AddEdge(source, i * n + j, num[i][j]); for (int k = 0; k < 4; k++) { int x = i + dir[k][0], y = j + dir[k][1]; if (x < 0 || x >= n || y < 1 || y > n) continue; dinic.AddEdge(i * n + j, x * n + y, INF); if (!mark[x][y]) dinic.AddEdge(x * n + y, sink, num[x][y]); mark[x][y] = true; } } printf("%d\n", Sum - dinic.Maxflow(source, sink)); } int main() { while (scanf("%d", &n) != EOF) { init(); } return 0; }