bzoj4042

比较好的树形dp，涉及到树上路径的题目，我们往往考虑对路径分类

当我们考虑以x为根的子树，有这样几类路径

1. 起点终点都在子树内

2. 一个点延伸到子树外

对于要选择另一个点在子树外的路径，要建立在不破坏儿子子树内的路径基础上

因为破坏会破坏多条，而只能多选择一条，不合适

因此我们考虑树dp，设f[x]为路径起点终点在子树内能选择的最多路径数目，d[x]表示还没用过的另一个点在x子树外的路径的集合

我们知道，对于x的每个孩子d，d[y]中最多只会有1条路径被选中，因此我们可以用状压dp解决

 type node=record
        po,next:longint;
      end;

 var a:array[..,..] of longint;
     w:array[..,..] of boolean;
     e:array[..] of node;
     g,f,q,b,s,p:array[..] of longint;
     len,i,tt,n,m,x,y:longint;

 function max(a,b:longint):longint;
   begin
     if a>b then exit(a) else exit(b);
   end;

 procedure add(x,y:longint);
   begin
     inc(len);
     e[len].po:=y;
     e[len].next:=p[x];
     p[x]:=len;
   end;

 function check(x,y:longint):boolean;
   var i,j:longint;
   begin
     for i:= to s[x] do
       for j:= to s[y] do
         if w[a[x,i],a[y,j]] then exit(true);
     exit(false);
   end;

 procedure dfs(x,fa:longint);
   var st,mx,i,y,t,j,k:longint;
   begin
     f[x]:=;
     s[x]:=;
     a[x,s[x]]:=x;
     i:=p[x];
     while i<> do
     begin
       y:=e[i].po;
       if y<>fa then
         dfs(y,x);
       i:=e[i].next;
     end;
     i:=p[x];
     t:=;
     while i<> do
     begin
       y:=e[i].po;
       if y<>fa then
       begin
         f[x]:=f[x]+f[y];
         if check(x,y) then inc(f[x])
         else begin
           inc(t);
           q[t]:=y;
         end;
       end;
       i:=e[i].next;
     end;
     k:=;
     for i:= to t- do
       for j:=i+ to t do
         if check(q[i],q[j]) then
         begin
           inc(k);
           b[k]:=( shl (i-)) or ( shl (j-));
         end;
     st:=;
     mx:=;
     for i:= to  shl t- do
     begin
       g[i]:=;
       for j:= to k do
         if b[j] and i=b[j] then
           g[i]:=max(g[i],g[i xor b[j]]+);
       if g[i]>mx then
       begin
         mx:=g[i];
         st:=;
       end;
       if g[i]=mx then st:=st or (( shl t-) xor i);
     end;
     f[x]:=f[x]+mx;
     for i:= to t do
       if ( shl (i-)) and st> then
       begin
         y:=q[i];
         for j:= to s[y] do
         begin
           inc(s[x]);
           a[x,s[x]]:=a[y,j];
         end;
       end;
   end;

 begin
   readln(tt);
   while tt> do
   begin
     dec(tt);
     len:=;
     fillchar(p,sizeof(p),);
     readln(n);
     for i:= to n- do
     begin
       readln(x,y);
       add(x,y);
       add(y,x);
     end;
     readln(m);
     fillchar(w,sizeof(w),false);
     for i:= to m do
     begin
       readln(x,y);
       w[x,y]:=true;
       w[y,x]:=true;
     end;
     dfs(,);
     writeln(f[]);
   end;
 end.