2014-2015 ACM-ICPC, Asia Xian Regional Contest(H - The Problem to Make You Happy )

Alice和Bob两人玩游戏。两人各有一个棋子，开始时在有向图$(n \le 100)$中不同的点上。Bob先手，两人轮流移动棋子，每次只能将棋子移到相邻的点上（沿着图上的有向边），不断重复这样的过程，当那人无法移动时，他就输了（规则1）。

还有两条附加规则：任何时候，当两个棋子到达同一点时，Alice胜出（规则2）。如果游戏永远无法结束，Bob胜出（规则3）。

两人都遵循最优策略，轮流移动棋子，问最终谁会赢。图中没有自环和重边。

可以列出方程，但本题有“无法结束”的情况 。
不妨只考虑能结束的情况
$f_{i,j,l}$表示 Bob 在i点 ，Alice在j点 ，轮到l的情况下，Bob是否必输。

根据规则2： $f_{i,i,l} =1 , f_{i,j,0} = 1(i无法移动)$
$显然f_{i,j,0}= 1 (i \ne j)$，当且仅当对于所有后继节点$v (i \to v)$,有$f_{v,j,0}=0$

#include<bits/stdc++.h>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define ForkD(i,k,n) for(int i=n;i>=k;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=Pre[x];p;p=Next[p])
#define Forpiter(x) for(int &p=iter[x];p;p=Next[p])  
#define Lson (o<<1)
#define Rson ((o<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define pb push_back
#define mp make_pair 
#define fi first
#define se second
#define MAXN (233+10)
#define MAXM (20000+10)
typedef long long ll;
typedef unsigned long long ull;
ll mul(ll a,ll b){return (a*b)%F;}
ll add(ll a,ll b){return (a+b)%F;}
ll sub(ll a,ll b){return (a-b+llabs(a-b)/F*F+F)%F;}
void upd(ll &a,ll b){a=(a%F+b%F)%F;}
int n,m,outdegree[MAXN];

int Next[MAXM],edge[MAXM],Pre[MAXM],siz=1;
void addedge(int u,int v){
    edge[++siz]=v;
    Next[siz]=Pre[u];
    Pre[u]=siz;
}

// 0 Bob 1 Alice
int kcase;

bool b[MAXN][MAXN][2]={0};
int f[MAXN][MAXN][2];
// f[x][y][0]=f[y][x][1]
int go[MAXN][MAXN][2];

#define pir pair< pair<int,int> ,int > 
queue<pir> q;
bool bfs() {
    while(!q.empty()) q.pop();
    For(i,n) 
    {
        f[i][i][0]=1; f[i][i][1]=1;
        q.push(mp(mp(i,i) ,0));
        q.push(mp(mp(i,i) ,1));
    }
    For(i,n)
    {
        if (!outdegree[i]) 
            For(j,n)
            {
                if (i==j) continue;
                f[i][j][0]=1;
                q.push(mp(mp(i,j) ,0));
            }   
    } 

    while(!q.empty()) {
        pir now = q.front();
        q.pop();
        int x=now.fi.fi,y=now.fi.se,l=now.se;
        if (l==1)
        {
            Forp(x) {
                int v=edge[p];
                if (f[v][y][0]>-1) continue;
                go[v][y][0]++;
                if (go[v][y][0] == outdegree[v] ) 
                {
                    f[v][y][0]= 1;
                    q.push(mp(mp(v,y) ,0) );
                }
            }
        }
        else {
            Forp(y) {
                int v=edge[p];
                if (f[x][v][1]>-1) continue;

                f[x][v][1]=1;
                q.push(mp(mp(x,v) ,1) );
            }       
        }       

    }
}



int main()
{
//  freopen("CF100548H.in","r",stdin);
//  freopen("CF100548H.out","w",stdout);

    int T;cin>>T;
    for(kcase=1;kcase<=T;kcase++) {

        MEM(Next) MEM(edge) MEM(Pre) siz=1;
        MEM(outdegree) MEMi(f) MEM(go)

        cin>>n>>m;
        For(i,m) {
            int u,v;
            scanf("%d%d",&u,&v);
            addedge(v,u);
            outdegree[u]++;
        }
        int x,y;
        cin>>x>>y;
        if (x==y) {
            printf("Case #%d: No\n",kcase);
            continue;
        } 

        bfs();      
        int flag=f[x][y][0]; 
        if (flag!=1) printf("Case #%d: Yes\n",kcase);
        else printf("Case #%d: No\n",kcase);    
    } 
    return 0;
}