[poi2009]gas

（ps：某牛说是这题陈题= =我太孤陋寡闻了。。竟然都不知道是陈题。。。）

 

 

这题这个s的限制非常要命。

如果没有s会很好做。

先按照深度排序。

按顺序枚举

此时选择最远的点（k-th祖先）必然最优。因为按照深度排完序了，所以最远的点可以覆盖所有它子孙节点的。比较显然，

对于这个节点 选择可以覆盖到的s个最深的覆盖。。

 

 

一个节点可以选择一个建好的消防站或者离他最远的节点新建一个消防站

因为已经按照深度排完序了。所以那么一个已经建好的消防站覆盖的节点没有新的覆盖得多。所以我们就利用已经建好的消防站而不新建。

大多数人是直接选择s个最深的覆盖，- -表示我不知道该如何这样实现，感觉复杂度很大。

 

于是我想了一个办法。。每个节点找可以覆盖到他的最深的节点。其实本质是一样的。只是实现手段有差

 

实现上

给每个节点标记。
f[i,k]
表示该节点i可以拓展k的有多少个位置。。（没次新加一个消防站就给它和它的k个祖先都加上s）

 

然后枚举节点向上找k个祖先。。看看哪个最深的。然后从那个点出发删掉在上下长度共为k的链上各删去1.

（无论删哪条都可以（因为已经排好序了，不会影响结果，深的可以选择的，浅的选择更好）这点上我纠结了很久）

 

 

由于我的思维有限，虽然过了全部24个data，但是还是不太敢证明这样是对的。。

 

 

 

 

（变量名还是严谨点。。k不要做为主要变量。还有，我一开始的深度判断出问题了）

 

program p1117; const inf = maxlongint shr 1; var q,d,dep,first,fa,num :array[0..200000]of longint; f :Array[0..200000,0..23]of longint; v :array[0..200000]of boolean; bian :array[0..500000]of record y,next:longint;end; tot,top,ans,n,s,k :longint; procedure connect(a,b:longint); begin inc(tot); with bian[tot] do begin y:=b; next:=first[a]; first[a]:=tot; end; end; procedure dfs(x,now:longint); var k:longint; begin v[x]:=true; dep[x]:=now; k:=first[x]; while k<>0 do with bian[k] do begin if not v[y] then begin fa[y]:=x; dfs(y,now+1); end; k:=next; end; end; procedure qsort(l,r:longint); var i,j,mid,tmp :longint; begin i:=l;j:=r; mid:=dep[(i+j) shr 1]; repeat while dep[i]>mid do inc(i); while dep[j]<mid do dec(j); if i<=j then begin tmp:=dep[i];dep[i]:=dep[j];dep[j]:=tmp; tmp:=num[i];num[i]:=num[j];num[j]:=tmp; inc(i);dec(j); end; until i>j; if l<j then qsort(l,j); if i<r then qsort(i,r); end; procedure init; var i,x,y,root :longint; begin readln(n,s,k); for i:=1 to n-1 do begin readln(x,y); num[i]:=i; inc(d[x]); inc(d[y]); connect(x,y); connect(y,x); end; num[n]:=n; for i:=1 to n do if d[i]=1 then begin root:=i; break; end; dfs(root,0); qsort(1,n); end; procedure find(x:longint); begin top:=0; while x<>0 do begin inc(top); q[top]:=x; if top>k then break; x:=fa[x]; end; end; procedure clear(x,yt:longint); var xx,yy,kk:longint; begin xx:=x; yy:=yt; while yy>=0 do begin dec(f[xx,yy]); xx:=fa[xx]; if xx=0 then break; dec(yy); end; xx:=x; yy:=yt+1; while yy<=k do begin kk:=first[xx]; while kk<>0 do with bian[kk] do begin if (fa[y]=xx)and(f[y,yy]>0) then begin dec(f[y,yy]); xx:=y; break; end; kk:=next; end; inc(yy); end; end; procedure mark(x,num:longint); begin while (x<>0)and(num>=0) do begin inc(f[x,num],s-1); x:=fa[x]; dec(num); end; end; procedure main; var ans,kk,i,j,tmp_i,tmp_j:longint; begin ans:=0; for kk:=1 to n do begin find(num[kk]); tmp_i:=inf; tmp_j:=-1; for i:=1 to top do begin for j:=i-1 to k+1 do if f[q[i],j]>0 then break; if j<=k then if (k-j-i>k-tmp_i-tmp_j)or(tmp_j=-1) then begin tmp_i:=i; tmp_j:=j; end; end; if tmp_j<>-1 then clear(q[tmp_i],tmp_j) else begin inc(ans); mark(q[top],k); end; end; writeln(ans); end; begin init; main; end.