[P3355骑士共存]

P3355

• 方格里面选若干个点放上骑士，骑士之间不能互相攻击到。问最多能放多少骑士。
• emm。可能是被费用流整自闭了，还以为这道题也是个神题，在想怎么用费用流跑。
• 看了一下自己之前居然写过这道题，代码居然是匈牙利？？好，再看一下题目。等等，最多放多少骑士，使得两两之间不冲突，那不就是最大独立集吗？
• 对于任意一个点，如果没有障碍，把它所能攻击到的点都向它连一条边。可以证明，这样建的图一定是没有奇环的，二分图染色一下，建图跑最大流就好了。
• 答案就是$n\ast n-m-maxflow$。

Coding

#include<bits/stdc++.h>
using namespace std;
const int N=4e4+100;
const int M=1e6+10;
const int inf=1e9;
const int dx[8]={-1,-2,1,2,-1,-2,1,2};
const int dy[8]={-2,-1,2,1,2,1,-2,-1};
int n,m,tot=1,maxflow,flow,s,t,d[N],ver[M],Next[M],edge[M],lin[N],color[N];
int a[210][210];
int cal(int i,int j){return (i-1)*n+j;}
void add(int x,int y){ver[++tot]=y;Next[tot]=lin[x];lin[x]=tot;}
void Add(int x,int y,int z){
    ver[++tot]=y;Next[tot]=lin[x];lin[x]=tot;edge[tot]=z;
    ver[++tot]=x;Next[tot]=lin[y];lin[y]=tot;edge[tot]=0;
}
void dfs(int x,int col){
    color[x]=col;
    for(int i=lin[x];i;i=Next[i]){
        int y=ver[i];
        if(!color[y]) dfs(y,3-col);
    }
}
bool bfs(){
    queue<int>q;
    memset(d,0,sizeof(d));
    d[s]=1,q.push(s);
    while(q.size()){
        int x=q.front();q.pop();
        for(int i=lin[x];i;i=Next[i]){
            int y=ver[i];
            if(edge[i]&&!d[y]){
                d[y]=d[x]+1;q.push(y);
                if(y==t) return 1;
            }
        }
    }return 0;
}
int dinic(int x,int flow){
    if(x==t) return flow;
    int rest=flow;
    for(int i=lin[x];i&&rest;i=Next[i]){
        int y=ver[i];
        if(edge[i]&&d[y]==d[x]+1){
            int k=dinic(y,min(edge[i],rest));
            if(!k) d[y]=0;
            rest-=k,edge[i]-=k,edge[i^1]+=k;
            if(!rest) return flow-rest;
        }
    }return flow-rest;
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;++i){
        int x,y;scanf("%d%d",&x,&y);
        a[x][y]=1;
    }s=0,t=n*n+1;
    for(int i=1;i<=n;++i) for(int j=1;j<=n;++j){
        if(a[i][j]) continue;
        for(int k=0;k<8;++k){
            int tx=i+dx[k],ty=j+dy[k];
            if(a[tx][ty]) continue;
            if(tx<=0||tx>n||ty<=0||ty>n) continue;
            add(cal(tx,ty),cal(i,j));add(cal(i,j),cal(tx,ty));
        }
    }
    for(int i=1;i<=n;++i) for(int j=1;j<=n;++j){
        if(!color[cal(i,j)]&&!a[i][j]) dfs(cal(i,j),1);
    }
    memset(lin,0,sizeof(lin));tot=1;
    for(int i=1;i<=n;++i) for(int j=1;j<=n;++j){
        if(a[i][j]) continue;
        if(color[cal(i,j)]==1) Add(s,cal(i,j),1);
        else Add(cal(i,j),t,1);
        for(int k=0;k<8;++k){
            int tx=i+dx[k],ty=j+dy[k];
            if(a[tx][ty]) continue;
            if(tx<=0||tx>n||ty<=0||ty>n) continue;
            int x=cal(i,j),y=cal(tx,ty);
            if(color[x]>color[y]) swap(x,y);
            Add(x,y,1);
        }
    }
    while(bfs()){
        while(flow=dinic(s,inf)) maxflow+=flow;
    }
    cout<<n*n-m-maxflow;
    return 0;
}