[BZOJ1880] [Sdoi2009] Elaxia的路线 (SPFA & 拓扑排序)
Description
Input
Output
Sample Input
1 6 7 8
1 2 1
2 5 2
2 3 3
3 4 2
3 9 5
4 5 3
4 6 4
4 7 2
5 8 1
7 9 1
Sample Output
HINT
对于30%的数据,N ≤ 100;
对于60%的数据,N ≤ 1000;
对于100%的数据,N ≤ 1500,输入数据保证没有重边和自环。
Source
Solution
补个条件:$m\leq 500000$
如果$dis_{s->u_i}+w_i=dis_{t->v_i}$,那么边$i$才可能成为答案
这些边组成的图一定是一个拓扑图,走一遍最长链即可。
其实主要的坑点在于因为是无向图,所以需要反着做一遍
也就是说,$x_1$->$y_1$和$y_2$->$x_2$的公共路径也可能是答案,也就是说,原题意是错的= =b
#include <bits/stdc++.h>
using namespace std;
struct edge
{
int v, w, nxt;
}e[];
int fst[][], dis[][], q[], indeg[];
int n, etot, sss1, ttt1, sss2, ttt2;
bool inq[]; void addedge(int *f, int u, int v, int w)
{
e[++etot] = (edge){v, w, f[u]}, f[u] = etot;
} bool check(int u, int i)
{
if(dis[][u] + e[i].w + dis[][e[i].v] != dis[][ttt1]) return false;
return dis[][u] + e[i].w + dis[][e[i].v] == dis[][ttt2];
} void SPFA(int sss, int *d)
{
int front = , back;
memset(d, , );
q[back = ] = sss, d[sss] = , inq[sss] = true;
while(front != back)
{
int u = q[++front & ];
front &= , inq[u] = false;
for(int i = fst[][u]; i; i = e[i].nxt)
if(d[e[i].v] > d[u] + e[i].w)
{
d[e[i].v] = d[u] + e[i].w;
if(!inq[e[i].v])
{
q[++back & ] = e[i].v;
back &= , inq[e[i].v] = true;
}
}
}
} int Topo_sort()
{
int front = , back = , ans = ;
for(int i = ; i <= n; ++i)
if(!indeg[i]) q[++back] = i;
while(front != back)
{
int u = q[++front];
for(int i = fst[][u]; i; i = e[i].nxt)
{
int v = e[i].v, w = e[i].w;
dis[][v] = max(dis[][v], dis[][u] + w);
if(!--indeg[e[i].v]) q[++back] = v;
}
}
for(int i = ; i <= n; ++i)
ans = max(ans, dis[][i]);
return ans;
} int main()
{
int m, u, v, w, ans;
scanf("%d%d", &n, &m);
scanf("%d%d%d%d", &sss1, &ttt1, &sss2, &ttt2);
for(int i = ; i <= m; ++i)
{
scanf("%d%d%d", &u, &v, &w);
addedge(fst[], u, v, w);
addedge(fst[], v, u, w);
}
SPFA(sss1, dis[]), SPFA(ttt1, dis[]);
SPFA(sss2, dis[]), SPFA(ttt2, dis[]);
for(int i = ; i <= n; ++i)
for(int j = fst[][i]; j; j = e[j].nxt)
if(check(i, j))
{
addedge(fst[], i, e[j].v, e[j].w);
++indeg[e[j].v];
}
ans = Topo_sort();
memset(fst[], , sizeof(fst[]));
memset(dis[], , sizeof(dis[]));
memset(indeg, , sizeof(indeg));
swap(sss2, ttt2);
SPFA(sss2, dis[]), SPFA(ttt2, dis[]);
for(int i = ; i <= n; ++i)
for(int j = fst[][i]; j; j = e[j].nxt)
if(check(i, j))
{
addedge(fst[], i, e[j].v, e[j].w);
++indeg[e[j].v];
}
ans = max(ans, Topo_sort());
printf("%d\n", ans);
return ;
}
[BZOJ1880] [Sdoi2009] Elaxia的路线 (SPFA & 拓扑排序)的更多相关文章
- bzoj1880: [Sdoi2009]Elaxia的路线(spfa,拓扑排序最长路)
1880: [Sdoi2009]Elaxia的路线 Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 1944 Solved: 759[Submit][St ...
- [luogu2149][bzoj1880][SDOI2009]Elaxia的路线【拓扑排序+最短路+DP】
题目描述 最近,Elaxia和w的关系特别好,他们很想整天在一起,但是大学的学习太紧张了,他们 必须合理地安排两个人在一起的时间. Elaxia和w每天都要奔波于宿舍和实验室之间,他们 希望在节约时间 ...
- BZOJ1880: [Sdoi2009]Elaxia的路线(最短路)
1880: [Sdoi2009]Elaxia的路线 Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 2049 Solved: 805 题目链接:https ...
- BZOJ1880:[SDOI2009]Elaxia的路线(最短路,拓扑排序)
Description 最近,Elaxia和w**的关系特别好,他们很想整天在一起,但是大学的学习太紧张了,他们 必须合理地安排两个人在一起的时间.Elaxia和w**每天都要奔波于宿舍和实验室之间, ...
- BZOJ1880: [Sdoi2009]Elaxia的路线
题意:求最短路最长公共距离. 考虑每一条边,如果满足dis(s1,u)+len+dis(v,t1)==dis(s1,t1) && dis(s2,u)+len+dis(v,t2)==di ...
- BZOJ1880 [Sdoi2009]Elaxia的路线 【最短路 + dp】
题目 最近,Elaxia和w的关系特别好,他们很想整天在一起,但是大学的学习太紧张了,他们 必须合理地安排两个人在一起的时间.Elaxia和w每天都要奔波于宿舍和实验室之间,他们 希望在节约时间的前提 ...
- 【BZOJ1880】[SDOI2009]Elaxia的路线 (最短路+拓扑排序)
[SDOI2009]Elaxia的路线 题目描述 最近,\(Elaxia\)和\(w**\)的关系特别好,他们很想整天在一起,但是大学的学习太紧张了,他们 必须合理地安排两个人在一起的时间. \(El ...
- BZOJ-1880 Elaxia的路线 SPFA+枚举
1880: [Sdoi2009]Elaxia的路线 Time Limit: 4 Sec Memory Limit: 64 MB Submit: 921 Solved: 354 [Submit][Sta ...
- 【BZOJ1880】[Sdoi2009]Elaxia的路线(最短路)
[BZOJ1880][Sdoi2009]Elaxia的路线(最短路) 题面 BZOJ 洛谷 题解 假装我们知道了任意两点间的最短路,那么我们怎么求解答案呢? 不难发现公共路径一定是一段连续的路径(如果 ...
随机推荐
- Linux下的压力测试工具:ab、http_load、webbench、siege
一.ab 1.1 介绍 ab是apache自带的一款功能强大的测试工具. 安装了apache一般就自带了. 1.2 下载 同apache. 1.3 安装 同apache. 1.4 安装结果 ...
- canvas-缩放
Canvas-图片缩放 由上一篇canvas-旋转的例子可以了解到canvas的一些特性,不熟悉的同学可以先去看看canvas-旋转. 我们在将图片引入canvas时,图片会一原始像素渲染.这样往往不 ...
- Jenkins系列——定时构建
1.环境说明 操作系统:win7旗舰版64bit jdk:sun JDK1.7.0_80 64bit tomcat:apache-tomcat-8.0.41 jenkins:2.32.3LST 本系列 ...
- 在kali安装中文输入法的教程
1终端下vi /etc/apt/sources.list 修改镜像元 (按E进行编辑 具体实例不同可能没有) 按 i进入编辑 擦除原有的几个官方源改为deb http://mirrors.ali ...
- Fantasia (Tarjan+树形DP)
Time Limit: 1000 ms Memory Limit: 256 MB Description 给定一张N个点.M条边的无向图 $G$ .每个点有个权值Wi. 我们定义 $G_i$ 为图 ...
- Docker mongodb 3.4 分片 一主 一副 一仲 鉴权集群部署.
非docker部署 为了避免过分冗余,并且在主节点挂了,还能顺利自动提升,所以加入仲裁节点 为什么要用docker部署,因为之前直接在虚拟机启动10个mongod 进程.多线程并发测试的时候,mong ...
- SPI知识总结
SPI知识总结 一.定义 SPI(Serial Peripheral Interface--串行外设接口)总线系统是一种同步串行外设接口,它可以使MCU与各种外围设备以串行方式进行通信以交换信息. 优 ...
- c#获取文件MD5算法
//获取文件MD5算法 private static string GetMD5FromFile(string fileName) { try { FileStream file = new File ...
- 关于 target="_blank"漏洞的分析
创建: 于 八月 30, 2016 关于 target="_blank"漏洞的分析 一.漏洞详情:首先攻击者能够将链接(指向攻击者自己控制的页面的,该被控页面的js脚本可以对母页 ...
- Linux CentOS7下安装python3
在CentOS7下,默认安装的就是python2.7,我现在来教大家如何安装python3: 1.首先安装python3.6可能使用的依赖 # yum -y install openssl-devel ...