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<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#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]<d[v])  
            {  
                  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<<ans-S.max_flow(s,t)<<endl;
    return 0;
}