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

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

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

Input

Output

Sample Input

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

Sample Output

1

对于60%的数据，n<=100。

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

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

#include <iostream>
#include <cstdio>
#include <memory.h>
#include <algorithm>
#include <queue>
using namespace std;
const int N=7e2+10;
const int P=76447;
struct Graph {int u,v,nx;
}g[2*N];
struct Cube {int u,v,nx,c,w;
}t[2*N*N+4*N];
int cnt,list[N],vis[2*N],dep[N],flow[2*N],dis[2*N],fa[2*N];
int g1,g2,sz[N],mxsize;
int h[N],xx[N],f[N][N],p[N];
int n,from[N],to[N],ans;void Add(int u,int v) {g[++cnt]=(Graph){u,v,list[u]};list[u]=cnt;g[++cnt]=(Graph){v,u,list[v]};list[v]=cnt;
}void Add_Cube(int u,int v,int w,int c) {t[++cnt]=(Cube){u,v,flow[u],c,w};flow[u]=cnt;t[++cnt]=(Cube){v,u,flow[v],0,-w};flow[v]=cnt;
}void Find_G(int u,int f) {int size=0;sz[u]=1;for (int i=list[u];i;i=g[i].nx)if (g[i].v!=f) {Find_G(g[i].v,u);sz[u]+=sz[g[i].v];size=max(size,sz[g[i].v]);}size=max(size,n-sz[u]);if (size==mxsize) g2=u;if (size<mxsize) {mxsize=size;g1=u;g2=0;}
}void Tree_Hash(int u,int f) {int xcnt=0;dep[u]=dep[f]+1;for (int i=list[u];i;i=g[i].nx) if (g[i].v!=f) Tree_Hash(g[i].v,u);for (int i=list[u];i;i=g[i].nx) if (g[i].v!=f) xx[++xcnt]=h[g[i].v];sort(xx+1,xx+xcnt+1);h[u]=17731;for (int i=1;i<=xcnt;i++) h[u]=(h[u]*7441%P+xx[i])%P;h[u]=h[u]*3119%P;
}bool CMP(int u,int v) {return dep[u]!=dep[v]?dep[u]>dep[v]:h[u]<h[v];
}bool SPFA() {queue<int> q;while (!q.empty()) q.pop();memset(dis,0x3f,sizeof dis);dis[0]=0;q.push(0);while (!q.empty()) {int u=q.front();q.pop();for (int i=flow[u];i;i=t[i].nx)if (dis[t[i].v]>dis[u]+t[i].w&&t[i].c) {dis[t[i].v]=dis[u]+t[i].w;fa[t[i].v]=i;if (!vis[t[i].v]) q.push(t[i].v);vis[t[i].v]=1;}vis[u]=0;}return dis[2*n+1]!=0x3f3f3f3f;
}void Get_Ans() {int x=2*n+1,mf=2147483647;while (fa[x]) {ans+=t[fa[x]].w;mf=min(mf,t[fa[x]].c);x=t[fa[x]].u;}x=2*n+1;while (fa[x]) {t[fa[x]].c-=mf;t[fa[x]^1].c+=mf;x=t[fa[x]].u;}
}int Min_Cost_Flow() {ans=0;while (SPFA()) Get_Ans();return ans;
}void DP() {for (int i=1;i<=n;i++) p[i]=i;sort(p+1,p+n+1,CMP);for (int i=1;i<=n;) {int j=i;for (;j<=n;j++)if (dep[p[i]]!=dep[p[j]]||h[p[i]]!=h[p[j]]) break;for (int k=i;k<j;k++)for (int l=i;l<j;l++) {int u=p[k],v=p[l];cnt=1;memset(flow,0,sizeof flow);for (int sv=list[v];sv;sv=g[sv].nx)if (dep[g[sv].v]==dep[v]+1)Add_Cube(n+g[sv].v,2*n+1,0,1);for (int su=list[u];su;su=g[su].nx)if (dep[g[su].v]==dep[u]+1) {Add_Cube(0,g[su].v,0,1);for (int sv=list[v];sv;sv=g[sv].nx)if (dep[g[sv].v]==dep[v]+1&&h[g[su].v]==h[g[sv].v])Add_Cube(g[su].v,n+g[sv].v,f[g[su].v][g[sv].v],1);}f[u][v]=Min_Cost_Flow()+(from[u]!=to[v]);}i=j;}
}int main() {scanf("%d",&n);for (int i=1,u,v;i<n;i++) scanf("%d%d",&u,&v),Add(u,v);for (int i=1;i<=n;i++) scanf("%d",&from[i]);for (int i=1;i<=n;i++) scanf("%d",&to[i]);mxsize=n;Find_G(1,0);if (g2) {cnt=0;memset(list,0,sizeof list);for (int i=1,j=1;i<=n-1;i++,j=i*2-1)if ((g[j].u!=g1||g[j].v!=g2)&&(g[j].v!=g1||g[j].u!=g2))Add(g[j].u,g[j].v);Add(g1,n+1),Add(n+1,g2);n++;g1=n;}Tree_Hash(g1,0);DP();printf("%d",f[g1][g1]);
}