NOI2003 逃学的小孩

这一题不会做啊……

我觉得真要比赛的话我可能会随机上几万次，然后再用LCA求距离，更新最优值，等到快超时的时候输出答案……

题解请看2007年陈瑜希论文

代码：

 const maxn=;
 type node=record
      w,go,next:longint;
      end;
 var  e:array[..maxn] of node;
      f:array[..maxn,..] of int64;
      first,u:array[..maxn] of longint;
      i,x,y,z,n,tot,m:longint;
      ans:int64;
 function max(x,y:int64):int64;
  begin
  if x>y then exit(x) else exit(y);
  end;
 function get(x,y,z:int64):int64;
  begin
  get:=x+*y+z;
  end;
 procedure insert(x,y,z:longint);
  begin
  inc(tot);e[tot].go:=y;e[tot].w:=z;e[tot].next:=first[x];first[x]:=tot;
  end;
 procedure init;
  begin
  readln(n,m);
  for i:= to m do
   begin
   readln(x,y,z);
   insert(x,y,z);
   insert(y,x,z);
   end;
  end;
 procedure update(x,y:int64);
  begin
   if y>f[x,] then
    begin
    f[x,]:=f[x,];
    f[x,]:=f[x,];
    f[x,]:=y;
    end
  else
   if y>f[x,] then
    begin
    f[x,]:=f[x,];
    f[x,]:=y;
    end
  else
    f[x,]:=max(f[x,],y);
  end;
 procedure dfs1(x,fa:longint);
  var i,y:longint;
  begin
  i:=first[x];
  while i<> do
   begin
   y:=e[i].go;
   if y<>fa then
    begin
    dfs1(y,x);
    u[y]:=e[i].w;
    update(x,f[y,]+u[y]);
    end;
   i:=e[i].next;
   end;
  end;
 procedure dfs2(x,fa:longint);
  var i,y:longint;
  begin
  if f[x,]+u[x]=f[fa,] then update(x,f[fa,]+u[x])
  else update(x,f[fa,]+u[x]);
  ans:=max(ans,get(f[x,],f[x,],f[x,]));
  i:=first[x];
  while i<> do
   begin
   y:=e[i].go;
   if y<>fa then dfs2(y,x);
   i:=e[i].next;
   end;
  end;
 procedure main;
  begin
  ans:=;
  dfs1(,);
  dfs2(,);
  writeln(ans);
  end;
 begin
  init;
  main;
 end.     

唉，这种神题，我什么时候才能自己想到呢？