Atcoder agc027F

先考虑至少有一个点未操作的情况。我们枚举一个未操作的点$rt$，令$A$和$B$以$rt$为根，那么需要操作的是在$A$和$B$中父亲不同的点，且合法必须满足未操作的点在两棵树分别形成与$rt$相连的连通块。如何判定是否可行呢？考虑我们操作一个点的时候它必须是叶子，于是对于一个需要操作的点，它在$A$中的儿子必须比它先操作，在$B$中的儿子必须比它后操作，用拓扑排序判定即可。
再考虑所有点都操作过的情况。枚举最后被操作的点$x$，它在$B$中是叶子，且它在操作前也是叶子，于是可以枚举它操作前所连的点的编号$y$，这样就知道了最后一次操作前树的形态，转化为$x$未操作的情况。
时间复杂度$\mathcal{O}(T{N}^3)$。

#include <bits/stdc++.h>
#define inf 0x3f3f3f3f

using namespace std;

vector <int> e1[55],e2[55],e3[55];
int d1[55],d2[55];

int p1[55],p2[55];

void dfs1(int x) {
  for(int i=0;i<e2[x].size();i++)
    if (e2[x][i]!=p1[x]) {
    	int u=e2[x][i];
    	p1[u]=x;
    	dfs1(u);
	}
}

void dfs2(int x) {
  for(int i=0;i<e3[x].size();i++)
    if (e3[x][i]!=p2[x]) {
    	int u=e3[x][i];
    	p2[u]=x;
    	dfs2(u);
	}
}

vector <int> ee[55];
int deg[55];
queue <int> q;

int solve(int n,int rt) {
  p1[rt]=p2[rt]=0;
  dfs1(rt);
  dfs2(rt);
  int s=0;
  for(int i=1;i<=n;i++) ee[i].clear();
  memset(deg,0,sizeof(deg));
  for(int i=1;i<=n;i++)
    if (p1[i]==p2[i]&&p1[p1[i]]!=p2[p1[i]]) return inf;
  for(int i=1;i<=n;i++)
    if (p1[i]!=p2[i]) {
    	s++;
    	for(int j=0;j<e2[i].size();j++)
    	  if (e2[i][j]!=p1[i]) {
		    ee[e2[i][j]].push_back(i);
		    deg[i]++;
		  }
		for(int j=0;j<e3[i].size();j++)
		  if (e3[i][j]!=p2[i]) {
		    ee[i].push_back(e3[i][j]); 
		    deg[e3[i][j]]++;
		  }
	}
  for(int i=1;i<=n;i++)
    if (!deg[i]) q.push(i);
  while (!q.empty()) {
  	int x=q.front();q.pop();
  	for(int i=0;i<ee[x].size();i++) {
  		int u=ee[x][i];
  		deg[u]--;
  		if (!deg[u]) q.push(u);
	  }
  }
  for(int i=1;i<=n;i++)
    if (deg[i]) return inf;
  return s;
}

int main() {
  int cases;
  scanf("%d",&cases);
  for(;cases;cases--) {
  	memset(d1,0,sizeof(d1));
  	memset(d2,0,sizeof(d2));
  	int n;
  	scanf("%d",&n);
  	for(int i=1;i<=n;i++) {
  		e1[i].clear();
  		e2[i].clear();
	  }
	for(int i=1;i<n;i++) {
		int x,y;
		scanf("%d%d",&x,&y);
		e2[x].push_back(y);
		e2[y].push_back(x);
		d1[x]++;d1[y]++;
	}
	for(int i=1;i<n;i++) {
		int x,y;
		scanf("%d%d",&x,&y);
		e1[x].push_back(y);
		e1[y].push_back(x);
		d2[x]++;d2[y]++;
	}
	for(int i=1;i<=n;i++) e3[i]=e1[i];
	int ans=inf;
	for(int i=1;i<=n;i++)
	  ans=min(ans,solve(n,i));
	if (ans==inf) {
		for(int i=1;i<=n;i++)
		  if (d2[i]==1) {
		  	for(int j=1;j<=n;j++)
		  	  if (j!=i) {
		  	  	  for(int k=1;k<=n;k++) 
					if (k!=i) {
		  	  	  	    e3[k].clear();
		  	  	  	    for(int l=0;l<e1[k].size();l++)
		  	  	  	      if (e1[k][l]!=i) e3[k].push_back(e1[k][l]);
					  }
				  e3[j].push_back(i);
				  e3[i].clear();
				  e3[i].push_back(j);
				  ans=min(ans,solve(n,i)+1);
				}
		  }
	}
	printf("%d\n",(ans<inf)?ans:-1);
  }
  return 0;
}