Submission #60190

# Submission time Handle Problem Language Result Execution time Memory
60190 2018-07-23T20:27:05 Z duality Duathlon (APIO18_duathlon) C++11
31 / 100
418 ms 42628 KB
#include <bits/stdc++.h>
using namespace std;
#define pb push_back
#define mp make_pair
typedef long long int LLI;
typedef pair<int,int> pii;
typedef vector<int> vi;
typedef vector<pii> vpii;

vi adjList[100000];
int disc[100000],low[100000],art[100000],num = 0;
vi com[100000];
vpii S;
int bcc = 0;
int doDFS(int u,int p) {
    int i,c = 0;
    disc[u] = low[u] = num++;
    for (i = 0; i < adjList[u].size(); i++) {
        int v = adjList[u][i];
        if (disc[v] == -1) {
            c++;
            S.pb(mp(u,v));
            doDFS(v,u);
            low[u] = min(low[u],low[v]);
            if (((p == -1) && (c > 1)) || ((p != -1) && (low[v] >= disc[u]))) {
                art[u] = 1;
                while (1) {
                    int uu = S.back().first,vv = S.back().second;
                    S.pop_back();
                    com[bcc].pb(uu),com[bcc].pb(vv);
                    if ((u == uu) && (v == vv)) break;
                }
                bcc++;
            }
        }
        else if ((v != p) && (disc[v] < disc[u])) low[u] = min(low[u],disc[v]),S.pb(mp(u,v));
    }
    return 0;
}
int size[100000];
vi adjList2[200000];
LLI ans = 0;
int visited[200000],sub[200000];
int doDFS2(int u,int p) {
    int i;
    sub[u] = (u < bcc) ? size[u]:0;
    visited[u] = 1;
    for (i = 0; i < adjList2[u].size(); i++) {
        int v = adjList2[u][i];
        if (v != p) sub[u] += doDFS2(v,u);
    }
    return sub[u];
}
int doDFS3(int u,int p,int r) {
    int i;
    if (u < bcc) ans += 2LL*size[u]*sub[u]*(sub[r]-sub[u]+size[u]);
    for (i = 0; i < adjList2[u].size(); i++) {
        int v = adjList2[u][i];
        if (v != p) doDFS3(v,u,r);
    }
    return 0;
}
int doDFS4(int u,int p,int s,int sum) {
    int i;
    if (u < bcc) {
        sum += com[u].size()-1;
        ans += (LLI) (com[s].size()-1)*(com[u].size()-1)*(sum);
    }
    for (i = 0; i < adjList2[u].size(); i++) {
        int v = adjList2[u][i];
        if (v != p) doDFS4(v,u,s,sum);
    }
    return 0;
}
int main() {
    int i;
    int n,m,u,v;
    scanf("%d %d",&n,&m);
    for (i = 0; i < m; i++) {
        scanf("%d %d",&u,&v);
        u--,v--;
        adjList[u].pb(v);
        adjList[v].pb(u);
    }

    int j;
    for (i = 0; i < n; i++) disc[i] = low[i] = -1;
    for (i = 0; i < n; i++) {
        if (disc[i] == -1) {
            doDFS(i,-1);
            if (!S.empty()) {
                while (!S.empty()) {
                    int u = S.back().first,v = S.back().second;
                    S.pop_back();
                    com[bcc].pb(u),com[bcc].pb(v);
                }
                bcc++;
            }
        }
    }
    for (i = 0; i < bcc; i++) {
        sort(com[i].begin(),com[i].end());
        com[i].resize(unique(com[i].begin(),com[i].end())-com[i].begin());
        for (j = 0; j < com[i].size(); j++) {
            if (art[com[i][j]]) {
                adjList2[i].pb(com[i][j]+bcc);
                adjList2[com[i][j]+bcc].pb(i);
            }
        }
        size[i] = com[i].size()-1;
        ans += (LLI) com[i].size()*(com[i].size()-1)*(com[i].size()-2);
        ans -= 2LL*size[i]*size[i]*size[i];
        ans += (LLI) size[i]*size[i];
    }
    for (i = 0; i < bcc; i++) {
        if (!visited[i]) {
            doDFS2(i,-1);
            doDFS3(i,-1,i);
            ans -= (LLI) sub[i]*sub[i];
        }
    }
    //for (i = 0; i < bcc; i++) doDFS4(i,-1,i,0);
    printf("%lld\n",ans);

    return 0;
}

Compilation message

count_triplets.cpp: In function 'int doDFS(int, int)':
count_triplets.cpp:18:19: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     for (i = 0; i < adjList[u].size(); i++) {
                 ~~^~~~~~~~~~~~~~~~~~~
count_triplets.cpp: In function 'int doDFS2(int, int)':
count_triplets.cpp:48:19: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     for (i = 0; i < adjList2[u].size(); i++) {
                 ~~^~~~~~~~~~~~~~~~~~~~
count_triplets.cpp: In function 'int doDFS3(int, int, int)':
count_triplets.cpp:57:19: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     for (i = 0; i < adjList2[u].size(); i++) {
                 ~~^~~~~~~~~~~~~~~~~~~~
count_triplets.cpp: In function 'int doDFS4(int, int, int, int)':
count_triplets.cpp:69:19: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
     for (i = 0; i < adjList2[u].size(); i++) {
                 ~~^~~~~~~~~~~~~~~~~~~~
count_triplets.cpp: In function 'int main()':
count_triplets.cpp:104:23: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
         for (j = 0; j < com[i].size(); j++) {
                     ~~^~~~~~~~~~~~~~~
count_triplets.cpp:78:10: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
     scanf("%d %d",&n,&m);
     ~~~~~^~~~~~~~~~~~~~~
count_triplets.cpp:80:14: warning: ignoring return value of 'int scanf(const char*, ...)', declared with attribute warn_unused_result [-Wunused-result]
         scanf("%d %d",&u,&v);
         ~~~~~^~~~~~~~~~~~~~~
# Verdict Execution time Memory Grader output
1 Correct 15 ms 9720 KB Output is correct
2 Correct 15 ms 9836 KB Output is correct
3 Correct 18 ms 9836 KB Output is correct
4 Correct 18 ms 9932 KB Output is correct
5 Correct 16 ms 10008 KB Output is correct
6 Correct 15 ms 10008 KB Output is correct
7 Incorrect 15 ms 10008 KB Output isn't correct
8 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 15 ms 9720 KB Output is correct
2 Correct 15 ms 9836 KB Output is correct
3 Correct 18 ms 9836 KB Output is correct
4 Correct 18 ms 9932 KB Output is correct
5 Correct 16 ms 10008 KB Output is correct
6 Correct 15 ms 10008 KB Output is correct
7 Incorrect 15 ms 10008 KB Output isn't correct
8 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 228 ms 28352 KB Output is correct
2 Correct 214 ms 29600 KB Output is correct
3 Correct 323 ms 30404 KB Output is correct
4 Correct 257 ms 30836 KB Output is correct
5 Correct 270 ms 30836 KB Output is correct
6 Correct 339 ms 30836 KB Output is correct
7 Correct 322 ms 30836 KB Output is correct
8 Correct 347 ms 30836 KB Output is correct
9 Correct 347 ms 30836 KB Output is correct
10 Correct 304 ms 30836 KB Output is correct
11 Correct 218 ms 30836 KB Output is correct
12 Correct 178 ms 30836 KB Output is correct
13 Correct 198 ms 30836 KB Output is correct
14 Correct 196 ms 30836 KB Output is correct
15 Correct 124 ms 30836 KB Output is correct
16 Correct 136 ms 30836 KB Output is correct
17 Correct 17 ms 30836 KB Output is correct
18 Correct 17 ms 30836 KB Output is correct
19 Correct 17 ms 30836 KB Output is correct
20 Correct 14 ms 30836 KB Output is correct
21 Correct 16 ms 30836 KB Output is correct
22 Correct 19 ms 30836 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 15 ms 30836 KB Output is correct
2 Correct 17 ms 30836 KB Output is correct
3 Correct 16 ms 30836 KB Output is correct
4 Correct 17 ms 30836 KB Output is correct
5 Correct 17 ms 30836 KB Output is correct
6 Correct 13 ms 30836 KB Output is correct
7 Correct 17 ms 30836 KB Output is correct
8 Correct 9 ms 30836 KB Output is correct
9 Correct 15 ms 30836 KB Output is correct
10 Correct 18 ms 30836 KB Output is correct
11 Correct 16 ms 30836 KB Output is correct
12 Correct 17 ms 30836 KB Output is correct
13 Correct 15 ms 30836 KB Output is correct
14 Correct 17 ms 30836 KB Output is correct
15 Correct 15 ms 30836 KB Output is correct
16 Correct 17 ms 30836 KB Output is correct
17 Correct 16 ms 30836 KB Output is correct
18 Correct 16 ms 30836 KB Output is correct
19 Correct 15 ms 30836 KB Output is correct
20 Correct 15 ms 30836 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 303 ms 30836 KB Output is correct
2 Correct 304 ms 30836 KB Output is correct
3 Correct 278 ms 31408 KB Output is correct
4 Correct 298 ms 31460 KB Output is correct
5 Correct 319 ms 31460 KB Output is correct
6 Correct 418 ms 42628 KB Output is correct
7 Correct 389 ms 42628 KB Output is correct
8 Correct 384 ms 42628 KB Output is correct
9 Correct 382 ms 42628 KB Output is correct
10 Correct 321 ms 42628 KB Output is correct
11 Correct 318 ms 42628 KB Output is correct
12 Correct 274 ms 42628 KB Output is correct
13 Correct 306 ms 42628 KB Output is correct
14 Correct 254 ms 42628 KB Output is correct
15 Correct 239 ms 42628 KB Output is correct
16 Correct 148 ms 42628 KB Output is correct
17 Correct 159 ms 42628 KB Output is correct
18 Correct 163 ms 42628 KB Output is correct
19 Correct 186 ms 42628 KB Output is correct
20 Correct 176 ms 42628 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 16 ms 42628 KB Output is correct
2 Correct 15 ms 42628 KB Output is correct
3 Incorrect 18 ms 42628 KB Output isn't correct
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 310 ms 42628 KB Output is correct
2 Correct 301 ms 42628 KB Output is correct
3 Incorrect 335 ms 42628 KB Output isn't correct
4 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 15 ms 9720 KB Output is correct
2 Correct 15 ms 9836 KB Output is correct
3 Correct 18 ms 9836 KB Output is correct
4 Correct 18 ms 9932 KB Output is correct
5 Correct 16 ms 10008 KB Output is correct
6 Correct 15 ms 10008 KB Output is correct
7 Incorrect 15 ms 10008 KB Output isn't correct
8 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 15 ms 9720 KB Output is correct
2 Correct 15 ms 9836 KB Output is correct
3 Correct 18 ms 9836 KB Output is correct
4 Correct 18 ms 9932 KB Output is correct
5 Correct 16 ms 10008 KB Output is correct
6 Correct 15 ms 10008 KB Output is correct
7 Incorrect 15 ms 10008 KB Output isn't correct
8 Halted 0 ms 0 KB -