BZOJ3641 : 货车运输
若一条边的v小于等于u,则贡献为l*w/v,否则贡献为l*w/u
将边按v从小到大排序,将询问按u从小到大排序
用树链剖分维护链上和,val[0]表示第一种情况下的贡献,val[1]表示第二种情况下的贡献
一开始val[0]都是0,val[1]=l*w,
然后每到一个询问(s,t,u),就把所有v小于等于u的边修改掉,val[0]改为l*w/v,val[1]改为0
因为是环套外向树,所以把额外的边(exa,exb,exl,exv)拿走后变成了一棵树。
查询时在下面三种情况中取min:
1.(s,t)链上val[0]的和+(s,t)链上val[1]的和/u
2.(s,exa)链上val[0]的和+(t,exb)链上val[0]的和+((s,exa)链上val[1]的和+(t,exb)链上val[1]的和)/u+exl*exw/min(u,exv)
3.(s,exb)链上val[0]的和+(t,exa)链上val[0]的和+((s,exb)链上val[1]的和+(t,exa)链上val[1]的和)/u+exl*exw/min(u,exv)
时间复杂度$O(n\log n)$
#include<cstdio>
#include<algorithm>
#define N 100010
using namespace std;
inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';}
struct edge{int l,x,v,w,id;}E[N];
struct que{int s,t,u,id;}Q[N];
int n,m,q,i,x,y,z,t,exa,exb,exl,exx,exv,exw,cnt,now,V[N],W[N],g[N<<1],nxt[N<<2],v[N<<2],ed,size[N<<1],d[N<<1],f[N<<1],son[N<<1],loc[N<<1],top[N<<1],dfn;
double ans[N],val[2][N<<3];
inline bool cmpe(edge a,edge b){return a.v<b.v;}
inline bool cmpq(que a,que b){return a.u<b.u;}
int F(int x){return f[x]==x?x:F(f[x]);}
inline void add(int x,int y){v[++ed]=y;nxt[ed]=g[x];g[x]=ed;}
void dfs1(int x,int pre){
size[x]=1;d[x]=d[f[x]=pre]+1;
int heavy=0,sizemax=0,i;
for(i=g[x];i;i=nxt[i])if(v[i]!=pre){
dfs1(v[i],x),size[x]+=size[v[i]];
if(size[v[i]]>sizemax)sizemax=size[v[i]],heavy=v[i];
}
if(heavy)son[x]=heavy;
}
void dfs2(int x,int pre,int t){
loc[x]=++dfn;top[x]=t;
if(son[x])dfs2(son[x],x,t);
for(int i=g[x];i;i=nxt[i])if(v[i]!=pre&&v[i]!=son[x])dfs2(v[i],x,v[i]);
}
void change(int p,int x,int a,int b,int c,double d){
if(a==b){val[p][x]=d;return;}
int mid=(a+b)>>1;
if(c<=mid)change(p,x<<1,a,mid,c,d);else change(p,x<<1|1,mid+1,b,c,d);
val[p][x]=val[p][x<<1]+val[p][x<<1|1];
}
double ask(int p,int x,int a,int b,int c,int d){
if(c<=a&&b<=d)return val[p][x];
int mid=(a+b)>>1;double t=0;
if(c<=mid)t=ask(p,x<<1,a,mid,c,d);
if(d>mid)t+=ask(p,x<<1|1,mid+1,b,c,d);
return t;
}
inline double query(int p,int x,int y){
double t=0;
while(top[x]!=top[y]){
if(d[top[x]]<d[top[y]])swap(x,y);
t+=ask(p,1,1,dfn,loc[top[x]],loc[x]);
x=f[top[x]];
}
if(d[x]>d[y])swap(x,y);
return t+ask(p,1,1,dfn,loc[x],loc[y]);
}
inline double getans(int x,int y,int u){
double t=query(0,x,y)+query(1,x,y)/(double)u,tmp;
tmp=query(0,x,exa)+query(0,y,exb)+(query(1,x,exa)+query(1,y,exb))/(double)u+(double)exl*(double)exw/(double)(exv<u?exv:u);
if(tmp<t)t=tmp;
tmp=query(0,x,exb)+query(0,y,exa)+(query(1,x,exb)+query(1,y,exa))/(double)u+(double)exl*(double)exw/(double)(exv<u?exv:u);
return tmp<t?tmp:t;
}
int main(){
read(n),read(m),read(q);
for(i=1;i<=n;i++)f[i]=i;
for(i=1;i<=n;i++){
read(x),read(y),read(z),read(t);
if(F(x)==F(y))exa=x,exb=y,exl=z,exx=t;
else{
f[f[x]]=f[y];
E[++cnt].l=z,E[cnt].x=t,E[cnt].id=n+cnt;
add(x,n+cnt),add(n+cnt,x),add(y,n+cnt),add(n+cnt,y);
}
}
for(i=1;i<=m;i++)read(V[i]),read(W[i]);
for(exv=V[exx],exw=W[exx],i=1;i<n;i++)E[i].v=V[E[i].x],E[i].w=W[E[i].x];
for(i=1;i<=q;i++)read(Q[i].s),read(Q[i].t),read(Q[i].u),Q[i].id=i;
sort(E+1,E+n,cmpe),sort(Q+1,Q+q+1,cmpq);
dfs1(1,0),dfs2(1,0,1);
for(i=1;i<n;i++)change(1,1,1,dfn,loc[E[i].id],(double)E[i].l*(double)E[i].w);
for(now=i=1;i<=q;i++){
while(now<n&&E[now].v<=Q[i].u)change(0,1,1,dfn,loc[E[now].id],(double)E[now].l*(double)E[now].w/(double)E[now].v),change(1,1,1,dfn,loc[E[now].id],0),now++;
ans[Q[i].id]=getans(Q[i].s,Q[i].t,Q[i].u);
}
for(i=1;i<=q;i++)printf("%.6f\n",ans[i]);
return 0;
}
BZOJ3641 : 货车运输的更多相关文章
- NOIP 2013 货车运输【Kruskal + 树链剖分 + 线段树 】【倍增】
NOIP 2013 货车运输[树链剖分] 树链剖分 题目描述 Description A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在 ...
- NOIP2013 货车运输 (最大生成树+树上倍增LCA)
死磕一道题,中间发现倍增还是掌握的不熟 ,而且深刻理解:SB错误毁一生,憋了近2个小时才调对,不过还好一遍AC省了更多的事,不然我一定会疯掉的... 3287 货车运输 2013年NOIP全国联赛提高 ...
- C++之路进阶——codevs3287(货车运输)
3287 货车运输 2013年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description A 国有 n ...
- NOIP2013 货车运输
3.货车运输 (truck.cpp/c/pas) [问题描述] A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有 q 辆货车在运输货 ...
- Codevs 3287 货车运输 2013年NOIP全国联赛提高组(带权LCA+并查集+最大生成树)
3287 货车运输 2013年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 传送门 题目描述 Description A 国有 n 座 ...
- 【NOIP2013提高组】货车运输
货车运输 (truck.cpp/c/pas) [问题描述] A国有n座城市,编号从1到n,城市之间有m条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有q辆货车在运输货物,司机们想知道每辆 ...
- Codevs3278[NOIP2013]货车运输
3287 货车运输 2013年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description A 国有 ...
- 洛谷 P1967 货车运输
洛谷 P1967 货车运输 题目描述 A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有 q 辆货车在运输货物, 司机们想知道每辆车在 ...
- P1967 货车运输
P1967 货车运输最大生成树+lca+并查集 #include<iostream> #include<cstdio> #include<queue> #inclu ...
随机推荐
- 安装 openSUSE Leap 42.1 之后要做的 8 件事
导读 openSUSE Leap 确实是个巨大的飞跃,它允许用户运行一个和 SUSE Linux 企业版拥有同样基因的发行版.和其它系统一样,为了实现最佳的使用效果,在使用它之前需要做些优化设置. 下 ...
- sql注入之你问我答小知识
/*每日更新,珍惜少年时博客*/ 1.问:为啥order by 是排序.而在注入当中后面加的却不是字段而是数字捏? 答:第N个字段嘛.order by 5就是第五个字段,如果5存在,6不存在.就说明只 ...
- 跟着百度学PHP[1]-if条件嵌套
权当自己的学习笔记.望大牛们切勿参考.如若发现错误,万望指出! 慕课任务 假设在发工资的时候,不仅判定性别,还要判定男性是否有房,没有房,可以发放住房补贴,对于女性,判定是否怀孕,怀孕还有怀孕补贴. ...
- 影像工作站的数据库安装错误之Win7系统下pg服务无法启动
1.关闭批处理 2.修改 PG安装路径下的Data文件下的pg_hba.conf文件中去掉IPv6的井号,如下图 3.结束pg进程 4.重启PG服务.
- django-cms 代码研究(一)djangocms是什么
首先用djangocms生成了一个站点(具体参考这里:http://www.cnblogs.com/Tommy-Yu/p/3878488.html),其文件结构如下: 本来以为会很有逼格,结果一看傻眼 ...
- ASP.NET - 视图状态概述
本文转载自dodream 视图状态是 ASP.NET 页框架用于在往返过程之间保留页和控件值的方法.在呈现页的 HTML 标记时,必须在回发过程中保留的页和值的当前状态将被序列化为Base64 编码字 ...
- Linux下提取IP至文件
ifconfig | grep 'inet[^6]' | sed 's/^\s*//g' | cut -d ' ' -f2 > ips.txt 排除127开头的IP: ifconfig | gr ...
- FastReport安装说明(中文版)
FastReport安装说明(中文版) 内容列表 I. IntroductionI. 介绍II. Manual installing of the FastReport packagesII. 手动安 ...
- Java for LeetCode 051 N-Queens
Given an integer n, return all distinct solutions to the n-queens puzzle. Each solution contains a d ...
- UVA11806-Cheerleaders(容斥原理+二进制)
In most professional sporting events, cheerleaders play a major role in entertaining the spectators. ...