nkoj 1313（noip2007）树网的核 与 nkoj 2650（SDOI2011）消防

1.树网的核

问题描述:

设T=(V, E, W) 是一个无圈且连通的无向图（也称为无根树），每条边到有正整数的权，我们称T为树网（treebetwork），其中V，E分别表示结点与边的集合，W表示各边长度的集合，并设T有n个结点。
路径：树网中任何两结点a，b都存在唯一的一条简单路径，用d(a, b)表示以a, b为端点的路径的长度，它是该路径上各边长度之和。我们称d(a, b)为a, b两结点间的距离。
D(v, P)=min{d(v, u), u为路径P上的结点}。
树网的直径：树网中最长的路径成为树网的直径。对于给定的树网T，直径不一定是唯一的，但可以证明：各直径的中点（不一定恰好是某个结点，可能在某条边的内部）是唯一的，我们称该点为树网的中心。
偏心距ECC(F)：树网T中距路径F最远的结点到路径F的距离，即 ECC(F)=max{d(v, F)，v∈V}
任务：对于给定的树网T=(V, E, W)和非负整数s，求一个路径F，他是某直径上的一段路径（该路径两端均为树网中的结点），其长度不超过s（可以等于s），使偏心距ECC(F)最小。我们称这个路径为树网T=(V, E, W)的核（Core）。必要时，F可以退化为某个结点。一般来说，在上述定义下，核不一定只有一个，但最小偏心距是唯一的。
下面的图给出了树网的一个实例。图中，A-B与A-C是两条直径，长度均为20。点W是树网的中心，EF边的长度为5。如果指定s=11，则树网的核为路径DEFG（也可以取为路径DEF），偏心距为8。如果指定s=0（或s=1、s=2），则树网的核为结点F，偏心距为12。

输入

第1行: 2个正整数n和s，中间用一个空格隔开。其中n为树网结点的个数，s为树网的核的长度的上界。设结点编号以此为1，2，……，n。 第2行..n行: 每行给出3个用空格隔开的正整数，依次表示每一条边的两个端点编号和长度。例如，“2 4 7”表示连接结点2与4的边的长度为7。 所给的数据都是正确的，不必检验。

输出

第1行：只有一个非负整数，为指定意义下的最小偏心距。

样例输入输出

输入样例1
5 2
1 2 5
2 3 2
2 4 4
2 5 3

输入样例2
8 6
1 3 2
2 3 2
3 4 6
4 5 3
4 6 4
4 7 2
7 8 3

输出样例1
5
输出样例2
5

水解：此题，数据贼小，所以可用弗洛伊德求各点最短路进行预处理，同时记录dis最大的两点，这两点就构成了树的直径的两端，这样就可以找出树的直径。
然后，枚举直径上的线段，通过判断（dis[i][l]+dis[i][r]==直径）来确定是否在直径上。
代码如下：

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
const int maxn=305,inf=1e9;
int dis[maxn][maxn];
int i,j,k,n,s,x,y,z,l,r,maxx=0;
inline void re(int &d)
{
    bool f=false;char t=getchar();
    while(t<'0'||t>'9'){if(t=='-')f=true;t=getchar();}
    for(d=0;t>='0'&&t<='9';t=getchar())d=(d<<1)+(d<<3)+t-'0';
    d=(f?(-d):(d));
}
int main()
{
    int tmp=0,ans=inf;
    scanf("%d%d",&n,&s); 
    for(i=1;i<=n;i++)
    for(j=1;j<=n;j++)
    {
        dis[i][j]=dis[j][i]=inf;
        dis[i][i]=0;
    }
    for(i=1;i<n;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        dis[x][y]=dis[y][x]=z;
    }
    for(k=1;k<=n;k++)
    for(i=1;i<=n;i++)
    {
        if(i==k)continue;
        for(j=1;j<=n;j++)
        {
            if(j==k||j==i)continue;
            if(dis[i][k]!=inf&&dis[j][k]!=inf)
            {
                dis[i][j]=min(dis[i][j],dis[i][k]+dis[j][k]);
                dis[j][i]=dis[i][j];
                if(dis[i][j]>maxx&&dis[i][j]!=inf)
                {maxx=dis[i][j];l=i;r=j;}
            }
        }
    }
    for(i=1;i<=n;i++)
    if(dis[i][l]+dis[i][r]==maxx)
    {
        for(j=1;j<=n;j++)
        if(dis[i][j]<=s&&dis[j][l]+dis[j][r]==maxx)
        {
            tmp=max(min(dis[i][l],dis[j][l]),min(dis[j][r],dis[i][r]));
            ans=min(tmp,ans);
        }
    }
    printf("%d",ans);
}

2.消防

时限2s
　　某个国家有n个城市，这n个城市中任意两个都连通且有唯一一条路径，每条连通两个城市的道路的长度为zi(zi<=1000)。
　　这个国家的人对火焰有超越宇宙的热情，所以这个国家最兴旺的行业是消防业。由于政府对国民的热情忍无可忍（大量的消防经费开销）可是却又无可奈何（总统竞选的国民支持率），所以只能想尽方法提高消防能力。
　　现在这个国家的经费足以在一条边长度和不超过s的路径（两端都是城市）上建立消防枢纽，为了尽量提高枢纽的利用率，要求其他所有城市到这条路径的距离的最大值最小。
　　你受命监管这个项目，你当然需要知道应该把枢纽建立在什么位置上。

输入输出样例都跟上题一样。
不过

【数据规模和约定】
对于20%的数据，n<=300。
对于50%的数据，n<=3000。
对于100%的数据，n<=300000，边长小等于1000。

倘若先做2题，就可轻松A1题，2题的规模达到了300000，我们需要更优秀的方法。
可以采用dfs+队列，先随便找个点开始搜，找到最深的点，再从这个点反过来搜，再找最深的点，这样就求出了树的直径。
我们要找的区间在直径上，对于我们要找的区间，枚举左端点，然后用队列维护一下最值即可。
代码如下：

#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
const int maxn=300050;
int n,s,x,y,z,mark=0,st=0,t=0,ans=1e9;
int f[maxn],road[maxn],l[maxn],r[maxn],vis[maxn],dis[maxn],p[maxn],father[maxn],fw[maxn];
struct node{int to,length;};
vector<node>v[maxn];queue<int>q;
inline void re(int &d)
{
    bool f=false;char t=getchar();
    while(t<'0'||t>'9'){if(t=='-')f=true;t=getchar();}
    for(d=0;t>='0'&&t<='9';t=getchar())d=(d<<1)+(d<<3)+t-'0';
    d=(f?(-d):(d));
}
int find(int u,int mark)  
{  
    dis[u]=0;vis[u]=mark;  
    int len=0;q.push(u);  
    while(!q.empty()){  
        int x=q.front();q.pop();  
        for(int i=0;i<v[x].size();i++) 
        {  
            int y=v[x][i].to;  
            if (vis[y]==mark)continue;
            father[y]=x;fw[y]=v[x][i].length;
            dis[y]=dis[x]+fw[y];len=max(len,dis[y]);
            vis[y]=mark;q.push(y);
        }  
    }  
    return len;  
} 
int main()
{
    int hea=1,tai=0,L,R=0;
    scanf("%d%d",&n,&s);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d%d",&x,&y,&z);
        v[x].push_back((node){y,z});
        v[y].push_back((node){x,z});
    }
    find(1,++mark);
    for(int i=1;i<=n;i++)if(dis[i]>dis[st])st=i;
    find(st,++mark);
    for(int i=1;i<=n;i++)if(dis[i]>dis[t])t=i;
    n=0;mark++;
    while(1)
    {road[++n]=t;vis[t]=mark;
    t=father[t];if(road[n]==st)break;}
    for(int i=n;i>=0;i--)r[i]=dis[road[i]]; 
    for(int i=1;i<=n;i++)l[i]=dis[road[1]]-dis[road[i]];
    for(int i=1;i<=n;i++)f[i]=find(road[i],mark);
    for(L=1;L<=n;L++)
    {
        while(p[hea]<L&&hea<=tai)hea++;
        while(R<n&&l[R+1]-l[L]<=s)
        {
            while(hea<=tai&&f[p[tai]]<=f[R+1])tai--;
            p[++tai]=++R;
        }
        int tmp=max(f[p[hea]],max(l[L],r[R]));
        ans=min(ans,tmp);
    }
    printf("%d",ans);
}