http://www.elijahqi.win/archives/1371
题目背景
公元 2044 年,人类进入了宇宙纪元。
题目描述
L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球。
小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如:有一艘物
流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道 是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之 间不会产生任何干扰。
为了鼓励科技创新,L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小 P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。
在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后, 这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的 物流公司的阶段性工作就完成了。
如果小 P 可以自由选择将哪一条航道改造成虫洞,试求出小 P 的物流公司完成阶段 性工作所需要的最短时间是多少?
输入输出格式
输入格式:
输入文件名为 transport.in。
第一行包括两个正整数 n、m,表示 L 国中星球的数量及小 P 公司预接的运输计划的数量,星球从 1 到 n 编号。
接下来 n-1 行描述航道的建设情况,其中第 i 行包含三个整数 ai, bi 和 ti,表示第
i 条双向航道修建在 ai 与 bi 两个星球之间,任意飞船驶过它所花费的时间为 ti。
接下来 m 行描述运输计划的情况,其中第 j 行包含两个正整数 uj 和 vj,表示第 j个 运输计划是从 uj 号星球飞往 vj 号星球。
输出格式:
输出 共1行,包含1个整数,表示小P的物流公司完成阶段性工作所需要的最短时间。
输入输出样例
输入样例#1: 复制
6 3
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5
输出样例#1: 复制
11
说明
所有测试数据的范围和特点如下表所示
请注意常数因子带来的程序效率上的影响。
真tmd一个卡常题 卡了半天没卡过去
二分答案 然后每次对所有的线路扫描一遍 然后 求出我交集 然后用线路中最大值–交集中的边的最大值 看能否小于我二分的答案
#include<cstdio>
#define N 300010
#include<algorithm>
#include<cstring>
using namespace std;
inline char gc(){
static char now[1<<16],*T,*S;
if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0;char ch=getchar();
while (ch<'0'||ch>'9') ch=getchar();
while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
return x;
}
struct node{
int y,next,z;
}data[N<<1];
struct node1{
int st,ed,value,lca;
}edge[N];
int num,n,m,h[N],Log[N],fa[N][20],sum[N][20],dep[N],flag[N];
void dfs(int x){
for (int i=h[x];i;i=data[i].next){
int y=data[i].y,z=data[i].z;
if (fa[x][0]==y) continue;fa[y][0]=x;dep[y]=dep[x]+1;sum[y][0]=z;
for (int j=1;j<=Log[dep[y]];++j){
fa[y][j]=fa[fa[y][j-1]][j-1];sum[y][j]=sum[y][j-1]+sum[fa[y][j-1]][j-1];
}dfs(y);
}
}
inline int calc(int x,int y,int id){
int sum1=0;
if (dep[x]<dep[y]) swap(x,y);
int dis=dep[x]-dep[y];
for (int i=0;i<=Log[dis];++i){
if ((1<<i)&dis) sum1+=sum[x][i],x=fa[x][i];
}edge[id].lca=x;
if (x==y)return sum1;
for (int i=Log[dep[y]];i>=0;--i){
if (fa[x][i]!=fa[y][i]){
sum1+=sum[x][i]+sum[y][i];x=fa[x][i];y=fa[y][i];
}
}edge[id].lca=fa[x][0];
sum1+=sum[x][0]+sum[y][0];return sum1;
}
void dfs1(int x){
for (int i=h[x];i;i=data[i].next){
int y=data[i].y;if (fa[x][0]==y) continue;
dfs1(y);flag[x]+=flag[y];
}
}
inline bool cmp(node1 a,node1 b){
return a.value>b.value;
}
inline bool check(int x){
memset(flag,0,sizeof(flag));
int max_dis=0,k=0;
for (int i=1;i<=m;++i){
if (edge[i].value>x){
flag[edge[i].st]++;flag[edge[i].ed]++;flag[edge[i].lca]-=2;
max_dis=max(max_dis,edge[i].value);++k;
}else break;
}dfs1(1);int maxans=0;
for (int i=2;i<=n;++i){
if (flag[i]==k){
maxans=max(maxans,sum[i][0]);
}
}return max_dis-maxans<=x;
}
int main(){
//freopen("2680.in","r",stdin);
n=read();m=read();
for (int i=1;i<n;++i){
int x=read(),y=read(),z=read();
data[++num].y=y;data[num].next=h[x];data[num].z=z;h[x]=num;
data[++num].y=x;data[num].next=h[y];data[num].z=z;h[y]=num;
}Log[0]=-1;num=0;
for (int i=1;i<=n;++i) Log[i]=Log[i>>1]+1;
dfs(1);
for (int i=1;i<=m;++i){
int x=read(),y=read();
edge[i].st=x;edge[i].ed=y;edge[i].value=calc(x,y,i);
}sort(edge+1,edge+m+1,cmp);
int l=0,r=edge[1].value;
while (l<=r){
int mid=l+r>>1;
if (check(mid)) r=mid-1;else l=mid+1;
}printf("%d",l);
return 0;
}