#include "bits/stdc++.h"
using namespace std;
#define int long long
const int MAXN = 3e5 + 10;
const int MOD = 1e9 + 7;
mt19937_64 rng((int)std::chrono::steady_clock::now().time_since_epoch().count());
int rnd(int x, int y) {
int u = uniform_int_distribution<int>(x, y)(rng); return u;
}
int bm(int b, int p) {
if(p==0) return 1 % MOD;
int r = bm(b, p >> 1);
if(p&1) return (((r*r) % MOD) * b) % MOD;
return (r*r) % MOD;
}
int inv(int b) {
return bm(b, MOD-2);
}
int fastlog(int x) {
return (x == 0 ? -1 : 64 - __builtin_clzll(x) - 1);
}
void printcase(int i) { cout << "Case #" << i << ": "; }
int freq[MAXN], ff[MAXN], n;
int block;
bool modui(pair<pair<int,int>,int> a, pair<pair<int,int>,int> b) {
if(a.first.first / block != b.first.first / block) return a.first.first / block < b.first.first / block;
return a.first.second < b.first.second;
}
bool byfreq(int x, int y) {
return freq[x] < freq[y];
}
void solve(int tc) {
int q;
cin >> n >> q;
int c = 300000;
block = ceil(sqrt(n));
int arr[n+1];
for(int i=1; i<=n; i++) cin >> arr[i];
pair<pair<int,int>,int> qry[q+1];
for(int i=1; i<=q; i++) {
cin >> qry[i].first.first >> qry[i].first.second;
qry[i].second = i;
}
sort(qry+1, qry+q+1, modui);
int res[q+1];
int l = qry[1].first.first, r = qry[1].first.first - 1;
vector<int> v;
// O(sqrt(n)) frequencies
set<int> all; // set of all frequencies
for(int i=1; i<=q; i++) {
while(l > qry[i].first.first) {
l--;
ff[freq[arr[l]]]--;
if(ff[freq[arr[l]]] == 0) all.erase(freq[arr[l]]);
freq[arr[l]]++;
ff[freq[arr[l]]]++;
if(ff[freq[arr[l]]] == 1) all.insert(freq[arr[l]]);
}
while(r > qry[i].first.second) {
ff[freq[arr[r]]]--;
if(ff[freq[arr[r]]] == 0) all.erase(freq[arr[r]]);
freq[arr[r]]--;
ff[freq[arr[r]]]++;
if(ff[freq[arr[r]]] == 1) all.insert(freq[arr[r]]);
r--;
}
while(r < qry[i].first.second) {
r++;
ff[freq[arr[r]]]--;
if(ff[freq[arr[r]]] == 0) all.erase(freq[arr[r]]);
freq[arr[r]]++;
ff[freq[arr[r]]]++;
if(ff[freq[arr[r]]] == 1) all.insert(freq[arr[r]]);
}
while(l < qry[i].first.first) {
ff[freq[arr[l]]]--;
if(ff[freq[arr[l]]] == 0) all.erase(freq[arr[l]]);
freq[arr[l]]--;
ff[freq[arr[l]]]++;
if(ff[freq[arr[l]]] == 1) all.insert(freq[arr[l]]);
l++;
}
queue<pair<int, int> > ln; stack<pair<int, int> > rn;
int m = 0;
for(int x: all) {
// {x, ff[x]}
if(ff[x] / 2 > 0) {
ln.push({x, ff[x] / 2});
rn.push({x, ff[x] / 2});
m += 2;
}
if(ff[x] % 2) {
if(ln.size() > rn.size()) rn.push({x, 1});
else ln.push({x, 1});
m++;
}
}
pair<int, int> b[m];
int it = 0;
while(ln.size()) {
b[it++] = ln.front(); ln.pop();
}
while(rn.size()) {
b[it++] = rn.top(); rn.pop();
}
int ans = 0;
int ws=0,wps=0,tps=0,ps=0;
for(int i=0; i<m; i++) {
int F = b[i].first;
int C = b[i].second;
int t1 = C * F * ps * ws;
int t2 = C * F * wps;
int t3 = F * C * (C+1) / 2 * ws - F * C * tps / 2 + F * C * ws;
ans += t1 - t2 + t3;
ws += C * F;
wps += C * F * ps;
tps += C * F * (C + 1);
ps += C;
}
for(int i=0; i<m; i++) ans += b[i].first * (b[i].first + 1) / 2 * b[i].second;
res[qry[i].second] = ans;
}
for(int i=1; i<=q; i++) cout << res[i] << "\n";
}
int32_t main() {
ios::sync_with_stdio(0); cin.tie(0);
int t = 1; //cin >> t;
for(int i=1; i<=t; i++) solve(i);
}
Compilation message
diversity.cpp: In function 'void solve(long long int)':
diversity.cpp:38:7: warning: unused variable 'c' [-Wunused-variable]
38 | int c = 300000;
| ^
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
0 ms |
336 KB |
Output is correct |
3 |
Correct |
0 ms |
340 KB |
Output is correct |
4 |
Correct |
0 ms |
340 KB |
Output is correct |
5 |
Correct |
0 ms |
332 KB |
Output is correct |
6 |
Incorrect |
0 ms |
212 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
2 ms |
340 KB |
Output is correct |
4 |
Correct |
14 ms |
1492 KB |
Output is correct |
5 |
Correct |
25 ms |
2768 KB |
Output is correct |
6 |
Correct |
42 ms |
3788 KB |
Output is correct |
7 |
Correct |
36 ms |
4452 KB |
Output is correct |
8 |
Correct |
42 ms |
3696 KB |
Output is correct |
9 |
Correct |
42 ms |
3700 KB |
Output is correct |
10 |
Correct |
39 ms |
4044 KB |
Output is correct |
11 |
Correct |
42 ms |
3820 KB |
Output is correct |
12 |
Correct |
41 ms |
3660 KB |
Output is correct |
13 |
Correct |
40 ms |
3704 KB |
Output is correct |
14 |
Correct |
40 ms |
3912 KB |
Output is correct |
15 |
Correct |
40 ms |
3808 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
2 ms |
340 KB |
Output is correct |
4 |
Correct |
14 ms |
1492 KB |
Output is correct |
5 |
Correct |
25 ms |
2768 KB |
Output is correct |
6 |
Correct |
42 ms |
3788 KB |
Output is correct |
7 |
Correct |
36 ms |
4452 KB |
Output is correct |
8 |
Correct |
42 ms |
3696 KB |
Output is correct |
9 |
Correct |
42 ms |
3700 KB |
Output is correct |
10 |
Correct |
39 ms |
4044 KB |
Output is correct |
11 |
Correct |
42 ms |
3820 KB |
Output is correct |
12 |
Correct |
41 ms |
3660 KB |
Output is correct |
13 |
Correct |
40 ms |
3704 KB |
Output is correct |
14 |
Correct |
40 ms |
3912 KB |
Output is correct |
15 |
Correct |
40 ms |
3808 KB |
Output is correct |
16 |
Correct |
0 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
328 KB |
Output is correct |
18 |
Correct |
2 ms |
340 KB |
Output is correct |
19 |
Correct |
15 ms |
1496 KB |
Output is correct |
20 |
Correct |
34 ms |
2636 KB |
Output is correct |
21 |
Correct |
48 ms |
3932 KB |
Output is correct |
22 |
Correct |
46 ms |
3688 KB |
Output is correct |
23 |
Correct |
46 ms |
3804 KB |
Output is correct |
24 |
Correct |
47 ms |
3664 KB |
Output is correct |
25 |
Correct |
48 ms |
3728 KB |
Output is correct |
26 |
Correct |
50 ms |
3672 KB |
Output is correct |
27 |
Correct |
47 ms |
3756 KB |
Output is correct |
28 |
Correct |
45 ms |
3896 KB |
Output is correct |
29 |
Correct |
46 ms |
3692 KB |
Output is correct |
30 |
Correct |
47 ms |
3724 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
2 ms |
340 KB |
Output is correct |
4 |
Correct |
14 ms |
1492 KB |
Output is correct |
5 |
Correct |
25 ms |
2768 KB |
Output is correct |
6 |
Correct |
42 ms |
3788 KB |
Output is correct |
7 |
Correct |
36 ms |
4452 KB |
Output is correct |
8 |
Correct |
42 ms |
3696 KB |
Output is correct |
9 |
Correct |
42 ms |
3700 KB |
Output is correct |
10 |
Correct |
39 ms |
4044 KB |
Output is correct |
11 |
Correct |
42 ms |
3820 KB |
Output is correct |
12 |
Correct |
41 ms |
3660 KB |
Output is correct |
13 |
Correct |
40 ms |
3704 KB |
Output is correct |
14 |
Correct |
40 ms |
3912 KB |
Output is correct |
15 |
Correct |
40 ms |
3808 KB |
Output is correct |
16 |
Correct |
0 ms |
340 KB |
Output is correct |
17 |
Correct |
1 ms |
328 KB |
Output is correct |
18 |
Correct |
2 ms |
340 KB |
Output is correct |
19 |
Correct |
15 ms |
1496 KB |
Output is correct |
20 |
Correct |
34 ms |
2636 KB |
Output is correct |
21 |
Correct |
48 ms |
3932 KB |
Output is correct |
22 |
Correct |
46 ms |
3688 KB |
Output is correct |
23 |
Correct |
46 ms |
3804 KB |
Output is correct |
24 |
Correct |
47 ms |
3664 KB |
Output is correct |
25 |
Correct |
48 ms |
3728 KB |
Output is correct |
26 |
Correct |
50 ms |
3672 KB |
Output is correct |
27 |
Correct |
47 ms |
3756 KB |
Output is correct |
28 |
Correct |
45 ms |
3896 KB |
Output is correct |
29 |
Correct |
46 ms |
3692 KB |
Output is correct |
30 |
Correct |
47 ms |
3724 KB |
Output is correct |
31 |
Incorrect |
1 ms |
340 KB |
Output isn't correct |
32 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
0 ms |
336 KB |
Output is correct |
3 |
Correct |
0 ms |
340 KB |
Output is correct |
4 |
Correct |
0 ms |
340 KB |
Output is correct |
5 |
Correct |
0 ms |
332 KB |
Output is correct |
6 |
Incorrect |
0 ms |
212 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
0 ms |
336 KB |
Output is correct |
3 |
Correct |
0 ms |
340 KB |
Output is correct |
4 |
Correct |
0 ms |
340 KB |
Output is correct |
5 |
Correct |
0 ms |
332 KB |
Output is correct |
6 |
Incorrect |
0 ms |
212 KB |
Output isn't correct |
7 |
Halted |
0 ms |
0 KB |
- |