[二分图最大独立集] BZOJ4808：马

棋盘01染色，然后把互相能打到的点连边，发现是个二分图。
二分图最大独立集 $=$ 点数 $-$ 最大匹配数
Dinic练手…

#include<cstdio>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=40005,maxe=400005,flag[8][2]={{-1,-2},{-1,2},{1,-2},{1,2},{2,1},{-2,1},{2,-1},{-2,-1}};
int n,m,ans,cnt,N,S,T,map[205][205],c[205][205];
bool check(int x,int y){ return 1<=x&&x<=n&&1<=y&&y<=m&&!map[x][y];  }
struct Edge{
    int from,to,cap,flow;
    Edge(int t1=0,int t2=0,int t3=0,int t4=0){from=t1;to=t2;cap=t3;flow=t4;}
} Es[maxe];
int fir[maxn],pos[maxn],nxt[maxe],tot=1; 
void add(int x,int y,int z){
    //printf("%d -> %d\n",x,y);
    Es[++tot]=Edge(x,y,z,0);
    nxt[tot]=fir[x]; fir[x]=tot;
    Es[++tot]=Edge(y,x,0,0);
    nxt[tot]=fir[y]; fir[y]=tot;
}
queue<int> Q;
int d[maxn],INF;
bool Bfs(){
    memset(d,63,sizeof(d)); INF=d[0];
    d[S]=0; Q.push(S);
    while(!Q.empty()){
        int x=Q.front(); Q.pop();
        for(int j=fir[x];j;j=nxt[j]) if(Es[j].cap>Es[j].flow&&d[Es[j].to]==INF){
            d[Es[j].to]=d[x]+1; Q.push(Es[j].to);
        }
    }
    return d[T]!=INF;
}
int Dfs_find(int x,int flw){
    if(x==T||flw==0) return flw;
    int res=0,t;
    for(int &j=pos[x];j;j=nxt[j])
     if(d[x]+1==d[Es[j].to]&&(t=Dfs_find(Es[j].to,min(flw,Es[j].cap-Es[j].flow)))>0){
        flw-=t; Es[j].flow+=t; Es[j^1].flow-=t; res+=t;
        if(flw==0) break;
     }
    return res;
}
int main(){
    freopen("bzoj4808.in","r",stdin);
    freopen("bzoj4808.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        c[i][1]=c[i-1][1]^1;
        for(int j=2;j<=m;j++) c[i][j]=c[i][j-1]^1;
        for(int j=1;j<=m;j++) scanf("%d",&map[i][j]);
    }
    N=n*m+2; S=n*m+1; T=n*m+2;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++) if(!map[i][j]){
        cnt++;
        if(c[i][j]) add(S,(i-1)*m+j,1); else add((i-1)*m+j,T,1);
     }
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++) if(c[i][j]&&!map[i][j])
      for(int k=0;k<=7;k++) if(check(i+flag[k][0],j+flag[k][1])){
        int x=i+flag[k][0], y=j+flag[k][1];
        add((i-1)*m+j,(x-1)*m+y,1);
      }
    while(Bfs()){
        for(int i=1;i<=N;i++) pos[i]=fir[i];
        ans+=Dfs_find(S,1e9);
    }
    printf("%d\n",cnt-ans);
    return 0;
}