#include<bits/stdc++.h>
#define pb emplace_back
using namespace std;
using ll = long long;
const int maxn = 300010;
int n, m;
vector<int> edge[maxn];
int low[maxn], st[maxn], stl, dep[maxn];
int cnt_b;
ll gsz;
vector<int> c[maxn];
bool is_ap[maxn];
ll res;
int h;
#define c2(i) ((ll)(i) * ((ll)(i)-1) / 2)
#define p3(i) ((ll)(i) * ((ll)(i)-1) * ((ll)(i)-2))
vector<int> tree_edge[maxn];
int sz[maxn];
void dfs2(int now, int last = 0){
sz[now] = c[now].size() + (now <= n);
for (int u : tree_edge[now]) if(u != last) {
dfs2(u, now);
sz[now] += sz[u] - 1;
}
if (now > n){
ll apc = 0;
//cout << "now : " ;
for (int u : c[now])
apc += is_ap[u];//, cout << u << ' ';
//cout << '\n';
if (!apc){
res += p3(c[now].size());
return ;
}
//cout << "apc " << apc << '\n';
auto add = [&](ll a_with_ap, ll b_with_ap, ll all) {
res += (a_with_ap - apc) * (b_with_ap - 1) * (a_with_ap - apc - 1) * 2;
//cout << a_with_ap - apc << ' ' << b_with_ap - 1 << ' ' << a_with_ap - apc -1 << '\n';
res += (apc - 1) * (b_with_ap - 1) * (a_with_ap - apc) * 2;
res += (b_with_ap - 1) * (all - a_with_ap - b_with_ap + 1) * (a_with_ap - apc);
};
for (int u : tree_edge[now]) if(u != last)
add(c[now].size(), sz[u], gsz);
if (last)
add(c[now].size(), gsz - sz[now]+1, gsz);
res += p3(c[now].size());
}
else {
auto add = [&](ll a, ll all){
res += (a-1) * (all-a);
};
for (int u : tree_edge[now]) if(u != last)
add(sz[u], gsz);
if (last)
add(gsz - sz[now]+1, gsz);
}
}
void cal(){
for (int i = h;i < cnt_b;++i){
for (int u : c[i]) if (is_ap[u]) {
tree_edge[u].pb(i);
tree_edge[i].pb(u);
}
}
dfs2(h);
gsz = 0, h = cnt_b;
}
void dfs(int now, int last = 0){
static int d;
++gsz;
low[now] = dep[now] = ++d;
st[stl++] = now;
int cnt_c = 0;
for (int u : edge[now]) if (u != last) {
if (!low[u]){
dfs(u, now);
++cnt_c;
low[now] = min(low[now], low[u]);
if (low[u] >= dep[now]){
do {
c[cnt_b].pb( st[--stl] );
} while (st[stl] != u);
c[cnt_b++].pb(now);
is_ap[now] = true;
}
}
else
low[now] = min(low[now], dep[u]);
}
if (last == 0 && cnt_c < 2)
is_ap[now] = false;
if (last == 0)
cal();
}
signed main(){
ios_base::sync_with_stdio(0), cin.tie(0);
cin >> n >> m;
for ( int a, b, i = 0 ; i < m ; ++i) {
cin >> a >> b;
edge[a].pb(b);
edge[b].pb(a);
}
h = cnt_b = n+1;
for (int i = 1;i <= n;++i) if (!low[i])
dfs(i);
cout << res << '\n';
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
19 ms |
21496 KB |
Output is correct |
2 |
Correct |
18 ms |
21496 KB |
Output is correct |
3 |
Correct |
18 ms |
21496 KB |
Output is correct |
4 |
Correct |
19 ms |
21496 KB |
Output is correct |
5 |
Correct |
19 ms |
21496 KB |
Output is correct |
6 |
Correct |
19 ms |
21496 KB |
Output is correct |
7 |
Incorrect |
19 ms |
21496 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
19 ms |
21496 KB |
Output is correct |
2 |
Correct |
18 ms |
21496 KB |
Output is correct |
3 |
Correct |
18 ms |
21496 KB |
Output is correct |
4 |
Correct |
19 ms |
21496 KB |
Output is correct |
5 |
Correct |
19 ms |
21496 KB |
Output is correct |
6 |
Correct |
19 ms |
21496 KB |
Output is correct |
7 |
Incorrect |
19 ms |
21496 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
98 ms |
34292 KB |
Output is correct |
2 |
Correct |
102 ms |
34292 KB |
Output is correct |
3 |
Correct |
126 ms |
40564 KB |
Output is correct |
4 |
Correct |
122 ms |
35316 KB |
Output is correct |
5 |
Correct |
114 ms |
35828 KB |
Output is correct |
6 |
Correct |
131 ms |
42100 KB |
Output is correct |
7 |
Correct |
131 ms |
35960 KB |
Output is correct |
8 |
Correct |
143 ms |
39032 KB |
Output is correct |
9 |
Correct |
138 ms |
33920 KB |
Output is correct |
10 |
Correct |
136 ms |
34272 KB |
Output is correct |
11 |
Correct |
101 ms |
29688 KB |
Output is correct |
12 |
Correct |
106 ms |
29560 KB |
Output is correct |
13 |
Correct |
101 ms |
29176 KB |
Output is correct |
14 |
Correct |
91 ms |
29048 KB |
Output is correct |
15 |
Correct |
69 ms |
27768 KB |
Output is correct |
16 |
Correct |
70 ms |
27768 KB |
Output is correct |
17 |
Correct |
21 ms |
22776 KB |
Output is correct |
18 |
Correct |
22 ms |
22776 KB |
Output is correct |
19 |
Correct |
21 ms |
22776 KB |
Output is correct |
20 |
Correct |
21 ms |
22780 KB |
Output is correct |
21 |
Correct |
21 ms |
22776 KB |
Output is correct |
22 |
Correct |
22 ms |
22776 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
19 ms |
21624 KB |
Output is correct |
2 |
Correct |
19 ms |
21624 KB |
Output is correct |
3 |
Correct |
19 ms |
21624 KB |
Output is correct |
4 |
Correct |
19 ms |
21880 KB |
Output is correct |
5 |
Correct |
19 ms |
21752 KB |
Output is correct |
6 |
Correct |
19 ms |
21756 KB |
Output is correct |
7 |
Correct |
19 ms |
21880 KB |
Output is correct |
8 |
Correct |
20 ms |
21628 KB |
Output is correct |
9 |
Correct |
20 ms |
21624 KB |
Output is correct |
10 |
Correct |
19 ms |
21624 KB |
Output is correct |
11 |
Correct |
20 ms |
21624 KB |
Output is correct |
12 |
Correct |
19 ms |
21624 KB |
Output is correct |
13 |
Correct |
19 ms |
21624 KB |
Output is correct |
14 |
Correct |
19 ms |
21624 KB |
Output is correct |
15 |
Correct |
19 ms |
21624 KB |
Output is correct |
16 |
Correct |
19 ms |
21624 KB |
Output is correct |
17 |
Correct |
19 ms |
21624 KB |
Output is correct |
18 |
Correct |
20 ms |
21624 KB |
Output is correct |
19 |
Correct |
19 ms |
21624 KB |
Output is correct |
20 |
Correct |
20 ms |
21624 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
132 ms |
34296 KB |
Output is correct |
2 |
Correct |
127 ms |
34296 KB |
Output is correct |
3 |
Correct |
146 ms |
34424 KB |
Output is correct |
4 |
Correct |
132 ms |
34296 KB |
Output is correct |
5 |
Correct |
138 ms |
34296 KB |
Output is correct |
6 |
Correct |
166 ms |
54648 KB |
Output is correct |
7 |
Correct |
155 ms |
41976 KB |
Output is correct |
8 |
Correct |
148 ms |
40056 KB |
Output is correct |
9 |
Correct |
152 ms |
38136 KB |
Output is correct |
10 |
Correct |
137 ms |
34300 KB |
Output is correct |
11 |
Correct |
142 ms |
34424 KB |
Output is correct |
12 |
Correct |
134 ms |
34428 KB |
Output is correct |
13 |
Correct |
136 ms |
34296 KB |
Output is correct |
14 |
Correct |
132 ms |
33528 KB |
Output is correct |
15 |
Correct |
123 ms |
32480 KB |
Output is correct |
16 |
Correct |
72 ms |
28920 KB |
Output is correct |
17 |
Correct |
78 ms |
32880 KB |
Output is correct |
18 |
Correct |
80 ms |
32880 KB |
Output is correct |
19 |
Correct |
82 ms |
33004 KB |
Output is correct |
20 |
Correct |
89 ms |
33016 KB |
Output is correct |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
19 ms |
21624 KB |
Output is correct |
2 |
Correct |
19 ms |
21624 KB |
Output is correct |
3 |
Incorrect |
20 ms |
21624 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
136 ms |
34296 KB |
Output is correct |
2 |
Correct |
143 ms |
34556 KB |
Output is correct |
3 |
Incorrect |
154 ms |
33316 KB |
Output isn't correct |
4 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
19 ms |
21496 KB |
Output is correct |
2 |
Correct |
18 ms |
21496 KB |
Output is correct |
3 |
Correct |
18 ms |
21496 KB |
Output is correct |
4 |
Correct |
19 ms |
21496 KB |
Output is correct |
5 |
Correct |
19 ms |
21496 KB |
Output is correct |
6 |
Correct |
19 ms |
21496 KB |
Output is correct |
7 |
Incorrect |
19 ms |
21496 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
19 ms |
21496 KB |
Output is correct |
2 |
Correct |
18 ms |
21496 KB |
Output is correct |
3 |
Correct |
18 ms |
21496 KB |
Output is correct |
4 |
Correct |
19 ms |
21496 KB |
Output is correct |
5 |
Correct |
19 ms |
21496 KB |
Output is correct |
6 |
Correct |
19 ms |
21496 KB |
Output is correct |
7 |
Incorrect |
19 ms |
21496 KB |
Output isn't correct |
8 |
Halted |
0 ms |
0 KB |
- |