BZOJ 4808(马-二分图最大独立集)
4808: 马
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 111 Solved: 46
[Submit][Status][Discuss]
Description
众所周知，马后炮是中国象棋中很厉害的一招必杀技。”马走日字”。本来，如果在要去的方向有别的棋子挡住（俗
称”蹩马腿”），则不允许走过去。为了简化问题，我们不考虑这一点。马跟马显然不能在一起打起来，于是rly在
一天再次借来了许多许多的马在棋盘上摆了起来……但这次，他实在没兴趣算方案数了，所以他只想知道在N×M的
矩形方格中摆马使其互不吃到的情况下的最多个数。但是，有一个很不幸的消息，rly由于玩得太Happy，质量本来
就不好的棋盘被rly弄坏了，不过幸好只是破了其中的一些格子（即不能再放子了），问题还是可以继续解决的。
Input
一行，两个正整数N和M。
接下来N行，每行M个数，要么为0，表示没坏，要么为1，表示坏了。
N<=200,M<=200
Output
一行，输出最多的个数。
Sample Input
2 3
0 1 0
0 1 0
Sample Output
2
HINT
Source
By FancyCoder
棋盘黑白染色，约束变为若干（黑，白）对点不能同时取
#include
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define MAXN (1000000+100)
#define MAXM (6000000+100)
long long mul(long long a,long long b){return (a*b)%F;}
long long add(long long a,long long b){return (a+b)%F;}
long long sub(long long a,long long b){return (a-b+(a-b)/F*F+F)%F;}
typedef long long ll;
int read()
{
int x=0,f=1; char ch=getchar();
while(!isdigit(ch)) {if (ch=='-') f=-1; ch=getchar();}
while(isdigit(ch)) { x=x*10+ch-'0'; ch=getchar();}
return x*f;
}
class Max_flow //dinic+当前弧优化
{
public:
int n,t;
int q[MAXN];
int edge[MAXM],Next[MAXM],Pre[MAXN],weight[MAXM],size;
void addedge(int u,int v,int w)
{
edge[++size]=v;
weight[size]=w;
Next[size]=Pre[u];
Pre[u]=size;
}
void addedge2(int u,int v,int w){addedge(u,v,w),addedge(v,u,0);}
bool b[MAXN];
int d[MAXN];
bool SPFA(int s,int t)
{
For(i,n) d[i]=INF;
MEM(b)
d[q[1]=s]=0;b[s]=1;
int head=1,tail=1;
while (head<=tail)
{
int now=q[head++];
Forp(now)
{
int &v=edge[p];
if (weight[p]&&!b[v])
{
d[v]=d[now]+1;
b[v]=1,q[++tail]=v;
}
}
}
return b[t];
}
int iter[MAXN];
int dfs(int x,int f)
{
if (x==t) return f;
Forpiter(x)
{
int v=edge[p];
if (weight[p]&&d[x]int nowflow=dfs(v,min(weight[p],f));
if (nowflow)
{
weight[p]-=nowflow;
weight[p^1]+=nowflow;
return nowflow;
}
}
}
return 0;
}
int max_flow(int s,int t)
{
(*this).t=t;
int flow=0;
while(SPFA(s,t))
{
For(i,n) iter[i]=Pre[i];
int f;
while (f=dfs(s,INF))
flow+=f;
}
return flow;
}
void mem(int n)
{
(*this).n=n;
size=1;
MEM(Pre)
}
}S;
int n,m;
int dx[8]={1,2,1,-2,-1,2,-1,-2};
int dy[8]={2,1,-2,1,2,-1,-2,-1};
bool inside(int i,int j) {
return 1<=i&&i<=n&&1<=j&&j<=m;
}
int id(int i,int j) {
return (i-1)*m+j;
}
int main()
{
// freopen("bzoj4808.in","r",stdin);
// freopen(".out","w",stdout);
n=read(),m=read();
int a[300][300];
int s=n*m+1,t=s+1;
S.mem(t);
int ans=0;
For(i,n) For(j,m) a[i][j]=read();
For(i,n) {
For(j,m) {
if (a[i][j]==0) {
++ans;
if ((i+j)&1) S.addedge2(s,id(i,j),1);
else S.addedge2(id(i,j),t,1);
Rep(k,8) {
int nx=i+dx[k];
int ny=j+dy[k];
if (inside(nx,ny)) if (!a[nx][ny]){
if ((i+j)&1) S.addedge2(id(i,j),id(nx,ny),1);
else S.addedge2(id(nx,ny),id(i,j),1);
}
}
}
}
}
cout<return 0;
}