【题意】给定n个点的图,正权无向边,正负权有向边,保证对有向边(u,v),v无法到达u,求起点出发到达所有点的最短距离。

【算法】拓扑排序+dijkstra

【题解】因为有负权边,直接对原图进行spfa,加slf优化后可过,但是这道题就没意思了。

理论上,最短路问题用spfa是不能保证复杂度的,但dijkstra的问题是不能处理负权边。

因为题目保证不能返回,实际上有向边将全图分成了几个部分。如果把仅由无向边连接的连通块看成点,则原图变成DAG。

对连通块内部进行dijkstra,在DAG上用拓扑序递推计算就可以保证O(n log n)出解。

具体实现:全图共用最短距离数组d[]。

对跨越连通块的有向边建新图,先用部分拓扑序删掉不从s出发的点。

开连通块个数的堆。

然后拓扑排序的过程中将到达别的连通块的点加入对应的堆,dijkstra时直接开始不用设置初始状态。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cctype>
#include<queue>
using namespace std;
const int maxn=,inf=0x3f3f3f3f;
int first[maxn],FIRST[maxn],d[maxn],col[maxn],in[maxn],A[maxn],B[maxn];
int tot,cnt,n,m,N,M,s,color;
struct edge{int v,w,from;}e[maxn*],E[maxn*];
struct cyc{
int x,d;
bool operator < (const cyc &a)const{
return d>a.d;//
}
};
priority_queue<cyc>q[maxn];
int read(){
char c;int s=,t=;
while(!isdigit(c=getchar()))if(c=='-')t=-;
do{s=s*+c-'';}while(isdigit(c=getchar()));
return s*t;
}
void insert(int u,int v,int w){tot++;e[tot].v=v;e[tot].w=w;e[tot].from=first[u];first[u]=tot;}
void INSERT(int u,int v,int w){cnt++;E[cnt].v=v;E[cnt].w=w;E[cnt].from=FIRST[u];FIRST[u]=cnt;in[v]++;}
void dfs(int x,int color){
col[x]=color;
for(int i=first[x];i;i=e[i].from)if(!col[e[i].v]){
dfs(e[i].v,color);
}
}
queue<int>Q;
void dijkstra(int k){
while(!q[k].empty()){
cyc y=q[k].top();q[k].pop();
if(y.d!=d[y.x])continue;
int x=y.x;
for(int i=first[x];i;i=e[i].from)if(d[e[i].v]>d[x]+e[i].w){
d[e[i].v]=d[x]+e[i].w;
q[k].push((cyc){e[i].v,d[e[i].v]});
}
}
}
int main(){
scanf("%d%d%d%d",&n,&m,&M,&s);
for(int i=;i<=m;i++){
int u=read(),v=read(),w=read();
insert(u,v,w);insert(v,u,w);
}
for(int i=;i<=n;i++)if(!col[i])dfs(i,++N);
for(int i=;i<=M;i++){
int u=read(),v=read(),w=read();
INSERT(col[u],col[v],w);
A[i]=u;B[i]=v;
}
for(int i=;i<=N;i++)if(!in[i]&&col[s]!=i)Q.push(i);
while(!Q.empty()){
int x=Q.front();Q.pop();
for(int i=FIRST[x];i;i=E[i].from){
in[E[i].v]--;
if(!in[E[i].v]&&col[s]!=E[i].v)Q.push(E[i].v);
}
}
Q.push(col[s]);
memset(d,0x3f,sizeof(d));d[s]=;q[col[s]].push((cyc){s,});
while(!Q.empty()){
int x=Q.front();Q.pop();
dijkstra(x);
for(int i=FIRST[x];i;i=E[i].from){
if(d[B[i]]>d[A[i]]+E[i].w){
d[B[i]]=d[A[i]]+E[i].w;
q[E[i].v].push((cyc){B[i],d[B[i]]});
}
in[E[i].v]--;
if(!in[E[i].v])Q.push(E[i].v);
}
}
for(int i=;i<=n;i++)if(d[i]<inf)printf("%d\n",d[i]);else printf("NO PATH\n");
return ;
}

dijkstra使用小根堆!每次加入距离最小的点。最致命的是写成大根堆也可以跑出答案,但是大数据就会很慢。

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

dijkstra使用小根堆!

【BZOJ】2200: [Usaco2011 Jan]道路和航线的更多相关文章

  1. [BZOJ 2200][Usaco2011 Jan]道路和航线 spfa+SLF优化

    Description Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查.他想把牛奶送到T个城镇 (1 <= T <= 25,000),编号为1T.这些城镇之间通过R条 ...

  2. bzoj 2200: [Usaco2011 Jan]道路和航线——拓扑+dijkstra

    Description Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查.他想把牛奶送到T个城镇 (1 <= T <= 25,000),编号为1T.这些城镇之间通过R条 ...

  3. BZOJ 2200: [Usaco2011 Jan]道路和航线

    Description Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查.他想把牛奶送到T个城镇 (1 <= T <= 25,000),编号为1T.这些城镇之间通过R条 ...

  4. bzoj 2200: [Usaco2011 Jan]道路和航线【spfa】

    直接跑最短路就行了--还不用判负环 #include<iostream> #include<cstdio> #include<queue> using namesp ...

  5. 2200: [Usaco2011 Jan]道路和航线 (拓扑排序+dijstra)

    Description Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查.他想把牛奶送到T个城镇 (1 <= T <= 25,000),编号为1T.这些城镇之间通过R条 ...

  6. [Usaco2011 Jan]道路和航线

    Description Farmer John正在一个新的销售区域对他的牛奶销售方案进行调查.他想把牛奶送到T个城镇 (1 <= T <= 25,000),编号为1T.这些城镇之间通过R条 ...

  7. bzoj2200: [Usaco2011 Jan]道路和航线

    先忽略航线,求出图中所有连通块,再用航线拓扑排序求出每个连通块的优先级 然后dijkstra时优先处理优先级高的块里的点就行了 ps:这题SPFA会TLE #include <iostream& ...

  8. BZOJ 2200--[Usaco2011 Jan]道路和航线(最短路&拓扑排序)

    2200: [Usaco2011 Jan]道路和航线 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 1128  Solved: 414[Submit] ...

  9. BZOJ 2199: [Usaco2011 Jan]奶牛议会

    2199: [Usaco2011 Jan]奶牛议会 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 375  Solved: 241[Submit][S ...

随机推荐

  1. 【IdentityServer4文档】- 贡献

    贡献 我们非常乐于接受社区贡献,但您应遵循一些指导原则,以便我们可以很方便的解决这个问题. 如何贡献? 最简单的方法是打开一个问题并开始讨论.然后,我们可以决定如何实现一个特性或一个变更.如果您即将提 ...

  2. 初入py

    1.下载工具sublime 我的网盘下载地址:https://pan.baidu.com/s/18-U1ZSg_zHoSAqUuvXj_PQ 直接解压即可 2.配置py27 在新建的文件里面编辑并保存 ...

  3. ICE checkbox 用法

    Hello everybody, I have a datable which contain multiple lines gotten from database, in the header o ...

  4. BZOJ 1834 网络扩容(最大流+费用流)

    对于第一问,直接求最大流. 对于第二问,建源点s和汇点t,s连1容量为INF,费用为0的边,n连t容量为最大流+k,费用为0的边.这样就把最大流限制为最多增加k了. 限制需要求扩充的最小费用,原图的边 ...

  5. HDOJ(HDU).1015 Safecracker (DFS)

    HDOJ(HDU).1015 Safecracker [从零开始DFS(2)] 从零开始DFS HDOJ.1342 Lotto [从零开始DFS(0)] - DFS思想与框架/双重DFS HDOJ.1 ...

  6. NOIP2017 列队——平衡树

    平衡树蒟蒻,敲了半天. 其实思路很简单,就是把许多个人合并成一个区间.必要的时候再拆开.(是不是和这个题的动态开点线段树有异曲同工之妙?) 每次操作最多多出来6个点. 理论上时间复杂度是nlogn,空 ...

  7. BZOJ1912 APIO2010 洛谷P3629 巡逻

    Description: 在一个地区中有 n 个村庄,编号为 1, 2, ..., n.有 n – 1 条道路连接着这些村 庄,每条道路刚好连接两个村庄,从任何一个村庄,都可以通过这些道路到达其 他任 ...

  8. eclipse里配置Android ndk环境,用eclipse编译.so文件

    做Android NDK开发时,c代码需要用ndk-build来进行编译,而java代码则需要用Android sdk编译. 编译c代码有两种方法: 一.写好c代码后,然后用cygwin搭建ndk-b ...

  9. [转]Android 如何根据网络地址获取网络图片方法

    http://blog.csdn.net/xiazdong/article/details/7724103 目录(?)[-] h2pre namecode classhtml stylefont-we ...

  10. 【DP】【P2340】奶牛会展

    传送门 Description 奶牛想证明它们是聪明而风趣的.为此,贝西筹备了一个奶牛博览会,她已经对N 头奶牛进行了面试,确定了每头奶牛的智商和情商. 贝西有权选择让哪些奶牛参加展览.由于负的智商或 ...