骑士共存问题【网络流24题】

题目链接:https://www.luogu.org/problem/P3355

首先,攻击关系不会出现奇环,所以可以对棋盘进行二分染色。

然后就和方格取数那个题一样了。

S连每个黑,边权为1.

每个黑连其能攻击的点,边权为inf。

每个白连T,边权为1。

#include <bits/stdc++.h>
#define rep(i, a, b) for(int i = (a); i <= (b); i++)
using namespace std;
const int N = 4e4+100;
int n,m,S,T,tot;
int head[N];
struct node{
    int v,cap,nxt;
}edge[int(3e6+100)];
int cur[N],deep[N];
void ae(int u,int v,int cap){  //前向星加边
    edge[++tot] = node{v,cap,head[u]};
    head[u] = tot;
    edge[++tot] = node{u,0,head[v]};
    head[v] = tot;
}
bool bfs(){  //构建分层图
    memset(deep,-1,sizeof(deep));  //-1表示未访问
    queue<int>q;
    q.push(S);
    deep[S] = 0;
    while(!q.empty()){
        int u = q.front();
        for(int i = head[u]; ~i;i = edge[i].nxt){
            int v = edge[i].v;
            if(deep[v]==-1&&edge[i].cap>0){
                q.push(v);
                deep[v] = deep[u]+1;
                if(v==T) return 1;   //提前退出
            }
        }
        q.pop();
    }
    return 0;
}
int dfs(int u,int f){  //寻找增广网
    int flow = 0,d;
    if(u==T||f==0) return f;
    for(int &i = cur[u]; ~i;i = edge[i].nxt){ //当前弧优化
        int v = edge[i].v;
        if(deep[v]>deep[u]&&edge[i].cap>0&&(d=dfs(v,min(f,edge[i].cap)))){
            edge[i].cap -= d;
            edge[i^1].cap += d;
            f -= d;
            flow += d;
            if(!f) break;
        }
    }
    if(flow==0) deep[u] = -1;  //删除无法增广的点
	return flow;
}
int dinic(){
    int ans = 0;
    while(bfs()){
        memcpy(cur, head, sizeof(head));  //重置cur数组
        ans += dfs(S,1e9);
    }
    return ans;
}

int dx[8]={-2,-1,1,2,2,1,-1,-2};
int dy[8]={1,2,2,1,-1,-2,-2,-1};
bool flag[220][220];
int col[220][220];
bool ok(int x,int y) {
    if(x<0||x>=n||y<0||y>=n) return 0;
    if(col[x][y]) return 0;
    return 1;
}
void Dfs(int x,int y,int now) {
    col[x][y] = now;
    rep(i, 0, 7) {
        int xx = x+dx[i];
        int yy = y+dy[i];
        if(ok(xx,yy)) Dfs(xx,yy,-now);
    }
}
int id(int x,int y) {
    return x*n+y;
}
int sum;
int main() {
    //freopen("a.txt","r",stdin);
    ios::sync_with_stdio(0);
    cin>>n>>m;
    sum = n*n-m;
    memset(head,-1,sizeof(head));  //head初始化-1
    tot = -1;

    rep(i, 1, m) {
        int x,y;
        cin>>x>>y;
        x--;y--;
        flag[x][y] = 1;
    }

    rep(i, 0, n-1)
        rep(j, 0, n-1) {
            if(col[i][j]) continue;
            Dfs(i,j,1);
        }
    S = n*n;
    T = n*n+1;
    rep(i, 0, n*n-1) {
        int x = i/n;
        int y = i%n;
        if(flag[x][y]) continue;
        if(col[x][y]==1) {
            ae(S,i,1);
            rep(j, 0, 7) {
                int xx = x+dx[j];
                int yy = y+dy[j];
                if((xx>=0&&xx<n&&yy>=0&&yy<n) && (!flag[xx][yy]))
                    ae(i,id(xx,yy),1e9);
            }
        }
        else ae(i,T,1);
    }

    cout<<sum-dinic()<<endl;
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值