| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 158211 | maruii | Duathlon (APIO18_duathlon) | C++14 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
int N, M, dfn[100005], dfnn, low[100005], cpn[100005], cpnn, sz[100005], T[100005];
vector<int> edge[100005];
vector<pii> bdge[100005];
ll ans;
void dfs(int x, int p) {
dfn[x] = low[x] = ++dfnn;
N++;
for (auto i : edge[x]) {
if (i == p) continue;
if (dfn[i]) {
if (dfn[i] < dfn[x]) low[x] = min(low[x], dfn[i]);
continue;
}
dfs(i, x);
low[x] = min(low[x], low[i]);
}
}
void bcc(int x, int c) {
cpn[x] = c;
sz[c]++;
for (auto i : edge[x]) {
if (cpn[i]) continue;
if (low[i] == dfn[i]) {
bdge[++cpnn].emplace_back(c, i);
bdge[c].emplace_back(cpnn, x);
bcc(i, cpnn);
}
else bcc(i, c);
}
}
void calc(int x, int p) {
ll n = sz[x], t = (N - n) * (N - n);
ans += n * (n - 1) * (n - 2) + 2 * (n - 1) * (n - 1) * (N - n);
for (auto i : bdge[x]) {
if (i.ff == p) continue;
calc(i.ff, x);
ans -= 2 * (n - 1) * sz[i.ff] * T[i.ss];
T[i.ss] += sz[i.ff];
sz[x] += sz[i.ff];
t -= 1ll * sz[i.ff] * sz[i.ff];
}
ans += n * (t - 1ll * (N - sz[x]) * (N - sz[x]));
}
int main() {
ios_base::sync_with_stdio(0), cin.tie(0);
int n; cin >> n >> M;
for (int i = 0; i < M; ++i) {
int x, y; cin >> x >> y;
edge[x].push_back(y);
edge[y].push_back(x);
}
for (int i = 1; i <= n; ++i) {
if (dfn[i]) continue;
N = 0;
dfs(i, i);
bcc(i, ++cpnn);
calc(cpnn, cpnn);
}
cout << ans;
return 0;
}