bzoj 1880: [Sdoi2009]Elaxia的路线(拓扑排序+spfa)
来源:互联网 发布:淘宝助理有什么好处 编辑:程序博客网 时间:2024/06/05 11:31
1880: [Sdoi2009]Elaxia的路线
Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 974 Solved: 382
[Submit][Status][Discuss]
Description
最近,Elaxia和w**的关系特别好,他们很想整天在一起,但是大学的学习太紧张了,他们 必须合理地安排两个人在一起的时间。Elaxia和w**每天都要奔波于宿舍和实验室之间,他们 希望在节约时间的前提下,一起走的时间尽可能的长。 现在已知的是Elaxia和w**所在的宿舍和实验室的编号以及学校的地图:地图上有N个路 口,M条路,经过每条路都需要一定的时间。 具体地说,就是要求无向图中,两对点间最短路的最长公共路径。
Input
第一行:两个整数N和M(含义如题目描述)。 第二行:四个整数x1、y1、x2、y2(1 ≤ x1 ≤ N,1 ≤ y1 ≤ N,1 ≤ x2 ≤ N,1 ≤ ≤ N),分别表示Elaxia的宿舍和实验室及w**的宿舍和实验室的标号(两对点分别 x1,y1和x2,y2)。 接下来M行:每行三个整数,u、v、l(1 ≤ u ≤ N,1 ≤ v ≤ N,1 ≤ l ≤ 10000),表 u和v之间有一条路,经过这条路所需要的时间为l。 出出出格格格式式式::: 一行,一个整数,表示每天两人在一起的时间(即最长公共路径的长度)。
Output
一行,一个整数,表示每天两人在一起的时间(即最长公共路径的长度)
Sample Input
9 10
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
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
3
HINT
对于30%的数据,N ≤ 100;
对于60%的数据,N ≤ 1000;
对于100%的数据,N ≤ 1500,输入数据保证没有重边和自环。
Source
Day2
以s1,t1,s2,t2为源点分别跑spfa ,记录下单源最短路,然后通过 f[u]+len[u,v]+g[v]=len的方式把所有在两个人最短路上的公共边加入有向无环图,然后利用拓扑排序求出最长链即为答案。
因为可能存在路径方向的问题,所有我们把s2,t2反过来再计算一边,就可以了 。
附样例一组:(可以卡掉某些AC程序)
8 6
1 7 3 8
1 2 1
1 4 1
4 8 1
2 3 1
2 7 3
4 7 3 output: 1
有些程序把1 7 3 8 改成3 8 1 7 出的答案会不一样
#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<queue>#include<cmath>#define N 2000using namespace std;int n,m,s1,t1,s2,t2,len1,len2,ans;int f1[N],g1[N],f2[N],g2[N],can[N],ins[N],q[N];int point[N],next[1000000],v[1000000],len[1000000],tot,cnt;int point1[N],next1[1000000],v1[1000000],c[1000000];int maxn[N],f[N];void insert(int x,int y,int z){tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; len[tot]=z;tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; len[tot]=z;}void add(int x,int y,int z){cnt++; next1[cnt]=point1[x]; point1[x]=cnt; v1[cnt]=y; c[cnt]=z;}void spfa(int *dis,int s){memset(can,0,sizeof(can));dis[s]=0; can[s]=1;queue<int> p;p.push(s);while (!p.empty()){int now=p.front(); p.pop();for (int i=point[now];i;i=next[i]) if (dis[v[i]]>dis[now]+len[i]) { dis[v[i]]=dis[now]+len[i]; if (can[v[i]]==0) { can[v[i]]=1; p.push(v[i]); } }can[now]=0;}}void solve(){memset(f,0,sizeof(f));int head=0,tail=0;for (int i=1;i<=n;i++) if (!ins[i]) q[++tail]=i;while (head<tail){int now=q[++head];for (int i=point1[now];i;i=next1[i]) { f[v1[i]]=f[now]+c[i]; if (!--ins[v1[i]]) q[++tail]=v1[i]; }}for (int i=1;i<=n;i++) ans=max(ans,f[i]);}int main(){ scanf("%d%d",&n,&m); scanf("%d%d%d%d",&s1,&t1,&s2,&t2); for (int i=1;i<=m;i++) { int x,y,z; scanf("%d%d%d",&x,&y,&z); insert(x,y,z); } memset(f1,0x7f,sizeof(f1)); memset(g1,0x7f,sizeof(g1)); memset(f2,0x7f,sizeof(f2)); memset(g2,0x7f,sizeof(g2)); spfa(f1,s1); len1=f1[t1]; spfa(g1,t1); spfa(f2,s2); len2=f2[t2]; spfa(g2,t2); for (int i=1;i<=n;i++) for (int j=point[i];j;j=next[j]) if (f1[i]+len[j]+g1[v[j]]==len1&&f2[i]+len[j]+g2[v[j]]==len2) add(i,v[j],len[j]),ins[v[j]]++; solve(); memset(f2,0x7f,sizeof(f2)); memset(g2,0x7f,sizeof(g2)); spfa(f2,t2); spfa(g2,s2); cnt=0; memset(ins,0,sizeof(ins)); memset(point1,0,sizeof(point1)); memset(next1,0,sizeof(next1)); for (int i=1;i<=n;i++) for (int j=point[i];j;j=next[j]) if (f1[i]+len[j]+g1[v[j]]==len1&&f2[i]+len[j]+g2[v[j]]==len2) add(i,v[j],len[j]),ins[v[j]]++; solve(); printf("%d\n",ans);}
0 0
- bzoj 1880: [Sdoi2009]Elaxia的路线(拓扑排序+spfa)
- BZOJ 1880 Sdoi2009 Elaxia的路线 SPFA+拓扑排序
- BZOJ-1880-Elaxia的路线-SDOI2009-SPFA+拓扑排序
- [BZOJ]1880: [Sdoi2009]Elaxia的路线 spfa+拓扑图DP
- [BZOJ1880][SDOI2009]Elaxia的路线(SPFA+拓扑排序+DP)
- SDOI2009 Elaxia的路线(最短路+拓扑排序)
- [BZOJ1880][Sdoi2009] Elaxia的路线 (最短路,拓扑排序)
- (bzoj 1880 [Sdoi2009]Elaxia的路线)<最短路+拓扑最长链>
- BZOJ1880[Sdoi2009]Elaxia的路线 spfa+拓扑序
- 【bzoj 1880】: [Sdoi2009]Elaxia的路线
- 【BZOJ 1880】 [Sdoi2009]Elaxia的路线
- 【39.87%】【BZOJ 1880】[Sdoi2009]Elaxia的路线
- BZOJ 1880 [Sdoi2009]Elaxia的路线
- BZOJ-1880 Elaxia的路线 SPFA+枚举
- bzoj1880: [Sdoi2009]Elaxia的路线 wikioi2309 SPFA
- [SDOI2009]Elaxia的路线 SPFA+Topo
- 洛谷P2149 [SDOI2009]Elaxia的路线(最短路,拓扑排序)
- 1880: [Sdoi2009]Elaxia的路线
- 必须知道的Android屏幕自适应解决方案
- 《Android 开发艺术探索》笔记——(3)View 的事件体系
- 生产者/消费者问题java实现
- Linux中后面带有波浪号的文件
- Logstash之日志多行合一行(日志错行)
- bzoj 1880: [Sdoi2009]Elaxia的路线(拓扑排序+spfa)
- 写matlab程序中用记录下来的小碎片
- C运行时库(C Run-time Library)详解
- 项目杂记(一)
- interest point VS descriptor
- 如何在shell中判断一个文件是否为空
- 设计模式 单例模式
- STL:二级空间配置器浅谈
- 一种简便的ios图片加密方法-对图片进行base64编码