Submission #1097552

# Submission time Handle Problem Language Result Execution time Memory
1097552 2024-10-07T14:49:58 Z alexander707070 Diversity (CEOI21_diversity) C++14
38 / 100
358 ms 61776 KB
#include<bits/stdc++.h>

#define MAXN 300007
using namespace std;

const int bucket_sz=500;

int n,a[MAXN],q;

int br[MAXN],maxs,l,r,p[MAXN],len,bucket[MAXN],pos[MAXN],nxt[MAXN][20],pr[MAXN][20],w[MAXN];
long long ans,mult,sum,s[MAXN],res[MAXN];

struct Fenwick{
    long long fenwick[MAXN];

    void update(int x,long long val){
        while(x<=maxs){
            fenwick[x]+=val;
            x+=(x & (-x));
        }
    }

    long long getsum(int x){
        long long res=0;
        
        while(x>=1){
            res+=fenwick[x];
            x^=(x & (-x));
        }

        return res;
    }
}pref,suff,sumpref,sumsuff;

struct query{
    int l,r,id;

    inline friend bool operator < (query fr,query sc){
        if(bucket[fr.l]!=bucket[sc.l])return bucket[fr.l]<bucket[sc.l];
        return fr.r<sc.r;
    }

}qs[MAXN];

long long calc(int x){
    return pref.getsum(x-1)-sumpref.getsum(x-1)*(maxs-x+1) + suff.getsum(maxs-x)-sumsuff.getsum(maxs-x)*x; 
}

void change(int x,long long val){
    if(val==0)return;

    long long mul=calc(x);

    ans-=s[x]*mul;
    ans-=s[x]*(s[x]-1)/2+s[x];

    s[x]+=val;

    sumpref.update(x,val);
    sumsuff.update(maxs-x+1,val);

    pref.update(x,(maxs-x+2)*val);
    suff.update(maxs-x+1,(x+1)*val);

    ans+=s[x]*mul;
    ans+=s[x]*(s[x]-1)/2+s[x];
}

void calcdp(){
    for(int j=1;j<20;j++){
        for(int i=1;i<=n;i++){
            nxt[i][j]=nxt[nxt[i][j-1]][j-1];
            pr[i][j]=pr[pr[i][j-1]][j-1];
        }
    }
}

void add(int x){
    br[x]++;

    int z=pos[x];
    for(int s=15;s>=0;s--){
        if(br[w[nxt[z][s]]]<br[x])z=nxt[z][s];
    }

    int y=w[z];

    if(x!=y){
        change(pos[x],br[y]-(br[x]-1));
        change(pos[y],br[x]-br[y]);

        swap(w[pos[x]],w[pos[y]]);
        swap(pos[x],pos[y]);
    }else change(pos[x],1);
}

void rem(int x){
    br[x]--;

    int z=pos[x];
    for(int s=15;s>=0;s--){
        if(pr[z][s]!=0 and br[w[pr[z][s]]]>br[x])z=pr[z][s];
    }

    int y=w[z];

    if(x!=y){
        change(pos[x],br[y]-(br[x]+1));
        change(pos[y],br[x]-br[y]);

        swap(w[pos[x]],w[pos[y]]);
        swap(pos[x],pos[y]);
    }else{
        change(pos[x],-1);
    }
}

int main(){

    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    cin>>n>>q;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        maxs=max(maxs,a[i]);

        bucket[i]=i/bucket_sz;
    }

    for(int i=1;i<=maxs;i++)pos[i]=w[i]=i;

    for(int i=1;i<=maxs;i++){
        if(i<=maxs/2)nxt[i][0]=maxs-i+1;
        else nxt[i][0]=maxs-i+2;
    }
    if(maxs%2==1)nxt[maxs/2+1][0]=maxs/2+1;

    for(int i=1;i<=maxs;i++)pr[nxt[i][0]][0]=i;

    if(maxs%2==1)pr[maxs/2+1][0]=maxs/2+2;
    else pr[maxs/2+1][0]=maxs/2;
    pr[1][0]=0;

    calcdp();

    for(int i=1;i<=q;i++){
        cin>>qs[i].l>>qs[i].r;
        qs[i].id=i;
    }

    sort(qs+1,qs+q+1);

    l=1; r=0;
    for(int i=1;i<=q;i++){
        while(r<qs[i].r){
            r++; add(a[r]);
        }

        while(l>qs[i].l){
            l--; add(a[l]);
        }

        while(r>qs[i].r){
            rem(a[r]); r--;
        }

        while(l<qs[i].l){
            rem(a[l]); l++;
        }

        res[qs[i].id]=ans;
    }

    for(int i=1;i<=q;i++){
        cout<<res[i]<<"\n";
    }
    
    return 0;
}
# Verdict Execution time Memory Grader output
1 Correct 1 ms 344 KB Output is correct
2 Correct 0 ms 344 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 1 ms 348 KB Output is correct
5 Correct 0 ms 348 KB Output is correct
6 Correct 0 ms 348 KB Output is correct
7 Correct 1 ms 344 KB Output is correct
8 Correct 1 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 1 ms 348 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 604 KB Output is correct
3 Correct 3 ms 2140 KB Output is correct
4 Correct 27 ms 16752 KB Output is correct
5 Correct 59 ms 33364 KB Output is correct
6 Correct 87 ms 49748 KB Output is correct
7 Correct 86 ms 49776 KB Output is correct
8 Correct 85 ms 49744 KB Output is correct
9 Correct 92 ms 49780 KB Output is correct
10 Correct 88 ms 49820 KB Output is correct
11 Correct 91 ms 49744 KB Output is correct
12 Correct 85 ms 49744 KB Output is correct
13 Correct 84 ms 49744 KB Output is correct
14 Correct 87 ms 49744 KB Output is correct
15 Correct 81 ms 49820 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 604 KB Output is correct
3 Correct 3 ms 2140 KB Output is correct
4 Correct 27 ms 16752 KB Output is correct
5 Correct 59 ms 33364 KB Output is correct
6 Correct 87 ms 49748 KB Output is correct
7 Correct 86 ms 49776 KB Output is correct
8 Correct 85 ms 49744 KB Output is correct
9 Correct 92 ms 49780 KB Output is correct
10 Correct 88 ms 49820 KB Output is correct
11 Correct 91 ms 49744 KB Output is correct
12 Correct 85 ms 49744 KB Output is correct
13 Correct 84 ms 49744 KB Output is correct
14 Correct 87 ms 49744 KB Output is correct
15 Correct 81 ms 49820 KB Output is correct
16 Correct 0 ms 344 KB Output is correct
17 Correct 1 ms 604 KB Output is correct
18 Correct 3 ms 2140 KB Output is correct
19 Correct 27 ms 16988 KB Output is correct
20 Correct 57 ms 33372 KB Output is correct
21 Correct 93 ms 49744 KB Output is correct
22 Correct 87 ms 49748 KB Output is correct
23 Correct 89 ms 49748 KB Output is correct
24 Correct 87 ms 49740 KB Output is correct
25 Correct 86 ms 49744 KB Output is correct
26 Correct 94 ms 49744 KB Output is correct
27 Correct 100 ms 49744 KB Output is correct
28 Correct 92 ms 49620 KB Output is correct
29 Correct 92 ms 49752 KB Output is correct
30 Correct 92 ms 49744 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 348 KB Output is correct
2 Correct 0 ms 604 KB Output is correct
3 Correct 3 ms 2140 KB Output is correct
4 Correct 27 ms 16752 KB Output is correct
5 Correct 59 ms 33364 KB Output is correct
6 Correct 87 ms 49748 KB Output is correct
7 Correct 86 ms 49776 KB Output is correct
8 Correct 85 ms 49744 KB Output is correct
9 Correct 92 ms 49780 KB Output is correct
10 Correct 88 ms 49820 KB Output is correct
11 Correct 91 ms 49744 KB Output is correct
12 Correct 85 ms 49744 KB Output is correct
13 Correct 84 ms 49744 KB Output is correct
14 Correct 87 ms 49744 KB Output is correct
15 Correct 81 ms 49820 KB Output is correct
16 Correct 0 ms 344 KB Output is correct
17 Correct 1 ms 604 KB Output is correct
18 Correct 3 ms 2140 KB Output is correct
19 Correct 27 ms 16988 KB Output is correct
20 Correct 57 ms 33372 KB Output is correct
21 Correct 93 ms 49744 KB Output is correct
22 Correct 87 ms 49748 KB Output is correct
23 Correct 89 ms 49748 KB Output is correct
24 Correct 87 ms 49740 KB Output is correct
25 Correct 86 ms 49744 KB Output is correct
26 Correct 94 ms 49744 KB Output is correct
27 Correct 100 ms 49744 KB Output is correct
28 Correct 92 ms 49620 KB Output is correct
29 Correct 92 ms 49752 KB Output is correct
30 Correct 92 ms 49744 KB Output is correct
31 Correct 1 ms 604 KB Output is correct
32 Correct 1 ms 604 KB Output is correct
33 Correct 1 ms 604 KB Output is correct
34 Correct 4 ms 2140 KB Output is correct
35 Correct 4 ms 2140 KB Output is correct
36 Correct 4 ms 2140 KB Output is correct
37 Correct 17 ms 8732 KB Output is correct
38 Correct 18 ms 8796 KB Output is correct
39 Correct 36 ms 17024 KB Output is correct
40 Correct 70 ms 33384 KB Output is correct
41 Correct 111 ms 49640 KB Output is correct
42 Correct 110 ms 49748 KB Output is correct
43 Correct 111 ms 49744 KB Output is correct
44 Correct 108 ms 49756 KB Output is correct
45 Correct 111 ms 49744 KB Output is correct
46 Correct 108 ms 49744 KB Output is correct
47 Correct 108 ms 49756 KB Output is correct
48 Correct 108 ms 49748 KB Output is correct
49 Correct 110 ms 49752 KB Output is correct
50 Correct 106 ms 49780 KB Output is correct
51 Correct 102 ms 49796 KB Output is correct
52 Correct 104 ms 49748 KB Output is correct
53 Correct 100 ms 49760 KB Output is correct
54 Correct 101 ms 49744 KB Output is correct
55 Correct 102 ms 49744 KB Output is correct
56 Correct 101 ms 49880 KB Output is correct
57 Correct 102 ms 49788 KB Output is correct
58 Correct 106 ms 49748 KB Output is correct
59 Correct 104 ms 49748 KB Output is correct
60 Correct 101 ms 49752 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 344 KB Output is correct
2 Correct 0 ms 344 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 1 ms 348 KB Output is correct
5 Correct 0 ms 348 KB Output is correct
6 Correct 0 ms 348 KB Output is correct
7 Correct 1 ms 344 KB Output is correct
8 Correct 1 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 1 ms 348 KB Output is correct
11 Correct 0 ms 348 KB Output is correct
12 Correct 0 ms 604 KB Output is correct
13 Correct 3 ms 2140 KB Output is correct
14 Correct 27 ms 16752 KB Output is correct
15 Correct 59 ms 33364 KB Output is correct
16 Correct 87 ms 49748 KB Output is correct
17 Correct 86 ms 49776 KB Output is correct
18 Correct 85 ms 49744 KB Output is correct
19 Correct 92 ms 49780 KB Output is correct
20 Correct 88 ms 49820 KB Output is correct
21 Correct 91 ms 49744 KB Output is correct
22 Correct 85 ms 49744 KB Output is correct
23 Correct 84 ms 49744 KB Output is correct
24 Correct 87 ms 49744 KB Output is correct
25 Correct 81 ms 49820 KB Output is correct
26 Correct 0 ms 344 KB Output is correct
27 Correct 1 ms 604 KB Output is correct
28 Correct 3 ms 2140 KB Output is correct
29 Correct 27 ms 16988 KB Output is correct
30 Correct 57 ms 33372 KB Output is correct
31 Correct 93 ms 49744 KB Output is correct
32 Correct 87 ms 49748 KB Output is correct
33 Correct 89 ms 49748 KB Output is correct
34 Correct 87 ms 49740 KB Output is correct
35 Correct 86 ms 49744 KB Output is correct
36 Correct 94 ms 49744 KB Output is correct
37 Correct 100 ms 49744 KB Output is correct
38 Correct 92 ms 49620 KB Output is correct
39 Correct 92 ms 49752 KB Output is correct
40 Correct 92 ms 49744 KB Output is correct
41 Correct 1 ms 604 KB Output is correct
42 Correct 1 ms 604 KB Output is correct
43 Correct 1 ms 604 KB Output is correct
44 Correct 4 ms 2140 KB Output is correct
45 Correct 4 ms 2140 KB Output is correct
46 Correct 4 ms 2140 KB Output is correct
47 Correct 17 ms 8732 KB Output is correct
48 Correct 18 ms 8796 KB Output is correct
49 Correct 36 ms 17024 KB Output is correct
50 Correct 70 ms 33384 KB Output is correct
51 Correct 111 ms 49640 KB Output is correct
52 Correct 110 ms 49748 KB Output is correct
53 Correct 111 ms 49744 KB Output is correct
54 Correct 108 ms 49756 KB Output is correct
55 Correct 111 ms 49744 KB Output is correct
56 Correct 108 ms 49744 KB Output is correct
57 Correct 108 ms 49756 KB Output is correct
58 Correct 108 ms 49748 KB Output is correct
59 Correct 110 ms 49752 KB Output is correct
60 Correct 106 ms 49780 KB Output is correct
61 Correct 102 ms 49796 KB Output is correct
62 Correct 104 ms 49748 KB Output is correct
63 Correct 100 ms 49760 KB Output is correct
64 Correct 101 ms 49744 KB Output is correct
65 Correct 102 ms 49744 KB Output is correct
66 Correct 101 ms 49880 KB Output is correct
67 Correct 102 ms 49788 KB Output is correct
68 Correct 106 ms 49748 KB Output is correct
69 Correct 104 ms 49748 KB Output is correct
70 Correct 101 ms 49752 KB Output is correct
71 Correct 18 ms 8796 KB Output is correct
72 Correct 18 ms 8792 KB Output is correct
73 Correct 21 ms 8844 KB Output is correct
74 Correct 22 ms 8740 KB Output is correct
75 Correct 19 ms 8792 KB Output is correct
76 Correct 43 ms 16988 KB Output is correct
77 Correct 44 ms 16976 KB Output is correct
78 Correct 41 ms 16988 KB Output is correct
79 Correct 42 ms 16988 KB Output is correct
80 Correct 41 ms 16988 KB Output is correct
81 Correct 94 ms 33876 KB Output is correct
82 Correct 100 ms 34832 KB Output is correct
83 Correct 97 ms 34672 KB Output is correct
84 Correct 96 ms 34640 KB Output is correct
85 Correct 94 ms 34648 KB Output is correct
86 Correct 112 ms 35276 KB Output is correct
87 Correct 108 ms 35412 KB Output is correct
88 Correct 111 ms 35412 KB Output is correct
89 Correct 111 ms 35408 KB Output is correct
90 Correct 114 ms 35408 KB Output is correct
91 Correct 252 ms 56400 KB Output is correct
92 Correct 253 ms 56596 KB Output is correct
93 Correct 253 ms 56404 KB Output is correct
94 Correct 253 ms 56488 KB Output is correct
95 Correct 236 ms 56400 KB Output is correct
96 Incorrect 358 ms 61776 KB Output isn't correct
97 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 1 ms 344 KB Output is correct
2 Correct 0 ms 344 KB Output is correct
3 Correct 0 ms 348 KB Output is correct
4 Correct 1 ms 348 KB Output is correct
5 Correct 0 ms 348 KB Output is correct
6 Correct 0 ms 348 KB Output is correct
7 Correct 1 ms 344 KB Output is correct
8 Correct 1 ms 348 KB Output is correct
9 Correct 0 ms 348 KB Output is correct
10 Correct 1 ms 348 KB Output is correct
11 Correct 0 ms 348 KB Output is correct
12 Correct 0 ms 604 KB Output is correct
13 Correct 3 ms 2140 KB Output is correct
14 Correct 27 ms 16752 KB Output is correct
15 Correct 59 ms 33364 KB Output is correct
16 Correct 87 ms 49748 KB Output is correct
17 Correct 86 ms 49776 KB Output is correct
18 Correct 85 ms 49744 KB Output is correct
19 Correct 92 ms 49780 KB Output is correct
20 Correct 88 ms 49820 KB Output is correct
21 Correct 91 ms 49744 KB Output is correct
22 Correct 85 ms 49744 KB Output is correct
23 Correct 84 ms 49744 KB Output is correct
24 Correct 87 ms 49744 KB Output is correct
25 Correct 81 ms 49820 KB Output is correct
26 Correct 0 ms 344 KB Output is correct
27 Correct 1 ms 604 KB Output is correct
28 Correct 3 ms 2140 KB Output is correct
29 Correct 27 ms 16988 KB Output is correct
30 Correct 57 ms 33372 KB Output is correct
31 Correct 93 ms 49744 KB Output is correct
32 Correct 87 ms 49748 KB Output is correct
33 Correct 89 ms 49748 KB Output is correct
34 Correct 87 ms 49740 KB Output is correct
35 Correct 86 ms 49744 KB Output is correct
36 Correct 94 ms 49744 KB Output is correct
37 Correct 100 ms 49744 KB Output is correct
38 Correct 92 ms 49620 KB Output is correct
39 Correct 92 ms 49752 KB Output is correct
40 Correct 92 ms 49744 KB Output is correct
41 Correct 1 ms 604 KB Output is correct
42 Correct 1 ms 604 KB Output is correct
43 Correct 1 ms 604 KB Output is correct
44 Correct 4 ms 2140 KB Output is correct
45 Correct 4 ms 2140 KB Output is correct
46 Correct 4 ms 2140 KB Output is correct
47 Correct 17 ms 8732 KB Output is correct
48 Correct 18 ms 8796 KB Output is correct
49 Correct 36 ms 17024 KB Output is correct
50 Correct 70 ms 33384 KB Output is correct
51 Correct 111 ms 49640 KB Output is correct
52 Correct 110 ms 49748 KB Output is correct
53 Correct 111 ms 49744 KB Output is correct
54 Correct 108 ms 49756 KB Output is correct
55 Correct 111 ms 49744 KB Output is correct
56 Correct 108 ms 49744 KB Output is correct
57 Correct 108 ms 49756 KB Output is correct
58 Correct 108 ms 49748 KB Output is correct
59 Correct 110 ms 49752 KB Output is correct
60 Correct 106 ms 49780 KB Output is correct
61 Correct 102 ms 49796 KB Output is correct
62 Correct 104 ms 49748 KB Output is correct
63 Correct 100 ms 49760 KB Output is correct
64 Correct 101 ms 49744 KB Output is correct
65 Correct 102 ms 49744 KB Output is correct
66 Correct 101 ms 49880 KB Output is correct
67 Correct 102 ms 49788 KB Output is correct
68 Correct 106 ms 49748 KB Output is correct
69 Correct 104 ms 49748 KB Output is correct
70 Correct 101 ms 49752 KB Output is correct
71 Correct 18 ms 8796 KB Output is correct
72 Correct 18 ms 8792 KB Output is correct
73 Correct 21 ms 8844 KB Output is correct
74 Correct 22 ms 8740 KB Output is correct
75 Correct 19 ms 8792 KB Output is correct
76 Correct 43 ms 16988 KB Output is correct
77 Correct 44 ms 16976 KB Output is correct
78 Correct 41 ms 16988 KB Output is correct
79 Correct 42 ms 16988 KB Output is correct
80 Correct 41 ms 16988 KB Output is correct
81 Correct 94 ms 33876 KB Output is correct
82 Correct 100 ms 34832 KB Output is correct
83 Correct 97 ms 34672 KB Output is correct
84 Correct 96 ms 34640 KB Output is correct
85 Correct 94 ms 34648 KB Output is correct
86 Correct 112 ms 35276 KB Output is correct
87 Correct 108 ms 35412 KB Output is correct
88 Correct 111 ms 35412 KB Output is correct
89 Correct 111 ms 35408 KB Output is correct
90 Correct 114 ms 35408 KB Output is correct
91 Correct 252 ms 56400 KB Output is correct
92 Correct 253 ms 56596 KB Output is correct
93 Correct 253 ms 56404 KB Output is correct
94 Correct 253 ms 56488 KB Output is correct
95 Correct 236 ms 56400 KB Output is correct
96 Incorrect 358 ms 61776 KB Output isn't correct
97 Halted 0 ms 0 KB -