#include <iostream>
#include <iomanip>
#include <cstdio>
#include <vector>
#include <bitset>
#include <string>
#include <cstring>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <deque>
#include <utility>
#include <algorithm>
#include <random>
#include <cmath>
#include <cassert>
#include <climits>
#include <ctime>
#include <chrono>
/*
#pragma GCC optimize("Ofast")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,sse4.1,sse4.2,popcnt,abm,mmx,avx,avx2,tune=native")
*/
#ifdef LOCAL
#define dbg(x) cout << #x << " : " << x << endl;
#else
#define dbg(x)
#endif
#define int long long
#define pb push_back
#define ppb pop_back()
#define mp make_pair
#define fi(a, b) for (int i = a; i < b; i++)
#define fj(a, b) for (int j = a; j < b; j++)
#define fk(a, b) for (int k = a; k < b; k++)
#define fi1(a, b) for (int i = a - 1; i >= b; i--)
#define fj1(a, b) for (int j = a - 1; j >= b; j--)
#define fk1(a, b) for (int k = a - 1; k >= b; k--)
#define fx(x, a) for (auto& x : a)
#define rep(i, a, b) for (int i = a; i < b; ++i)
#define rep1(i, a, b) for (int i = a - 1; i >= b; --i)
#define siz(x) (int)x.size()
#define lb lower_bound
#define ub upper_bound
#define all(x) x.begin(), x.end()
using namespace std;
template<typename T1, typename T2>inline void mine(T1 &x, const T2 &y) { if (y < x) x = y; }
template<typename T1, typename T2>inline void maxe(T1 &x, const T2 &y) { if (x < y) x = y; }
ostream& operator << (ostream &out, const vector<int> &b) {
for (auto k : b) out << k << ' ';
return out;
}
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef char ch;
typedef string str;
typedef pair<int, int> pii;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<pii> vpii;
typedef vector<vpii> vvpii;
typedef vector<ch> vch;
typedef vector<vch> vvch;
typedef vector<str> vs;
const int MOD = 1000000007;
const int INF = 1000000050;
const long long BIG = (long long)2e18 + 50;
const int MX = 200010;
const double EPS = 1e-9;
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
int n, m;
int c[MX];
int from[MX];
int to[MX];
vi bandij(int s, int ban, const vvi &g) {
vi dist(n, BIG);
vch used(n, 0);
dist[s] = 0;
for (int i = 0; i < n - 1; ++i) {
int v = -1;
for (int j = 0; j < n; ++j) {
if (!used[j] && (v == -1 || dist[j] < dist[v])) {
v = j;
}
}
if (dist[v] == BIG) break;
fx(id, g[v]) {
if (id == ban) continue;
mine(dist[to[id]], dist[v] + c[id]);
}
used[v] = 1;
}
return dist;
}
pair<vvi, vi> dij(int s, const vvi & g) {
vi dist(n, BIG);
vch used(n, 0);
vvpii gb(n);
dist[s] = 0;
for (int i = 0; i < n - 1; ++i) {
int v = -1;
for (int j = 0; j < n; ++j) {
if (!used[j] && (v == -1 || dist[j] < dist[v])) {
v = j;
}
}
if (dist[v] == BIG) break;
fx(id, g[v]) {
if (dist[to[id]] > dist[v] + c[id]) {
gb[to[id]].clear();
dist[to[id]] = dist[v] + c[id];
}
if (dist[to[id]] == dist[v] + c[id])
gb[to[id]].pb({from[id], id});
}
used[v] = 1;
}
vvpii gb2(n);
for (int i = 0; i < n; ++i) {
used[i] = 0;
for (auto[x, id] : gb[i]) {
gb2[x].pb({i, id});
}
}
swap(gb2, gb);
vi ord;
function<void(int)> dfs = [&](int v) {
used[v] = 1;
for (auto[x, id] : gb[v]) {
if (!used[x]) {
dfs(x);
}
}
ord.pb(v);
};
dfs(s);
reverse(all(ord));
vector<pair<ull, int>> ways(n, {0, 0});
ways[s] = {1, 1};
vch was(n, 0);
fx(v, ord) {
was[v] = 1;
ways[v].second %= MOD;
for (auto[x, id] : gb[v]) {
if (was[x]) continue;
ways[x].first += ways[v].first;
ways[x].second += ways[v].second;
}
}
vi bridges;
fx(v, ord) {
for (auto[to, id] : gb[v]) {
if (ways[to] == ways[v]) {
bridges.pb(id);
}
}
}
sort(all(bridges));
vvi dists;
fx(id, bridges)
dists.pb(bandij(s, id, g));
dists.pb(dist);
fx(x, bridges)
x %= m;
return mp(dists, bridges);
}
int d[MX];
int get(vi &bridges, int id) {
auto it = lb(all(bridges), id);
if (it == bridges.end() || *it != id)
return bridges.size();
return it - bridges.begin();
}
int32_t main() {
ios_base::sync_with_stdio(0);
cin.tie(0);
cout.tie(0);
cin >> n >> m;
vvi g(n);
vvi gc(n);
for (int i = 0; i < m; ++i) {
cin >> from[i] >> to[i] >> c[i] >> d[i];
--from[i], --to[i];
from[i + m] = to[i];
to[i + m] = from[i];
c[i + m] = c[i];
g[from[i]].pb(i);
gc[from[i + m]].pb(i + m);
}
const int s = 0;
const int t = n - 1;
auto[dsv, bsv] = dij(s, g);
auto[dvs, bvs] = dij(s, gc);
auto[dtv, btv] = dij(t, g);
auto[dvt, bvt] = dij(t, gc);
/*dbg(bsv);
dbg(bvs);
dbg(btv);
dbg(bvt);*/
int ans = dsv.back()[t] + dtv.back()[s];
for (int id = 0; id < m; ++id) {
int u = from[id];
int v = to[id];
int bid = get(bsv, id);
int st = dsv[bid][t];
int sv = dsv[bid][v];
bid = get(bvs, id);
int us = dvs[bid][u];
bid = get(btv, id);
int ts = dtv[bid][s];
int tv = dtv[bid][v];
bid = get(bvt, id);
int ut = dvt[bid][u];
int cur = d[id];
cur += min(st, sv + c[id] + ut);
cur += min(ts, tv + c[id] + us);
/*dbg(id);
dbg(cur);*/
mine(ans, cur);
}
cout << (ans >= BIG ? -1 : ans) << '\n';
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
26 ms |
1132 KB |
Output is correct |
2 |
Correct |
4 ms |
620 KB |
Output is correct |
3 |
Correct |
68 ms |
2156 KB |
Output is correct |
4 |
Correct |
70 ms |
2156 KB |
Output is correct |
5 |
Correct |
2 ms |
492 KB |
Output is correct |
6 |
Correct |
8 ms |
876 KB |
Output is correct |
7 |
Correct |
1 ms |
364 KB |
Output is correct |
8 |
Correct |
1 ms |
364 KB |
Output is correct |
9 |
Correct |
1 ms |
492 KB |
Output is correct |
10 |
Correct |
73 ms |
2156 KB |
Output is correct |
11 |
Correct |
73 ms |
2156 KB |
Output is correct |
12 |
Correct |
71 ms |
2156 KB |
Output is correct |
13 |
Correct |
21 ms |
1388 KB |
Output is correct |
14 |
Correct |
45 ms |
1772 KB |
Output is correct |
15 |
Correct |
45 ms |
1644 KB |
Output is correct |
16 |
Correct |
46 ms |
2028 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
24 ms |
4716 KB |
Output is correct |
2 |
Correct |
25 ms |
4544 KB |
Output is correct |
3 |
Correct |
26 ms |
4716 KB |
Output is correct |
4 |
Correct |
2 ms |
620 KB |
Output is correct |
5 |
Correct |
2 ms |
492 KB |
Output is correct |
6 |
Correct |
1 ms |
492 KB |
Output is correct |
7 |
Correct |
1 ms |
492 KB |
Output is correct |
8 |
Correct |
1 ms |
364 KB |
Output is correct |
9 |
Correct |
27 ms |
4588 KB |
Output is correct |
10 |
Correct |
25 ms |
4460 KB |
Output is correct |
11 |
Correct |
26 ms |
4588 KB |
Output is correct |
12 |
Correct |
25 ms |
4588 KB |
Output is correct |
13 |
Correct |
26 ms |
4588 KB |
Output is correct |
14 |
Correct |
25 ms |
4716 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
23 ms |
1004 KB |
Output is correct |
2 |
Correct |
13 ms |
1024 KB |
Output is correct |
3 |
Correct |
33 ms |
5440 KB |
Output is correct |
4 |
Correct |
8 ms |
876 KB |
Output is correct |
5 |
Correct |
42 ms |
6892 KB |
Output is correct |
6 |
Correct |
1 ms |
364 KB |
Output is correct |
7 |
Correct |
1 ms |
364 KB |
Output is correct |
8 |
Correct |
29 ms |
5472 KB |
Output is correct |
9 |
Correct |
30 ms |
5524 KB |
Output is correct |
10 |
Correct |
35 ms |
6824 KB |
Output is correct |
11 |
Correct |
37 ms |
7736 KB |
Output is correct |
12 |
Correct |
42 ms |
7856 KB |
Output is correct |
13 |
Correct |
1 ms |
364 KB |
Output is correct |
14 |
Correct |
1 ms |
364 KB |
Output is correct |
15 |
Correct |
1 ms |
364 KB |
Output is correct |
16 |
Correct |
1 ms |
364 KB |
Output is correct |
17 |
Correct |
1 ms |
364 KB |
Output is correct |
18 |
Correct |
1 ms |
364 KB |
Output is correct |
19 |
Correct |
38 ms |
8048 KB |
Output is correct |
20 |
Correct |
36 ms |
7744 KB |
Output is correct |
21 |
Correct |
39 ms |
7744 KB |
Output is correct |
22 |
Correct |
36 ms |
7676 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
26 ms |
1132 KB |
Output is correct |
2 |
Correct |
4 ms |
620 KB |
Output is correct |
3 |
Correct |
68 ms |
2156 KB |
Output is correct |
4 |
Correct |
70 ms |
2156 KB |
Output is correct |
5 |
Correct |
2 ms |
492 KB |
Output is correct |
6 |
Correct |
8 ms |
876 KB |
Output is correct |
7 |
Correct |
1 ms |
364 KB |
Output is correct |
8 |
Correct |
1 ms |
364 KB |
Output is correct |
9 |
Correct |
1 ms |
492 KB |
Output is correct |
10 |
Correct |
73 ms |
2156 KB |
Output is correct |
11 |
Correct |
73 ms |
2156 KB |
Output is correct |
12 |
Correct |
71 ms |
2156 KB |
Output is correct |
13 |
Correct |
21 ms |
1388 KB |
Output is correct |
14 |
Correct |
45 ms |
1772 KB |
Output is correct |
15 |
Correct |
45 ms |
1644 KB |
Output is correct |
16 |
Correct |
46 ms |
2028 KB |
Output is correct |
17 |
Correct |
24 ms |
4716 KB |
Output is correct |
18 |
Correct |
25 ms |
4544 KB |
Output is correct |
19 |
Correct |
26 ms |
4716 KB |
Output is correct |
20 |
Correct |
2 ms |
620 KB |
Output is correct |
21 |
Correct |
2 ms |
492 KB |
Output is correct |
22 |
Correct |
1 ms |
492 KB |
Output is correct |
23 |
Correct |
1 ms |
492 KB |
Output is correct |
24 |
Correct |
1 ms |
364 KB |
Output is correct |
25 |
Correct |
27 ms |
4588 KB |
Output is correct |
26 |
Correct |
25 ms |
4460 KB |
Output is correct |
27 |
Correct |
26 ms |
4588 KB |
Output is correct |
28 |
Correct |
25 ms |
4588 KB |
Output is correct |
29 |
Correct |
26 ms |
4588 KB |
Output is correct |
30 |
Correct |
25 ms |
4716 KB |
Output is correct |
31 |
Correct |
23 ms |
1004 KB |
Output is correct |
32 |
Correct |
13 ms |
1024 KB |
Output is correct |
33 |
Correct |
33 ms |
5440 KB |
Output is correct |
34 |
Correct |
8 ms |
876 KB |
Output is correct |
35 |
Correct |
42 ms |
6892 KB |
Output is correct |
36 |
Correct |
1 ms |
364 KB |
Output is correct |
37 |
Correct |
1 ms |
364 KB |
Output is correct |
38 |
Correct |
29 ms |
5472 KB |
Output is correct |
39 |
Correct |
30 ms |
5524 KB |
Output is correct |
40 |
Correct |
35 ms |
6824 KB |
Output is correct |
41 |
Correct |
37 ms |
7736 KB |
Output is correct |
42 |
Correct |
42 ms |
7856 KB |
Output is correct |
43 |
Correct |
1 ms |
364 KB |
Output is correct |
44 |
Correct |
1 ms |
364 KB |
Output is correct |
45 |
Correct |
1 ms |
364 KB |
Output is correct |
46 |
Correct |
1 ms |
364 KB |
Output is correct |
47 |
Correct |
1 ms |
364 KB |
Output is correct |
48 |
Correct |
1 ms |
364 KB |
Output is correct |
49 |
Correct |
38 ms |
8048 KB |
Output is correct |
50 |
Correct |
36 ms |
7744 KB |
Output is correct |
51 |
Correct |
39 ms |
7744 KB |
Output is correct |
52 |
Correct |
36 ms |
7676 KB |
Output is correct |
53 |
Correct |
235 ms |
7032 KB |
Output is correct |
54 |
Correct |
234 ms |
7020 KB |
Output is correct |
55 |
Correct |
237 ms |
7020 KB |
Output is correct |
56 |
Correct |
70 ms |
2156 KB |
Output is correct |
57 |
Correct |
70 ms |
2176 KB |
Output is correct |
58 |
Correct |
200 ms |
5940 KB |
Output is correct |
59 |
Correct |
200 ms |
5996 KB |
Output is correct |
60 |
Correct |
199 ms |
6124 KB |
Output is correct |
61 |
Correct |
197 ms |
5996 KB |
Output is correct |
62 |
Correct |
198 ms |
5996 KB |
Output is correct |
63 |
Correct |
199 ms |
5996 KB |
Output is correct |
64 |
Correct |
332 ms |
6252 KB |
Output is correct |
65 |
Correct |
325 ms |
6124 KB |
Output is correct |
66 |
Correct |
287 ms |
6124 KB |
Output is correct |
67 |
Correct |
18 ms |
4192 KB |
Output is correct |
68 |
Correct |
73 ms |
6508 KB |
Output is correct |
69 |
Correct |
72 ms |
6380 KB |
Output is correct |
70 |
Correct |
160 ms |
6892 KB |
Output is correct |
71 |
Correct |
159 ms |
6892 KB |
Output is correct |
72 |
Correct |
162 ms |
6636 KB |
Output is correct |
73 |
Correct |
160 ms |
6764 KB |
Output is correct |
74 |
Correct |
162 ms |
7020 KB |
Output is correct |