Submission #1195319

#	Time	Username	Problem	Language	Result	Execution time	Memory
1195319		Zbyszek99	Diversity (CEOI21_diversity)	C++20	100 / 100	6258 ms	241696 KiB

diversity

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#define ll long long
#define ld long double
#define ull unsigned long long
#define ff first
#define ss second
#define pii pair<int,int>
#define pll pair<long long, long long>
#define vi vector<int>
#define vl vector<long long>
#define pb push_back
#define rep(i, b) for(int i = 0; i < (b); ++i)
#define rep2(i,a,b) for(int i = a; i <= (b); ++i)
#define rep3(i,a,b,c) for(int i = a; i <= (b); i+=c)
#define count_bits(x) __builtin_popcountll((x))
#define all(x) (x).begin(),(x).end()
#define siz(x) (int)(x).size()
#define forall(it,x) for(auto& it:(x))
using namespace __gnu_pbds;
using namespace std;
typedef tree<int, null_type, less<int>, rb_tree_tag,tree_order_statistics_node_update> ordered_set;
//mt19937 mt;void random_start(){mt.seed(chrono::time_point_cast<chrono::milliseconds>(chrono::high_resolution_clock::now()).time_since_epoch().count());}
//ll los(ll a, ll b) {return a + (mt() % (b-a+1));}
const int INF = 1e9+50;
const ll INF_L = 1e18+40;
const ll MOD = 1e9+7;

int sort_type = 0;

struct query
{
    int l,r,ind;
    bool operator<(const query& other) const
    {
        if(sort_type == 0) return l < other.l;
        else return r < other.r;
    }
};

const int sqr2 = 100;
const int sqr = 520;

int elm_cnt[300'003];
ll oc_cnt[300'003];
unordered_set<int> non_zero;
int n,q;
int arr[300'001];
ll sqr_sum[300'001][sqr2+1];
ll query_ans[50'001];
vector<query> block[300'001/sqr+10];

void change_num(int x, int p)
{
    oc_cnt[elm_cnt[x]]--;
    if(oc_cnt[elm_cnt[x]] == 0) non_zero.erase(non_zero.find(elm_cnt[x]));
    elm_cnt[x]+=p;
    oc_cnt[elm_cnt[x]]++;
    if(oc_cnt[elm_cnt[x]] == 1) non_zero.insert(elm_cnt[x]);
}

vl oc;

ll f(ll n2, ll it, ll pref, ll ile)
{
    if(ile == 0) return 0;
    return (pref - it >= 0 ? -sqr_sum[pref-it][it] : 0) + sqr_sum[pref+(ile-1)*it][it] + (n2-it) * (ile*pref + it * ((ile * (ile-1))/2));
}

ll get_ans(int n2)
{
    oc = {};
    forall(it,non_zero) oc.pb(it);
    sort(all(oc));
    ll ans = 0;
    ll pref = 0;
    ll suf = 0;
    forall(it,oc)
    {
        ans += oc_cnt[it] * (it * (n2-it) + (it * (it-1)) / 2 + it);
        if(it <= sqr2)
        {
            if(pref > suf) swap(pref,suf);
            ll ile = min(oc_cnt[it],(suf-pref+it-1)/it);
            ans += f(n2,it,pref,ile+(oc_cnt[it]-ile)/2);
            pref += ((oc_cnt[it]-ile)/2+ile)*it;
            ans += f(n2,it,suf,(oc_cnt[it]-ile+1)/2);
            suf += (oc_cnt[it]-ile+1)/2 * it;
        }
        else
        {
            rep(k,oc_cnt[it])
            {
                if(pref < suf)
                {
                    ans += pref * (n2-pref-it);
                    pref += it;
                }
                else
                {
                    ans += suf * (n2-suf-it);
                    suf += it;
                }
            }
        }
    }
    return ans;
}

int main()
{
    ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    //random_start();
    cin >> n >> q;
    rep(i,n)
    {
        cin >> arr[i];
    }
    rep2(i,1,300'000)
    {
        rep2(d,1,sqr2)
        {
            sqr_sum[i][d] = -(ll)i*(ll)i + (i-d >= 0 ? sqr_sum[i-d][d] : 0);
        }
    }
    oc_cnt[0] = 10000000;
    vector<query> queries;
    rep(i,q)
    {
        int l,r;
        cin >> l >> r;
        l--;
        r--;
        queries.pb({l,r,i});
    }
    sort_type = 0;
    sort(all(queries));
    forall(it,queries)
    {
        block[it.l/sqr].pb(it);
    }
    rep(b,300'001/sqr+10)
    {
        sort_type = 1;
        sort(all(block[b]));
        if(siz(block[b]) == 0) continue;
        rep2(k,block[b][0].l,block[b][0].r) change_num(arr[k],1);
        query_ans[block[b][0].ind] = get_ans(block[b][0].r-block[b][0].l+1);
        int cur_l = block[b][0].l;
        int cur_r = block[b][0].r;
        rep2(i,1,siz(block[b])-1)
        {
            while(cur_r < block[b][i].r)
            {
                cur_r++;
                change_num(arr[cur_r],1);
            }
            while(cur_l < block[b][i].l)
            {
                cur_l++;
                change_num(arr[cur_l-1],-1);
            }
            while(cur_l > block[b][i].l)
            {
                cur_l--;
                change_num(arr[cur_l],1);
            }
            query_ans[block[b][i].ind] = get_ans(block[b][i].r-block[b][i].l+1);
        }
        rep2(k,block[b].back().l,block[b].back().r) change_num(arr[k],-1);
    }
    rep(i,q) cout << query_ans[i] << "\n";
}

• Subtask 1: 4 / 4
• Subtask 2: 10 / 10
• Subtask 3: 8 / 8
• Subtask 4: 16 / 16
• Subtask 5: 26 / 26
• Subtask 6: 36 / 36