BZOJ2282 SDOI2011消防/NOIP2007树网的核（二分答案+树形dp）

　　要求最大值最小容易想到二分答案。首先对每个点求出子树中与其最远的距离是多少，二分答案后就可以标记上一些必须在所选择路径中的点，并且这些点是不应存在祖先关系的。那么如果剩下的点数量>=3，显然该答案不可行；=0，显然可行；=1，由该点沿其到根的路径往上爬，并计算最远距离判断是否合法；=2，求出两点lca后与1的做法类似。

　　noip原题是要求该路径在一条直径上，事实上这条最优路径一定是在直径上的，不过并不太懂这之间的关系。

　　写的太丑于是常数巨大，bzoj上愉快的T掉了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 500010
int n,m,p[N],deep[N],up[N],f[N],fa[N],t=,ans=,cnt,q[];
struct data{int to,nxt,len;
}edge[N<<];
void addedge(int x,int y,int z){t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].len=z,p[x]=t;}
void dfs(int k)
{
    for (int i=p[k];i;i=edge[i].nxt)
    if (edge[i].to!=fa[k])
    {
        fa[edge[i].to]=k;
        deep[edge[i].to]=deep[k]+;
        up[edge[i].to]=edge[i].len;
        dfs(edge[i].to);
        f[k]=max(f[k],f[edge[i].to]+edge[i].len);
    }
}
void paint(int k,int x)
{
    if (cnt>) return;
    for (int i=p[k];i;i=edge[i].nxt)
    if (edge[i].to!=fa[k])
    {
        paint(edge[i].to,x);
        if (cnt>) break;
        if (f[edge[i].to]+edge[i].len>x&&f[edge[i].to]<=x)
        q[++cnt]=edge[i].to;
        if (cnt>) break;
    }
}
bool check(int k)
{
    cnt=;
    paint(,k);
    if (cnt>) return ;
    if (cnt==)
    {
        int x=q[],y=q[],len=;
        while (x!=y)
        {
            if (deep[x]<deep[y]) swap(x,y);
            len+=up[x],x=fa[x];
        }
        if (len>m) return ;
        len=;
        while (fa[x])
        {
            y=x,len+=up[x],x=fa[x];
            if (len>k) return ;
            for (int i=p[x];i;i=edge[i].nxt)
            if (edge[i].to!=y&&edge[i].to!=fa[x]&&f[edge[i].to]+edge[i].len+len>k) return ;
        }
    }
    if (cnt==)
    {
        int x=q[],y,len=,flag=;
        while (fa[x])
        {
            y=x,len+=up[x],x=fa[x];
            if (!flag&&len>m) flag=,len=up[y];
            if (len*flag>k) return ;
            for (int i=p[x];i;i=edge[i].nxt)
            if (edge[i].to!=y&&edge[i].to!=fa[x]&&f[edge[i].to]+edge[i].len+len*flag>k) return ;
        }
    }
    return ;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj2282.in","r",stdin);
    freopen("bzoj2282.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    int l=,r=;
    for (int i=;i<n;i++)
    {
        int x=read(),y=read(),z=read();r+=z;
        addedge(x,y,z),addedge(y,x,z);
    }
    dfs();
    while (l<=r)
    {
        int mid=l+r>>;
        if (check(mid)) ans=mid,r=mid-;
        else l=mid+;
    }
    cout<<ans;
    return ;
}