LOJ#3087. 「GXOI / GZOI2019」旅行者

正着求一遍dij,反着求一遍,然后枚举每条边,从u到v,如果到u最近的点和v能到的最近的点不同,那么可以更新答案

没了

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 100005
#define ba 47
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
struct node {
int to,next;
int64 val;
}E[500005];
int N,M,K;
int head[MAXN],sumE,u[500005],v[500005],c[500005];
vector<int> st;
int64 dis[2][MAXN],from[2][MAXN];
bool vis[MAXN];
void add(int u,int v,int64 c) {
E[++sumE].to = v;
E[sumE].next = head[u];
E[sumE].val = c;
head[u] = sumE;
}
void Init() {
read(N);read(M);read(K);
int a,b;
for(int i = 1 ; i <= M ; ++i) {
read(u[i]);read(v[i]);read(c[i]);
}
st.clear();
for(int i = 1 ; i <= K ; ++i) {
read(a);st.pb(a);
}
}
priority_queue<pair<int64,int> > Q;
void Dijkstra(int id) {
for(int i = 1 ; i <= N ; ++i) {
dis[id][i] = 1e18;
}
for(auto t : st) {
dis[id][t] = 0;from[id][t] = t;
Q.push(mp(-dis[id][t],t));
}
memset(vis,0,sizeof(vis));
while(!Q.empty()) {
auto now = Q.top();Q.pop();
if(vis[now.se]) continue;
int t = now.se;
vis[t] = 1;
for(int i = head[t] ; i ; i = E[i].next) {
int v = E[i].to;
if(dis[id][v] > dis[id][t] + E[i].val) {
dis[id][v] = dis[id][t] + E[i].val;
from[id][v] = from[id][t];
Q.push(mp(-dis[id][v],v));
}
}
}
}
void Solve() {
sumE = 0;memset(head,0,sizeof(head));
for(int i = 1 ; i <= M ; ++i) {
add(u[i],v[i],c[i]);
} Dijkstra(0);
sumE = 0;memset(head,0,sizeof(head));
for(int i = 1 ; i <= M ; ++i) {
add(v[i],u[i],c[i]);
}
Dijkstra(1);
int64 ans = 1e18;
for(int i = 1 ; i <= M ; ++i) {
if(from[0][u[i]] != from[1][v[i]]) {
ans = min(ans,c[i] + dis[0][u[i]] + dis[1][v[i]]);
}
}
out(ans);enter;
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
int T;
read(T);
while(T--) {
Init();
Solve();
}
}

最新文章

  1. iOS沙盒处理
  2. spark1.5引进内置函数
  3. 安装sql server提示挂起报错
  4. Sliding Window Maximum 解答
  5. Linux MySql安装步骤
  6. 浅谈访问控制列表(ACL)
  7. [LeetCode] Maximum Average Subarray II 子数组的最大平均值之二
  8. 【STM32H7教程】第4章 STM32H7工程模板建立(MDK5)
  9. .NET Core微服务之基于App.Metrics+InfluxDB+Grafana实现统一性能监控
  10. [Localization] SSD - Single Shot MultiBoxDetector
  11. tf.nn.embedding_lookup函数【转载】
  12. 架构师素养及从小菜进阶架构(CTO)的书籍【转】
  13. PAT 1021 Deepest Root[并查集、dfs][难]
  14. Informatica 常用组件Filter之一 概述
  15. codeforces#512 Div2
  16. 跟着百度学PHP[9]-session与cookie的异同
  17. yii2使用 db log
  18. 中文 bootstrapValidator
  19. Controller与View数据传递 多Model传递
  20. Crontab 学习

热门文章

  1. 「BZOJ 4565」「HAOI 2016」字符合并「区间状压DP」
  2. ie中兼容性问题
  3. 2016 Multi-University Training Contest 2 部分题解
  4. cs231n lecture2 image classification
  5. arcgis python desc.dataType
  6. 上有传参下传json的接口调用
  7. jQuery常用Method-API
  8. .netcore centos配置systemctl自动启动
  9. SSM批量添加数据
  10. Jmeter使用实践- 基础介绍