結果
問題 | No.899 γatheree |
ユーザー | t9unkubj |
提出日時 | 2024-06-01 07:38:14 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,015 ms / 2,000 ms |
コード長 | 9,835 bytes |
コンパイル時間 | 3,895 ms |
コンパイル使用メモリ | 248,776 KB |
実行使用メモリ | 27,424 KB |
最終ジャッジ日時 | 2024-06-01 07:38:37 |
合計ジャッジ時間 | 22,778 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,812 KB |
testcase_01 | AC | 3 ms
6,940 KB |
testcase_02 | AC | 2 ms
6,940 KB |
testcase_03 | AC | 2 ms
6,944 KB |
testcase_04 | AC | 2 ms
6,940 KB |
testcase_05 | AC | 2 ms
6,944 KB |
testcase_06 | AC | 991 ms
26,464 KB |
testcase_07 | AC | 924 ms
26,528 KB |
testcase_08 | AC | 907 ms
27,404 KB |
testcase_09 | AC | 940 ms
26,624 KB |
testcase_10 | AC | 928 ms
26,532 KB |
testcase_11 | AC | 901 ms
26,496 KB |
testcase_12 | AC | 898 ms
26,460 KB |
testcase_13 | AC | 937 ms
26,496 KB |
testcase_14 | AC | 917 ms
26,472 KB |
testcase_15 | AC | 910 ms
27,424 KB |
testcase_16 | AC | 893 ms
26,480 KB |
testcase_17 | AC | 1,015 ms
26,664 KB |
testcase_18 | AC | 924 ms
26,624 KB |
testcase_19 | AC | 909 ms
26,456 KB |
testcase_20 | AC | 928 ms
27,420 KB |
testcase_21 | AC | 685 ms
27,092 KB |
testcase_22 | AC | 725 ms
27,240 KB |
testcase_23 | AC | 721 ms
27,252 KB |
ソースコード
#ifdef t9unkubj #include"template.h" #else #undef _GLIBCXX_DEBUG #pragma GCC optimize("O3") #define dbg(...) 199958 using namespace std; #include<bits/stdc++.h> using uint=unsigned; using ll=long long; using ull=unsigned long long; using ld=long double; using pii=pair<int,int>; using pll=pair<ll,ll>; template<class T>using vc=vector<T>; template<class T>using vvc=vc<vc<T>>; template<class T>using vvvc=vvc<vc<T>>; using vi=vc<int>; using vvi=vc<vi>; using vvvi=vc<vvi>; using vl=vc<ll>; using vvl=vc<vl>; using vvvl=vc<vvl>; template<class T>using smpq=priority_queue<T,vector<T>,greater<T>>; template<class T>using bipq=priority_queue<T>; #define rep(i,n) for(ll (i)=0;i<(ll)(n);i++) #define REP(i,j,n) for(ll (i)=(j);i<(ll)(n);i++) #define DREP(i,n,m) for(ll (i)=(n);i>=(m);i--) #define drep(i,n) for(ll i=((n)-1);i>=0;i--) #define all(x) x.begin(),x.end() #define rall(x) x.rbegin(),x.rend() #define mp make_pair #define mt make_tuple #define pb push_back #define eb emplace_back #define fi first #define se second #define is insert #define bg begin() #define ed end() void scan(int&a) { cin >> a; } void scan(ll&a) { cin >> a; } void scan(string&a) { cin >> a; } void scan(char&a) { cin >> a; } void scan(uint&a) { cin >> a; } void scan(ull&a) { cin >> a; } void scan(bool&a) { cin >> a; } void scan(ld&a){ cin>> a;} template<class T> void scan(vector<T>&a) { for(auto&x:a) scan(x); } void read() {} template<class Head, class... Tail> void read(Head&head, Tail&... tail) { scan(head); read(tail...); } #define INT(...) int __VA_ARGS__; read(__VA_ARGS__); #define LL(...) ll __VA_ARGS__; read(__VA_ARGS__); #define ULL(...) ull __VA_ARGS__; read(__VA_ARGS__); #define STR(...) string __VA_ARGS__; read(__VA_ARGS__); #define CHR(...) char __VA_ARGS__; read(__VA_ARGS__); #define DBL(...) double __VA_ARGS__; read(__VA_ARGS__); #define LD(...) ld __VA_ARGS__; read(__VA_ARGS__); #define VC(type, name, ...) vector<type> name(__VA_ARGS__); read(name); #define VVC(type, name, size, ...) vector<vector<type>> name(size, vector<type>(__VA_ARGS__)); read(name); void print(int a) { cout << a; } void print(ll a) { cout << a; } void print(string a) { cout << a; } void print(char a) { cout << a; } void print(uint a) { cout << a; } void print(bool a) { cout << a; } void print(ull a) { cout << a; } void print(ld a){ cout<< a; } template<class T> void print(vector<T>a) { for(int i=0;i<(int)a.size();i++){if(i)cout<<" ";print(a[i]);}cout<<endl;} void PRT() { cout <<endl; return ; } template<class T> void PRT(T a) { print(a); cout <<endl; return; } template<class Head, class... Tail> void PRT(Head head, Tail ... tail) { print(head); cout << " "; PRT(tail...); return; } template<class T,class F> bool chmin(T &x, F y){ if(x>y){ x=y; return true; } return false; } template<class T, class F> bool chmax(T &x, F y){ if(x<y){ x=y; return true; } return false; } void YesNo(bool b){ cout<<(b?"Yes":"No")<<endl; } void Yes(){ cout<<"Yes"<<endl; } void No(){ cout<<"No"<<endl; } template<class T> int popcount(T n){ return __builtin_popcount(n); } template<class T> T sum(vc<T>&a){ return accumulate(all(a),T(0)); } template<class T> T max(vc<T>&a){ return *max_element(all(a)); } template<class T> T min(vc<T>&a){ return *min_element(all(a)); } template<class T> void unique(vc<T>&a){ a.erase(unique(all(a)),a.end()); } vvi readgraph(int n,int m,int off = -1){ vvi g(n); rep(i, m){ int u,v; cin>>u>>v; u+=off,v+=off; g[u].push_back(v); g[v].push_back(u); } return g; } vvi readtree(int n,int off=-1){ return readgraph(n,n-1,off); } template<class T> vc<T> presum(vc<T> &a){ vc<T> ret(a.size()+1); rep(i,a.size())ret[i+1]=ret[i]+a[i]; return ret; } template<class T, class F> vc<T> &operator+=(vc<T> &a,F b){ for (auto&v:a)v += b; return a; } template<class T, class F> vc<T> &operator-=(vc<T>&a,F b){ for (auto&v:a)v-=b; return a; } template<class T, class F> vc<T> &operator*=(vc<T>&a,F b){ for (auto&v:a)v*=b; return a; } #endif double pass_time=0; template<class T,class F> struct lazy_segtree{ int n; vector<T>data; vector<F>lazy; using TT=function<T(T,T)>; using TF=function<T(T,F)>; using FF=function<F(F,F)>; TT op; T e; TF mapping; FF merge; F id; int log; lazy_segtree(int n1,TT op,T e,TF mapping,FF merge,F id):n(n1),op(op),e(e),mapping(mapping),merge(merge),id(id){ int n_=1; log=0; while(n_<=n)n_*=2,log++; data=vector<T>(n_*2,e); lazy=vector<F>(n_,id); n=n_; for(int i=n_-1;i>0;i--)update(i); } void update(int k){ data[k]=op(data[k*2],data[k*2+1]); } void apply(int k,F x){// data[k]=mapping(data[k],x); if(k<n){ lazy[k]=merge(lazy[k],x); if(lazy[k]==id)return; //segment tree beats! if(data[k].fail)push(k),update(k); } } void push(int k){//上から遅延素を伝播させる apply(k*2,lazy[k]); apply(k*2+1,lazy[k]); lazy[k]=id; } void set(int p,T x){ p+=n; for(int i=log;i>=1;i--){ push(p>>i); } data[p]=x; for(int i=1;i<=log;i++){ update(p>>i); } } void apply(int l,int r,F x){ if(l==r)return; l+=n,r+=n; for(int i=log;i>=1;i--){ if(((l>>i)<<i)!=l)push(l>>i); if(((r>>i)<<i)!=r)push((r-1)>>i); } int l2{l},r2{r}; while(l<r){ if(l&1)apply(l++,x); if(r&1)apply(--r,x); l>>=1,r>>=1; } l=l2,r=r2; for(int i=1;i<=log;i++){ if(((l>>i)<<i)!=l)update(l>>i); if(((r>>i)<<i)!=r)update((r-1)>>i); } } T prod(int l,int r){ l+=n,r+=n; for(int i=log;i>=1;i--){ if(((l>>i)<<i)!=l)push(l>>i); if(((r>>i)<<i)!=r)push((r-1)>>i); } T sml=e,smr=e; while(l<r){ if(l&1)sml=op(sml,data[l++]); if(r&1)smr=op(data[--r],smr); l>>=1,r>>=1; } return op(sml,smr); } T get(int p){ p+=n; for(int i=log;i>=1;i--)push(p>>i); return data[p]; } }; struct bfs_euler_tour{ vector<vector<int>>g; vector<int>ord; vector<int>par; vector<vector<pair<int,int>>>segment; int dmax; bfs_euler_tour(vector<vector<int>>g,int dmax):g(g),ord(g.size(),-1),par(g.size(),-1),segment(g.size(),vector<pair<int,int>>(dmax+1)),dmax(dmax){ queue<int>que; que.push(0); int t=0; while(que.size()){ auto p=que.front();que.pop(); ord[p]=t++; for(auto x:g[p]){ if(ord[x]==-1){ par[x]=p; que.push(x); } } } for(int i=0;i<int(g.size());i++)segment[i][0]=make_pair(ord[i],ord[i]+1); for(int i=1;i<=dmax;i++){ for(int j=0;j<int(g.size());j++){ segment[j][i]=make_pair(int(1e9),-1); for(auto x:g[j]){ if(x==par[j])continue; segment[j][i].first=min(segment[j][i].first,segment[x][i-1].first); segment[j][i].second=max(segment[j][i].second,segment[x][i-1].second); } } } } vector<pair<int,int>>merge_segments(vector<pair<int,int>>&segments){ sort(segments.begin(),segments.end(),[](auto a,auto b){ return a.second<b.second; }); vector<pair<int,int>>res; pair<int,int>now{1e9,-1};//[a,b) for(auto [a,b]:segments){ if(now.second>=a){//[a,b) [c,d) now.first=min(now.first,a); now.second=max(now.second,b); }else{ if(now!=make_pair(int(1e9),-1))res.push_back(now); now=make_pair(a,b); } } res.push_back(now); return res; } pair<int,int>vertex(int u,int d){//下にuからの距離がdである頂点集合を返す(半開) if(u==-1)return make_pair(0,1); return segment[u][d]; } //上に距離がd,上に距離がd-1,..下に距離がdの順に返す(半開) vector<pair<int,int>>query(int u,int d){ assert(d<=dmax); vector<pair<int,int>>res; res.push_back(make_pair(ord[u],ord[u]+1)); auto u_=u; for(int d_=1;d_<=d;d_++){ if(u_!=-1)u_=par[u_]; res.push_back(vertex(u_,d-d_)); res.push_back(vertex(u_,0)); } reverse(res.begin(),res.end()); for(int d_=1;d_<=d;d_++){ res.push_back(vertex(u,d_)); } return merge_segments(res); } }; void solve(){ INT(n); auto g=readtree(n,0); VC(int,a,n); bfs_euler_tour t(g,2); auto op=[](ll a,ll b){return a+b;}; auto op2=[](ll a,ll b){if(b)return 0ll;return a;}; auto op3=[](ll a,ll b){return a|b;}; lazy_segtree<ll,ll>seg(n,op,0,op2,op3,0); rep(i,n)seg.set(t.ord[i],a[i]); INT(q); while(q--){ INT(x); ll sum=0; for(auto [a,b]:t.query(x,2)){ sum+=seg.prod(a,b); } for(auto [a,b]:t.query(x,2)){ seg.apply(a,b,1); } seg.set(t.ord[x],sum); PRT(sum); } } signed main(){ #ifdef t9unkubj freopen("input.txt", "r", stdin); freopen("output.txt", "w", stdout); #endif cin.tie(0)->sync_with_stdio(0); pass_time=clock(); int t=1; while(t--)solve(); pass_time=clock()-pass_time; dbg(pass_time/CLOCKS_PER_SEC); }