HDU 1569 체크 수 (2) (최대 점 권 독립 집합)

2801 단어
/*
  :
         ,          。
           ,i + j      s    ,              ,    t    ,              ,                 ,     INF。                  。
  :
1、        = sum -        。
2、        =     =    
  :dinic  
  :
http://wenku.baidu.com/view/98deaf06b52acfc789ebc91a.html###
http://www.cnblogs.com/ltang/archive/2010/11/17/1879573.html
*/

#include <iostream>
//#define re(i, n) for(int i = 0; i < n; ++ i)
using namespace std;

const int nMax = 2505;
const int INF = 0x7fffffff;
int queue[nMax];//            
int dis[nMax];//              
struct Edge
	//   ,  :    、    、         
{
	int v, w, next;
	Edge(){}
	Edge(int v, int w, int next):v(v), w(w), next(next){}
}adj[8 * nMax];
int V[nMax];//V[u]     u     , Edge    ,          
int cnt;
int s, t;

int min(int a, int b)
{
	return a < b ? a : b;
}

void add(int u, int v, int w)//        u - > v   
{
	adj[cnt] = Edge(v, w, V[u]);
	V[u] = cnt ++;
	adj[cnt] = Edge(u, 0, V[v]);
	V[v] = cnt ++;
}

int bfs()//    
{
	int front, rear;
	int v;
	memset(dis, 0, sizeof(dis));
	front = rear = 0;
	dis[s] = 1;
	queue[front ++] = s;
	while(rear < front)
	{
		int u = queue[rear ++];
		for(int i = V[u]; i != -1; i = adj[i].next)// u    
			if(adj[i].w && dis[v = adj[i].v] == 0)//      v         
			{
				dis[v] = dis[u] + 1;
				if(v == t) return 1;
				queue[front ++] = v;
			}
	}
	return 0;
}

int dfs(int u, int limit = INF)//   u   t,        
{
	if(u == t) return limit;
	int count = 0;
	for(int i = V[u]; i != -1; i = adj[i].next)// u     
	{
		int v = adj[i].v;
		if((dis[v] == dis[u] + 1) && adj[i].w)//       ,           
		{
			int z = dfs(v, min(limit - count, adj[i].w));
			if(z > 0)//          0, v t   
			{
				count += z;
				adj[i].w -= z;
				adj[i ^ 1].w += z;//  adj[i + 1] += z  ,    !
			}
			else
				dis[v] = -1;//        v      
		}
	}
	return count;
}

int dinic()
{
	int ans = 0;
	while(bfs())//            
		ans += dfs(s);
	return ans;
}

int main()
{
	//freopen("f://data.in", "r", stdin);
	int m, n;
	int sum;
	while(scanf("%d %d", &m, &n) != EOF)
	{
		cnt = 0;
		int x;
		sum = 0;
		s = 0;
		t = m * n + 1;
		memset(V, -1, sizeof(V));
		for(int i = 1; i <= m; ++ i)
			for(int j = 1; j <= n; ++ j)
			{
				scanf("%d", &x);
				sum += x;
				if((i + j) & 1)
				{
					add(s, (i - 1) * n + j, x);
					// 
					if(i > 1) add((i - 1) * n + j, (i - 2) * n + j, INF);
					// 
					if(i < m) add((i - 1) * n + j, i * n + j, INF);
					// 
					if(j > 1) add((i - 1) * n + j, (i -1) * n + j - 1, INF);
					// 
					if(j < n) add((i - 1) * n + j, (i - 1) * n + j + 1, INF);
				}
				else
					add((i - 1) * n + j, t, x);
			}
		printf("%d
",sum - dinic()); } return 0; }

좋은 웹페이지 즐겨찾기