4326: NOIP2015 运输计划
Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 1841 Solved: 1173
[ Submit][ Status][ Discuss]
Description
公元 2044 年,人类进入了宇宙纪元。L 国有 n 个星球,还有 n−1 条双向航道,每条航道建立在两个星球之间,这 n−1 条航道连通了 L 国的所有星球。小 P 掌管一家物流公司, 该公司有很多个运输计划,每个运输计划形如:有一艘物流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之间不会产生任何干扰。为了鼓励科技创新, L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后,这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的物流公司的阶段性工作就完成了。如果小 P 可以自由选择将哪一条航道改造成虫洞, 试求出小 P 的物流公司完成阶段性工作所需要的最短时间是多少?
Input
第一行包括两个正整数 n,m,表示 L 国中星球的数量及小 P 公司预接的运输计划的数量,星球从 1 到 n 编号。接下来 n−1 行描述航道的建设情况,其中第 i 行包含三个整数 ai,bi 和 ti,表示第 i 条双向航道修建在 ai 与 bi 两个星球之间,任意飞船驶过它所花费的时间为 ti。数据保证 1≤ai,bi≤n 且 0≤ti≤1000。接下来 m 行描述运输计划的情况,其中第 j 行包含两个正整数 uj 和 vj,表示第 j 个运输计划是从 uj 号星球飞往 vj号星球。数据保证 1≤ui,vi≤n
Output
输出文件只包含一个整数,表示小 P 的物流公司完成阶段性工作所需要的最短时间。
Sample Input
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5
Sample Output
HINT
Source
刚开始看到题,最大值最小,树上路径问题又不带修改,直接树上差分+二分。然后还要在树上做下路径以及lca,尽管知识点很多,很杂,但是不难想,算是noip中比较宽容的一个d2t3吧。。。
先二分一个答案mid,然后把路径长度大于mid的运输计划取出来做树上差分,判断所有路径都有交集的点(因为只有这样才能把所有路径长度减到mid以下),if(f[j]==i&&mid>=p[1].dis-dist[j][0])如果最长路径可行说明下面都可行(事先按照降序排好)(dist[i][j]表示i起点向上2^1长度的路径总长)。如果这个答案可行,更新ans,然后r=mid,否则l=mid,一直更新答案即可。
具体细节见代码
#include <cstdio>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
vector<pair<int,int> >e[300005];
int n,m,ai,bi,ti,l=0,r,f[300005],anc[524289][20],dist[524289][20],deep[300005]={0,1},cnt,ans=2147483647;
struct trans
{
int u,v,dis,Lca;
}p[300005];
int read()
{
int ans=0,f=1;
char c=getchar();
while(c<'0'||c>'9')
{
if(c=='-') f=-1;
c=getchar();
}
while(c>='0'&&c<='9')
{
ans=ans*10+c-'0';
c=getchar();
}
return ans*f;
}
bool cmp(trans a,trans b)
{
return a.dis>b.dis;
}
void dfs1(int now)
{
for(int i=1;i<=19;i++)
{
anc[now][i]=anc[anc[now][i-1]][i-1];
dist[now][i]=dist[anc[now][i-1]][i-1]+dist[now][i-1];
}
for(int i=0;i<e[now].size();i++)
if(e[now][i].first!=anc[now][0])
{
deep[e[now][i].first]=deep[now]+1;
dist[e[now][i].first][0]=e[now][i].second;
anc[e[now][i].first][0]=now;
dfs1(e[now][i].first);
}
}
void lca(int cur,int u,int v)
{
int ret=0;
if(deep[u]<deep[v]) swap(u,v);
for(int i=19;i>=0;i--)
if(deep[v]<=deep[anc[u][i]])
{
ret+=dist[u][i];
u=anc[u][i];
}
if(u==v)
{
p[cur].Lca=u;
p[cur].dis=ret;
return;
}
for(int i=19;i>=0;i--)
if(anc[u][i]!=anc[v][i])
{
ret+=dist[u][i]+dist[v][i];
u=anc[u][i];
v=anc[v][i];
}
ret+=dist[u][0]+dist[v][0];
p[cur].Lca=anc[u][0];
p[cur].dis=ret;
}
void dfs2(int now)
{
for(int i=0;i<e[now].size();i++)
if(e[now][i].first!=anc[now][0])
{
dfs2(e[now][i].first);
f[now]+=f[e[now][i].first];
}
}
int main()
{
n=read();m=read();
for(int i=1;i<n;i++)
{
ai=read();bi=read();ti=read();
e[ai].push_back(make_pair(bi,ti));
e[bi].push_back(make_pair(ai,ti));
}
dfs1(1);
for(int i=1;i<=m;i++)
{
p[i].u=read();p[i].v=read();
lca(i,p[i].u,p[i].v);//把lca和dis预处理好,避免在线计算
r=max(r,p[i].dis);
}
sort(p+1,p+m+1,cmp);
while(l<r)
{
int mid=cnt?r:(l+r)>>1,i=1,ok=0;//防止死循环以及漏掉答案
memset(f,0,sizeof(f));
for(;i<=m;i++)
{
if(p[i].dis<=mid) break;
f[p[i].Lca]-=2;
f[p[i].u]++;
f[p[i].v]++;//树上差分
}
i--;
if(!i)
{
ans=min(ans,r=mid);
continue;
}
dfs2(1);
for(int j=2;j<=n;j++)
if(f[j]==i&&mid>=p[1].dis-dist[j][0]) ok=1;
if(ok&&!cnt) ans=min(ans,r=mid);
else if(!ok&&!cnt) l=mid;
if(r-l==1) cnt++;//防止死循环以及漏掉答案
if(cnt==2) break;//防止死循环以及漏掉答案
}
printf("%d\n",ans==2147483647?0:ans);//细节要注意
}