【树形DP】爱心蜗牛

 

爱心蜗牛 

提交文件：badnews.pas/c/cpp

输入文件：badnews.in

输出文件：badnews.out 

猫猫把嘴伸进池子里，正准备“吸”鱼吃，却听到门铃响了。猫猫擦了擦脸上的水，打开门一看，那人正是她的好朋友——川川。

川川手里拿着一辆玩具汽车，对猫猫说：“这是我的新汽车！”接着，伴随一阵塑料叩击声，玩具汽车的车门竟开了，一只蜗牛慢慢吞吞爬了出来。

“哇！这么大的蜗牛……”猫猫惊讶道。

“这是我的宠物蜗牛，他叫点点。”川川介绍道。

“把他送给我好吗？”猫猫央求道。

“可以让他陪你几天，但是不能送给你……”

点点沿着川川的身体，爬到了地上，又移到了猫猫的池子旁边，只听见猫猫向川川介绍她的“创意吃鱼法”，心里不禁起了一丝凉意：“这个女生太毒了……吃鱼前还要玩鱼……”转眼一看，池中的鱼依旧畅快地游来游去。

“或许这些鱼听不懂猫语吧……好在我会一点儿猫语，也会一点鱼语……阿弥陀佛，善哉善哉。我还是救救这些鱼吧……”点点自言自语，一边费力地移动着身躯。他认识到——单凭自己的力量，把猫猫的阴谋告诉每一条鱼，似乎不太可能——自己底盘太低，走不快，看来只得想其他办法来传达信息。一翻认真思考之后，点点想到，如果把猫猫的计划告诉其中一条鱼，再让鱼们互相传达消息，那么在相对较短的时间内，每条鱼都会得知猫猫的计划。

鱼们的社会等级森严，除了国王菜鱼之外，每条鱼均有且只有一个直接上级，菜鱼则没有上级。如果A是B的上级，B是C的上级，那么A就是C的上级。

绝对不会出现这样两条鱼A、B：A是B的上级，B也是A的上级。

　　最开始的时刻是0，点点要做的，就只是用1单位的时间把猫猫的阴谋告诉某一条“信息源鱼”，让鱼们自行散布消息。在任意一个时间单位中，任何一条已经接到通知的鱼，都可以把消息告诉他的一个直接上级或者直接下属。

现在，点点想知道：

1.到底需要多长时间，消息才能传遍池子里的鱼？

2.使消息传递过程消耗的时间最短，可供选择的“信息源鱼”有那些？ 

输入数据：

第一行有一个数N（N≤1000），表示池中的鱼数，池鱼按照1到n编上了号码，国王菜鱼的标号是1。

第2行到第N行（共N-1行），每一行一个数，第I行的数表示鱼I的直接上级的标号。 

输出数据：

　　对于每组输入数据：

　　第一行有一个数，表示最后一条鱼接到通知的最早时间。

　　第二行有若干个数，表示可供选择的鱼F的标号，按照标号从小到大的顺序输出，中间用空格分开。 

输入样例：

8

1

1

3

4

4

4

3 

输出样例：

5

3 4 5 6 7

==========================

===================

type
  pnode=^node;
  node=record
         x:longint;
         next:pnode;
       end;
var
  n:longint;
  map:array[1..1000]of pnode;
  ans:longint;
  ansi:array[1..1000]of longint;
  f_bo:array[1..1000]of boolean;
  c:array[1..1000,1..1000]of longint;
  cs:array[1..1000]of longint;
  anss:longint;
procedure init;
begin
  assign(input,'badnews.in');
  assign(output,'badnews.out');
  reset(input); rewrite(output);
end;

procedure terminate;
begin
  close(input); close(output);
  halt;
end;

procedure insert(a,b:longint);
var
  p:pnode;
begin
  new(p);
  p^.next:=map[a];
  p^.x:=b;
  map[a]:=p;
end;

procedure qsort(s,t,n:longint);
var
  i,j:longint;
  tem:longint;
  x:longint;
begin
  i:=s; j:=t;
  x:=c[n,(s+t)shr 1];
  repeat
    while x<c[n,j] do dec(j);
    while c[n,i]<x do inc(i);
    if i<=j then
      begin
        tem:=c[n,j];
        c[n,j]:=c[n,i];
        c[n,i]:=tem;
        inc(i); dec(j);
      end;
  until i>j;
  if i<t then qsort(i,t,n);
  if s<j then qsort(s,j,n);
end;

function dfs(t:longint):longint;
var
  p:pnode;
  i:longint;
begin
  dfs:=0;
  f_bo[t]:=false; //标记当前点
  p:=map[t];
  cs[t]:=0;
  //表示t号节点有多少个元素..
  while p<>nil do
    begin
      if f_bo[p^.x] then
        begin
          inc(cs[t]);
          c[t,cs[t]]:=dfs(p^.x);
        end;
      p:=p^.next;
    end;
  if cs[t]=0 then exit(0);
  qsort(1,cs[t],t);
  for i:=1 to cs[t] do
    if dfs<c[t,i]+cs[t]+1-i then dfs:=c[t,i]+cs[t]+1-i;
  exit(dfs);
end;

procedure main;
var
  i:longint;
  x:longint;
  t:longint;
begin
  readln(n);
  for i:=1 to n do map[i]:=nil;
  for i:=2 to n do
    begin
      read(x);
      insert(i,x);
      insert(x,i);
    end;
  ans:=maxlongint;
  for i:=1 to n do
    begin
      fillchar(f_bo,sizeof(f_bo),true);
      t:=dfs(i);
      if t<ans then
        begin
          ans:=t;
          anss:=1;
          ansi[anss]:=i;
        end
        else
        if t=ans then
          begin
            inc(anss);
            ansi[anss]:=i;
          end;
    end;
  writeln(ans+1);
  for i:=1 to anss do write(ansi[i],' ');
end;

begin
  init;
  main;
  terminate;
end.