#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using lint = long long; using pint = pair; using plint = pair; struct fast_ios { fast_ios(){ cin.tie(nullptr), ios::sync_with_stdio(false), cout << fixed << setprecision(20); }; } fast_ios_; #define ALL(x) (x).begin(), (x).end() #define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i=i##_begin_;i--) #define REP(i, n) FOR(i,0,n) #define IREP(i, n) IFOR(i,0,n) template void ndarray(vector& vec, const V& val, int len) { vec.assign(len, val); } template void ndarray(vector& vec, const V& val, int len, Args... args) { vec.resize(len), for_each(begin(vec), end(vec), [&](T& v) { ndarray(v, val, args...); }); } template bool chmax(T &m, const T q) { return m < q ? (m = q, true) : false; } template bool chmin(T &m, const T q) { return m > q ? (m = q, true) : false; } int floor_lg(long long x) { return x <= 0 ? -1 : 63 - __builtin_clzll(x); } template pair operator+(const pair &l, const pair &r) { return make_pair(l.first + r.first, l.second + r.second); } template pair operator-(const pair &l, const pair &r) { return make_pair(l.first - r.first, l.second - r.second); } template vector sort_unique(vector vec) { sort(vec.begin(), vec.end()), vec.erase(unique(vec.begin(), vec.end()), vec.end()); return vec; } template int arglb(const std::vector &v, const T &x) { return std::distance(v.begin(), std::lower_bound(v.begin(), v.end(), x)); } template int argub(const std::vector &v, const T &x) { return std::distance(v.begin(), std::upper_bound(v.begin(), v.end(), x)); } template istream &operator>>(istream &is, vector &vec) { for (auto &v : vec) is >> v; return is; } template ostream &operator<<(ostream &os, const vector &vec) { os << '['; for (auto v : vec) os << v << ','; os << ']'; return os; } template ostream &operator<<(ostream &os, const array &arr) { os << '['; for (auto v : arr) os << v << ','; os << ']'; return os; } #if __cplusplus >= 201703L template istream &operator>>(istream &is, tuple &tpl) { std::apply([&is](auto &&... args) { ((is >> args), ...);}, tpl); return is; } template ostream &operator<<(ostream &os, const tuple &tpl) { os << '('; std::apply([&os](auto &&... args) { ((os << args << ','), ...);}, tpl); return os << ')'; } #endif template ostream &operator<<(ostream &os, const deque &vec) { os << "deq["; for (auto v : vec) os << v << ','; os << ']'; return os; } template ostream &operator<<(ostream &os, const set &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; } template ostream &operator<<(ostream &os, const unordered_set &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; } template ostream &operator<<(ostream &os, const multiset &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; } template ostream &operator<<(ostream &os, const unordered_multiset &vec) { os << '{'; for (auto v : vec) os << v << ','; os << '}'; return os; } template ostream &operator<<(ostream &os, const pair &pa) { os << '(' << pa.first << ',' << pa.second << ')'; return os; } template ostream &operator<<(ostream &os, const map &mp) { os << '{'; for (auto v : mp) os << v.first << "=>" << v.second << ','; os << '}'; return os; } template ostream &operator<<(ostream &os, const unordered_map &mp) { os << '{'; for (auto v : mp) os << v.first << "=>" << v.second << ','; os << '}'; return os; } #ifdef HITONANODE_LOCAL const string COLOR_RESET = "\033[0m", BRIGHT_GREEN = "\033[1;32m", BRIGHT_RED = "\033[1;31m", BRIGHT_CYAN = "\033[1;36m", NORMAL_CROSSED = "\033[0;9;37m", RED_BACKGROUND = "\033[1;41m", NORMAL_FAINT = "\033[0;2m"; #define dbg(x) cerr << BRIGHT_CYAN << #x << COLOR_RESET << " = " << (x) << NORMAL_FAINT << " (L" << __LINE__ << ") " << __FILE__ << COLOR_RESET << endl #define dbgif(cond, x) ((cond) ? cerr << BRIGHT_CYAN << #x << COLOR_RESET << " = " << (x) << NORMAL_FAINT << " (L" << __LINE__ << ") " << __FILE__ << COLOR_RESET << endl : cerr) #else #define dbg(x) (x) #define dbgif(cond, x) 0 #endif #include #include #include #include // This implementation is based on: https://gist.github.com/wata-orz/d3037bd0b919c76dd9ddc0379e1e3192 template struct SSSUP { int V; std::vector>> to; SSSUP(int n) : V(n), to(n) { static_assert(INF > 0, "INF must be positive"); } void add_edge(int u, int v, T len, G g) { assert(u >= 0 and u < V); assert(v >= 0 and v < V); assert(len >= 0); to[u].emplace_back(v, len, g); to[v].emplace_back(u, len, g); } int s; std::vector dist_sp, dist; std::vector parent, depth; std::vector psi; std::vector uf_ps; int find(int x) { if (uf_ps[x] == -1) { return x; } else { return uf_ps[x] = find(uf_ps[x]); } } void unite(int r, int c) { uf_ps[c] = r; } void solve(int s) { assert(s >= 0 and s < V); // Solve SSSP { dist_sp.assign(V, INF); depth.assign(V, -1), parent.assign(V, -1); psi.assign(V, e()); std::priority_queue, std::vector>, std::greater<>> que; dist_sp[s] = 0, depth[s] = 0; que.emplace(0, s); while (que.size()) { T d, l; int u, v; G g; std::tie(d, u) = que.top(); que.pop(); if (dist_sp[u] != d) continue; for (const auto &p : to[u]) { std::tie(v, l, g) = p; const auto d2 = d + l; if (dist_sp[v] > d2) { dist_sp[v] = d2, depth[v] = depth[u] + 1, parent[v] = u, psi[v] = op(psi[u], g); que.emplace(d2, v); } } } } // sssup() uf_ps.assign(V, -1); std::priority_queue, std::vector>, std::greater<>> que; for (int u = 0; u < V; u++) { for (int i = 0; i < int(to[u].size()); i++) { int v; T l; G g; std::tie(v, l, g) = to[u][i]; if (u < v and op(psi[u], g) != psi[v]) que.emplace(dist_sp[u] + dist_sp[v] + l, u, i); } } dist.assign(V, INF); while (que.size()) { T h; int u0, i; std::tie(h, u0, i) = que.top(); que.pop(); const int v0 = std::get<0>(to[u0][i]); int u = find(u0), v = find(v0); std::vector bs; while (u != v) { if (depth[u] > depth[v]) { bs.push_back(u), u = find(parent[u]); } else { bs.push_back(v), v = find(parent[v]); } } for (const int x : bs) { unite(u, x); dist[x] = h - dist_sp[x]; for (int i = 0; i < int(to[x].size()); i++) { int y; T l; G g; std::tie(y, l, g) = to[x][i]; if (op(psi[x], g) != psi[y]) { que.emplace(dist[x] + dist_sp[y] + l, x, i); } } } } for (int i = 0; i < V; i++) { if (psi[i] != e() and dist_sp[i] < dist[i]) dist[i] = dist_sp[i]; } } }; // https://gist.github.com/wata-orz/d3037bd0b919c76dd9ddc0379e1e3192 constexpr lint INF = 1LL << 60; struct SSSP { int s; vector dist; vector parent, depth; }; SSSP sssp(const vector> &g, int s) { vector dist(g.size(), INF); vector depth(g.size()), parent(depth.size()); priority_queue que; dist[s] = 0; depth[s] = 0; que.emplace(0, s); while (que.size()) { auto [d, u] = que.top(); que.pop(); d = -d; if (dist[u] != d) { continue; } for (auto [v, l] : g[u]) { lint d2 = d + l; if (chmin(dist[v], d2)) { depth[v] = depth[u] + 1; parent[v] = u; que.emplace(-d2, v); } } } return SSSP{s, dist, parent, depth}; } // UnionFind Tree (0-indexed), based on size of each disjoint set struct UnionFind { vector ps; UnionFind(int N = 0) : ps(N, -1) { } int find(int x) { if (ps[x] == -1) return x; else { int r = find(ps[x]); ps[x] = r; return r; } } void unite(int r, int c) { ps[c] = r; } }; vector sssup(const vector> &g, const SSSP &sp, const vector> is_consistent) { REP(i, g.size()) assert(g[i].size() == is_consistent[i].size()); vector dist(g.size(), INF); UnionFind uf(g.size()); using P = tuple; priority_queue

que; REP(u, g.size()) { REP(i, g[u].size()) { auto [v, l] = g[u][i]; if (u < v and !is_consistent[u][i]) { que.emplace(-(sp.dist[u] + sp.dist[v] + l), u, i); } } } while (que.size()) { auto [h, u0, i] = que.top(); que.pop(); h = -h; const int v0 = g[u0][i].first; int u = uf.find(u0), v = uf.find(v0); vector bs; while (u != v) { if (sp.depth[u] > sp.depth[v]) { bs.push_back(u); u = uf.find(sp.parent[u]); } else { bs.push_back(v); v = uf.find(sp.parent[v]); } } for (auto x : bs) { uf.unite(u, x); dist[x] = h - sp.dist[x]; REP(i, g[x].size()) { if (is_consistent[x][i]) { const auto [y, l] = g[x][i]; que.emplace(-(dist[x] + sp.dist[y] + l), x, i); } } } } return dist; } // Weighted UnionFind template struct WeightedUnionFind { std::vector par, sz; std::vector pot; WeightedUnionFind(int N = 0) : par(N), sz(N, 1), pot(N) { std::iota(par.begin(), par.end(), 0); } int find(int x) { if (par[x] != x) { int r = find(par[x]); pot[x] = pot[x] ^ pot[par[x]], par[x] = r; } return par[x]; } bool unite(int s, int t, T rel_diff) { // Relate s and t by t = s + rel_diff // Return false iff contradiction happens. rel_diff = rel_diff ^ weight(s) ^ weight(t); if ((s = find(s)) == (t = find(t))) return rel_diff == 0; if (sz[s] < sz[t]) std::swap(s, t); par[t] = s, sz[s] += sz[t], pot[t] = rel_diff; return true; } T weight(int x) { find(x); return pot[x]; } T diff(int s, int t) { return weight(t) ^ weight(s); } int count(int x) { return sz[find(x)]; } bool same(int s, int t) { return find(s) == find(t); } }; unsigned op(unsigned x, unsigned y) { return x ^ y; } unsigned e() { return 0; } int main() { int N, M, K; cin >> N >> M >> K; // vector> to(N); // vector>> to2(N); constexpr lint INF = 1LL << 60; SSSUP graph(N); vector> edges; REP(e, M) { int a, b, c; string x; cin >> a >> b >> c >> x; unsigned m = 0; for (auto c : x) m = m * 2 + c - '0'; a--, b--; graph.add_edge(a, b, c, m); graph.add_edge(b, a, c, m); // to[a].emplace_back(b, c); // to[b].emplace_back(a, c); // to2[a].emplace_back(b, c, m); // to2[b].emplace_back(a, c, m); } graph.solve(N - 1); REP(i, N - 1) { cout << (graph.dist[i] == INF ? -1 : graph.dist[i]) << '\n'; } // auto sp = sssp(to, N - 1); // WeightedUnionFind wuf(N); // REP(i, N) { // for (auto [j, w, g] : to2[i]) { // if (sp.parent[j] == i and sp.dist[i] + w == sp.dist[j] and !wuf.same(i, j)) wuf.unite(i, j, g); // } // } // vector> is_consistent(N); // REP(i, N) { // for (auto [j, w, g] : to2[i]) { // auto diff = wuf.diff(i, j); // is_consistent[i].push_back(diff == g); // } // } // to2.clear(); // vector ret = sssup(to, sp, is_consistent); // REP(i, N - 1) { // lint x = ret[i]; // if (wuf.diff(N - 1, i)) chmin(x, sp.dist[i]); // cout << (x == INF ? -1 : x) << '\n'; // } }