제출 #966747

#제출 시각아이디문제언어결과실행 시간메모리
966747PM1Sprinkler (JOI22_sprinkler)C++17
53 / 100
4077 ms129884 KiB
#include <bits/stdc++.h> using namespace std; #pragma GCC optimize("O3,unroll-loops") #pragma GCC target("avx2") #define ll long long const int mxn=2e5+5; int n,M,q,h[mxn],st[mxn],fn[mxn],cnt,dis[mxn],par[mxn]; int fnd[mxn][42][2]; vector<int>v[mxn],d[mxn]; struct segment{ vector<ll>val; int sz=0; void build(){ while(__builtin_popcount(sz)!=1){ sz++; } val.resize(sz*2,1); } void add(int i,int l,int r,int tl,int tr,int w){ if(l>r||l>tr||r<tl||tl>tr){ return ; } if(l>=tl&&r<=tr){ val[i]*=w; val[i]%=M; return ; } int m=(l+r)>>1; add((i<<1),l,m,tl,tr,w); add((i<<1)^1,m+1,r,tl,tr,w); return ; } int get(int i){ if(i==0){ return 1; } int ret=val[i]*get((i>>1))%M; return ret; } }seg[mxn]; void dfs(int z){ st[z]=++cnt; d[dis[z]].push_back(st[z]); for(auto i:v[z]){ if(i!=par[z]){ par[i]=z; dis[i]=dis[z]+1; dfs(i); } } fn[z]=cnt; } void add(int x,int y,ll w){ if(d[x].size()==0)return; int l=fnd[y][x-dis[y]][0]; int r=fnd[y][x-dis[y]][1]; seg[x].add(1,0,seg[x].sz-1,l,r,w); } void up(int x,int dd,int w){ for(;dd>=1 && x!=1;dd--){ int z=dis[x]+dd,zz=dis[x]+dd-1; add(z,x,w); add(zz,x,w); x=par[x]; } for(;dd>=0;dd--){ add(dis[x]+dd,x,w); } } void get(int z,int x){ int l=fnd[x][0][0]; cout<<seg[z].get(l+seg[z].sz)<<'\n'; } int main(){ ios::sync_with_stdio(false); cin.tie(0); cout.tie(0); cin>>n>>M; for(int i=1;i<n;i++){ int x,y; cin>>x>>y; v[x].push_back(y); v[y].push_back(x); } for(int i=1;i<=n;i++) cin>>h[i]; dfs(1); for(int i=0;i<=n;i++){ seg[i].sz=d[i].size(); seg[i].build(); } for(int i=1;i<=n;i++){ fnd[i][0][1]=fnd[i][0][0]=lower_bound(d[dis[i]].begin(),d[dis[i]].end(),st[i])-d[dis[i]].begin(); for(int j=1;j<=40 && dis[i]+j<mxn;j++){ int dd=dis[i]+j; fnd[i][j][0]=lower_bound(d[dd].begin(),d[dd].end(),st[i])-d[dd].begin(); fnd[i][j][1]=upper_bound(d[dd].begin(),d[dd].end(),fn[i])-d[dd].begin()-1; } } for(int i=1;i<=n;i++){ int l=fnd[i][0][0]; seg[dis[i]].add(1,0,seg[dis[i]].sz-1,l,l,h[i]); } cin>>q; while(q--){ int ty,x,dd,w; cin>>ty; if(ty==1){ cin>>x>>dd>>w; up(x,dd,w); } else{ cin>>x; get(dis[x],x); } } return 0; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...