CQOI2014 危桥

危桥

 

这题告诉我们人有多大胆，乱搞是标算。

乱猜一下交换一对源汇后仍然可以成立才是Yes

然后。。。。就这样反直觉的证出来了，

对于每一个a1->b2,b1->a2的路径

可以交换一对源汇（a1,a2），然后之前每一个a1->b2,b1->a2的路径就会对应变成了a2->a1,b1->b2的一条路径

新产生的a2->b2,b1->a1的路径对应之前的a1->a2,b1->b2的路径。

两个最大流相互佐证，都满足条件就是yes

最大流无向图建边就是把有向图的建边函数的反向边流量设为和正向边一样。

不要自己乱搞调不出来。。。。。。

#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#define maxn 55
#define maxm 55*55*8
#define inf 0x3f3f3f3f
using namespace std;

int n,a[3],b[3];
int info[maxn],Prev[maxm],to[maxm],cap[maxm],cnt_e;
int buff[maxm];
inline void Node(int u,int v,int c)
{
    Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v;
    cap[cnt_e]=c;
}
inline void Line(int u,int v,int c)
{
    Node(u,v,c),Node(v,u,c);
}

int dis[maxn],stm,S,T;

int Q[maxn],L,R;
bool BFS()
{
	memset(dis,0x3f,sizeof dis);
	L=R=0;
	dis[Q[R++]=T]=0;
	for(int now;L<R;)
	{
		now=Q[L++];
		for(int i=info[now];i;i=Prev[i])
			if(cap[i^1] && dis[to[i]]>dis[now]+1)
			{
				dis[to[i]]=dis[now]+1;
				Q[R++]=to[i];
			}
	}
	return dis[S]<0x3f3f3f3f;
}

int aug(int now,int Max)
{
    if(now==T)
        return Max;
    int inc,st=Max;
    for(int i=info[now];i && st;i=Prev[i])
        if(cap[i]&&dis[to[i]]+1==dis[now])
		{
			inc=aug(to[i],min(st,cap[i]));
			st-=inc,cap[i]-=inc,cap[i^1]+=inc;
		}
    return Max-st;
}

int sap()
{
    for(stm=0;BFS();) 
		stm+=aug(S,inf);
    return stm;
}

int main()
{
	
	//freopen("1.in","r",stdin);
	
    char ch[100];
    while(~scanf("%d%d%d%d%d%d%d",&n,&a[0],&a[1],&a[2],&b[0],&b[1],&b[2]))
    {
		memset(cap,0,sizeof cap);
        memset(info,0,sizeof info); cnt_e=1;
        S=n+1,T=n+2;
        a[0]++,a[1]++,b[0]++,b[1]++;a[2]<<=1,b[2]<<=1;
        for(int i=1;i<=n;i++)
        {
            scanf("%s",ch+1);
            for(int j=1;j<=i;j++)
			{
				if(ch[j]=='O')  Line(i,j,2);
				if(ch[j]=='N')	Line(i,j,inf);
			}
        }
                
        memcpy(buff,cap,sizeof cap);
        Line(S,a[0],a[2]),Line(a[1],T,a[2]),Line(S,b[0],b[2]),Line(b[1],T,b[2]);
        if(sap()<a[2]+b[2]) puts("No");
        else 
        {
            memcpy(cap,buff,sizeof buff);
            swap(a[0],a[1]);
            Line(S,a[0],a[2]),Line(a[1],T,a[2]),Line(S,b[0],b[2]),Line(b[1],T,b[2]);
            if(sap()<a[2]+b[2]) puts("No");
            else puts("Yes");
        }
    }
}