题意:一个方格n*m,取出一些点,要求两两不相邻,求最大和。思路:建图,相邻的点有一条边,则建立了一个二分图,求最大点权独立集(所取点两两无公共边,权值和最大),问题转化为求总权和-最小点权覆盖集(点集I覆盖所有边,点权之和最小),(对应于原题,就是求拿掉最小点集,这些点覆盖所有边,拿掉后,每个点必然两两不相邻,否则:假设u,v相邻,则u->v这条边未被覆盖,矛盾),在建立超级源汇点s,t,s连向所有X中的点(设二分图G(X,Y)),Y联向t,,权值为点权,原来X->Y的所有边权值改为inf,问题转化为:求s->t最小割(一组权值和最小的割边集,去掉后s->t不连通),而每去掉一条割边,相当于去掉原图一个点,这个点必然牵着下面X->Y的边,故最小割即为最小点权覆盖集!(部分证明:摘自某大牛:可以这样理解:X到Y的边权为INF,自然不会成为最小割中的边,那就只有可能

是S到X和Y到T中的边,而:S到X中点x的边e1, 权为点x的点权,点x和Y中的所有临边e2,都需要受

到e1的流量的限制,同样,X到Y中点y的所有边也会受到点y到T的容量限制。这样求得割就能保证覆

盖掉所有的边。

我们可以用反证法证明一下:假设有边<x, y>没有被覆盖掉,则边<S, x>流量为0且边<y, T>流量为0,

而<x, y>流量为INF,自然可以找到一条S到T的增流路径<S, x, y, T>,与以求得流为最大流相矛盾,

则可以说明,在最大流的情况下,所有的边都已经被覆盖掉。)

建图这里注意一下:二分图,只有X->Y有连边!!!

结论(二分图):最小点权覆盖集=最小割=最大流; 最大点权覆盖集=总权和-最小点权覆盖集

PS:网络流博大精深。。。。。

#include<iostream>  //15ms
#include<queue>
#include<cstdio>
using namespace std;
int n,m,nume;const int inf=0x3f3f3f3f; int e[40000][3];
int level[2509];int vis[2509]; int head[2509];
bool bfs() //dinic
{
for(int i=0;i<=n*m+1;i++)
vis[i]=level[i]=0;
queue<int>q; q.push(0); vis[0]=1;
while(!q.empty())
{
int cur=q.front(); q.pop();
for(int i=head[cur];i!=-1;i=e[i][1])
{ int v=e[i][0];
if(!vis[v]&&e[i][2]>0)
{
level[v]=level[cur]+1;
if(v==m*n+1)return 1;
vis[v]=1;
q.push(v);
}
}
}
return vis[n*m+1];
}
int dfs(int u,int minf)
{
if(u==n*m+1||minf==0)return minf;
int sumf=0,f;
for(int i=head[u];i!=-1&&minf;i=e[i][1])
{ int v=e[i][0];
if(level[v]==level[u]+1&&e[i][2]>0)
{
f=dfs(v,minf<e[i][2]?minf:e[i][2]);
e[i][2]-=f; e[i^1][2]+=f;
minf-=f; sumf+=f;
}
}
return sumf;
}
int dinic()
{
int sum=0;
while(bfs())
{
sum+=dfs(0,inf);
}
return sum;
}
int a[52][52];
void addegde(int u,int v,int w)
{
e[nume][0]=v;e[nume][1]=head[u];head[u]=nume;
e[nume++][2]=w;
e[nume][0]=u;e[nume][1]=head[v];head[v]=nume;
e[nume++][2]=0;
}
int main()
{
while(scanf("%d",&n,&m)!=EOF)
{
for(int i=0;i<=n*m+1;i++)
head[i]=-1;
int temp,sum=0; nume=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
scanf("%d",&temp);
a[i][j]=temp; sum+=temp;
if((i+j)%2)addegde(0,(i-1)*m+j,temp);
else addegde((i-1)*m+j,n*m+1,temp);
}
for(int i=1;i<=n;i++) //建图这里注意一下:二分图,只有X->Y有连边!!!
for(int j=1;j<=m;j++)
{
if((i+j)%2==0)continue;
if(i-1>=1)
addegde((i-1)*m+j,(i-2)*m+j,inf);
if(j+1<=m)
addegde((i-1)*m+j,(i-1)*m+j+1,inf);
if(i+1<=n)
addegde((i-1)*m+j,i*m+j,inf);
if(j-1>=1)
addegde((i-1)*m+j,(i-1)*m+j-1,inf);
}
int ans=dinic();
printf("%d\n",sum-ans);
}
return 0;
}

hdu1569 方格取数 求最大点权独立集的更多相关文章

  1. hdu1569 方格取数(2) 最大点权独立集=总权和-最小点权覆盖集 (最小点权覆盖集=最小割=最大流)

    /** 转自:http://blog.csdn.net/u011498819/article/details/20772147 题目:hdu1569 方格取数(2) 链接:https://vjudge ...

  2. HDU 1565 1569 方格取数(最大点权独立集)

    HDU 1565 1569 方格取数(最大点权独立集) 题目链接 题意:中文题 思路:最大点权独立集 = 总权值 - 最小割 = 总权值 - 最大流 那么原图周围不能连边,那么就能够分成黑白棋盘.源点 ...

  3. HDU 1569 - 方格取数(2) - [最大点权独立集与最小点权覆盖集]

    嗯,这是关于最大点权独立集与最小点权覆盖集的姿势,很简单对吧,然后开始看题. 题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1569 Time Limi ...

  4. TZOJ 3665 方格取数(2)(最大点权独立集)

    描述 给你一个m*n的格子的棋盘,每个格子里面有一个非负数. 从中取出若干个数,使得任意的两个数所在的格子没有公共边,就是说所取数所在的2个格子不能相邻,并且取出的数的和最大. 输入 包括多个测试实例 ...

  5. hdu - 1565 方格取数(1) && 1569 方格取数(2) (最大点权独立集)

    http://acm.hdu.edu.cn/showproblem.php?pid=1565 两道题只是数据范围不同,都是求的最大点权独立集. 我们可以把下标之和为奇数的分成一个集合,把下标之和为偶数 ...

  6. HDU1569 方格取数(2) —— 二分图点带权最大独立集、最小割最大流

    题目链接:https://vjudge.net/problem/HDU-1569 方格取数(2) Time Limit: 10000/5000 MS (Java/Others)    Memory L ...

  7. hdu1565+hdu1569(最大点权独立集)

    传送门:hdu1565 方格取数(1) 传送门:hdu1569 方格取数(2) 定理:1. 最小点权覆盖集=最小割=最大流2. 最大点权独立集=总权-最小点权覆盖集 步骤: 1. 先染色,取一个点染白 ...

  8. HDU1565 方格取数(1) —— 状压DP or 插头DP(轮廓线更新) or 二分图点带权最大独立集(最小割最大流)

    题目链接:https://vjudge.net/problem/HDU-1565 方格取数(1) Time Limit: 10000/5000 MS (Java/Others)    Memory L ...

  9. HDU 1565 最大点权独立集

    首先要明白图论的几个定义: 点覆盖.最小点覆盖: 点覆盖集即一个点集,使得所有边至少有一个端点在集合里.或者说是“点” 覆盖了所有“边”.. 最小点覆盖(minimum vertex covering ...

随机推荐

  1. No package python-pip available

    因为没有此rpm包,此包包含在epel源里面 输入rpm -ivh http://dl.fedoraproject.org/pub/epel/6/i386/epel-release-6-8.noarc ...

  2. Alfred的配置和使用

    http://www.jianshu.com/p/f77ad047f7b0   Alfred:mac上的神兵利器,提升工作效率*n,快捷键:option + 空格.鉴于是看了池老师的<人生元编程 ...

  3. css--背景和列表

    背景 背景样式: background-color    设置元素的背景颜色 background-image  把图像设置为背景 background-repeat     设置背景图像是否重复及如 ...

  4. Codeforces 727C Guess the Array

    题目传送门 长度为\(n\)的数组,询问\(n\)次来求出数组中每一个数的值 我们可以先询问三次 \(a[1]+a[2]\) \(a[1]+a[3]\) \(a[2]+a[3]\) 然后根据这三次询问 ...

  5. 使用一位数组解决 1 1 2 3 5 8 13 数列问题 斐波纳契数列 Fibonacci

    斐波纳契数列 Fibonacci 输出这个数列的前20个数是什么? 1 1 2 3 5 8 13 21 34 55 89 144 233 377 610 987 1597 使用数组实现输出数列的前30 ...

  6. 【贪心】「poj1328」Radar Installation

    建模:二维转一维:贪心 Description Assume the coasting is an infinite straight line. Land is in one side of coa ...

  7. 【树状数组 思维题】luoguP3616 富金森林公园

    树状数组.差分.前缀和.离散化 题目描述 博艾的富金森林公园里有一个长长的富金山脉,山脉是由一块块巨石并列构成的,编号从1到N.每一个巨石有一个海拔高度.而这个山脉又在一个盆地中,盆地里可能会积水,积 ...

  8. (58)zabbix网络拓扑图配置network map

    zabbix网络地图介绍 “zabbix network map”可以简单的理解为动态网络拓扑图,可以针对业务来配置zabbix map, 通过map可以了解应用的整体状况:服务器是否异常.网络是否有 ...

  9. 自动化运维工具Ansible

    一.简介 当下有许多的运维自动化工具( 配置管理 ),例如:Ansible.SaltStack.Puppet.Fabric 等. Ansible 一种集成 IT 系统的配置管理.应用部署.执行特定任务 ...

  10. '>>' should be '> >' within a nested template argument list

    在编译关于opencv相机标定的工程的时候出现了这个问题 vector<vector<Point3f>>  objectPoints;  error: 'objectPoint ...