Submission #385062

# Submission time Handle Problem Language Result Execution time Memory
385062 2021-04-03T06:47:30 Z qpwoeirut Horses (IOI15_horses) C++17
100 / 100
1122 ms 73836 KB
#include "horses.h"
#include <algorithm>
#include <cassert>
#include <iostream>
#include <set>

using namespace std;

typedef long long ll;

const int MN = 500005;
const ll MOD = 1e9 + 7;

ll binpow(ll x, ll e, const ll mod) {
    x %= mod;
    ll ret = 1;
    for (; e>0; e>>=1) {
        if (e&1) ret = (ret * x) % mod;
        x = (x * x) % mod;
    }
    return ret;
}

struct Node {
    int lo, hi;
    ll mx, pd;
};

int N;
ll A[MN], B[MN];
Node st[MN << 2];

void build() {
    int PN = N-1; PN |= PN >> 1; PN |= PN >> 2; PN |= PN >> 4; PN |= PN >> 8; PN |= PN >> 16; ++PN;
    for (int i=0; i<PN; ++i) {
        st[PN+i].lo = i;
        st[PN+i].hi = i+1;
        st[PN+i].mx = (i < N ? B[i] : 0);
        st[PN+i].pd = (i < N ? A[i] : 0);
    }
    for (int i=PN-1; i>0; --i) {
        st[i].lo = st[i << 1].lo;
        st[i].hi = st[(i << 1) + 1].hi;
        st[i].mx = max(st[i << 1].mx, st[(i << 1) + 1].mx);
        st[i].pd = (st[i << 1].pd * st[(i << 1) + 1].pd) % MOD;
    }
}

void upd_pd(const int node, const int pos, const ll val) {
    if (pos >= st[node].hi || st[node].lo >= pos + 1) return;
    if (pos <= st[node].lo && st[node].hi <= pos + 1) {
        st[node].pd = val;
        return;
    }

    upd_pd(node << 1, pos, val);
    upd_pd((node << 1) + 1, pos, val);

    st[node].pd = (st[node << 1].pd * st[(node << 1) + 1].pd) % MOD;
}

void upd_mx(const int node, const int pos, const ll val) {
    if (pos >= st[node].hi || st[node].lo >= pos + 1) return;
    if (pos <= st[node].lo && st[node].hi <= pos + 1) {
        st[node].mx = val;
        return;
    }

    upd_mx(node << 1, pos, val);
    upd_mx((node << 1) + 1, pos, val);

    st[node].mx = max(st[node << 1].mx, st[(node << 1) + 1].mx);
}

ll query_pd(const int node, const int lo, const int hi) {
    if (lo >= st[node].hi || st[node].lo >= hi) return 1;
    if (lo <= st[node].lo && st[node].hi <= hi) return st[node].pd;
    return (query_pd(node << 1, lo, hi) * query_pd((node << 1) + 1, lo, hi)) % MOD;
}

ll query_mx(const int node, const int lo, const int hi) {
    if (lo >= st[node].hi || st[node].lo >= hi) return 0;
    if (lo <= st[node].lo && st[node].hi <= hi) return st[node].mx;
    return max(query_mx(node << 1, lo, hi), query_mx((node << 1) + 1, lo, hi));
}

set<int, greater<int>> idxs;
int get_ans() {
    __int128 suff = 1;
    int idx = 0;
    for (const int i: idxs) {
        suff *= A[i];
        idx = i;
        if (suff >= MOD) break;
    }

    if (suff < MOD) idx = 0;

    ll pref = query_pd(1, 0, idx);
    __int128 best = query_mx(1, 0, N);
    for (const int i: idxs) {
        best = max(best, suff * query_mx(1, i, N));
        suff /= A[i];
        if (idx == i) break;
    }
    assert(suff == 1);
    best %= MOD;

    return (pref * best) % MOD;
}

int init(int n, int X[], int Y[]) {
    N = n;
    copy(X, X+N, A);
    copy(Y, Y+N, B);

    build();
    
    for (int i=0; i<N; ++i) {
        if (A[i] > 1) idxs.insert(i);
    }

    return get_ans();
}

int updateX(int pos, int val) {
    upd_pd(1, pos, val);
    if (A[pos] == 1 && val > 1) {
        idxs.insert(pos);
    } else if (A[pos] > 1 && val == 1) {
        idxs.erase(pos);
    }
    A[pos] = val;
    return get_ans();
}

int updateY(int pos, int val) {
    upd_mx(1, pos, val);
    B[pos] = val;
	return get_ans();
}

/*
pref/suff of X (point update pduct)
max pref[i] * Y[i] (lazy multiply max segtree)

inverse mod p^2 = pow(x, (p*(p-1)) - 1, p^2)
*/

Compilation message

horses.cpp: In function 'int get_ans()':
horses.cpp:109:26: warning: conversion from '__int128' to 'int' may change value [-Wconversion]
  109 |     return (pref * best) % MOD;
      |            ~~~~~~~~~~~~~~^~~~~
# Verdict Execution time Memory Grader output
1 Correct 1 ms 364 KB Output is correct
2 Correct 1 ms 364 KB Output is correct
3 Correct 1 ms 364 KB Output is correct
4 Correct 1 ms 364 KB Output is correct
5 Correct 1 ms 364 KB Output is correct
6 Correct 1 ms 364 KB Output is correct
7 Correct 1 ms 364 KB Output is correct
8 Correct 2 ms 364 KB Output is correct
9 Correct 1 ms 364 KB Output is correct
10 Correct 1 ms 364 KB Output is correct
11 Correct 1 ms 364 KB Output is correct
12 Correct 1 ms 364 KB Output is correct
13 Correct 1 ms 364 KB Output is correct
14 Correct 1 ms 384 KB Output is correct
15 Correct 1 ms 364 KB Output is correct
16 Correct 1 ms 364 KB Output is correct
17 Correct 1 ms 384 KB Output is correct
18 Correct 1 ms 364 KB Output is correct
19 Correct 1 ms 364 KB Output is correct
20 Correct 1 ms 364 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 364 KB Output is correct
2 Correct 1 ms 364 KB Output is correct
3 Correct 1 ms 364 KB Output is correct
4 Correct 1 ms 364 KB Output is correct
5 Correct 1 ms 364 KB Output is correct
6 Correct 1 ms 364 KB Output is correct
7 Correct 1 ms 364 KB Output is correct
8 Correct 1 ms 364 KB Output is correct
9 Correct 1 ms 364 KB Output is correct
10 Correct 1 ms 364 KB Output is correct
11 Correct 1 ms 364 KB Output is correct
12 Correct 1 ms 364 KB Output is correct
13 Correct 1 ms 364 KB Output is correct
14 Correct 1 ms 364 KB Output is correct
15 Correct 2 ms 364 KB Output is correct
16 Correct 1 ms 364 KB Output is correct
17 Correct 1 ms 512 KB Output is correct
18 Correct 1 ms 364 KB Output is correct
19 Correct 1 ms 364 KB Output is correct
20 Correct 1 ms 364 KB Output is correct
21 Correct 1 ms 364 KB Output is correct
22 Correct 1 ms 364 KB Output is correct
23 Correct 2 ms 492 KB Output is correct
24 Correct 2 ms 492 KB Output is correct
25 Correct 3 ms 492 KB Output is correct
26 Correct 2 ms 492 KB Output is correct
27 Correct 6 ms 492 KB Output is correct
28 Correct 3 ms 492 KB Output is correct
29 Correct 2 ms 364 KB Output is correct
30 Correct 2 ms 492 KB Output is correct
31 Correct 4 ms 492 KB Output is correct
32 Correct 6 ms 492 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 943 ms 61964 KB Output is correct
2 Correct 373 ms 61420 KB Output is correct
3 Correct 400 ms 61548 KB Output is correct
4 Correct 436 ms 61400 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 364 KB Output is correct
2 Correct 1 ms 364 KB Output is correct
3 Correct 1 ms 364 KB Output is correct
4 Correct 1 ms 364 KB Output is correct
5 Correct 2 ms 364 KB Output is correct
6 Correct 1 ms 364 KB Output is correct
7 Correct 2 ms 364 KB Output is correct
8 Correct 1 ms 384 KB Output is correct
9 Correct 1 ms 364 KB Output is correct
10 Correct 1 ms 364 KB Output is correct
11 Correct 1 ms 364 KB Output is correct
12 Correct 1 ms 364 KB Output is correct
13 Correct 2 ms 364 KB Output is correct
14 Correct 1 ms 364 KB Output is correct
15 Correct 1 ms 364 KB Output is correct
16 Correct 1 ms 384 KB Output is correct
17 Correct 1 ms 364 KB Output is correct
18 Correct 1 ms 364 KB Output is correct
19 Correct 1 ms 364 KB Output is correct
20 Correct 1 ms 364 KB Output is correct
21 Correct 1 ms 364 KB Output is correct
22 Correct 1 ms 364 KB Output is correct
23 Correct 2 ms 492 KB Output is correct
24 Correct 2 ms 492 KB Output is correct
25 Correct 3 ms 492 KB Output is correct
26 Correct 2 ms 492 KB Output is correct
27 Correct 7 ms 492 KB Output is correct
28 Correct 3 ms 512 KB Output is correct
29 Correct 2 ms 364 KB Output is correct
30 Correct 2 ms 492 KB Output is correct
31 Correct 4 ms 492 KB Output is correct
32 Correct 7 ms 492 KB Output is correct
33 Correct 72 ms 41068 KB Output is correct
34 Correct 66 ms 40940 KB Output is correct
35 Correct 207 ms 71148 KB Output is correct
36 Correct 202 ms 71148 KB Output is correct
37 Correct 156 ms 39020 KB Output is correct
38 Correct 129 ms 51820 KB Output is correct
39 Correct 47 ms 38764 KB Output is correct
40 Correct 181 ms 66284 KB Output is correct
41 Correct 86 ms 38892 KB Output is correct
42 Correct 112 ms 38892 KB Output is correct
43 Correct 171 ms 66668 KB Output is correct
44 Correct 167 ms 66796 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 1 ms 364 KB Output is correct
2 Correct 1 ms 364 KB Output is correct
3 Correct 1 ms 364 KB Output is correct
4 Correct 2 ms 364 KB Output is correct
5 Correct 1 ms 512 KB Output is correct
6 Correct 1 ms 364 KB Output is correct
7 Correct 1 ms 364 KB Output is correct
8 Correct 1 ms 364 KB Output is correct
9 Correct 1 ms 364 KB Output is correct
10 Correct 1 ms 364 KB Output is correct
11 Correct 1 ms 364 KB Output is correct
12 Correct 1 ms 364 KB Output is correct
13 Correct 2 ms 364 KB Output is correct
14 Correct 1 ms 364 KB Output is correct
15 Correct 1 ms 364 KB Output is correct
16 Correct 1 ms 364 KB Output is correct
17 Correct 1 ms 364 KB Output is correct
18 Correct 1 ms 364 KB Output is correct
19 Correct 1 ms 364 KB Output is correct
20 Correct 2 ms 364 KB Output is correct
21 Correct 1 ms 364 KB Output is correct
22 Correct 1 ms 364 KB Output is correct
23 Correct 2 ms 492 KB Output is correct
24 Correct 2 ms 528 KB Output is correct
25 Correct 4 ms 492 KB Output is correct
26 Correct 2 ms 492 KB Output is correct
27 Correct 6 ms 492 KB Output is correct
28 Correct 3 ms 492 KB Output is correct
29 Correct 2 ms 364 KB Output is correct
30 Correct 2 ms 492 KB Output is correct
31 Correct 4 ms 492 KB Output is correct
32 Correct 5 ms 492 KB Output is correct
33 Correct 966 ms 65392 KB Output is correct
34 Correct 415 ms 73836 KB Output is correct
35 Correct 403 ms 65260 KB Output is correct
36 Correct 471 ms 68820 KB Output is correct
37 Correct 66 ms 40812 KB Output is correct
38 Correct 63 ms 41100 KB Output is correct
39 Correct 207 ms 71148 KB Output is correct
40 Correct 200 ms 71388 KB Output is correct
41 Correct 152 ms 39148 KB Output is correct
42 Correct 129 ms 51820 KB Output is correct
43 Correct 56 ms 38764 KB Output is correct
44 Correct 184 ms 66156 KB Output is correct
45 Correct 86 ms 38892 KB Output is correct
46 Correct 109 ms 38892 KB Output is correct
47 Correct 168 ms 66668 KB Output is correct
48 Correct 169 ms 66540 KB Output is correct
49 Correct 286 ms 43884 KB Output is correct
50 Correct 246 ms 43884 KB Output is correct
51 Correct 425 ms 73068 KB Output is correct
52 Correct 322 ms 72684 KB Output is correct
53 Correct 1122 ms 42252 KB Output is correct
54 Correct 460 ms 55788 KB Output is correct
55 Correct 162 ms 39916 KB Output is correct
56 Correct 340 ms 68076 KB Output is correct
57 Correct 555 ms 40628 KB Output is correct
58 Correct 784 ms 40940 KB Output is correct
59 Correct 169 ms 66668 KB Output is correct