Submission #978656

# Submission time Handle Problem Language Result Execution time Memory
978656 2024-05-09T12:31:41 Z vjudge1 Sprinkler (JOI22_sprinkler) C++17
41 / 100
4000 ms 123144 KB
#include <bits/stdc++.h>
#pragma GCC optimize("unroll-loops")
#pragma GCC target("avx2")

using namespace std;

int n,L,q;
vector<int> g[200005],f[200005];
int t[200005];
int h[200005],tin[200005],curtin;
int v[200005],pos[200005];
int firstson[200005][45],lastson[200005][45];

void dfs(int nod)
{
    tin[nod] = ++curtin;
    for (auto vecin : g[nod])
    {
        if (vecin != t[nod])
        {
            t[vecin] = nod;
            f[nod].push_back(vecin);
            h[vecin] = 1 + h[nod];
            dfs(vecin);
        }
    }
}

bool cmp(int x,int y)
{
    if (h[x] != h[y])
        return h[x] < h[y];
    return tin[x] < tin[y];
}

int aint[800005];

int vkk[200005];

int bulan = 30;

void update(int nod,int l,int r,int st,int dr,int val)
{
    if (r < st or dr < l)
        return;
    st = max(st,l);
    dr = min(dr,r);
    if (st == l and r == dr)
    {
        aint[nod] = 1ll * aint[nod] * val % L;
        //cout << nod << ' ' << l << ' ' << r << ' ' << aint[nod] << endl;
        return;
    }
    if (dr - st + 1 <= bulan)
    {
        for (int i = st; i <= dr; i++)
            vkk[i] = 1ll * vkk[i] * val % L;
        return;
    }
    int mij = (l + r) / 2;
    update(2 * nod,l,mij,st,dr,val);
    update(2 * nod + 1,mij + 1,r,st,dr,val);
}

int paint[200005];

void build(int nod,int l,int r)
{
    if (l == r)
        paint[l] = nod;
    else
    {
        int mij = (l + r) / 2;
        build(2 * nod,l,mij);
        build(2 * nod + 1,mij + 1,r);
    }
}

signed main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    cout.tie(NULL);
    cin >> n >> L;
    for (int i = 1; i < n; i++)
    {
        int x,y;
        cin >> x >> y;
        g[x].push_back(y);
        g[y].push_back(x);
    }
    dfs(1);
    for (int i = 1; i <= n; i++)
        v[i] = i;
    sort(v + 1,v + n + 1,cmp);
    for (int i = 1; i <= n; i++)
        pos[v[i]] = i;
    for (int i = 1; i <= 4 * n; i++)
        aint[i] = 1;
    build(1,1,n);
    for (int i = 1; i <= n; i++)
        vkk[i] = 1;
    for (int i = 1; i <= n; i++)
    {
        int xx;
        cin >> xx;
        update(1,1,n,pos[i],pos[i],xx);
    }
    for (int i = 1; i <= n; i++)
    {
        int nod = i;
        for (int j = 0; j <= 40; j++)
        {
            if (firstson[nod][j] == 0 or pos[i] < pos[firstson[nod][j]])
                firstson[nod][j] = i;
            if (lastson[nod][j] == 0 or pos[i] > pos[lastson[nod][j]])
                lastson[nod][j] = i;
            nod = t[nod];
            if (nod == 0)
                break;
        }
    }
    cin >> q;
    for (int i = 1; i <= q; i++)
    {
        int tip;
        cin >> tip;
        if (tip == 1)
        {
            int x,d,w;
            cin >> x >> d >> w;
            vector<int> ruta;
            int nod = x;
            for (int j = 0; j <= min(d,h[x]); j++)
            {
                ruta.push_back(nod);
                nod = t[nod];
            }
            nod = ruta.back();
            ruta.pop_back();
            int dist = d - (h[x] - h[nod]);
            for (int j = 0; j <= dist; j++)
            {
                if (firstson[nod][j] == 0)
                    break;
                int lu = pos[firstson[nod][j]],ru = pos[lastson[nod][j]];
                update(1,1,n,lu,ru,w);
            }
            while (!ruta.empty())
            {
                nod = ruta.back();
                ruta.pop_back();
                int dist = d - (h[x] - h[nod]);
                for (int j = dist - 1; j <= dist; j++)
                {
                    if (firstson[nod][j] == 0)
                        break;
                    int lu = pos[firstson[nod][j]],ru = pos[lastson[nod][j]];
                    if (ru - lu + 1 <= bulan)
                    {
                        for (int jj = lu; jj <= ru; jj++)
                            vkk[jj] = 1ll * vkk[jj] * w % L;
                    }
                    else
                        update(1,1,n,lu,ru,w);
                }
            }
        }
        else
        {
            int x;
            cin >> x;
            int vjj = vkk[pos[x]];
            x = paint[pos[x]];
            int rsp = 1ll * aint[1] * vjj % L;
            while (x != 1)
            {
                rsp = 1ll * rsp * aint[x] % L;
                x /= 2;
            }
            //cout << pos[x] << ' ';
            cout << rsp << '\n';
        }
    }
    return 0;
}

/**
Rezolvarea se bazeaza pe faptul ca D mic
Dupa ce fac ordinea dfs, pun nodurile intr-un sir sortate dupa inaltime si, la egalitate, dupa tin
La un update, iau al max(h[i],D)-lea stramos al lui x, fie el y. Updatez, pe inaltimile h[y] + k pentru k de la 0 la D - (h[x] - h[y]), subarborele
Apoi, merg nod dupa nod cu y spre x, la fiecare pas updatez ultimele doua nivele la care ajunge nodul curent y
Un update este de range prod si query de point value, deci pot face cu un aib care tine chestiile modulo L
Complexitate update: D * log, query: log
Pentru a afla care sunt primul si ultimul nod de la o inaltime care sunt in subarborele meu, inaltime <= D + h[nod], pot doar sa precalculez in N*D

Complexitate: O(N * D + N * logN + Q * D * logN)

Exista o micuta problema, aib-ul nu cred ca isi face treaba deoarece L nu e prim, incercam aint si fingers crossed sa nu dea TLE
**/
# Verdict Execution time Memory Grader output
1 Correct 4 ms 14684 KB Output is correct
2 Correct 3 ms 14832 KB Output is correct
3 Correct 3 ms 14684 KB Output is correct
4 Correct 4 ms 15196 KB Output is correct
5 Correct 5 ms 15196 KB Output is correct
6 Correct 6 ms 15196 KB Output is correct
7 Correct 5 ms 15192 KB Output is correct
8 Correct 4 ms 15192 KB Output is correct
9 Correct 3 ms 14684 KB Output is correct
10 Correct 3 ms 14684 KB Output is correct
11 Correct 4 ms 14684 KB Output is correct
12 Correct 3 ms 14684 KB Output is correct
13 Correct 5 ms 14684 KB Output is correct
14 Correct 3 ms 14684 KB Output is correct
15 Correct 4 ms 14680 KB Output is correct
16 Correct 3 ms 14684 KB Output is correct
17 Correct 3 ms 14684 KB Output is correct
18 Correct 3 ms 14684 KB Output is correct
19 Correct 3 ms 14684 KB Output is correct
20 Correct 3 ms 14684 KB Output is correct
21 Correct 3 ms 14684 KB Output is correct
22 Correct 3 ms 14684 KB Output is correct
23 Correct 4 ms 14936 KB Output is correct
24 Correct 3 ms 14684 KB Output is correct
25 Correct 3 ms 14684 KB Output is correct
26 Correct 3 ms 14684 KB Output is correct
27 Correct 3 ms 14680 KB Output is correct
28 Correct 3 ms 14684 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 14684 KB Output is correct
2 Correct 567 ms 102948 KB Output is correct
3 Correct 629 ms 99548 KB Output is correct
4 Correct 663 ms 115364 KB Output is correct
5 Correct 586 ms 101204 KB Output is correct
6 Correct 405 ms 100440 KB Output is correct
7 Correct 412 ms 100172 KB Output is correct
8 Correct 288 ms 98812 KB Output is correct
9 Correct 643 ms 123144 KB Output is correct
10 Correct 737 ms 118740 KB Output is correct
11 Correct 570 ms 103112 KB Output is correct
12 Correct 653 ms 99664 KB Output is correct
13 Correct 387 ms 97736 KB Output is correct
14 Correct 465 ms 97688 KB Output is correct
15 Correct 399 ms 97608 KB Output is correct
16 Correct 381 ms 97632 KB Output is correct
17 Correct 438 ms 98388 KB Output is correct
18 Correct 4 ms 14684 KB Output is correct
19 Correct 3 ms 14684 KB Output is correct
20 Correct 3 ms 14712 KB Output is correct
21 Correct 3 ms 14684 KB Output is correct
22 Correct 3 ms 14684 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 14684 KB Output is correct
2 Correct 567 ms 102948 KB Output is correct
3 Correct 629 ms 99548 KB Output is correct
4 Correct 663 ms 115364 KB Output is correct
5 Correct 586 ms 101204 KB Output is correct
6 Correct 405 ms 100440 KB Output is correct
7 Correct 412 ms 100172 KB Output is correct
8 Correct 288 ms 98812 KB Output is correct
9 Correct 643 ms 123144 KB Output is correct
10 Correct 737 ms 118740 KB Output is correct
11 Correct 570 ms 103112 KB Output is correct
12 Correct 653 ms 99664 KB Output is correct
13 Correct 387 ms 97736 KB Output is correct
14 Correct 465 ms 97688 KB Output is correct
15 Correct 399 ms 97608 KB Output is correct
16 Correct 381 ms 97632 KB Output is correct
17 Correct 438 ms 98388 KB Output is correct
18 Correct 4 ms 14684 KB Output is correct
19 Correct 3 ms 14684 KB Output is correct
20 Correct 3 ms 14712 KB Output is correct
21 Correct 3 ms 14684 KB Output is correct
22 Correct 3 ms 14684 KB Output is correct
23 Correct 3 ms 14680 KB Output is correct
24 Correct 599 ms 102840 KB Output is correct
25 Correct 751 ms 100040 KB Output is correct
26 Correct 690 ms 120740 KB Output is correct
27 Correct 614 ms 101200 KB Output is correct
28 Correct 426 ms 100592 KB Output is correct
29 Correct 414 ms 99924 KB Output is correct
30 Correct 309 ms 98628 KB Output is correct
31 Correct 674 ms 117412 KB Output is correct
32 Correct 784 ms 118580 KB Output is correct
33 Correct 596 ms 103032 KB Output is correct
34 Correct 720 ms 99428 KB Output is correct
35 Correct 3 ms 14680 KB Output is correct
36 Correct 3 ms 14684 KB Output is correct
37 Correct 3 ms 14684 KB Output is correct
38 Correct 3 ms 14684 KB Output is correct
39 Correct 3 ms 14684 KB Output is correct
40 Correct 3 ms 14684 KB Output is correct
41 Correct 3 ms 14684 KB Output is correct
42 Correct 3 ms 14684 KB Output is correct
43 Correct 3 ms 14680 KB Output is correct
44 Correct 4 ms 14684 KB Output is correct
45 Correct 3 ms 14684 KB Output is correct
46 Correct 4 ms 14684 KB Output is correct
47 Correct 3 ms 14684 KB Output is correct
48 Correct 3 ms 14680 KB Output is correct
49 Correct 3 ms 14680 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 14684 KB Output is correct
2 Correct 905 ms 120312 KB Output is correct
3 Correct 2585 ms 115440 KB Output is correct
4 Correct 1189 ms 116700 KB Output is correct
5 Correct 1418 ms 99648 KB Output is correct
6 Correct 1400 ms 98916 KB Output is correct
7 Correct 1014 ms 98900 KB Output is correct
8 Correct 323 ms 97472 KB Output is correct
9 Correct 903 ms 115360 KB Output is correct
10 Correct 2573 ms 118160 KB Output is correct
11 Correct 839 ms 99920 KB Output is correct
12 Correct 3209 ms 99908 KB Output is correct
13 Execution timed out 4057 ms 97672 KB Time limit exceeded
14 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 3 ms 14680 KB Output is correct
2 Correct 870 ms 116760 KB Output is correct
3 Correct 2617 ms 111328 KB Output is correct
4 Correct 1245 ms 114436 KB Output is correct
5 Correct 1295 ms 101048 KB Output is correct
6 Correct 1317 ms 100520 KB Output is correct
7 Correct 1156 ms 100168 KB Output is correct
8 Correct 319 ms 98628 KB Output is correct
9 Correct 872 ms 121816 KB Output is correct
10 Correct 2546 ms 119076 KB Output is correct
11 Correct 897 ms 102608 KB Output is correct
12 Correct 3335 ms 99780 KB Output is correct
13 Execution timed out 4097 ms 97580 KB Time limit exceeded
14 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 4 ms 14684 KB Output is correct
2 Correct 3 ms 14832 KB Output is correct
3 Correct 3 ms 14684 KB Output is correct
4 Correct 4 ms 15196 KB Output is correct
5 Correct 5 ms 15196 KB Output is correct
6 Correct 6 ms 15196 KB Output is correct
7 Correct 5 ms 15192 KB Output is correct
8 Correct 4 ms 15192 KB Output is correct
9 Correct 3 ms 14684 KB Output is correct
10 Correct 3 ms 14684 KB Output is correct
11 Correct 4 ms 14684 KB Output is correct
12 Correct 3 ms 14684 KB Output is correct
13 Correct 5 ms 14684 KB Output is correct
14 Correct 3 ms 14684 KB Output is correct
15 Correct 4 ms 14680 KB Output is correct
16 Correct 3 ms 14684 KB Output is correct
17 Correct 3 ms 14684 KB Output is correct
18 Correct 3 ms 14684 KB Output is correct
19 Correct 3 ms 14684 KB Output is correct
20 Correct 3 ms 14684 KB Output is correct
21 Correct 3 ms 14684 KB Output is correct
22 Correct 3 ms 14684 KB Output is correct
23 Correct 4 ms 14936 KB Output is correct
24 Correct 3 ms 14684 KB Output is correct
25 Correct 3 ms 14684 KB Output is correct
26 Correct 3 ms 14684 KB Output is correct
27 Correct 3 ms 14680 KB Output is correct
28 Correct 3 ms 14684 KB Output is correct
29 Correct 3 ms 14684 KB Output is correct
30 Correct 567 ms 102948 KB Output is correct
31 Correct 629 ms 99548 KB Output is correct
32 Correct 663 ms 115364 KB Output is correct
33 Correct 586 ms 101204 KB Output is correct
34 Correct 405 ms 100440 KB Output is correct
35 Correct 412 ms 100172 KB Output is correct
36 Correct 288 ms 98812 KB Output is correct
37 Correct 643 ms 123144 KB Output is correct
38 Correct 737 ms 118740 KB Output is correct
39 Correct 570 ms 103112 KB Output is correct
40 Correct 653 ms 99664 KB Output is correct
41 Correct 387 ms 97736 KB Output is correct
42 Correct 465 ms 97688 KB Output is correct
43 Correct 399 ms 97608 KB Output is correct
44 Correct 381 ms 97632 KB Output is correct
45 Correct 438 ms 98388 KB Output is correct
46 Correct 4 ms 14684 KB Output is correct
47 Correct 3 ms 14684 KB Output is correct
48 Correct 3 ms 14712 KB Output is correct
49 Correct 3 ms 14684 KB Output is correct
50 Correct 3 ms 14684 KB Output is correct
51 Correct 3 ms 14680 KB Output is correct
52 Correct 599 ms 102840 KB Output is correct
53 Correct 751 ms 100040 KB Output is correct
54 Correct 690 ms 120740 KB Output is correct
55 Correct 614 ms 101200 KB Output is correct
56 Correct 426 ms 100592 KB Output is correct
57 Correct 414 ms 99924 KB Output is correct
58 Correct 309 ms 98628 KB Output is correct
59 Correct 674 ms 117412 KB Output is correct
60 Correct 784 ms 118580 KB Output is correct
61 Correct 596 ms 103032 KB Output is correct
62 Correct 720 ms 99428 KB Output is correct
63 Correct 3 ms 14680 KB Output is correct
64 Correct 3 ms 14684 KB Output is correct
65 Correct 3 ms 14684 KB Output is correct
66 Correct 3 ms 14684 KB Output is correct
67 Correct 3 ms 14684 KB Output is correct
68 Correct 3 ms 14684 KB Output is correct
69 Correct 3 ms 14684 KB Output is correct
70 Correct 3 ms 14684 KB Output is correct
71 Correct 3 ms 14680 KB Output is correct
72 Correct 4 ms 14684 KB Output is correct
73 Correct 3 ms 14684 KB Output is correct
74 Correct 4 ms 14684 KB Output is correct
75 Correct 3 ms 14684 KB Output is correct
76 Correct 3 ms 14680 KB Output is correct
77 Correct 3 ms 14680 KB Output is correct
78 Correct 3 ms 14684 KB Output is correct
79 Correct 905 ms 120312 KB Output is correct
80 Correct 2585 ms 115440 KB Output is correct
81 Correct 1189 ms 116700 KB Output is correct
82 Correct 1418 ms 99648 KB Output is correct
83 Correct 1400 ms 98916 KB Output is correct
84 Correct 1014 ms 98900 KB Output is correct
85 Correct 323 ms 97472 KB Output is correct
86 Correct 903 ms 115360 KB Output is correct
87 Correct 2573 ms 118160 KB Output is correct
88 Correct 839 ms 99920 KB Output is correct
89 Correct 3209 ms 99908 KB Output is correct
90 Execution timed out 4057 ms 97672 KB Time limit exceeded
91 Halted 0 ms 0 KB -