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[0..1200,0..1200] of longint;
    w:array[0..1010,0..1010] of boolean;
    e:array[0..2400] of node;
    g,f,q,b,s,p:array[0..1200] 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:=1 to s[x] do
      for j:=1 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]:=0;
    s[x]:=1;
    a[x,s[x]]:=x;
    i:=p[x];
    while i<>0 do
    begin
      y:=e[i].po;
      if y<>fa then
        dfs(y,x);
      i:=e[i].next;
    end;
    i:=p[x];
    t:=0;
    while i<>0 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:=0;
    for i:=1 to t-1 do
      for j:=i+1 to t do
        if check(q[i],q[j]) then
        begin
          inc(k);
          b[k]:=(1 shl (i-1)) or (1 shl (j-1));
        end;
    st:=0;
    mx:=0;
    for i:=0 to 1 shl t-1 do
    begin
      g[i]:=0;
      for j:=1 to k do
        if b[j] and i=b[j] then
          g[i]:=max(g[i],g[i xor b[j]]+1);
      if g[i]>mx then
      begin
        mx:=g[i];
        st:=0;
      end;
      if g[i]=mx then st:=st or ((1 shl t-1) xor i);
    end;
    f[x]:=f[x]+mx;
    for i:=1 to t do
      if (1 shl (i-1)) and st>0 then
      begin
        y:=q[i];
        for j:=1 to s[y] do
        begin
          inc(s[x]);
          a[x,s[x]]:=a[y,j];
        end;
      end;
  end;

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