답안 #735517

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
735517 2023-05-04T09:33:38 Z GrindMachine Zamjene (COCI16_zamjene) C++17
140 / 140
2475 ms 121120 KB
// Om Namah Shivaya

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

using namespace std;
using namespace __gnu_pbds;

template<typename T> using Tree = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
typedef long long int ll;
typedef long double ld;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

#define fastio ios_base::sync_with_stdio(false); cin.tie(NULL)
#define pb push_back
#define endl '\n'
#define sz(a) a.size()
#define setbits(x) __builtin_popcountll(x)
#define ff first
#define ss second
#define conts continue
#define ceil2(x, y) ((x + y - 1) / (y))
#define all(a) a.begin(), a.end()
#define rall(a) a.rbegin(), a.rend()
#define yes cout << "Yes" << endl
#define no cout << "No" << endl

#define rep(i, n) for(int i = 0; i < n; ++i)
#define rep1(i, n) for(int i = 1; i <= n; ++i)
#define rev(i, s, e) for(int i = s; i >= e; --i)
#define trav(i, a) for(auto &i : a)

template<typename T>
void amin(T &a, T b) {
    a = min(a, b);
}

template<typename T>
void amax(T &a, T b) {
    a = max(a, b);
}

#ifdef LOCAL
#include "debug.h"
#else
#define debug(x) 42
#endif

/*

refs:
https://bigprimes.org/ (for generating big prime MOD)


sol idea:
dsu + multiset hashing

*/

const ll MOD = 964863244005699737;
const int N = 1e5 + 5;
const int inf1 = int(1e9) + 5;
const ll inf2 = ll(1e18) + 5;

struct DSU {
    vector<int> par, rankk, siz;
    vector<ll> hsh;
    unordered_map<ll, ll> hashes;
    ll cnt = 0;

    DSU() {

    }

    DSU(int n) {
        init(n);
    }

    void init(int n) {
        par = vector<int>(n + 1);
        rankk = vector<int>(n + 1);
        siz = vector<int>(n + 1);
        hsh = vector<ll>(n + 1);
        rep(i, n + 1) create(i);
    }

    void create(int u) {
        par[u] = u;
        rankk[u] = 0;
        siz[u] = 1;
    }

    int find(int u) {
        if (u == par[u]) return u;
        else return par[u] = find(par[u]);
    }

    bool same(int u, int v) {
        return find(u) == find(v);
    }

    void merge(int u, int v) {
        u = find(u), v = find(v);
        if (u == v) return;

        if (rankk[u] == rankk[v]) rankk[u]++;
        if (rankk[u] < rankk[v]) swap(u, v);

        par[v] = u;

        cnt -= siz[u] * hashes[MOD - hsh[u]];
        hashes[hsh[u]] -= siz[u];

        cnt -= siz[v] * hashes[MOD - hsh[v]];
        hashes[hsh[v]] -= siz[v];

        hsh[u] = (hsh[u] + hsh[v]) % MOD;
        siz[u] += siz[v];

        cnt += siz[u] * hashes[MOD - hsh[u]];
        hashes[hsh[u]] += siz[u];
    }
};

void solve(int test_case)
{
    ll n, q; cin >> n >> q;
    vector<ll> a(n + 5);
    rep1(i, n) cin >> a[i];

    mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
    uniform_int_distribution<ll> dist(1, MOD - 1);

    auto b = a;
    sort(b.begin() + 1, b.begin() + n + 1);

    map<ll, ll> mp;
    rep1(i, n) {
        if (mp.count(a[i])) conts;
        mp[a[i]] = dist(rng);
    }

    DSU dsu(n + 5);

    rep1(i, n) {
        ll v = (mp[a[i]] - mp[b[i]] + MOD) % MOD;
        dsu.hsh[i] = v;
        dsu.cnt += dsu.hashes[MOD - v];
        dsu.hashes[v]++;
    }

    while (q--) {
        ll t; cin >> t;
        if (t == 1) {
            ll u, v; cin >> u >> v;
            ll pu = dsu.find(u), pv = dsu.find(v);


            dsu.cnt -= dsu.siz[pu] * dsu.hashes[MOD - dsu.hsh[pu]];
            dsu.hashes[dsu.hsh[pu]] -= dsu.siz[pu];

            dsu.hsh[pu] += mp[a[v]] - mp[a[u]];
            dsu.hsh[pu] = (dsu.hsh[pu] % MOD + MOD) % MOD;

            dsu.cnt += dsu.siz[pu] * dsu.hashes[MOD - dsu.hsh[pu]];
            dsu.hashes[dsu.hsh[pu]] += dsu.siz[pu];



            dsu.cnt -= dsu.siz[pv] * dsu.hashes[MOD - dsu.hsh[pv]];
            dsu.hashes[dsu.hsh[pv]] -= dsu.siz[pv];

            dsu.hsh[pv] += mp[a[u]] - mp[a[v]];
            dsu.hsh[pv] = (dsu.hsh[pv] % MOD + MOD) % MOD;

            dsu.cnt += dsu.siz[pv] * dsu.hashes[MOD - dsu.hsh[pv]];
            dsu.hashes[dsu.hsh[pv]] += dsu.siz[pv];


            swap(a[u], a[v]);
        }
        else if (t == 2) {
            ll u, v; cin >> u >> v;
            dsu.merge(u, v);
        }
        else if (t == 3) {
            if (dsu.hashes[0] == n) cout << "DA" << endl;
            else cout << "NE" << endl;
        }
        else {
            ll ans = dsu.cnt;
            cout << ans << endl;
        }
    }
}

int main()
{
    fastio;

    int t = 1;
    // cin >> t;

    rep1(i, t) {
        solve(i);
    }

    return 0;
}
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 1 ms 212 KB Output is correct
3 Correct 1 ms 212 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 212 KB Output is correct
2 Correct 1 ms 212 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 324 KB Output is correct
2 Correct 1 ms 340 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 340 KB Output is correct
2 Correct 2 ms 340 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 1 ms 328 KB Output is correct
2 Correct 1 ms 340 KB Output is correct
3 Correct 1 ms 340 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 4 ms 724 KB Output is correct
2 Correct 5 ms 760 KB Output is correct
3 Correct 5 ms 724 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 45 ms 4188 KB Output is correct
2 Correct 68 ms 5284 KB Output is correct
3 Correct 86 ms 7220 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 841 ms 30436 KB Output is correct
2 Correct 2173 ms 99524 KB Output is correct
3 Correct 2475 ms 121120 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 1164 ms 63580 KB Output is correct
2 Correct 1707 ms 79608 KB Output is correct
3 Correct 1303 ms 77644 KB Output is correct
# 결과 실행 시간 메모리 Grader output
1 Correct 752 ms 41808 KB Output is correct
2 Correct 1166 ms 65916 KB Output is correct
3 Correct 1260 ms 77524 KB Output is correct