POI 2001 跳舞蝇

02年的论文题，论文里的方法看不懂QAQ。

然后07年的论文里有一个有向图的同构判定，于是拿来水一水。

结果……

呵呵呵呵呵呵呵呵呵

让我静一静TAT

我觉得肯定是评测机太慢了，我的复杂度是d*n*k的

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef unsigned long long ll;
const int a=566;
const int b=346;
const int c=489;
ll f1[2005][55],f2[2005][55],s[2005];
int p1[2005],p2[2005],n;
void solve(ll f[2005][55],int p[2005],int k){
	for(int i=1;i<=n;i++)s[i]=0;
	for(int i=1;i<=n;i++){
		f[i][k]=f[i][k-1]*a+f[p[i]][k-1]*b;
		s[p[i]]+=f[i][k-1];
	}
	for(int i=1;i<=n;i++)
	f[i][k]+=s[i]*c;
}
int t1[2005],t2[2005];
bool check(){
	for(int i=1;i<=n;i++)
	f1[i][0]=f2[i][0]=1;
	for(int i=1;i<=33;i++){
		solve(f1,p1,i);
		solve(f2,p2,i);
		for(int j=1;j<=n;j++){
			t1[j]=f1[j][i];
			t2[j]=f2[j][i];
		}
		sort(t1+1,t1+1+n);
		sort(t2+1,t2+1+n);
		for(int j=1;j<=n;j++)
		if(t1[j]!=t2[j])return 0;
	}
	return 1;
}
int main(){
	//freopen("a.in","r",stdin);
	int T;scanf("%d",&T);
	while(T--){
		 scanf("%d",&n);
		 for(int i=1;i<=n;i++)scanf("%d",&p1[i]);
		 for(int i=1;i<=n;i++)scanf("%d",&p2[i]); 
		 puts(check()?"T":"N");
	}
	return 0;
}

昨天睡觉的时候突然想起来这个是基环内向树哎。

既然是这样的话我们当然不能用普通的图的hash了。

首先考虑树的部分，当然是有根树的hash了，这个很好搞。

然后设有根树的根在环上。

当然我们知道每条链的尽头都是环（自环也是环）

于是我们想一下环可以怎么hash。

环？循环？字符串的循环同构？（不要问我怎么想到的）

于是我们可以把环看成一个串，串上每个字符就是这个点所连的树的hash值。

然后对每个环求hash值。

最后排序一下比较所有环的hash值（一一对应）

啦啦啦啦啦啦啦啦啦了~

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
typedef unsigned long long ll;
const int a=233;
const int b=213;
const int c=417;
const int N=2000+5;
vector<ll>v[N];
ll w[N],xa[N];
int p[N],du[N],n;
void toposort(){
	queue<int>q;
	for(int i=1;i<=n;i++)
	if(!du[i])q.push(i);
	while(!q.empty()){
		int u=q.front();q.pop();
		ll ans=0;
		sort(v[u].begin(),v[u].end());
		for(int i=0;i<v[u].size();i++)
		ans=ans*a+v[u][i];
		w[u]=ans+xa[v[u].size()+1];
		du[p[u]]--;
		v[p[u]].push_back(w[u]);
		if(!du[p[u]])q.push(p[u]);
	}
	for(int i=1;i<=n;i++)
	if(du[i]){
		ll ans=0;
		sort(v[i].begin(),v[i].end());
		for(int j=0;j<v[i].size();j++)
		ans=ans*a+v[i][j];
		w[i]=ans+xa[v[i].size()+1];
		w[i]+=c;
	}
}
void findcircle(vector<ll>& g){
	g.clear();
	for(int i=1;i<=n;i++)
	if(du[i]){
		vector<ll>f;
		int j=i;
		do{
			f.push_back(w[j]);
			j=p[j];
		}while(j!=i);
		sort(f.begin(),f.end());
		ll ans=0;
		for(j=0;j<f.size();j++)
		ans=ans*b+f[j];
		g.push_back(ans);
	}
	sort(g.begin(),g.end());
}
vector<ll>f,g;
void init(){
	for(int i=1;i<=n;i++){
		v[i].clear();
		du[i]=0;
	}
}
int main(){
	//freopen("a.in","r",stdin);
	int T;scanf("%d",&T);
	xa[0]=1;
	for(int i=1;i<=2000;i++)xa[i]=xa[i-1]*a;
	while(T--){
		scanf("%d",&n);
		init();
		for(int i=1;i<=n;i++){
			scanf("%d",&p[i]);
			du[p[i]]++;
		}
		toposort();findcircle(f);
		init();
		for(int i=1;i<=n;i++){
			scanf("%d",&p[i]);
			du[p[i]]++;
		}
		toposort();findcircle(g);
		puts(f==g?"T":"N");
	}
	return 0;
}