NOIP2015 运输计划 (树上差分+二分答案)
---恢复内容开始---
题目大意:给你一颗树,你可以把其中一条边的边权改成0,使给定的一些树链的权值和的最大值最小
把lenth定义为未修改边权时的答案
考虑二分答案,如果二分的答案成立,设修改成0的边边权为k,那么所有原链长>mid的链都要被这条边影响,显然这些链存在边权为k的公共边
那么我们二分出一个答案,然后把所有原链长>mid的链在树上打差分来记录每条边被覆盖多少次
如果某条边的满足lenth-边权<=mid 且 这条边被覆盖次数==原链长>mid的链数,说明这条边作为0边可行。
会卡常,据说树剖不会被卡但我用了结果还是被卡了。懒得优化了
#include <cstdio>
#include <algorithm>
#include <cstring>
#define N 300100
using namespace std; int n,m,cte,mxl,cnt,lth,ct,idd;
int head[N],sz[N],tp[N],hvy[N],fa[N],dep[N],dif[N],sum[N],d[N],dis[N];
struct EDGE{
int to,nxt,val;
}edge[N*];
struct FLY{
int x,y,ff,len;
}fly[N];
void edge_add(int u,int v,int w)
{
cte++;
edge[cte].to = v;
edge[cte].nxt= head[u];
edge[cte].val= w;
head[u]=cte;
}
int gc()
{
int rett=,fh=;char p=getchar();
while(p<''||p>'') {p=getchar();}
while(p>=''&&p<='') {rett=rett*+p-'';p=getchar();}
return rett*fh;
}
void dfs_order1(int u,int dad)
{
int h=,id=;
for(int j=head[u];j!=-;j=edge[j].nxt)
{
int v=edge[j].to;
if(v==dad) continue;
dis[v]=dis[u]+edge[j].val;
dep[v]=dep[u]+;
dfs_order1(v,u);
d[v]=edge[j].val;
fa[v]=u;
if(sz[v]>h) h=sz[v],id=v;
sz[u]+=sz[v];
}
sz[u]++;
if(head[u]==-) return;
hvy[u]=id;
}
void dfs_order2(int u,int dad,int tpp)
{
tp[u]=tpp;
if(hvy[u]) tp[hvy[u]]=tp[u],dfs_order2(hvy[u],u,tpp);
for(int j=head[u];j!=-;j=edge[j].nxt)
{
int v=edge[j].to;
if(v==dad||v==hvy[u]) continue;
dfs_order2(v,u,v);
}
}
void dfs_maxlenth(int u,int dad)
{
sum[u]+=dif[u];
for(int j=head[u];j!=-;j=edge[j].nxt)
{
int v=edge[j].to;
if(v==dad) continue;
dfs_maxlenth(v,u);
sum[u]+=sum[v];
}
if(sum[u]>=ct&&d[u]>mxl) idd=u,mxl=d[u];
}
int LCA(int x,int y)
{
while(tp[x]!=tp[y])
{
if(dep[tp[x]]<dep[tp[y]]) swap(x,y);
x=fa[tp[x]];
}
return dep[x]<=dep[y]?x:y;
}
bool check(int mid)
{
memset(sum,,sizeof(sum));
mxl=,ct=;;
for(int i=;i<=m;i++)
if(fly[i].len>mid)
{
sum[fly[i].x]++,sum[fly[i].y]++;
sum[fly[i].ff]-=;
ct++;
}
dfs_maxlenth(,-);
if(lth-mxl<=mid) return ;
else return ;
} int main()
{
//freopen("testdata.in","r",stdin);
scanf("%d%d",&n,&m);
memset(head,-,sizeof(head));
int x,y,z;
for(int i=;i<n;i++)
{
x=gc(),y=gc(),z=gc();
edge_add(x,y,z);
edge_add(y,x,z);
}
dfs_order1(,-);
tp[]=;
dfs_order2(,-,);
lth=;
for(int i=;i<=m;i++)
{
fly[i].x=gc(),fly[i].y=gc();
fly[i].ff=LCA(fly[i].x,fly[i].y);
fly[i].len=dis[fly[i].x]+dis[fly[i].y]-*dis[fly[i].ff];
lth=max(lth,fly[i].len);
}
int l=,r=;
r=lth;
int ret=;
while(l<r)
{
int mid=(l+r)>>;
if(check(mid)) r=mid,ret=mid;
else l=mid+;
}
printf("%d",ret);
return ;
}
---恢复内容结束---
NOIP2015 运输计划 (树上差分+二分答案)的更多相关文章
- [NOIP 2015]运输计划-[树上差分+二分答案]-解题报告
[NOIP 2015]运输计划 题面: A[NOIP2015 Day2]运输计划 时间限制 : 20000 MS 空间限制 : 262144 KB 问题描述 公元 2044 年,人类进入了宇宙纪元. ...
- 洛谷P2680 运输计划(倍增LCA + 树上差分 + 二分答案)
[题目链接] [思路]: 根据题意可以明显看出,当所有任务都完成时的时间是最终的结果,也就是说本题要求,求出最小的最大值. 那这样的话就暗示了将答案二分,进行check. [check方法]: 如果说 ...
- 洛谷P2680 运输计划——树上差分
题目:https://www.luogu.org/problemnew/show/P2680 久违地1A了好高兴啊! 首先,要最大值最小,很容易想到二分: 判断当前的 mid 是否可行,需要看看有没有 ...
- [NOIP2015]运输计划 线段树or差分二分
目录 [NOIP2015]运输计划 链接 思路1 暴力数据结构 思路2 二分树上差分 总的 代码1 代码2 [NOIP2015]运输计划 链接 luogu 好久没写博客了,水一篇波. 思路1 暴力数据 ...
- BZOJ 4326 NOIP2015 运输计划(树上差分+LCA+二分答案)
4326: NOIP2015 运输计划 Time Limit: 30 Sec Memory Limit: 128 MB Submit: 1388 Solved: 860 [Submit][Stat ...
- BZOJ 4326 NOIP2015 运输计划 (二分+树上差分)
4326: NOIP2015 运输计划 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 1930 Solved: 1231[Submit][Statu ...
- [NOIP2015]运输计划 D2 T3 LCA+二分答案+差分数组
[NOIP2015]运输计划 D2 T3 Description 公元2044年,人类进入了宇宙纪元. L国有n个星球,还有n-1条双向航道,每条航道建立在两个星球之间,这n-1条航道连通了L国的所有 ...
- LOJ2425 NOIP2015 运输计划 【二分+LCA+树上差分】*
LOJ2425 NOIP2015 运输计划 LINK 题意:给你一颗树,可以将任意一条边的权值变成0,然后求m条路径的长度的最小值 思路: 先二分最后的距离ans,然后我们把路程大于ans的所有路径拿 ...
- 【BZOJ-4326】运输计划 树链剖分 + 树上差分 + 二分
4326: NOIP2015 运输计划 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 703 Solved: 461[Submit][Status] ...
随机推荐
- IOS开发:使用lipo合并armv7,i386,armv7s库文件
假设多个版本的lib分别是 libxxx.armv7.a , libxxx.armv7s.a, libxxx.i386.a我们的目标是 把他们合并成超级通用版的libxxx.a 打开命令行 Term ...
- node中exports和module.exports的关系及使用
在node中,需要记住,在使用exports和module.exports的时候,实际输出的是module.exports. exports指向module.exports,是module.expor ...
- 常用HTML小部件
- Java多线程中Sleep与Wait的区别
Java中的多线程是一种抢占式的机制 而不是分时机制.抢占式机制指的是有多个线程处于可运行状态,但是只有一个线程在运行. 共同点: 1. 他们都是在多线程的环境下,都可以在程序的调用处阻塞指定的毫秒数 ...
- nyoj256-C小加之级数求和
C小加 之 级数求和 时间限制:3000 ms | 内存限制:65535 KB 难度:2 描述 最近,C小加 又遇到难题了,正寻求你的帮助. 已知:Sn= 1+1/2+1/3+-+1/n. 显然对 ...
- Linux 字符设备驱动简单总结(转)
http://my.oschina.net/u/1169027/blog/191538
- 洛谷—— P1126 机器人搬重物
https://www.luogu.org/problem/show?pid=1126 题目描述 机器人移动学会(RMI)现在正尝试用机器人搬运物品.机器人的形状是一个直径1.6米的球.在试验阶段,机 ...
- [React Router] Create a ProtectedRoute Component in React Router (setState callback to force update)
In this lesson we'll create a protected route just for logged in users. We'll combine a Route with a ...
- SVN打tag
SVN打tag是一个非经常常使用的功能,要谈打tag,还得从SVN官方推荐的文件夹结构说起.SVN官方推荐在一个版本号库的根文件夹下先建立trunk.branches.tags这三个文件夹 ...
- 四、基于HTTPS协议的12306抢票软件设计与实现--水平DNS并发查询分享
一.基于HTTPS协议的12306抢票软件设计与实现--实现效果 二.基于HTTPS协议的12306抢票软件设计与实现--相关接口以及数据格式 三.基于HTTPS协议的12306抢票软件设计与实现-- ...