codevs 1922 骑士共存问题 最小割

10 篇文章 0 订阅
1 篇文章 0 订阅

分析:比较经典的最小割题目,我们把横纵坐标之和%2=0的位置标记为二分图的左半部分,另外的作为右半部份,当然删除掉的不能用的。
然后s向左半部分连边,右半部分向T连边。对于能够互相到达的点,连长度为inf的边(保证最小割不会包含可以互相到达的点),然后可以证明互相到达的两个点,肯定属于不同的部分,然后跑最小割,答案就是n^2-m-ans.
手残打错了一些小细节,,非常不应该啊。。还是打得不是很熟练,要多打。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define inf 0x7fffffff
using namespace std;
typedef long long ll;
const int N=5e5+5;
int n,m,part1,part2;
const int M=505;
int a[M][M],del[M][M];
const int dx[8]={-2,-2,-1,1,2,2,1,-1};
const int dy[8]={-1,1,2,2,1,-1,-2,-2};
int tot=1;
int head[N],next[N],go[N],q[N],dis[50001],val[N];
inline void add(int x,int y,int z)
{
    go[++tot]=y;
    val[tot]=z;
    next[tot]=head[x];
    head[x]=tot;
}
inline void ins(int x,int y,int z)
{
    add(x,y,z);
    add(y,x,0);
}
ll ans=0;
inline bool bfs()
{
    int t=0,w=1;
    memset(dis,-1,sizeof(dis));
    dis[0]=0;
    q[1]=0;
    while (t<w)
    {
        int x=q[++t];
        for(int i=head[x];i;i=next[i])
        {
            int v=go[i];
            if (val[i]&&dis[v]==-1)
            {
                dis[v]=dis[x]+1;
                q[++w]=v;
            }
        }
    }
    return dis[part2]!=-1;
}
inline int dfs(int x,int f)
{
    if (x==part2)return f;
    int w,used=0;
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (val[i]&&dis[v]==dis[x]+1)
        {
            w=f-used;
            w=dfs(v,min(w,val[i]));
            val[i]-=w;
            val[i^1]+=w;
            used+=w;
            if (used==f)return f;
        }
    }
    if (!used)dis[x]=-1;
    return used;
}
inline void dinic()
{
    while (bfs())ans+=dfs(0,inf);
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,m)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        del[x][y]=1;
    }
    part1=1,part2=(n*n+1)/2+1;
    fo(i,1,n)
    {
        fo(j,1,n)
        {
            if ((i+j)%2==0)
            a[i][j]=part1++;
            else a[i][j]=part2++;
        }
    }
    fo(i,1,n)
    fo(j,1,n)
    if (del[i][j])continue;
    else if (a[i][j]<part1)
    {
        fo(k,0,7)
        {
            int x=i+dx[k];
            int y=j+dy[k];
            if (x<1||x>n||y<1||y>n||del[x][y])continue;
            ins(a[i][j],a[x][y],inf);
        }
        ins(0,a[i][j],1);
    }
    else ins(a[i][j],part2,1);
    dinic();
    printf("%d\n",n*n-m-ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值