1. 程式人生 > >bzoj3504:[Cqoi2014]危橋

bzoj3504:[Cqoi2014]危橋

clu class line 可能 zoj3504 fin set tchar fine

傳送門

好像還是個比較簡單的網絡流
說實話我不知道網上雙向邊建成容量為\(2\)的為什麽是對的,我感覺容量為\(1\)才是對的啊
然後唯一的坑點就是直接建圖\(a1\)可能會流到\(b2\)\(b1\)可能會流到\(a2\),這就會誤判
然後交換一組起點和終點就好了,再跑就行了(這也是我看了題解的地方)
代碼:

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
#define rg register
const int maxn=3e4+10,inf=1e9;char c[61][61];queue<int>q;
int dis[maxn],n,a1,s,t,a2,an,b1,b2,bn,cnt=1,ans1,ans2,pre[maxn*2],nxt[maxn*2],h[maxn],v[maxn],cur[maxn];
void add(int x,int y,int z)
{
    pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt,v[cnt]=z;
    pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt,v[cnt]=0;
}
bool bfs()
{
    memset(dis,0,sizeof dis);
    q.push(s),dis[s]=1;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(rg int i=h[x];i;i=nxt[i])
            if(!dis[pre[i]]&&v[i])dis[pre[i]]=dis[x]+1,q.push(pre[i]);
    }
    return dis[t];
}
int dfs(int x,int flow)
{
    if(x==t||!flow)return flow;
    int f=flow;
    for(rg int &i=cur[x];i;i=nxt[i])
        if(dis[pre[i]]>dis[x]&&v[i])
        {
            int y=dfs(pre[i],min(f,v[i]));
            f-=y,v[i]-=y,v[i^1]+=y;
            if(!f)return flow;
        }
    if(f==flow)dis[x]=-1;
    return flow-f;
}
int main()
{
    while(scanf("%d%d%d%d%d%d%d",&n,&a1,&a2,&an,&b1,&b2,&bn)!=EOF)
    {
        a1++,a2++,b1++,b2++,cnt=1,ans1=0,memset(h,0,sizeof h);
        for(rg int i=1;i<=n;i++)scanf("%s",c[i]+1);
        for(rg int i=1;i<=n;i++)
            for(rg int j=1;j<=n;j++)
                if(c[i][j]=='O')add(i,j,1);
                else if(c[i][j]=='N')add(i,j,inf);
        s=0,t=n+1,add(s,a1,an),add(a2,t,an),add(s,b1,bn),add(b2,t,bn);
        for(;bfs();ans1+=dfs(s,inf))memcpy(cur,h,sizeof h);
        if(ans1<an+bn){printf("No\n");continue;}
        ans1=0,memset(h,0,sizeof h),cnt=1;
        for(rg int i=1;i<=n;i++)
            for(rg int j=1;j<=n;j++)
                if(c[i][j]=='O')add(i,j,1);
                else if(c[i][j]=='N')add(i,j,inf);
        add(s,a1,an),add(a2,t,an),add(s,b2,bn),add(b1,t,bn);
        for(;bfs();ans1+=dfs(s,inf))memcpy(cur,h,sizeof h);
        if(ans1<an+bn){printf("No\n");continue;}
        printf("Yes\n");
    }
}

bzoj3504:[Cqoi2014]危橋