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;
}