P3163 [CQOI2014]危桥

49 篇文章 0 订阅
41 篇文章 0 订阅

类型:网络最大流
建模分析:
普通的边按图中说的建立就好。
由于是来回 我们需要 从源点向a1,b1分别建立一条容量为 2 * an 和 2 * bn的边。
从a2 ,b2 向汇点分别建立一条容量为 2 * an和2 * bn的边。
然后跑一下网络最大流 看看最大流是不是 2*(an+bn)
你以为这就完了?
不存在的,上述的建图方式如果在a1跑到 b2的话,是不对的,但是我们判断出了合法。
那该砸处理?
我也不知道为啥傻里吧唧的先把图重建一遍 然后源点连a2,b2。汇点连a1,b2。
这样搞跟上面没区别啊。。。。
我真的菜。
后来想了5分钟,忽然发现我们只需让B反过来建就好了。。。。
Luogu
BZOJ

#include <cstdio>
#include <iostream>
#include <queue>
#include <cstring> 
#define il inline
using namespace std;
const int inf=0x7fffffff;
const int maxm=11000;
int head[maxm],to[maxm*2],cap[maxm*2],net[maxm*2],deep[maxm],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 n,a1,a2,an,b1,b2,bn;
char bri[90][90];
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;
}
il void init()
{
    memset(head,0,sizeof(head));
    cnt=1;
}
il bool work()
{
    for(int i=0;i<n;i++)
     for(int j=0;j<n;j++)
      {
        cin>>bri[i][j];
        if(bri[i][j]=='O') adx(i,j,2);
        if(bri[i][j]=='N') adx(i,j,inf);
      }
    int s=1001,t=1000;
    adx(s,a1,2*an),adx(s,b1,2*bn);
    adx(a2,t,2*an),adx(b2,t,2*bn);
    if(2*(an+bn)!=dinic(s,t)) return 0;
    memset(head,0,sizeof(head));
    cnt=1;
    for(int i=0;i<n;i++)
     for(int j=0;j<n;j++)
      {
        if(bri[i][j]=='O') adx(i,j,2);
        if(bri[i][j]=='N') adx(i,j,inf);
      }
    adx(s,a1,2*an),adx(s,b2,2*bn);
    adx(a2,t,2*an),adx(b1,t,2*bn);
    if(2*(an+bn)!=dinic(s,t)) return 0;
    return 1;
}
int main()
{
    while((scanf("%d%d%d%d%d%d%d",&n,&a1,&a2,&an,&b1,&b2,&bn))!=EOF)
    {
        init();
        if(work()) printf("Yes\n");
        else printf("No\n");
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值