从下午坑到网上。。noip的数据太弱,若干的地方写挂结果还随便过= =

最坑的就是网上有些题解没考虑周全。。。

第一步是找直径,用两次bfs(或者dfs,Linux下系统栈挺大的。。)解决。找出其中一条直径就可以了,虽然蒟蒻不会证明但是看起来似乎挺有道理的

要看证明的话可以看这个题解:http://trinklee.blog.163.com/blog/static/238158060201411175015709/

直径上的路径本来有n^2条,但是我们发现,首先对于同一起点/终点的路径,路径长度越长(合法的前提下)偏心距肯定越小,但是路径又有长度的约束,所以可以用单调队列来优化。

假设直径两端点为x,y,对于一条路径(a,b),离路径最远的点要么是x和y中的一个,要么是 a到b的路径上的节点中 不走直径到达的最远的节点。

假设x离a比较近,b离y比较近,如果有不在直径上的节点k使得dist(k,a)>dist(x,a)的话,那么(x,y)就不是直径了((k,a)(a,y)可以拼成更长的一条路径);对于b同理

但是需要注意的是,以上的k节点必须满足 k到路径最近的点是路径的端点。。。

不然的话假设k到路径最近的点为c,那么dis(c,k)也有可能是偏心距。。。最直观的例子是给定的s比直径长度还大,那样偏心距就全得靠k了。。

当然了noip数据足够水noip数据足够水noip数据足够水所以。。。

单调队列的话蒟蒻弄了两个单调队列,一个是递增地储存x到路径上的节点的距离,另一个递减地储存路径上的点到 其他节点(具体看上面)最远的距离。。。

当然了因为第一个单调队列每次删掉非法的队头后只要直接在队尾入队,所以实际上可以不用开队列,只要记录队头队尾的位置就好了= =

 1 #include<cstdio>
 2 #include<math.h>
 3 #include<algorithm>
 4 #include<iostream>
 5 #include<cstring>
 6 using namespace std;
 7 const int maxn=500233;
 8 struct edge{
 9     int too,pre,dis;
10 }e[maxn<<1];
11 struct zs{
12     int id,num;
13 }dll[maxn],q[maxn];
14 int last[maxn],num[maxn],dl[maxn],f[maxn],p[maxn];
15 int dist[4][maxn];
16 int i,j,k,a,b;
17 int n,m,l,r,now,tot,posa,posb,s,c,ans,l1,r1;
18 bool u[maxn],arr[maxn],cant[maxn];
19 
20 
21 inline void insert(int a,int b,int c){e[++tot].too=b;e[tot].dis=c;e[tot].pre=last[a];last[a]=tot;
22 e[++tot].too=a;e[tot].dis=c;e[tot].pre=last[b];last[b]=tot;
23 }
24 void bfs(int st,int ed,int id){
25     int l=0,r=1,now,i,j;dl[1]=st;u[st]=1;dist[id][st]=0;
26     while(l<r){
27         now=dl[++l];
28         for(i=last[now];i;i=e[i].pre)if(!u[e[i].too]){
29             dl[++r]=e[i].too,dist[id][e[i].too]=dist[id][now]+e[i].dis,u[e[i].too]=1;
30             if(id==3)if(dl[r]==ed){l=r;break;
31             }
32         }
33     }
34     memset(u,0,sizeof(u));
35     if(id==3){
36         arr[ed]=1;num[1]=ed;num[0]=1;
37         for(i=r;i;i--)if(dist[2][dl[i]]+dist[3][dl[i]]==dist[2][st])num[++num[0]]=dl[i],arr[dl[i]]=1;
38         
39     }
40 }
41 void bii(int x){
42     int i;u[x]=1;
43     for(i=last[x];i;i=e[i].pre)if(!arr[e[i].too]&&!u[e[i].too]){
44         bii(e[i].too);
45         f[x]=max(f[x],f[e[i].too]+e[i].dis);
46     }
47 }
48 int main(){
49     scanf("%d%d",&n,&s);
50     for(i=1;i<n;i++)scanf("%d%d%d",&a,&b,&c),insert(a,b,c);
51     bfs(1,0,1);posa=0;dist[1][0]=-1;for(i=1;i<=n;i++)if(dist[1][i]>dist[1][posa])posa=i;
52     bfs(posa,0,2);posb=0;dist[2][0]=-1;for(i=1;i<=n;i++)if(dist[2][i]>dist[2][posb])posb=i;
53     bfs(posb,posa,3);
54     for(i=1;i<=num[0];i++)bii(num[i]);
55     ans=1233333333;
56     for(i=1;i<=num[0];i++)p[i]=dist[2][num[i]];
57     l=1;r=0;l1=1;r1=0;
58     for(i=1;i<=num[0];i++){
59         while(l<=r&&dist[2][num[i]]-dist[2][num[dll[l].id]]>s)cant[dll[l].id]=1,l++;
60         while(l1<=r1&&cant[q[l1].id])l1++;
61         while(l1<=r1&&q[r1].num<=f[num[i]])r1--;
62         q[++r1].num=f[num[i]];q[r1].id=i;
63         dll[++r].num=p[i];dll[r].id=i;
64         ans=min(ans,max(max(dll[l].num,dist[3][num[i]]),q[l1].num));
65     }
66     printf("%d\n",ans);
67     return 0;
68 }
View Code

相关文章: