QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#386694 | #7177. Many Many Cycles | iee | WA | 1ms | 3816kb | C++14 | 2.2kb | 2024-04-11 19:25:35 | 2024-04-11 19:25:35 |
Judging History
answer
#include <bits/stdc++.h>
#define int long long
using namespace std;
constexpr int N = 5005, Log = 13;
int n, m;
vector<array<int, 2>> e[N];
int dis[N];
vector<array<int, 3>> q;
int col[N], blk, dfc;
int dfn[N], st[N][Log];
bool vis[N];
void dfs(int u, int fa) {
dfn[u] = ++dfc, st[dfc][0] = fa;
col[u] = blk;
vis[u] = 1;
for (auto [v, w] : e[u]) {
if (vis[v]) {
if (v != fa) q.push_back({min(u, v), max(u, v), w});
continue;
}
dis[v] = dis[u] + w;
dfs(v, u);
}
}
int amin(int u, int v) {
return dfn[u] < dfn[v] ? u : v;
}
int LCA(int u, int v) {
if (u == v) return u;
if (dfn[u] > dfn[v]) swap(u, v);
int k = __lg(dfn[v] - dfn[u]);
return amin(st[dfn[u] + 1][k], st[dfn[v] - (1 << k) + 1][k]);
}
int dist(int u, int v) {
return dis[u] + dis[v] - dis[LCA(u, v)] * 2;
}
int cross(int u1, int v1, int u2, int v2) {
if (dist(u1, LCA(u2, v2)) + dist(v1, LCA(u2, v2)) != dist(u1, v1)) {
swap(u1, u2), swap(v1, v2);
}
if (dist(u1, LCA(u2, v2)) + dist(v1, LCA(u2, v2)) != dist(u1, v1)) {
return 0;
}
int l = LCA(u2, v2), res = 0;
if (dist(u1, l) + dist(l, LCA(u1, v1)) == dist(u1, LCA(u1, v1))) {
res += dis[LCA(u2, u1) ^ LCA(v2, u1) ^ l] - dis[LCA(u1, v1)];
}
if (dist(v1, l) + dist(l, LCA(u1, v1)) == dist(v1, LCA(u1, v1))) {
res += dis[LCA(u2, v1) ^ LCA(v2, v1) ^ l] - dis[LCA(u1, v1)];
}
return res;
}
signed main() {
cin.tie(0)->sync_with_stdio(0);
cin >> n >> m;
for (int i = 1; i <= m; i++) {
int u, v, w;
cin >> u >> v >> w;
e[u].push_back({v, w});
e[v].push_back({u, w});
}
for (int i = 1; i <= n; i++) {
if (!vis[i]) {
++blk;
dfs(i, 0);
}
}
for (int j = 1; j < Log; j++) {
for (int i = 1; i + (1 << j) - 1 <= n; i++) {
st[i][j] = amin(st[i][j - 1], st[i + (1 << (j - 1))][j - 1]);
}
}
sort(q.begin(), q.end()), q.erase(unique(q.begin(), q.end()), q.end());
int d = 0;
for (auto [u, v, w] : q) {
d = __gcd(d, dist(u, v) + w);
}
for (int i = 0; i < q.size(); i++) {
for (int j = 0; j < i; j++) {
auto [u1, v1, w1] = q[i];
auto [u2, v2, w2] = q[j];
if (col[u1] == col[u2]) {
d = __gcd(d, dist(u1, v1) + w1 + dist(u2, v2) + w2 - cross(u1, v1, u2, v2) * 2);
}
}
}
cout << d << "\n";
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 1ms
memory: 3740kb
input:
4 4 1 2 1 2 3 1 3 4 1 4 1 1
output:
4
result:
ok answer is '4'
Test #2:
score: 0
Accepted
time: 1ms
memory: 3816kb
input:
4 5 1 2 1 1 3 2 1 4 1 2 3 1 3 4 1
output:
4
result:
ok answer is '4'
Test #3:
score: 0
Accepted
time: 0ms
memory: 3804kb
input:
20 50 1 2 8 1 3 1 3 4 5 3 5 9 3 6 5 6 7 6 7 8 8 2 9 2 8 10 3 8 11 7 8 12 5 3 13 4 7 14 3 6 15 7 9 16 6 8 17 7 16 18 9 16 19 3 18 20 10 11 3 2 17 1 1 16 2 2 15 1 1 10 3 2 9 1 2 19 2 1 6 1 2 7 3 1 17 3 2 15 3 2 8 6 2 5 1 2 8 1 2 12 1 1 12 7 1 4 1 2 18 2 1 11 7 1 14 1 1 18 1 1 18 9 1 10 6 1 14 3 2 20 2...
output:
2
result:
ok answer is '2'
Test #4:
score: 0
Accepted
time: 1ms
memory: 3808kb
input:
20 50 1 2 18468 1 3 26501 3 4 15725 3 5 29359 3 6 24465 6 7 28146 7 8 16828 2 9 492 8 10 11943 8 11 5437 8 12 14605 3 13 154 7 14 12383 6 15 18717 9 16 19896 8 17 21727 16 18 11539 16 19 19913 18 20 26300 11 3 2 17 1 1 16 2 2 15 1 1 10 3 2 9 1 2 19 2 1 6 1 2 7 3 1 17 3 2 15 3 2 8 6 2 5 1 2 8 1 2 12 ...
output:
1
result:
ok answer is '1'
Test #5:
score: -100
Wrong Answer
time: 1ms
memory: 3756kb
input:
100 150 1 2 184676335 1 3 191705725 1 4 293606963 1 5 57078146 2 6 168279962 6 7 29961943 5 8 54392392 5 9 39020154 5 10 123837422 7 11 197199896 3 12 217274772 7 13 18709913 6 14 263007036 11 15 287053812 3 16 303347674 9 17 151417712 17 18 68705548 15 19 326652758 12 20 128598724 2 21 275290779 11...
output:
1
result:
wrong answer expected '3', found '1'