#pragma GCC target("avx512f") #pragma GCC optimize("Ofast") #pragma GCC optimize("unroll-loops") // ループの部分展開 #include // #include // using namespace atcoder; #ifdef LOCAL #include #define dbg(...) debug::dbg(#__VA_ARGS__, __VA_ARGS__) #else #define dbg(...) void(0) #endif using namespace std; // #define int ll #define rep(i, a, b) for(auto i=static_cast(a), i##_end__=static_cast(b); i < i##_end__; i++) #define rep_r(i, a, b) for(auto i=static_cast(a), i##_end__=static_cast(b); i >= i##_end__; i--) #define fore(i, a) for(auto& i: a) #define forei(i, it, a) for(auto [i, it] = std::pair{0, std::begin(a)}; it != std::end(a); i++, it++) #define all(x) std::begin(x), std::end(x) using ll = long long; // __int128; using ull = unsigned long long; template using priority_queue_asc = priority_queue, greater>; template ostream& operator<<(ostream& os, const vector& v){ rep(i, 0, v.size()){ os << (i == 0 ? "" : " ") << v[i]; } return os; } template istream& operator>>(istream& is, pair& p) { is >> p.first >> p.second; return is; } template istream& operator>>(istream& is, vector& v) { fore(e, v){ is >> e; } return is; } template int siz(const C& c) { return static_cast(c.size()); } template constexpr int siz(const T (&)[N]) { return static_cast(N); } template constexpr void fil(A& a, const V &v){ if constexpr(is_array_v) for(auto& s: a) fil(s, v); else a = v; }; template constexpr auto vct(size_t n, const V& v){ return vector(n, v); } template constexpr auto vct(size_t n, A... a){ return vector(n, vct(a...)); } template constexpr bool chmax(T& a, S b) { if (a < b) { a = b; return 1; } return 0; } template constexpr bool chmin(T& a, S b) { if (a > b) { a = b; return 1; } return 0; } template constexpr T pow2(T x){ return x * x; } template constexpr T divceil(T x, S div) { return (x % div == 0) ? x / div : (x >= 0) ? (x / div) + 1 : -((-x) / div); } template constexpr T divfloor(T x, S div){ return (x % div == 0) ? x / div : (x >= 0) ? (x / div) : -((-x) / div) - 1; } template inline void outif(bool c, T t, S f){ if(c) cout << t << endl; else cout << f << endl; } inline void yesno(bool x){ cout << (x ? "Yes" : "No") << endl; } constexpr long long INF = 1LL << 60; // 1.15e18 constexpr int MOD = (int)1e9 + 7; using Type = ll; constexpr int MAX_N = 2048; int N; Type g[MAX_N][MAX_N]; // 最短路木 - 密グラフバージョン O(N^2) // 戻り値: d[v]=最短距離(重み), p[v]=親頂点(rootの場合、-1), l[v]=通る経路の色分け。最初の頂点(root直下の頂点番号) tuple, vector, vector> shortest_path_tree(int root, Type inf = numeric_limits::max()){ vector d(N, inf); vector p(N, -1), l(N, -1); vector done(N); d[root] = 0; while(true){ int v = -1; Type mn = inf; rep(i, 0, N) if(!done[i] && d[i] < mn) v = i, mn = d[i]; if(v == -1) break; done[v] = true; rep(i, 0, N){ bool f = done[i] || i == v || g[v][i] >= inf || d[i] <= mn + g[v][i]; d[i] = f ? d[i] : mn + g[v][i]; p[i] = f ? p[i] : v; l[i] = f ? l[i] : v == root ? i : l[v]; } } return {d, p, l}; } // 無向グラフ 最小サイクル - 密グラフバージョン O(N^2) pair> shortest_cycle_undirected(int root, Type inf = numeric_limits::max()){ auto [d, p, l] = shortest_path_tree(root, inf); Type mn = inf; int mu = -1, mv = -1; rep(u, 0, N){ if(u == root || d[u] >= inf) continue; rep(v, u + 1, N){ bool f = v == root || l[u] == l[v] || d[v] >= inf || g[u][v] >= inf || mn <= d[u] + d[v] + g[u][v]; mn = f ? mn : d[u] + d[v] + g[u][v]; mu = f ? mu : u; mv = f ? mv : v; } } if(mn == inf) return {inf, {}}; vector path; int v = mu; while(v != root) path.push_back(v), v = p[v]; path.push_back(root); reverse(all(path)); v = mv; while(v != root) path.push_back(v), v = p[v]; return {mn, path}; } pair> shortest_cycle_undirected_all(Type inf = numeric_limits::max()){ Type mn = inf; vector path; rep(i, 0, N){ auto [m, p] = shortest_cycle_undirected(i, inf); if(mn > m) mn = m, path = p; } return {mn, path}; } // 有向グラフ 最小サイクル - 密グラフバージョン O(N^2) pair> shortest_cycle_directed(int root, Type inf = numeric_limits::max()){ auto [d, p, l] = shortest_path_tree(root, inf); Type mn = inf; int v = -1; rep(i, 0, N){ bool f = i == root || d[i] >= inf || g[i][root] >= inf || mn <= d[i] + g[i][root]; mn = f ? mn : d[i] + g[i][root]; v = f ? v : i; } if(mn == inf) return {inf, {}}; vector path; while(v != root) path.push_back(v), v = p[v]; path.push_back(root); reverse(all(path)); return {mn, path}; } pair> shortest_cycle_directed_all(Type inf = numeric_limits::max()){ Type mn = inf; vector path; rep(i, 0, N){ auto [m, p] = shortest_cycle_directed(i, inf); if(mn > m) mn = m, path = p; } return {mn, path}; } void _main() { int T, M; cin >> T >> N >> M; constexpr auto inf = numeric_limits::max(); fil(g, inf); rep(i, 0, M){ int u, v, w; cin >> u >> v >> w; u--; v--; g[u][v] = w; if(T == 0) g[v][u] = w; } if(T == 0){ auto [mn, path] = shortest_cycle_undirected_all(inf); outif(mn < inf, mn, -1); } else{ auto [mn, path] = shortest_cycle_directed_all(inf); outif(mn < inf, mn, -1); } } signed main() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); cerr << fixed << setprecision(15); _main(); }