oj.uz

Submission #809235

# Submission time Handle Problem Language Result Execution time Memory
809235 2023-08-05T23:06:28 Z ashcomp Sumtree (INOI20_sumtree) C++14
50 / 100
 820 ms 275292 KB 
#include<bits/stdc++.h>

using namespace std;
 
#define all(x) x.begin() , x.end()
#define pii pair<int, int>
#define piii pair<int, pii>
#define pll pair<ll, ll>
#define plll pair<ll, pll>
#define pb push_back
#define qb pop_back
#define F first
#define S second
#define wall cerr<<"--------------------------------------"<<endl
 
mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#pragma GCC optimize("-O2")
 
typedef long long ll;
typedef double db;

const ll INF = 1e9;
const ll maxn = 3e5 + 10;
const ll N = 2e6 + 10;
const ll mod = 1e9 + 7;

ll n , q , k , sz[maxn] , st[maxn] , en[maxn] , h[maxn] , ti=1 , ans;
vector<ll> g[maxn] , ANS; 

ll fac[N] , inv[N];
ll pw(ll a , ll b)
{
        ll res = 1;
        while(b){
                if(b&1)res = (res*a)%mod;
                b>>=1; a=(a*a)%mod;
        }
        return res;
}
ll choose(ll a , ll b)
{
        if(b>a || b<0)return 0;

        ll A = fac[a];
        ll B = (inv[b] * inv[a-b])%mod;
        ll C = (A*B)%mod;
        return C;
}
void preps()
{
        fac[0] = 1;
        for(ll i=1; i<N; i++){
                fac[i] = (fac[i-1] * i)%mod;
        }
        inv[N-1] = pw(fac[N-1],mod-2);
        for(ll i=N-1; i>0; i--){
                inv[i-1] = (inv[i]*i)%mod;
        }
        for(ll i=1; i<maxn; i++){
                sz[i] = 1;
        }
}//ok

void dfs(ll v , ll p)
{
        st[v] = ti++;
        en[v] = st[v];
        for(auto u : g[v]){
                if(u==p)continue;
                h[u] = h[v] + 1;
                dfs(u,v);
                sz[v] += sz[u];
                en[v] = max(en[v],en[u]);
        }
}//ok

struct PAR_segment{
        set<pll> seg[maxn<<2];
        void build(ll tl=1 , ll tr=n , ll v=1)
        {
                seg[v].insert({INF,INF});
                if(tl==tr)return;
                ll mid=(tl+tr)>>1 , ml=(v<<1) , mr=(v<<1)|1;
                build(tl,mid,ml);
                build(mid+1,tr,mr);
        }
        void upd(pll val , ll l , ll r , ll tl=1 , ll tr=n , ll v=1)
        {
                if(l>r)return;
                if(l==tl && r==tr){
                        seg[v].insert(val);
                        return;
                }
                ll mid=(tl+tr)>>1 , ml=(v<<1) , mr=(v<<1)|1;
                upd(val,l,min(r,mid),tl,mid,ml);
                upd(val,max(l,mid+1),r,mid+1,tr,mr);
        }
        void del(pll val , ll l , ll r , ll tl=1 , ll tr=n , ll v=1)
        {
                if(l>r)return;
                if(l==tl && r==tr){
                        seg[v].erase(val);
                        return;
                }
                ll mid=(tl+tr)>>1 , ml=(v<<1) , mr=(v<<1)|1;
                del(val,l,min(r,mid),tl,mid,ml);
                del(val,max(l,mid+1),r,mid+1,tr,mr);
        }
        pll ask(ll pos , ll tl=1 , ll tr=n , ll v=1){
                if(tl==tr){
                        pll res = *seg[v].begin();
                        return res;
                }
                ll mid=(tl+tr)>>1 , ml=(v<<1) , mr=(v<<1)|1;
                pll res = *seg[v].begin() , e;
                if(pos<=mid)e = ask(pos,tl,mid,ml);
                else e = ask(pos,mid+1,tr,mr);
                if(e.F < res.F){
                        return e;
                }else{
                        return res;
                }
        }
};

struct SZ_segment{
        ll fen[maxn];
        void upd(ll pos , ll val)
        {
                for(; pos<=n; pos+=(pos&-pos)){
                        fen[pos]+=val;
                }
        }
        ll find(ll pos)
        {
                ll res = 0;
                for(; pos>0; pos-=(pos&-pos)){
                        res += fen[pos];
                }
                return res;
        }
        ll ask(ll l , ll r)
        {
                return find(r)-find(l-1);
        }

};//ok

PAR_segment Par;
SZ_segment Sz , Val;
ll cnt[maxn] , temp[maxn];

int main()
{
        ios_base::sync_with_stdio(0);cin.tie(NULL);cout.tie(NULL);
        preps();
        cin>>n>>k;
        for(int i=1; i<n; i++){
                ll v,u; cin>>v>>u;
                g[v].push_back(u);
                g[u].push_back(v);
        }
        dfs(1,1);
        /**/
        ans = choose(k+sz[1]-1 , sz[1]-1);
        ANS.push_back(ans);
        Par.build();
        Par.upd({-h[1],1},1,n);
        Sz.upd(st[1],n);
        Val.upd(st[1],k);
        cnt[1] = n;
        temp[1] = k;
        /**/
        cin>>q;
        for(int i=1; i<=q; i++){
                int tp; cin>>tp;
                if(tp==1){
                        ll v,X; cin>>v>>X;
                        X -= Val.ask(st[v],en[v]);
                        Val.upd(st[v],X);
                        ll siz=sz[v] - Sz.ask(st[v],en[v]);
                        Sz.upd(st[v],siz);
                        pll e = Par.ask(st[v]);ll p = e.S;
                        Par.upd({-h[v],v},st[v],en[v]);
                        Val.upd(st[p],-X);
                        Sz.upd(st[p],-siz);
                        //cout<<p<<" , "<<siz<<" , "<<X<<"\n";
                        /**/
                        ll x = choose(temp[p]+cnt[p]-1 , cnt[p]-1);
                        ans = (ans * pw(x,mod-2))%mod;
                        cnt[p] -= siz;
                        cnt[v] += siz;
                        temp[p] -= X;
                        temp[v] += X;
                        x = choose(temp[p]+cnt[p]-1 , cnt[p]-1);
                        ans = (ans * x)%mod;
                        x = choose(temp[v]+cnt[v]-1 , cnt[v]-1);
                        ans = (ans * x)%mod;
                        ANS.push_back(ans);
                        /**/
                }else{
                        ll v; cin>>v;
                        Par.del({-h[v],v},st[v],en[v]);
                        pll e = Par.ask(st[v]);ll p = e.S;
                        ll siz = cnt[v];
                        Val.upd(st[v],-temp[v]);
                        Val.upd(st[p],temp[v]);
                        Sz.upd(st[v],-siz);
                        Sz.upd(st[p],siz);
                        /**/
                        ll x = choose(temp[p]+cnt[p]-1 , cnt[p]-1);
                        ll y = choose(temp[v]+cnt[v]-1 , cnt[v]-1);
                        x = (y*x)%mod;
                        ans = (ans * pw(x,mod-2))%mod;
                        cnt[p] += siz;
                        cnt[v] -= siz;
                        temp[p] += temp[v];
                        temp[v] -= temp[v];
                        x = choose(temp[p]+cnt[p]-1 , cnt[p]-1);
                        ans = (ans * x)%mod;
                        ANS.push_back(ans);
                        /**/
                }
        }
        for(auto u : ANS){
                cout<<u<<"\n";
        }
        return 0;
}
/*
4 3
1 2
2 3
2 4
2
1 3 1
1 2 1

*/
/*

*/
/*
by
          _____                    _____                    _____                    _____                  
         /\    \                  /\    \                  /\    \                  /\    \                 
        /::\    \                /::\    \                /::\____\                /::\    \                
       /::::\    \              /::::\    \              /:::/    /                \:::\    \               
      /::::::\    \            /::::::\    \            /:::/    /                  \:::\    \              
     /:::/\:::\    \          /:::/\:::\    \          /:::/    /                    \:::\    \             
    /:::/__\:::\    \        /:::/__\:::\    \        /:::/____/                      \:::\    \            
   /::::\   \:::\    \       \:::\   \:::\    \      /::::\    \                      /::::\    \           
  /::::::\   \:::\    \    ___\:::\   \:::\    \    /::::::\    \   _____    ____    /::::::\    \          
 /:::/\:::\   \:::\    \  /\   \:::\   \:::\    \  /:::/\:::\    \ /\    \  /\   \  /:::/\:::\    \         
/:::/  \:::\   \:::\____\/::\   \:::\   \:::\____\/:::/  \:::\    /::\____\/::\   \/:::/  \:::\____\        
\::/    \:::\  /:::/    /\:::\   \:::\   \::/    /\::/    \:::\  /:::/    /\:::\  /:::/    \::/    /        
 \/____/ \:::\/:::/    /  \:::\   \:::\   \/____/  \/____/ \:::\/:::/    /  \:::\/:::/    / \/____/         
          \::::::/    /    \:::\   \:::\    \               \::::::/    /    \::::::/    /                  
           \::::/    /      \:::\   \:::\____\               \::::/    /      \::::/____/                   
           /:::/    /        \:::\  /:::/    /               /:::/    /        \:::\    \                   
          /:::/    /          \:::\/:::/    /               /:::/    /          \:::\    \                  
         /:::/    /            \::::::/    /               /:::/    /            \:::\    \                 
        /:::/    /              \::::/    /               /:::/    /              \:::\____\                
        \::/    /                \::/    /                \::/    /                \::/    /                
         \/____/                  \/____/                  \/____/                  \/____/                 
                                                                                                            
*/
//NOGHRE SHODANAM BAD NIST :_)

Subtask #1
10.0 / 10.0

# Verdict Execution time Memory Grader output
1 Correct 168 ms 143924 KB Output is correct
2 Correct 214 ms 143968 KB Output is correct
3 Correct 178 ms 143884 KB Output is correct
4 Correct 197 ms 143960 KB Output is correct
5 Correct 184 ms 138196 KB Output is correct
6 Correct 64 ms 98200 KB Output is correct
7 Correct 64 ms 98008 KB Output is correct
8 Correct 63 ms 98028 KB Output is correct
9 Correct 174 ms 134692 KB Output is correct
10 Correct 204 ms 134520 KB Output is correct
11 Correct 175 ms 134684 KB Output is correct
12 Correct 162 ms 132788 KB Output is correct
13 Correct 156 ms 140136 KB Output is correct

Subtask #2
20.0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 61 ms 97388 KB Output is correct
2 Correct 65 ms 97308 KB Output is correct
3 Correct 61 ms 97308 KB Output is correct
4 Correct 61 ms 97356 KB Output is correct
5 Correct 61 ms 97492 KB Output is correct
6 Correct 63 ms 98004 KB Output is correct
7 Correct 64 ms 97964 KB Output is correct
8 Correct 64 ms 97916 KB Output is correct
9 Correct 65 ms 98216 KB Output is correct
10 Correct 65 ms 98288 KB Output is correct
11 Correct 67 ms 98476 KB Output is correct
12 Correct 63 ms 98236 KB Output is correct
13 Correct 73 ms 98284 KB Output is correct
14 Correct 67 ms 98336 KB Output is correct
15 Correct 66 ms 98508 KB Output is correct
16 Correct 65 ms 98124 KB Output is correct
17 Correct 66 ms 98388 KB Output is correct
18 Correct 64 ms 97908 KB Output is correct
19 Correct 64 ms 98236 KB Output is correct
20 Correct 65 ms 97960 KB Output is correct
21 Correct 69 ms 97884 KB Output is correct
22 Correct 62 ms 97532 KB Output is correct

Subtask #3
20.0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 211 ms 152936 KB Output is correct
2 Correct 252 ms 160336 KB Output is correct
3 Correct 263 ms 154156 KB Output is correct
4 Correct 374 ms 173092 KB Output is correct
5 Correct 809 ms 224640 KB Output is correct
6 Correct 68 ms 98732 KB Output is correct
7 Correct 65 ms 98176 KB Output is correct
8 Correct 66 ms 98508 KB Output is correct
9 Correct 358 ms 154700 KB Output is correct
10 Correct 352 ms 151460 KB Output is correct
11 Correct 321 ms 150172 KB Output is correct
12 Correct 381 ms 152052 KB Output is correct
13 Correct 820 ms 275164 KB Output is correct
14 Correct 705 ms 275156 KB Output is correct
15 Correct 697 ms 275228 KB Output is correct
16 Correct 706 ms 275156 KB Output is correct
17 Correct 697 ms 275292 KB Output is correct

Subtask #4
0 / 20.0

# Verdict Execution time Memory Grader output
1 Correct 485 ms 144724 KB Output is correct
2 Correct 448 ms 144420 KB Output is correct
3 Correct 441 ms 144596 KB Output is correct
4 Correct 453 ms 144664 KB Output is correct
5 Correct 419 ms 142160 KB Output is correct
6 Correct 441 ms 144624 KB Output is correct
7 Correct 325 ms 122848 KB Output is correct
8 Correct 335 ms 122964 KB Output is correct
9 Correct 446 ms 144656 KB Output is correct
10 Correct 435 ms 144596 KB Output is correct
11 Correct 435 ms 144504 KB Output is correct
12 Correct 333 ms 122788 KB Output is correct
13 Incorrect 232 ms 118804 KB Output isn't correct
14 Halted 0 ms 0 KB -

Subtask #5
0 / 30.0

# Verdict Execution time Memory Grader output
1 Correct 168 ms 143924 KB Output is correct
2 Correct 214 ms 143968 KB Output is correct
3 Correct 178 ms 143884 KB Output is correct
4 Correct 197 ms 143960 KB Output is correct
5 Correct 184 ms 138196 KB Output is correct
6 Correct 64 ms 98200 KB Output is correct
7 Correct 64 ms 98008 KB Output is correct
8 Correct 63 ms 98028 KB Output is correct
9 Correct 174 ms 134692 KB Output is correct
10 Correct 204 ms 134520 KB Output is correct
11 Correct 175 ms 134684 KB Output is correct
12 Correct 162 ms 132788 KB Output is correct
13 Correct 156 ms 140136 KB Output is correct
14 Correct 61 ms 97388 KB Output is correct
15 Correct 65 ms 97308 KB Output is correct
16 Correct 61 ms 97308 KB Output is correct
17 Correct 61 ms 97356 KB Output is correct
18 Correct 61 ms 97492 KB Output is correct
19 Correct 63 ms 98004 KB Output is correct
20 Correct 64 ms 97964 KB Output is correct
21 Correct 64 ms 97916 KB Output is correct
22 Correct 65 ms 98216 KB Output is correct
23 Correct 65 ms 98288 KB Output is correct
24 Correct 67 ms 98476 KB Output is correct
25 Correct 63 ms 98236 KB Output is correct
26 Correct 73 ms 98284 KB Output is correct
27 Correct 67 ms 98336 KB Output is correct
28 Correct 66 ms 98508 KB Output is correct
29 Correct 65 ms 98124 KB Output is correct
30 Correct 66 ms 98388 KB Output is correct
31 Correct 64 ms 97908 KB Output is correct
32 Correct 64 ms 98236 KB Output is correct
33 Correct 65 ms 97960 KB Output is correct
34 Correct 69 ms 97884 KB Output is correct
35 Correct 62 ms 97532 KB Output is correct
36 Correct 211 ms 152936 KB Output is correct
37 Correct 252 ms 160336 KB Output is correct
38 Correct 263 ms 154156 KB Output is correct
39 Correct 374 ms 173092 KB Output is correct
40 Correct 809 ms 224640 KB Output is correct
41 Correct 68 ms 98732 KB Output is correct
42 Correct 65 ms 98176 KB Output is correct
43 Correct 66 ms 98508 KB Output is correct
44 Correct 358 ms 154700 KB Output is correct
45 Correct 352 ms 151460 KB Output is correct
46 Correct 321 ms 150172 KB Output is correct
47 Correct 381 ms 152052 KB Output is correct
48 Correct 820 ms 275164 KB Output is correct
49 Correct 705 ms 275156 KB Output is correct
50 Correct 697 ms 275228 KB Output is correct
51 Correct 706 ms 275156 KB Output is correct
52 Correct 697 ms 275292 KB Output is correct
53 Correct 485 ms 144724 KB Output is correct
54 Correct 448 ms 144420 KB Output is correct
55 Correct 441 ms 144596 KB Output is correct
56 Correct 453 ms 144664 KB Output is correct
57 Correct 419 ms 142160 KB Output is correct
58 Correct 441 ms 144624 KB Output is correct
59 Correct 325 ms 122848 KB Output is correct
60 Correct 335 ms 122964 KB Output is correct
61 Correct 446 ms 144656 KB Output is correct
62 Correct 435 ms 144596 KB Output is correct
63 Correct 435 ms 144504 KB Output is correct
64 Correct 333 ms 122788 KB Output is correct
65 Incorrect 232 ms 118804 KB Output isn't correct
66 Halted 0 ms 0 KB -