Submission #859909

# Submission time Handle Problem Language Result Execution time Memory
859909 2023-10-11T06:49:56 Z c2zi6 Duathlon (APIO18_duathlon) C++14
0 / 100
92 ms 27332 KB
#include <bits/stdc++.h>
#define ff first
#define ss second
#define pb push_back
#define all(a) (a).begin(), (a).end()
#define rep(i, a, b) for (int i = int(a); i <= int(b); ++i)
#define repn(i, n) for (int i = 0; i < int(n); ++i)
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int, int> PII;
typedef vector<int> VI;
typedef vector<PII> VPI;
typedef vector<VI> VVI;
typedef vector<VPI> VVPI;
template<class T> T setmax(T& a, T& b) {if (a < b) return a = b; return a;}
template<class T> T setmin(T& a, T& b) {if (a < b) return a; return a = b;}

const int MAXN = 2e5;   
int n, m;
VVI gp, gp2, comp;
VI stk;
int low[MAXN], num[MAXN], vis[MAXN], ap[MAXN], id[MAXN], cnt[MAXN];
int tc;
int itsz[MAXN];
int N;

void tarjan(int u, int p = -1) {
    vis[u] = 1;
    stk.pb(u);
    num[u] = low[u] = tc++;
    for (int v : gp[u]) {
        if (v == p) continue;
        if (vis[v] == 1) {
            setmin(low[u], num[v]);
        } else if (vis[v] == 0) {
            tarjan(v, u);
            setmin(low[u], low[v]);
            if (low[v] >= num[u]) {
                ap[u] = (num[u] > 0 || num[v] > 1);
                comp.pb(VI{u});
                while (comp.back().back() != v) {
                    comp.back().pb(stk.back());
                    stk.pop_back();
                }
            }
        }
    }
    vis[u] = 2;
}

ll ans = 0;

int dfs(int u, int p = -1) {
    vis[u] = true;
    int ret = 0;
    // cout << "SKSVAV " << u << endl;
    for (int v : gp2[u]) {
        if (v == p) continue;
        int r = dfs(v, u); 
        ret += r - 1;
    }
    // cout << "PRCAV " << u << ", ret = " << ret << endl;
    if (cnt[u] == -1) {
        int retrn = ret + 1;
        // cout << u << " " << retrn << endl;
    // cout << "MYUS@ " << u << " " << retrn << endl;
        return retrn;//ret+1;
    }
    // cout << u+1 << " " << cnt[u] << " " << N << endl;
    ans -= 1ll*cnt[u]*(cnt[u]-1)*(N-cnt[u]);
    // cout << u << " " << cnt[u] << " " << ret << endl;
    ans -= 2ll*cnt[u]*ret*(N-cnt[u]-ret);


    int retrn = cnt[u] + ret;
    // cout << "MYUS@ " << u << " " << retrn << endl;
    return retrn;
}

void solve() {
    cin >> n >> m;
    gp = VVI(n);
    repn(i, m) {
        int u, v;
        cin >> u >> v;
        u--, v--;
        gp[u].pb(v);
        gp[v].pb(u);
    }
    int iter = 0;
    repn(u, n) {
        if (!vis[u]) {
            tc = 0;
            tarjan(u);
            itsz[iter++] = tc;
        }
    }
    int last = 0;
    gp2 = VVI(2*n);
    repn(u, n) if (ap[u]) cnt[last] = -1, id[u] = last++;
    for (VI v : comp) {
        for (int x : v) {
            if (ap[x]) {
                gp2[id[x]].pb(last);
                gp2[last].pb(id[x]);
            } else id[x] = last;
        }
        cnt[last] = v.size();
        last++;
    }
    while (gp2.size() > last) gp2.pop_back();
    n = last;


    fill(vis, vis+n, 0);
    iter = 0;
    repn(u, n) {
        if (!vis[u]) {
            N = itsz[iter];
            ans += N*(N-1)*(N-2);
            // cout << ans << endl;
            dfs(u);
            iter++;
        }
    }
    cout << ans << endl;
}

int main() {
    // ios_base::sync_with_stdio(false);
    // cin.tie(null); cout.tie(null);
    ll t = 1;
    while (t--) solve();
}

Compilation message

count_triplets.cpp: In function 'void solve()':
count_triplets.cpp:112:23: warning: comparison of integer expressions of different signedness: 'std::vector<std::vector<int> >::size_type' {aka 'long unsigned int'} and 'int' [-Wsign-compare]
  112 |     while (gp2.size() > last) gp2.pop_back();
      |            ~~~~~~~~~~~^~~~~~
# Verdict Execution time Memory Grader output
1 Correct 1 ms 2392 KB Output is correct
2 Correct 1 ms 2392 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2648 KB Output is correct
5 Correct 1 ms 2396 KB Output is correct
6 Correct 1 ms 2488 KB Output is correct
7 Incorrect 1 ms 2396 KB Output isn't correct
8 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 2392 KB Output is correct
2 Correct 1 ms 2392 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2648 KB Output is correct
5 Correct 1 ms 2396 KB Output is correct
6 Correct 1 ms 2488 KB Output is correct
7 Incorrect 1 ms 2396 KB Output isn't correct
8 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 81 ms 27332 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 1 ms 2648 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 88 ms 26372 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 1 ms 2652 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Incorrect 92 ms 26448 KB Output isn't correct
2 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 2392 KB Output is correct
2 Correct 1 ms 2392 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2648 KB Output is correct
5 Correct 1 ms 2396 KB Output is correct
6 Correct 1 ms 2488 KB Output is correct
7 Incorrect 1 ms 2396 KB Output isn't correct
8 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 2392 KB Output is correct
2 Correct 1 ms 2392 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2648 KB Output is correct
5 Correct 1 ms 2396 KB Output is correct
6 Correct 1 ms 2488 KB Output is correct
7 Incorrect 1 ms 2396 KB Output isn't correct
8 Halted 0 ms 0 KB -