[codevs1022]覆盖(染色+二分图最大匹配)

传送门
将格子染色之后就可以将问题转化成二分图最大匹配的问题了。
代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<queue>
#define ll long long
using namespace std;
inline int read(){
    int x=0;char ch=' ';int f=1;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int inf=0x3f3f3f3f;
struct edge{
    int to,next,cap,flow;
}e[200001];
int n,m,s,t,tot=-1;
int head[10002];
inline void addedge(int x,int y,int l){
    e[++tot].to=y;e[tot].next=head[x];e[tot].cap=l;e[tot].flow=0;head[x]=tot;
    e[++tot].to=x;e[tot].next=head[y];e[tot].cap=0;e[tot].flow=0;head[y]=tot;
}
int dis[10002];
inline bool bfs(){
    memset(dis,-1,sizeof(dis));
    queue<int> q;
    q.push(s);
    dis[s]=0;
    while(!q.empty()){
        int x=q.front();
        q.pop();
        for(int i=head[x];i!=-1;i=e[i].next){
            int u=e[i].to;
            if(dis[u]==-1&&e[i].flow<e[i].cap){
                dis[u]=dis[x]+1;
                q.push(u);
            }
        }
    }
    return dis[t]!=-1;
}
inline int dfs(int x,int a){
    if(x==t||a==0)return a;
    int flow=0,f;
    for(int i=head[x];i!=-1;i=e[i].next){
        int u=e[i].to;
        if(e[i].flow<e[i].cap&&dis[u]==dis[x]+1&&(f=dfs(u,min(a,e[i].cap-e[i].flow)))>0){
            e[i].flow+=f;
            e[i^1].flow-=f;
            flow+=f;
            a-=f;
            if(a==0)break;
        }
    }
    if(!flow)dis[x]=0;
    return flow;
}
inline ll dinic(){
    ll flow=0;
    while(bfs()){
        flow+=dfs(s,inf);
    }
    return flow;
}
int k;
int color[101][101];
int main(){
    memset(head,-1,sizeof(head));
    n=read();m=read();s=0;t=n*m+1;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if((i+j)&1)color[i][j]=1;
        }
    }
    k=read();
    for(int i=1;i<=k;i++){
        int x=read(),y=read();
        color[x][y]=-1;
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(color[i][j]==-1)continue;
            if(color[i][j]==1)addedge(s,(i-1)*m+j,1);
            else{
                addedge((i-1)*m+j,t,1);continue;
            }
            if(i>1&&color[i-1][j]!=-1)addedge((i-1)*m+j,(i-2)*m+j,1);
            if(i<n&&color[i+1][j]!=-1)addedge((i-1)*m+j,i*m+j,1);
            if(j>1&&color[i][j-1]!=-1)addedge((i-1)*m+j,(i-1)*m+j-1,1);
            if(j<n&&color[i][j+1]!=-1)addedge((i-1)*m+j,(i-1)*m+j+1,1);
        }
    }
    printf("%lld",dinic());
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值