[bzoj2140][强联通]稳定婚姻

Description

我国的离婚率连续7年上升，今年的头两季，平均每天有近5000对夫妇离婚，大城市的离婚率上升最快，有研究婚
姻问题的专家认为，是与简化离婚手续有关。25岁的姗姗和男友谈恋爱半年就结婚，结婚不到两个月就离婚，是典
型的“闪婚闪离”例子，而离婚的导火线是两个人争玩电脑游戏，丈夫一气之下，把电脑炸烂。有社会工作者就表
示，80后求助个案越来越多，有些是与父母过多干预有关。而根据民政部的统计，中国离婚五大城市首位是北京，
其次是上海、深圳，广州和厦门，那么到底是什么原因导致我国成为离婚大国呢?有专家分析说，中国经济急速发
展，加上女性越来越来越独立，另外，近年来简化离婚手续是其中一大原因。——以上内容摘自第一视频门户现代
生活给人们施加的压力越来越大，离婚率的不断升高已成为现代社会的一大问题。而其中有许许多多的个案是由婚
姻中的“不安定因素”引起的。妻子与丈夫吵架后，心如绞痛，于是寻求前男友的安慰，进而夫妻矛盾激化，最终
以离婚收场，类似上述的案例数不胜数。我们已知n对夫妻的婚姻状况，称第i对夫妻的男方为Bi，女方为Gi。若某
男Bi与某女Gj曾经交往过（无论是大学，高中，亦或是幼儿园阶段，i≠j），则当某方与其配偶（即Bi与Gi或Bj与
Gj）感情出现问题时，他们有私奔的可能性。不妨设Bi和其配偶Gi感情不和，于是Bi和Gj旧情复燃，进而Bj因被戴
绿帽而感到不爽，联系上了他的初恋情人Gk……一串串的离婚事件像多米诺骨牌一般接踵而至。若在Bi和Gi离婚的
前提下，这2n个人最终依然能够结合成n对情侣，那么我们称婚姻i为不安全的，否则婚姻i就是安全的。给定所需
信息，你的任务是判断每对婚姻是否安全。 Input

第一行为一个正整数n，表示夫妻的对数； 以下n行，每行包含两个字符串，表示这n对夫妻的姓名（先女后男），由一个空格隔开；
第n+2行包含一个正整数m，表示曾经相互喜欢过的情侣对数；
以下m行，每行包含两个字符串，表示这m对相互喜欢过的情侣姓名（先女后男），由一个空格隔开。
所有姓名字符串中只包含英文大小写字母，大小写敏感，长度不大于8，
保证每对关系只在输入文件中出现一次，输入文件的最后m行不会出现未在之前出现过的姓名， 这2n个人的姓名各不相同，
1≤n≤4000，0≤m≤20000。

Output

输出文件共包含n行，第i行为“Safe”（如果婚姻i是安全的）或“Unsafe”（如果婚姻i是不安全的）。

Sample Input

【样例输入1】

2

Melanie Ashley

Scarlett Charles

1

Scarlett Ashley

【样例输入2】

2

Melanie Ashley

Scarlett Charles

2

Scarlett Ashley

Melanie Charles

Sample Output

【样例输出1】

Safe

Safe

【样例输出2】

Unsafe

Unsafe

题解

裸跑二分图匹配的话，很明显是mn^2的。
考虑一个这样的建图
现任的时候，男连向女
前任的时候，女连向男
跑一边tarjan缩点看看两任在不在一个连通块里，在的话就一定是Unsafe
证明：
由于有向边，同时连通块中的边一定是男->女->男…，且为现任->前任->现任…，类似于二分图匹配模型。那么假如在一个连通块时，可以保证至少有一个人能选择到前任啦

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node{int x,y,next;}a[211000];int len,last[8100];
void ins(int x,int y){len++;a[len].x=x;a[len].y=y;a[len].next=last[x];last[x]=len;}
struct Tire{int son[75],s;Tire(){memset(son,-1,sizeof(son));s=0;}}tr[75000];int trlen,cnt;
char ch[25];int root;
void add()
{
    int len=strlen(ch+1);int p=root;
    for(int i=1;i<=len;i++)
    {
        int y=ch[i]-'A'+1;
        if(tr[p].son[y]==-1)tr[p].son[y]=++trlen;
        p=tr[p].son[y];
    }
    tr[p].s=++cnt;
}
int fd(){int len=strlen(ch+1);int p=root;for(int i=1;i<=len;i++)p=tr[p].son[ch[i]-'A'+1];return tr[p].s;}
int belong[8100],dfn[8100],low[8100],id,block,tp;
int sta[8100];bool v[8100];
void dfs(int x)
{
    dfn[x]=low[x]=++id;sta[++tp]=x;
    v[x]=true;
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(dfn[y]==-1)
        {
            dfs(y);
            low[x]=min(low[x],low[y]);
        }
        else
        {
            if(v[y])low[x]=min(low[x],dfn[y]);
        }
    }
    if(dfn[x]==low[x])
    {
        block++;int i;
        do
        {
            i=sta[tp--];
            v[i]=false;
            belong[i]=block;
        }while(i!=x);
    }
}
int n,m;
int main()
{
    scanf("%d",&n);
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=n;i++)
    {
        scanf("%s",ch+1);add();int tmpx=cnt;
        scanf("%s",ch+1);add();int tmpy=cnt;
        ins(tmpx,tmpy);
    }
    scanf("%d",&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%s",ch+1);int tmpx=fd();
        scanf("%s",ch+1);int tmpy=fd();
        ins(tmpy,tmpx);
    }
    id=block=tp=0;
    memset(dfn,-1,sizeof(dfn));
    memset(v,false,sizeof(v));
    for(int i=1;i<=cnt;i++)if(dfn[i]==-1)dfs(i);
    for(int i=1;i<=n;i++)
    {
        if(belong[2*i-1]==belong[2*i])printf("Unsafe\n");
        else printf("Safe\n");
    }
    return 0;
}