【Luogu】P3355骑士共存问题(最小割)

  题目链接

  像题面那样把棋盘染成红黄点。发现骑士迈一步能到达的点的颜色一定是跟他所在的格子的颜色不同的。于是(woc哪来的于是?这个性质有这么明显吗?)从源点向所有红点连边,从所有黄点向汇点连边,红点向从它一步能迈到的黄点连边,然后求出最小割(其实就是最大流啦)用可以站骑士的总点数减去。

  网络流构图好喵啊     一脸懵逼

  

#include<cstdio>
#include<cctype>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<queue>
#define maxn 80000
using namespace std;

inline long long read(){
    long long num=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-1;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*10+ch-'0';
        ch=getchar();
    }
    return num*f;
}

int n,m;
inline int count(int i){    return i&1?i+1:i-1;    }
inline int calc(int x,int y){    return (x-1)*n+y;    }

int u[9]={0,1,2,2,1,-1,-2,-2,-1};
int w[9]={0,-2,-1,1,2,2,1,-1,-2};

struct Edge{
    int next,to,val;
}edge[5*maxn];
int head[maxn*3],num;
inline void addedge(int from,int to,int val){
    edge[++num]=(Edge){head[from],to,val};
    head[from]=num;
}
inline void add(int from,int to,int val){
    addedge(from,to,val);
    addedge(to,from,0);
}

bool vis[maxn*3];
bool ext[500][800];
int list[maxn*3];
int dfn[maxn*3];
int Start,End;

bool bfs(){
    memset(vis,0,sizeof(vis));
    queue<int> q;
    vis[Start]=1;    dfn[Start]=1;    q.push(Start);
    while(!q.empty()){
        int from=q.front();q.pop();
        for(int i=head[from];i;i=edge[i].next){
            int to=edge[i].to;
            if(edge[i].val<=0||vis[to])    continue;
            vis[to]=1;
            dfn[to]=dfn[from]+1;
            q.push(to);
        }
    }
    return vis[End];
}

int dfs(int x,int val){
    if(val==0||x==End)    return val;
    int flow=0;    vis[x]=1;
    for(int &i=list[x];i;i=edge[i].next){
        int to=edge[i].to;
        if(vis[to]||dfn[to]!=dfn[x]+1||edge[i].val==0)    continue;
        int now=dfs(to,min(val,edge[i].val));
        val-=now;edge[i].val-=now;edge[count(i)].val+=now;flow+=now;
        if(val==0)    break;
    }
    if(val!=flow)    dfn[x]=-1;
    return flow;
}

int maxflow(){
    int ans=0;
    while(bfs()){
        memset(vis,0,sizeof(vis));
        for(int i=Start;i<=End;++i)    list[i]=head[i];
        int now=dfs(Start,0x7fffffff);
        if(!now)    break;
        ans+=now;
    }
    return ans;
}


int main(){
    n=read(),m=read();    End=n*n+1;
    for(int i=1;i<=m;++i){
        int x=read(),y=read();
        ext[x][y]=1;
    }
    for(int i=1;i<=n;++i)
        for(int j=1;j<=n;++j){
            if(ext[i][j])    continue;
            if(!((i+j)&1)){
                add(Start,calc(i,j),1);
                for(int k=1;k<=8;++k){
                    int x=i+u[k];int y=j+w[k];
                    if(x<1||x>n||y<1||y>n)    continue;
                    if(ext[x][y])    continue;
                    add(calc(i,j),calc(x,y),0x7fffffff);
                }
            }
            else add(calc(i,j),End,1);
        }
    int ans=maxflow();
    printf("%d",n*n-m-ans);
}

 

posted @ 2018-01-05 08:51  Konoset  阅读(145)  评论(0编辑  收藏  举报