【BZOJ2140】稳定婚姻 Tarjan

【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对相互喜欢过的情侣姓名（先女后男），由一个空格隔开。

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
【数据规模和约定】
对于100%的数据，所有姓名字符串中只包含英文大小写字母，大小写敏感，长度不大于8，保证每对关系只在输入文件中出现一次，输入文件的最后m行不会出现未在之前出现过的姓名，这2n个人的姓名各不相同，1≤n≤4000，0≤m≤20000。

HINT

题解：看到数据范围想要直接Tarjan缩点看一下每对夫妻是否在一个边双里？too naive，下面有一个反例

虽然B1G1,B3G3在同一个边双里，但是显然他们是安全的。所以我们不能将问题转化成一个无向图中的模型，而是人为的为这些关系定向。

令原配中的边从女方指向男方，旧情中的边从男方指向女方，此时再找强联通分量就行了。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <map>
using namespace std;
typedef unsigned long long ull;
map<ull,int> mp;
int n,m,cnt,sum,tot,top;
int to[50000],next[50000],head[10000],dep[10000],low[10000],ins[10000],sta[10000],bel[10000];
char str[20];
ull gethash()
{
	scanf("%s",str);
	ull hash=0;
	int i,len=strlen(str);
	for(i=0;i<len;i++)	hash=hash*131+str[i];
	return hash;
}
void tarjan(int x)
{
	dep[x]=low[x]=++tot,ins[x]=1,sta[++top]=x;
	for(int i=head[x];i!=-1;i=next[i])
	{
		if(!dep[to[i]])	tarjan(to[i]),low[x]=min(low[x],low[to[i]]);
		else	if(ins[to[i]])	low[x]=min(low[x],dep[to[i]]);
	}
	if(dep[x]==low[x])
	{
		sum++;
		int t;
		do
		{
			t=sta[top--],ins[t]=0,bel[t]=sum;
		}while(t!=x);
	}
}
void add(int a,int b)
{
	to[cnt]=b,next[cnt]=head[a],head[a]=cnt++;
}
int main()
{
	scanf("%d",&n);
	int i,a,b;
	memset(head,-1,sizeof(head));
	for(i=1;i<=n;i++)
	{
		mp[gethash()]=i*2-1;
		mp[gethash()]=i*2;
		add(i*2,i*2-1);
	}
	scanf("%d",&m);
	for(i=1;i<=m;i++)
	{
		a=mp[gethash()],b=mp[gethash()];
		if(a&1)	add(a,b);
		else	add(b,a);
	}
	for(i=1;i<=2*n;i++)	if(!dep[i])	tarjan(i);
	for(i=1;i<=n;i++)
	{
		if(bel[i*2-1]==bel[i*2])	printf("Unsafe\n");
		else	printf("Safe\n");
	}
	return 0;
}