骑士共存问题

类型：最小割

建模分析：

建模同方格取数问题。
不同的是本题有位置不能放骑士，那么此位置不能加入图中。

题目链接：

COGS
Luogu

Ac Code：

#include 
#include 
#include 
#include  
#define il inline
using namespace std;
const int inf=0x7fffffff;
const int maxm=1100000;
int head[maxm],to[maxm*2],cap[maxm*2],net[maxm*2],deep[210*210],cnt=1;
il void add(int x,int y,int c){cnt++,to[cnt]=y,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
queue <int> dl;
int id[300][300],map[300][300];
int dx[]={0,1,-1,1,-1,2,2,-2,-2};
int dy[]={0,2,2,-2,-2,1,-1,1,-1};
il bool BFS(int s,int t)
{
    while(!dl.empty()) dl.pop();
    memset(deep,-1,sizeof(deep));
    dl.push(s),deep[s]=0;
    while(!dl.empty())
    {
        int x=dl.front();dl.pop();
        for(int i=head[x];i;i=net[i])
         if(cap[i]>0&&deep[to[i]]==-1)
          dl.push(to[i]),deep[to[i]]=deep[x]+1;
    }
    return deep[t]==-1?0:1;
}
int dfs(int now,int flow,int t)
{
    if(now==t) return flow;
    int w,used=0;
    for(int i=head[now];i;i=net[i])
    {
        int v=to[i];
        if(deep[v]==deep[now]+1&&cap[i])
        {
            w=dfs(v,min(flow-used,cap[i]),t);
            cap[i]-=w;
            cap[i^1]+=w;
            used+=w;
            if(used==flow) return flow;
        }
    }
    if(!used) deep[now]=-1;
    return used;
}
il int dinic(int s,int t)
{
    int maxflow=0;
    while(BFS(s,t)) maxflow+=dfs(s,inf,t);
    return maxflow;
}
inline void adx(int x,int y,int cax)
{
    add(x,y,cax),add(y,x,0);
}
il int read()
{
    int x=0,w=1;
    char ch=0;
    while(ch<'0'||ch>'9') 
    {
        if(ch=='-') w=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
     x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
    return x*w;
}
int main()
{
    //freopen("knight.in","r",stdin);
    //freopen("knight.out","w",stdout);
    int n=read(),m=read();
    int s=0,t=n*n+3,sum=n*n-m;
    int num=0;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      id[i][j]=++num;
    for(int i=1,x,y;i<=m;i++)
     x=read(),y=read(),map[x][y]=1;
    //puts("Yes\n");
    for(int i=1;i<=n;i++)
     for(int j=1;j<=n;j++)
      if(!map[i][j])
      {
        if((i+j)%2)
        {
            adx(s,id[i][j],1);
            for(int k=1;k<=8;k++)
            {
                int bx=i+dx[k],by=j+dy[k];
                if(bx<1||by<1||bx>n||by>n) continue;
                if(map[bx][by]) continue;
                adx(id[i][j],id[bx][by],inf);
            }
        }   
        else adx(id[i][j],t,1);
      } 
    return printf("%d\n",sum-dinic(s,t))*0;
}