[BZOJ3124]直径

Description

小Q最近学习了一些图论知识。根据课本，有如下定义。树：无回路且连通的无向图，每条边都有正整数的权值来表示其长度。如果一棵树有N个节点，可以证明其有且仅有N-1 条边。 路径：一棵树上，任意两个节点之间最多有一条简单路径。我们用 dis(a,b)
表示点a和点b的路径上各边长度之和。称dis(a,b)为a、b两个节点间的距离。  
 直径：一棵树上，最长的路径为树的直径。树的直径可能不是唯一的。 
现在小Q想知道，对于给定的一棵树，其直径的长度是多少，以及有多少条边满足所有的直径都经过该边。

Input

第一行包含一个整数N，表示节点数。 
接下来N-1行，每行三个整数a, b, c ，表示点 a和点b之间有一条长度为c
的无向边。

Output

共两行。第一行一个整数，表示直径的长度。第二行一个整数，表示被所有
直径经过的边的数量。

Sample Input

6
3 1 1000
1 4 10
4 2 100
4 5 50
4 6 100

Sample Output

1110
2

【样例说明】
直径共有两条，3 到2的路径和3到6的路径。这两条直径都经过边(3, 1)和边(1, 4)。

HINT

对于100%的测试数据：2≤N≤200000，所有点的编号都在1..N的范围内，

边的权值≤10^9。

就是让你求直径的长和直径并的数量。

直径当然好求，而直径并，一定是在一条直径上。

所以我们可以先求出一条最长链。而所有直径的并一定是最长链上连续的一段。

证明很简单：如果中间有分开而最后又和在一起，显然会形成一个环。

然后我们对于最长链上的每个点，dfs出其子树中理他最远的点，若两点之间的距离等于该点到直径一个端点的距离，那么显然这个点到端点之间这一段就不能用来统计答案了。

然后我们可以从左往右做一便这个操作，反向再做一遍，中间部分即为直径的并。

代码：

#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
vector<int> a[],b[];
int last[],u,v,next[];
long long dis[],mmm[],op;
bool vv[];
void dfs1(int o,long long p,int q)
{
    if(p>op){op=p;u=o;}
    for(int i=;i<a[o].size();i++)
        if((!vv[a[o][i]])&&(a[o][i]!=q))
        {
            vv[a[o][i]]=true;
            dfs1(a[o][i],p+b[o][i],o);
        }
}
void dfs2(int o,long long p,int q)
{
    last[o]=q;
    dis[o]=p;
    if(p>op){op=p;v=o;}
    for(int i=;i<a[o].size();i++)
        if((!vv[a[o][i]])&&(a[o][i]!=q))
        {
            vv[a[o][i]]=true;
            dfs2(a[o][i],p+b[o][i],o);
        }
}
int main()
{
    int n;
    cin>>n;
    for(int i=;i<n;i++)
    {
        int x,y,z;
        cin>>x>>y>>z;
        a[x].push_back(y);
        b[x].push_back(z);
        a[y].push_back(x);
        b[y].push_back(z);
    }
    memset(vv,,sizeof(vv));op=;
    dfs1(,,);
    memset(vv,,sizeof(vv));op=;
    dfs2(u,,);
    int distance=dis[v];
    cout<<dis[v]<<endl;
    memset(vv,,sizeof(vv));
    for(int i=v;i!=;i=last[i]) vv[i]=true;
    for(int i=v;i!=;i=last[i])
    {
        op=;
        dfs1(i,,);
        mmm[i]=op;
    }
    int j=v;
    for(int i=last[v];i!=;i=last[i]) next[i]=j,j=i;
    int ans=;
    int i;
    for(i=j;i!=;i=next[i])
        if(dis[v]-dis[i]==mmm[i]) break;
    for(;i!=;i=last[i])
    {
        if(dis[i]==mmm[i]) break;
        ans++;
    }
    cout<<ans<<endl;
    return ;
}