博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ 4808 二分图最大独立集
阅读量:5223 次
发布时间:2019-06-14

本文共 1924 字,大约阅读时间需要 6 分钟。

思路:

棋盘是个二分图

那就把一个可以走的白点  向所有可以走的黑点连边

跑一个最大匹配   (匹配上了就代表这两个点不能共存)

最大独立集=sum-最大匹配

//By SiriusRen#include 
#include
#include
#include
using namespace std;const int N=222,M=N*N*21;int n,m,a[N][N],color[N][N],id[N][N],S=0,jy,ans,sum;int first[M],next[M],v[M],w[M],vis[M/16],tot,T,cnt;char xx[]={
1,-1,2,-2,1,-1,2,-2};char yy[]={
2,-2,1,-1,-2,2,-1,1};void Add(int x,int y,int z){w[tot]=z,v[tot]=y,next[tot]=first[x],first[x]=tot++;}void add(int x,int y,int z){Add(x,y,z),Add(y,x,0);}bool tell(){ memset(vis,-1,sizeof(vis));vis[S]=0; queue
q;q.push(S); while(!q.empty()){ int t=q.front();q.pop(); for(int i=first[t];~i;i=next[i]) if(w[i]&&vis[v[i]]==-1) q.push(v[i]),vis[v[i]]=vis[t]+1; } return vis[T]!=-1;}int zeng(int x,int y){ if(x==T)return y; int r=0; for(int i=first[x];y>r&&~i;i=next[i]) if(w[i]&&vis[v[i]]==vis[x]+1){ int t=zeng(v[i],min(y-r,w[i])); w[i]-=t,w[i^1]+=t,r+=t; } if(!r)vis[x]=-1; return r;}bool check(int x,int y){ if(x<1||x>n||y<1||y>m||a[x][y])return false; return true;}int main(){ memset(first,-1,sizeof(first)); scanf("%d%d",&n,&m);T=n*m+1; for(int i=1;i<=n;i++) for(int j=1;j<=m;j++){ id[i][j]=++cnt; scanf("%d",&a[i][j]); if((i+j)&1){color[i][j]=1,add(S,id[i][j],1);} else add(id[i][j],T,1); if(!a[i][j])sum++; } for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(color[i][j]&&!a[i][j]) for(int k=0;k<8;k++){ int tx=i+xx[k],ty=j+yy[k]; if(check(tx,ty))add(id[i][j],id[tx][ty],1); } while(tell())while(jy=zeng(S,0x3f3f3f3f))ans+=jy; printf("%d\n",sum-ans);}

 

转载于:https://www.cnblogs.com/SiriusRen/p/6654369.html

你可能感兴趣的文章
web富客户端应用下,前端架构、系列(二)。
查看>>
CentOS yum 安装node.js
查看>>
[HEOI2016/TJOI2016]游戏 解题报告
查看>>
selenium模块用法详解
查看>>
class的继承,从基类开始
查看>>
一个简单的ssm项目
查看>>
webapi+entityframework分享
查看>>
Python—模块介绍
查看>>
各浏览器下 scrollTop的差异
查看>>
html-id,class
查看>>
DOM_radio
查看>>
如何更高效的使用谷歌解决问题
查看>>
Linux使用wget仿站
查看>>
C#学习-子类的初始化顺序
查看>>
内存泄露 内存溢出
查看>>
commons-lang3-3.4类库
查看>>
使用 windows 计划任务播放音乐文件
查看>>
java 堆、栈
查看>>
Xdebug配置
查看>>
ES6常用知识点整理(一)
查看>>