BZOJ 3694 最短路
233333想简单了。。。。
题解参见http://hzwer.com/3710.html
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxv 4050
#define maxe 200500
#define inf 0x7f7f7f7f7f7f7f7fLL
using namespace std;
struct edge
{
long long u,v,w,flag,nxt;
}e[maxe];
long long n,m,x,y,a,b,g[maxv],nume=,dis[maxv];
long long top[maxv],fath[maxv],dfn[maxv],dep[maxv],size[maxv],son[maxv],times=;
long long root,tot=,ls[maxv<<],rs[maxv<<],lazy[maxv<<];
void addedge(long long u,long long v,long long w,long long flag)
{
e[++nume].u=u;
e[nume].v=v;
e[nume].w=w;
e[nume].flag=flag;
e[nume].nxt=g[u];
g[u]=nume;
}
void dfs1(long long x)
{
size[x]=;son[x]=;
for (long long i=g[x];i;i=e[i].nxt)
{
long long v=e[i].v;
if ((e[i].flag) && (v!=fath[x]))
{
fath[v]=x;dep[v]=dep[x]+;dis[v]=dis[x]+e[i].w;
dfs1(v);
size[x]+=size[v];
if (size[v]>size[son[x]]) son[x]=v;
}
}
}
void dfs2(long long x,long long father)
{
dfn[x]=++times;top[x]=father;
if (son[x]) dfs2(son[x],father);
for (long long i=g[x];i;i=e[i].nxt)
{
long long v=e[i].v;
if ((e[i].flag) && (v!=fath[x]) && (v!=son[x]))
dfs2(v,v);
}
}
void build(long long &now,long long left,long long right)
{
now=++tot;lazy[now]=inf;
if (left==right) return;
long long mid=(left+right)>>;
build(ls[now],left,mid);
build(rs[now],mid+,right);
}
void modify(long long now,long long left,long long right,long long l,long long r,long long x)
{
if ((left==l) && (right==r))
{
lazy[now]=min(lazy[now],x);
return;
}
long long mid=(left+right)>>;
if (r<=mid) modify(ls[now],left,mid,l,r,x);
else if (l>=mid+) modify(rs[now],mid+,right,l,r,x);
else {modify(ls[now],left,mid,l,mid,x);modify(rs[now],mid+,right,mid+,r,x);}
}
long long lca(long long u,long long v)
{
long long f1=top[u],f2=top[v];
while (f1!=f2)
{
if (dep[f1]<dep[f2]) {swap(f1,f2);swap(u,v);}
u=fath[f1];f1=top[u];
}
if (dep[u]<dep[v]) return u;
else return v;
}
void work(long long x)
{
long long u=e[x].u,v=e[x].v,ret=dis[u]+e[x].w+dis[v];
long long t=lca(u,v);
if (t==v) return;
u=t;
long long f1=top[u],f2=top[v];
while (f1!=f2)
{
if (dep[f1]<dep[f2]) {swap(f1,f2);swap(u,v);}
modify(root,,n,dfn[f1],dfn[u],ret);
u=fath[f1];f1=top[u];
}
if (u!=v)
{
if (dep[u]>dep[v]) swap(u,v);
modify(root,,n,dfn[u]+,dfn[v],ret);
}
}
long long ask(long long now,long long left,long long right,long long pos)
{
if (left==right)
return lazy[now];
long long mid=(left+right)>>;
if (pos<=mid) return min(ask(ls[now],left,mid,pos),lazy[now]);
else return min(ask(rs[now],mid+,right,pos),lazy[now]);
}
int main()
{
scanf("%lld%lld",&n,&m);
for (long long i=;i<=m;i++)
{
scanf("%lld%lld%lld%lld",&x,&y,&a,&b);
addedge(x,y,a,b);
addedge(y,x,a,b);
}
dfs1();
dfs2(,);
build(root,,n);
for (long long i=;i<=nume;i++)
{
long long u=e[i].u,v=e[i].v,ret=dis[u]+e[i].w+dis[v];
if (e[i].flag) continue;
work(i);
}
for (long long i=;i<=n;i++)
{
long long ret=ask(root,,n,dfn[i]);
if (ret==inf) printf("-1 ");
else printf("%lld ",ret-dis[i]);
}
printf("\n");
return ;
}
BZOJ 3694 最短路的更多相关文章
- Bzoj 3694: 最短路 树链剖分
3694: 最短路 Time Limit: 5 Sec Memory Limit: 256 MBSubmit: 67 Solved: 34[Submit][Status][Discuss] Des ...
- BZOJ 3694&&DTOJ 1972: 最短路
题目描述 给出一个n个点m条边的无向图,n个点的编号从1~n,定义源点为1.定义最短路树如下:从源点1经过边集T到任意一点i有且仅有一条路径,且这条路径是整个图1到i的最短路径,边集T构成最短路树. ...
- BZOJ 3040 最短路 (堆优化dijkstra)
这题不是裸的最短路么?但是一看数据范围就傻了.点数10^6,边数10^7.这个spfa就别想了(本来spfa就是相当不靠谱的玩意),看来是要用堆优化dijkstra了.但是,平时写dijkstra时为 ...
- BZOJ 3040: 最短路(road) ( 最短路 )
本来想学一下配对堆的...结果学着学着就偏了... 之前 kpm 写过这道题 , 前面的边不理它都能 AC .. 我也懒得去写前面的加边了... 用 C++ pb_ds 库里的 pairing_hea ...
- BZOJ 3040: 最短路(road) [Dijkstra + pb_ds]
3040: 最短路(road) Time Limit: 60 Sec Memory Limit: 200 MBSubmit: 2476 Solved: 814[Submit][Status][Di ...
- 【刷题】BZOJ 2125 最短路
Description 给一个N个点M条边的连通无向图,满足每条边最多属于一个环,有Q组询问,每次询问两点之间的最短路径. Input 输入的第一行包含三个整数,分别表示N和M和Q 下接M行,每行三个 ...
- BZOJ 2125: 最短路
2125: 最短路 Time Limit: 1 Sec Memory Limit: 259 MBSubmit: 756 Solved: 331[Submit][Status][Discuss] D ...
- BZOJ.2125.最短路(仙人掌 最短路Dijkstra)
题目链接 多次询问求仙人掌上两点间的最短路径. 如果是在树上,那么求LCA就可以了. 先做着,看看能不能把它弄成树. 把仙人掌看作一个图(实际上就是),求一遍根节点到每个点的最短路dis[i]. 对于 ...
- BZOJ 3040最短路
题目描述 给定一个 NN 个点, MM 条有向边的带权图,请你计算从 SS 出发,到每个点的距离. 数据保证你能从 SS 出发到任意点. 输入输出格式 输入格式: 第一行两个整数 NN . MM ,表 ...
随机推荐
- php调试工具firephp
PHP调试有很多中方法.今天给大家推荐的是 FirePHP 这个工具,是基于firefox的一个组件.Firephp组件安装地址.使用之前需要先安装下这个组件.完了还需要down一个FirePHP的包 ...
- Java 8怎么了之二:函数和原语
[编者按]本文作者为专注于自然语言处理多年的 Pierre-Yves Saumont,Pierre-Yves 著有30多本主讲 Java 软件开发的书籍,自2008开始供职于 Alcatel-Luce ...
- 导出含有图片的项目成jar文件后运行,图片不显示
在编写完Java程序后,打包成Jar时发布,会发现找不到Jar文件中的图片和文本文件,其原因是程序中载入图片或文本文件时,使用了以当前工作路径为基准的方式来指定文件和路径.这与用户运行Jar包时的当前 ...
- 详解 $_SERVER 函数中QUERY_STRING和REQUEST_URI区别
详解 $_SERVER 函数中QUERY_STRING和REQUEST_URI区别 http://blog.sina.com.cn/s/blog_686999de0100jgda.html 实例: ...
- VS2012简单的使用感受+插件推荐
VS2012简单的使用感受+插件推荐http://www.cnblogs.com/tangge/archive/2013/03/12/2955367.html
- Linux下Sublime Text 2的安装
安装方法1: 通过apt-get install来安装,可以如下来做: sudo add-apt-repository ppa:webupd8team/sublime-text-2 sudo apt- ...
- 重载和覆盖的区别?(overload vs override)
override与overload的区别? override 表示重写,overload 表示重载.override是子类和父类之间的关系,是垂直关系:overload是同一个类中方法之间的关系,是水 ...
- 安装DirectX SDK (June 2010) 失败(Error Code S1023)(转)
解决: 在控制面板里面搜索“Microsoft Visual C++*”,然后卸载,这个工具的名字大概是: Microsoft Visual C++ 2010 x86/x64 redistribuab ...
- 2013 ACM/ICPC Asia Regional Changsha Online - C Color Representation Conversion
这个纯粹是一个细节题啊!!! 由于某个地方的浮点数比较写错了,WA了无数次啊…… 代码如下: #include<iostream> #include<cstdio> #incl ...
- js正则表达式用法大全
匹配中文字符的正则表达式: [u4e00-u9fa5] 评注:匹配中文还真是个头疼的事,有了这个表达式就好办了 匹配双字节字符(包括汉字在内):[^x00-xff] 评注:可以用来计算字符串的长度(一 ...