[BZOJ3197][SDOI2013]刺客信条assassin

bzoj
luogu

Description

故事发生在1486 年的意大利,Ezio原本只是一个文艺复兴时期的贵族,后来因为家族成员受到圣殿骑士的杀害,决心成为一名刺客。最终,凭借着他的努力和出众的天赋,成为了杰出的刺客大师,他不仅是个身手敏捷的武林高手,飞檐走壁擅长各种暗杀术。刺客组织在他的带领下,为被剥削的平民声张正义,赶跑了原本统治意大利的圣殿骑士首领-教皇亚历山大六世。在他的一生中,经历了无数次惊心动魄、扣人心弦的探险和刺杀。

曾经有一次,为了寻找Altair留下的线索和装备,Ezio在佛罗伦萨中的刺客墓穴进行探索。这个刺客墓穴中有许多密室,且任何两个密室之间只存在一条唯一的路径。这些密室里都有一个刺客标记,他可以启动或者关闭该刺客标记。为了打开储存着线索和装备的储藏室,Ezio必须操作刺客标记来揭开古老的封印。要想解开这个封印,他需要通过改变某些刺客标记的启动情况,使得所有刺客标记与封印密码“看起来一样”。

在这里,“看起来一样”的定义是:存在一种“标记”密室与“密码”密室之间一一对应的关系,使得密室间的连接情况和启动情况相同(提示中有更详细解释)。幸运的是,在Ezio来到刺客墓穴之前,在Da Vinci 的帮助下,Ezio 已经得知了打开储藏室所需要的密码。

而你的任务则是帮助Ezio 找出达成目标所需要最少的改动标记次数。

Input

第一行给出一个整数n, 表示密室的个数,
第二行至第 n 行, 每行绐出两个整数 a 和 b, 表示第a个密室和第b个密室之间存在一条通道。
第 n+1行给出 n 个整数,分别表示当时每个密室的启动情况 (0表示关闭, 1表示启动)。
第 n+2行给出 n 个整数, 分别表示密码中每个密室的启动情况。

Output

输出只有一行,即输出最少改动标记次数

Sample Input

4
1 2
2 3
3 4
0 0 1 1
1 0 0 0

Sample Output

1

Hint

密室的编号是可以变的!将第三个密室关闭后,在当前标记和密码之间,存在1->4,2->3,3->2,4->1 的对应关系,重新编号后连接情况没有改变,且标记与密码对应。对于更一般的情况,存在一个1 到n 的置换P,使得对于任意密室之间的道路u-v,都一定有密码密室中的道路P(u)-P(v);如果不存在密室之间的道路u-v,则一定没有密码密室中的道路P(u)-P(v)。

对于100%的数据,n<=700,且每个密室至多与11个密室相通

sol

首先判断树同构可以使用树\(hash\)的方法,为了避免出错可以把很多东西\(hash\)进去,比如说子树\(size\),儿子个数等等。
然后就是转移了。我们设\(f_{i,j}\)表示\(i\)号点去与\(j\)号点对应匹配的最小代价,这里要求\(i,j\)同构,那么\(i,j\)的各个儿子也应该会分别对应同构。我们现在最关键的问题就是处理\(i,j\)的各个儿子的匹配关系,使得总代价最小。
我们按照深度由大到小处理,这样在处理\(i,j\)时,它们的儿子的答案就已经算好了。在双方同构的儿子之间连权值为答案的边,这样就成了二分图最小权匹配,直接使用费用流即可。

code

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
int gi(){
	int x=0,w=1;char ch=getchar();
	while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
	if (ch=='-') w=0,ch=getchar();
	while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
	return w?x:-x;
}
#define ull unsigned long long
const int N = 1005;
int n,col1[N],col2[N],id[N],f[N][N];
struct MCMF{
	struct edge{int to,nxt,w,cost;}a[N*N];
	int S,T,head[N],cnt,dis[N],vis[N],pe[N],ans;
	queue<int>Q;
	void init(){memset(head,0,sizeof(head));cnt=1;}
	void link(int u,int v,int w,int cost){
		a[++cnt]=(edge){v,head[u],w,cost};head[u]=cnt;
		a[++cnt]=(edge){u,head[v],0,-cost};head[v]=cnt;
	}
	bool spfa(){
		memset(dis,63,sizeof(dis));
		dis[S]=0;Q.push(S);
		while (!Q.empty()){
			int u=Q.front();Q.pop();
			for (int e=head[u];e;e=a[e].nxt){
				int v=a[e].to;
				if (a[e].w&&dis[v]>dis[u]+a[e].cost){
					dis[v]=dis[u]+a[e].cost;pe[v]=e;
					if (!vis[v]) vis[v]=1,Q.push(v);
				}
			}
			vis[u]=0;
		}
		if (dis[T]==dis[0]) return false;
		ans+=dis[T];
		for (int i=T;i!=S;i=a[pe[i]^1].to)
			--a[pe[i]].w,++a[pe[i]^1].w;
		return true;
	}
	int work(){ans=0;while(spfa());return ans;}
}G;
const ull base1 = 20020415;
const ull base2 = 20011118;
struct TREE{
	int to[N<<1],nxt[N<<1],head[N<<1],cnt;
	int root,sz[N],w[N],fa[N],dep[N];ull Hash[N];
	void link(int u,int v){
		to[++cnt]=v;nxt[cnt]=head[u];head[u]=cnt;
	}
	void getroot(int u,int f){
		sz[u]=1;
		for (int e=head[u];e;e=nxt[e])
			if (to[e]!=f){
				getroot(to[e],u);
				sz[u]+=sz[to[e]];
				w[u]=max(w[u],sz[to[e]]);
			}
		w[u]=max(w[u],n-sz[u]);
		if (w[u]<w[root]) root=u;
	}
	ull tmp[N];
	void dfs(int u,int f){
		fa[u]=f;dep[u]=dep[f]+1;sz[u]=1;
		for (int e=head[u];e;e=nxt[e])
			if (to[e]!=f) dfs(to[e],u),sz[u]+=sz[to[e]];
		int len=0;
		for (int e=head[u];e;e=nxt[e])
			if (to[e]!=f) tmp[++len]=Hash[to[e]];
		sort(tmp+1,tmp+len+1);
		Hash[u]=(base1*len)^sz[u];
		for (int i=1;i<=len;++i) Hash[u]=((Hash[u]*base2)^tmp[i])+tmp[i];
	}
	void work(){
		w[0]=n;
		getroot(1,0);getroot(root,0);
		for (int e=head[root],lst=0;e;lst=e,e=nxt[e])
			if (sz[to[e]]*2==n){
				++n;
				if (e==head[root]) head[root]=nxt[e];
				else nxt[lst]=nxt[e];
				for (int i=head[to[e]],Lst=0;i;Lst=i,i=nxt[i])
					if (to[i]==root){
						if (i==head[to[e]]) head[to[e]]=nxt[i];
						else nxt[Lst]=nxt[i];
						break;
					}
				link(n,root);link(root,n);link(n,to[e]);link(to[e],n);
				root=n;break;
			}
		dfs(root,0);
	}
	int cal(int u,int v){
		int s1[N],s2[N],t1=0,t2=0;
		for (int e=head[u];e;e=nxt[e])
			if (to[e]!=fa[u]) s1[++t1]=to[e];
		for (int e=head[v];e;e=nxt[e])
			if (to[e]!=fa[v]) s2[++t2]=to[e];
		G.init();G.S=t1+t2+1;G.T=G.S+1;
		for (int i=1;i<=t1;++i) G.link(G.S,i,1,0);
		for (int i=1;i<=t2;++i) G.link(i+t1,G.T,1,0);
		for (int i=1;i<=t1;++i)
			for (int j=1;j<=t2;++j)
				if (Hash[s1[i]]==Hash[s2[j]])
					G.link(i,j+t1,1,f[s1[i]][s2[j]]);
		return G.work()+(col1[u]!=col2[v]);
	}
}T;
bool cmp(int i,int j){
	return T.dep[i]==T.dep[j]?T.Hash[i]<T.Hash[j]:T.dep[i]>T.dep[j];
}
int main(){
	n=gi();
	for (int i=1;i<n;++i){
		int u=gi(),v=gi();
		T.link(u,v);T.link(v,u);
	}
	for (int i=1;i<=n;++i) col1[i]=gi();
	for (int i=1;i<=n;++i) col2[i]=gi();
	T.work();
	for (int i=1;i<=n;++i) id[i]=i;
	sort(id+1,id+n+1,cmp);
	for (int i=1,j;i<=n;i=j){
		j=i+1;
		while (j<=n&&T.dep[id[i]]==T.dep[id[j]]&&T.Hash[id[i]]==T.Hash[id[j]]) ++j;
		for (int k=i;k<j;++k)
			for (int l=i;l<j;++l)
				f[id[k]][id[l]]=T.cal(id[k],id[l]);
	}
	printf("%d\n",f[T.root][T.root]);return 0;
}
posted @ 2018-05-28 21:14  租酥雨  阅读(182)  评论(0编辑  收藏