线段树优化建图 || CF786B Legacy
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#define ll long long
using namespace std;
inline ll rd(){
ll x=,f=;char c=getchar();
while(c<''||c>''){if(c=='-')f=-; c=getchar();}
while(c>=''&&c<=''){x=x*+c-''; c=getchar();}
return f*x;
}
const int maxn=(1e5)+,maxq=maxn;
const ll inf=1ll<<;
int N,S,Q,o,v,u,ql,qr,num_edge=,edge_head[maxn*];
int tr[maxn*][],cnt;
ll w,Dis[maxn*];
bool vis[maxn*];
struct Num_{
int id;ll dis;
bool operator < (const Num_&a) const{
return a.dis<dis;
}
}e;
priority_queue<Num_>pri;
struct Edge{ int to,nx; ll dis; }edge[maxn*];
inline void Add_edge(int from,int to,int dis){
edge[++num_edge].nx=edge_head[from];
edge[num_edge].to=to;
edge[num_edge].dis=dis;
edge_head[from]=num_edge;
return;
}
void Build(int x,int l,int r){
if(l==r){
tr[x][]=tr[x][]=l;
return;
}
int m=(l+r)>>;
tr[x][]=++cnt; tr[x][]=++cnt;
Build(x<<,l,m); Build(x<<|,m+,r);
Add_edge(tr[x][],tr[x<<][],);
Add_edge(tr[x][],tr[x<<|][],);
Add_edge(tr[x<<][],tr[x][],);
Add_edge(tr[x<<|][],tr[x][],);
return;
}
void Update(int x,int l,int r,int ql,int qr,int u,int val,int op){
if(ql<=l&&r<=qr){
if(op==) Add_edge(u,tr[x][],val);
else Add_edge(tr[x][],u,val);
return;
}
int m=(l+r)>>;
if(ql<=m)Update(x<<,l,m,ql,qr,u,val,op);
if(qr>m)Update(x<<|,m+,r,ql,qr,u,val,op);
return;
}
inline void Dijkstra(int s){
for(int i=;i<=cnt;i++)Dis[i]=inf;
Dis[s]=;
e.id=s;e.dis=;
pri.push(e);
while(!pri.empty()){
int x=(pri.top()).id;
pri.pop();
if(vis[x])continue;
vis[x]=;
for(int i=edge_head[x];i;i=edge[i].nx){
int y=edge[i].to;
if(Dis[x]+edge[i].dis<Dis[y]){
Dis[y]=Dis[x]+edge[i].dis;
e.id=y;e.dis=Dis[y];
pri.push(e);
}
}
}
return;
}
int main(){
N=rd();Q=rd();S=rd();
cnt=N;
Build(,,N);
while(Q--){
o=rd();
if(o==){
v=rd();u=rd();w=rd();
Add_edge(v,u,w);
}
else {
v=rd();ql=rd();qr=rd();w=rd();
if(o==)
Update(,,N,ql,qr,v,w,);
else//o==3
Update(,,N,ql,qr,v,w,);
}
}
Dijkstra(S);
for(int i=;i<=N;i++)
if(Dis[i]!=inf)printf("%lld ",Dis[i]);
else printf("-1 ");
return ;
}
邻接表建图版
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#define make(a,b) (make_pair(a,b))
#define ll long long
using namespace std;
inline ll rd(){
ll x=,f=;char c=getchar();
while(c<''||c>''){if(c=='-')f=-; c=getchar();}
while(c>=''&&c<=''){x=x*+c-''; c=getchar();}
return f*x;
}
const int maxn=(1e5)+,maxq=maxn;
const ll inf=1ll<<;
int N,Q,S,trIn[maxn<<],trOut[maxn<<];
int u,v,ql,qr,o,cnt;
bool vis[maxn<<];
vector<pair<int,ll> >edge[maxn<<];//to,val
ll w,Dis[maxn<<];
struct Num_{
ll dis;
int id;
bool operator < (const Num_&a) const {
return a.dis<dis;
}
}b;
priority_queue<Num_>pri;
void Build(int x,int l,int r){
if(l==r){
trIn[x]=l;
trOut[x]=l;
return;
}
int mid=(l+r)>>;
trIn[x]=++cnt;trOut[x]=++cnt;
Build(x<<,l,mid);Build(x<<|,mid+,r);
edge[trIn[x]].push_back(make(trIn[x<<],));
edge[trIn[x]].push_back(make(trIn[x<<|],));
edge[trOut[x<<]].push_back(make(trOut[x],));
edge[trOut[x<<|]].push_back(make(trOut[x],));
return;
}
void Update(int x,int l,int r,int ql,int qr,int u,int val,int o){
if(ql<=l&&r<=qr){
if(o==)//In
edge[u].push_back(make(trIn[x],val));
else//Out
edge[trOut[x]].push_back(make(u,val));
return;
}
int mid=(l+r)>>;
if(ql<=mid)Update(x<<,l,mid,ql,qr,u,val,o);
if(qr>mid)Update(x<<|,mid+,r,ql,qr,u,val,o);
return;
}
inline void Dijkstra(int s){
for(int i=;i<=cnt;i++)Dis[i]=inf;
Dis[s]=;
b.dis=;b.id=s;
pri.push(b);
while(!pri.empty()){
int x=(pri.top()).id;
pri.pop();
if(vis[x])continue;
vis[x]=;
int toi=edge[x].size();
for(int i=;i<toi;i++){
int y=edge[x][i].first;
if(Dis[y]>Dis[x]+edge[x][i].second){
Dis[y]=Dis[x]+edge[x][i].second;
b.id=y;b.dis=Dis[y];
pri.push(b);
}
}
}
return;
}
int main(){
N=rd();Q=rd();S=rd();
cnt=N;
Build(,,N);
while(Q--){
o=rd();
if(o==){
v=rd();u=rd();w=rd();
edge[v].push_back(make(u,w));
}
else if(o==){
v=rd();ql=rd();
qr=rd();w=rd();
Update(,,N,ql,qr,v,w,);
}
else {//o==3
v=rd();ql=rd();
qr=rd();w=rd();
Update(,,N,ql,qr,v,w,);
}
}
Dijkstra(S);
for(int i=;i<=N;i++){
if(Dis[i]!=inf) printf("%lld ",Dis[i]);
else printf("-1 ");
}
return ;
}
Vector建图版
By:AlenaNuna
线段树优化建图 || CF786B Legacy的更多相关文章
- CF786B Legacy 线段树优化建图 + spfa
CodeForces 786B Rick和他的同事们做出了一种新的带放射性的婴儿食品(???根据图片和原文的确如此...),与此同时很多坏人正追赶着他们.因此Rick想在坏人们捉到他之前把他的遗产留给 ...
- G. 神圣的 F2 连接着我们 线段树优化建图+最短路
这个题目和之前写的一个线段树优化建图是一样的. B - Legacy CodeForces - 787D 线段树优化建图+dij最短路 基本套路 之前这个题目可以相当于一个模板,直接套用就可以了. 不 ...
- BZOJ5017 [SNOI2017]炸弹 - 线段树优化建图+Tarjan
Solution 一个点向一个区间内的所有点连边, 可以用线段树优化建图来优化 : 前置技能传送门 然后就得到一个有向图, 一个联通块内的炸弹可以互相引爆, 所以进行缩点变成$DAG$ 然后拓扑排序. ...
- 【BZOJ3681】Arietta 树链剖分+可持久化线段树优化建图+网络流
[BZOJ3681]Arietta Description Arietta 的命运与她的妹妹不同,在她的妹妹已经走进学院的时候,她仍然留在山村中.但是她从未停止过和恋人 Velding 的书信往来.一 ...
- 【ARC069F】Flags 2-sat+线段树优化建图+二分
Description 数轴上有 n 个旗子,第 ii 个可以插在坐标 xi或者 yi,最大化两两旗子之间的最小距离. Input 第一行一个整数 N. 接下来 N 行每行两个整数 xi, ...
- 【bzoj5017】[Snoi2017]炸弹 线段树优化建图+Tarjan+拓扑排序
题目描述 在一条直线上有 N 个炸弹,每个炸弹的坐标是 Xi,爆炸半径是 Ri,当一个炸弹爆炸时,如果另一个炸弹所在位置 Xj 满足: Xi−Ri≤Xj≤Xi+Ri,那么,该炸弹也会被引爆. 现在 ...
- 【bzoj4699】树上的最短路(树剖+线段树优化建图)
题意 给你一棵 $n$ 个点 $n-1$ 条边的树,每条边有一个通过时间.此外有 $m$ 个传送条件 $(x_1,y_1,x_2,y_2,c)$,表示从 $x_1$ 到 $x_2$ 的简单路径上的点可 ...
- 【BZOJ4276】[ONTAK2015]Bajtman i Okrągły Robin 线段树优化建图+费用流
[BZOJ4276][ONTAK2015]Bajtman i Okrągły Robin Description 有n个强盗,其中第i个强盗会在[a[i],a[i]+1],[a[i]+1,a[i]+2 ...
- 【bzoj3073】[Pa2011]Journeys 线段树优化建图+堆优化Dijkstra
题目描述 Seter建造了一个很大的星球,他准备建造N个国家和无数双向道路.N个国家很快建造好了,用1..N编号,但是他发现道路实在太多了,他要一条条建简直是不可能的!于是他以如下方式建造道路:(a, ...
随机推荐
- N2RR
一.背景 氨(NH 3)是一种新兴的能量载体,在液氨中含有17.6%(重量)的氢,而在甲醇中的含量为12.5%(重量),很可能成为未来氢能经济的有希望的候选者.然而,如果基于NH 3的肥料不足以养活世 ...
- golang 通过reflect反射修改值
不是所有的反射值都可以修改.对于一个反射值是否可以修改,可以通过CanSet()进行检查. 要修改值,必须满足: 可以寻址 可寻址的类型: 指针指向的具体元素 slice的元素 可寻址的结构体的字段( ...
- Jmeter(一) - 调用数据的参数化
1. 做性能测试, 不可避免的一点一定会有使用不同的用户密码进行登陆. 如何使登陆用户参数化呢?
- CSS3 —— 盒子模型
盒子模型 主要的属性就5个:width.height.padding.border.margin.如下: width和height:内容的宽度.高度(不是盒子的宽度.高度). padding:内边距 ...
- 使用UI Automation实现自动化测试--1-4
Introduction UI Automation是Microsoft .NET 3.0框架下提供的一种用于自动化测试的技术,是在MSAA基础上建立的,MSAA就是Microsoft Active ...
- 100+ Python挑战性编程练习(1)
目前,这个项目已经获得了7.1k Stars,4.1k Forks. 初级水平是指刚刚通过Python入门课程的人.他可以用1或2个Python类或函数来解决一些问题.通常,答案可以直接在教科书中找 ...
- ansible-playbook 案例
nginx的安装 编写nginx的自动部署文件nginx.yml hosts主机更改为自己定义的 访问目标主机组的IP地址,查看测试页面 测试页面:显示的是本机ip 1 <h1> ...
- MyBatis按时间排序
测试代码 ActivityReadExample readExample = new ActivityReadExample(); readExample.setOrderByClause(" ...
- PHP 数组辅助函数
/** * 取多维数据中某字段的值 * @param array $array 数据源数组 * @param string|array $field 要获取的字段 * @return array 结果 ...
- javascript异步延时加载及判断是否已加载js/css文件
<html> <head> <script type="text/javascript"> /**======================= ...