NOI2003 逃学的小孩

这一题不会做啊……

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

题解请看2007年陈瑜希论文

代码：

const maxn=400100;
type node=record
     w,go,next:longint;
     end;
var  e:array[0..maxn] of node;
     f:array[0..maxn,1..3] of int64;
     first,u:array[0..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+2*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:=1 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,1] then
   begin
   f[x,3]:=f[x,2];
   f[x,2]:=f[x,1];
   f[x,1]:=y;
   end
 else
  if y>f[x,2] then
   begin
   f[x,3]:=f[x,2];
   f[x,2]:=y;
   end
 else
   f[x,3]:=max(f[x,3],y);
 end;
procedure dfs1(x,fa:longint);
 var i,y:longint;
 begin
 i:=first[x];
 while i<>0 do
  begin
  y:=e[i].go;
  if y<>fa then
   begin
   dfs1(y,x);
   u[y]:=e[i].w;
   update(x,f[y,1]+u[y]);
   end;
  i:=e[i].next;
  end;
 end;
procedure dfs2(x,fa:longint);
 var i,y:longint;
 begin
 if f[x,1]+u[x]=f[fa,1] then update(x,f[fa,2]+u[x])
 else update(x,f[fa,1]+u[x]);
 ans:=max(ans,get(f[x,1],f[x,2],f[x,3]));
 i:=first[x];
 while i<>0 do
  begin
  y:=e[i].go;
  if y<>fa then dfs2(y,x);
  i:=e[i].next;
  end;
 end;
procedure main;
 begin
 ans:=0;
 dfs1(1,0);
 dfs2(1,0);
 writeln(ans);
 end;
begin
 init;
 main;
end.     

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