답안 #696862

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
696862 2023-02-07T12:41:05 Z cig32 Diversity (CEOI21_diversity) C++17
18 / 100
50 ms 4452 KB
#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 -