結果
問題 | No.235 めぐるはめぐる (5) |
ユーザー | fumofumofuni |
提出日時 | 2021-04-03 14:46:03 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2,416 ms / 10,000 ms |
コード長 | 7,418 bytes |
コンパイル時間 | 2,887 ms |
コンパイル使用メモリ | 221,484 KB |
実行使用メモリ | 42,148 KB |
最終ジャッジ日時 | 2023-08-26 05:36:05 |
合計ジャッジ時間 | 12,988 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2,416 ms
41,480 KB |
testcase_01 | AC | 1,411 ms
42,148 KB |
testcase_02 | AC | 2,204 ms
41,448 KB |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,n) for(ll i=0;i<n;i++) #define repl(i,l,r) for(ll i=(l);i<(r);i++) #define per(i,n) for(ll i=(n)-1;i>=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define pb push_back #define ins insert #define pqueue(x) priority_queue<x,vector<x>,greater<x>> #define all(x) (x).begin(),(x).end() #define CST(x) cout<<fixed<<setprecision(x) #define rev(x) reverse(x); using ll=long long; using vl=vector<ll>; using vvl=vector<vector<ll>>; using pl=pair<ll,ll>; using vpl=vector<pl>; using vvpl=vector<vpl>; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; const ll dy[8]={-1,0,1,0,1,1,-1,-1}; const ll dx[8]={0,-1,0,1,1,-1,1,-1}; template <typename T> inline bool chmax(T &a, T b) { return ((a < b) ? (a = b, true) : (false)); } template <typename T> inline bool chmin(T &a, T b) { return ((a > b) ? (a = b, true) : (false)); } template <typename X, typename M> struct SegTreeLazy {//遅延セグ木 単位元に注意(updateなら選ばれない数、affineなら(1,0)) using FX = function<X(X, X)>; using FA = function<X(X, M, int)>; using FM = function<M(M, M)>; int n; FX fx; FA fa; FM fm; const X ex; const M em; vector<X> dat; vector<M> lazy; SegTreeLazy(int n_, FX fx_, FA fa_, FM fm_, X ex_, M em_) : n(), fx(fx_), fa(fa_), fm(fm_), ex(ex_), em(em_), dat(n_ * 4, ex), lazy(n_ * 4, em) { int x = 1; while (n_ > x) x *= 2; n = x; } void set(int i, X x) { dat[i + n - 1] = x; } void build() { for (int k = n - 2; k >= 0; k--) dat[k] = fx(dat[2 * k + 1], dat[2 * k + 2]); } /* lazy eval */ void eval(int k, int len) { if (lazy[k] == em) return; // 更新するものが無ければ終了 if (k < n - 1) { // 葉でなければ子に伝搬 lazy[k * 2 + 1] = fm(lazy[k * 2 + 1], lazy[k]); lazy[k * 2 + 2] = fm(lazy[k * 2 + 2], lazy[k]); } // 自身を更新 dat[k] = fa(dat[k],lazy[k],len);//fa(dat[k], fp(lazy[k], len)); lazy[k] = em; } void update(int a, int b, M x, int k, int l, int r) { eval(k, r - l); if (a <= l && r <= b) { // 完全に内側の時 lazy[k] = fm(lazy[k], x); eval(k, r - l); } else if (a < r && l < b) { // 一部区間が被る時 update(a, b, x, k * 2 + 1, l, (l + r) / 2); // 左の子 update(a, b, x, k * 2 + 2, (l + r) / 2, r); // 右の子 dat[k] = fx(dat[k * 2 + 1], dat[k * 2 + 2]); } } void update(int a, int b, M x) { update(a, b, x, 0, 0, n); } X query_sub(int a, int b, int k, int l, int r) { eval(k, r - l); if (r <= a || b <= l) { // 完全に外側の時 return ex; } else if (a <= l && r <= b) { // 完全に内側の時 return dat[k]; } else { // 一部区間が被る時 X vl = query_sub(a, b, k * 2 + 1, l, (l + r) / 2); X vr = query_sub(a, b, k * 2 + 2, (l + r) / 2, r); return fx(vl, vr); } } X query(int a, int b) { return query_sub(a, b, 0, 0, n); } X operator[](int i){ return query(i,i+1); } }; struct HLD { int n,pos; vector<vector<int> > G; vector<int> vid, head, sub, par, dep, inv, type; HLD(){} HLD(int n): n(n),pos(0),G(n),vid(n,-1),head(n),sub(n,1), par(n,-1),dep(n,0),inv(n),type(n){} void add_edge(int u, int v) { G[u].push_back(v); G[v].push_back(u); } void build(vector<int> rs={0}) { int c=0; for(int r:rs){ dfs_sz(r); head[r]=r; dfs_hld(r,c++); } } void dfs_sz(int v) { for(int &u:G[v]){ if(u==par[v]) continue; par[u]=v; dep[u]=dep[v]+1; dfs_sz(u); sub[v]+=sub[u]; if(sub[u]>sub[G[v][0]]) swap(u,G[v][0]); } } void dfs_hld(int v,int c) { vid[v]=pos++; inv[vid[v]]=v; type[v]=c; for(int u:G[v]){ if(u==par[v]) continue; head[u]=(u==G[v][0]?head[v]:u); dfs_hld(u,c); } } // for_each(vertex) // [l,r] <- attention!! template<typename F> void for_each(int u, int v, const F& f) { while(1){ if(vid[u]>vid[v]) swap(u,v); f(max(vid[head[v]],vid[u]),vid[v]); if(head[u]!=head[v]) v=par[head[v]]; else break; } } template<typename T,typename Q,typename F> T for_each(int u,int v,T ti,const Q &q,const F &f){ T l=ti,r=ti; while(1){ if(vid[u]>vid[v]){ swap(u,v); swap(l,r); } l=f(l,q(max(vid[head[v]],vid[u]),vid[v])); if(head[u]!=head[v]) v=par[head[v]]; else break; } return f(l,r); } // for_each(edge) // [l,r] <- attention!! template<typename F> void for_each_edge(int u, int v,const F& f) { while(1){ if(vid[u]>vid[v]) swap(u,v); if(head[u]!=head[v]){ f(vid[head[v]],vid[v]); v=par[head[v]]; }else{ if(u!=v) f(vid[u]+1,vid[v]); break; } } } int lca(int u,int v){ while(1){ if(vid[u]>vid[v]) swap(u,v); if(head[u]==head[v]) return u; v=par[head[v]]; } } int distance(int u,int v){ return dep[u]+dep[v]-2*dep[lca(u,v)]; } }; const int mod = 1000000007; const int max_n = 200005; struct mint { ll x; // typedef long long ll; mint(ll x=0):x((x%mod+mod)%mod){} mint operator-() const { return mint(-x);} mint& operator+=(const mint a) { if ((x += a.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint a) { if ((x += mod-a.x) >= mod) x -= mod; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= mod; return *this;} mint operator+(const mint a) const { return mint(*this) += a;} mint operator-(const mint a) const { return mint(*this) -= a;} mint operator*(const mint a) const { return mint(*this) *= a;} mint pow(ll t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } bool operator==(const mint &p) const { return x == p.x; } bool operator!=(const mint &p) const { return x != p.x; } // for prime mod mint inv() const { return pow(mod-2);} mint& operator/=(const mint a) { return *this *= a.inv();} mint operator/(const mint a) const { return mint(*this) /= a;} }; istream& operator>>(istream& is, mint& a) { return is >> a.x;} ostream& operator<<(ostream& os, const mint& a) { return os << a.x;} int main(){ ll n;cin >> n; auto fx=[](pair<mint,mint> a,pair<mint,mint> b){return make_pair(a.fi+b.fi,a.se+b.se);}; auto fm=[](pair<mint,mint> a,mint b,int len){return make_pair(a.fi+a.se*b,a.se);}; auto fa=[](mint a,mint b){return a+b;}; SegTreeLazy<pair<mint,mint>,mint> st(n,fx,fm,fa,{0,0},0); vl s(n),c(n);rep(i,n)cin >> s[i];rep(i,n)cin >> c[i]; HLD hld(n); rep(i,n-1){ ll a,b;cin >> a >> b;a--;b--; hld.add_edge(a,b); } hld.build(); rep(i,n)st.set(hld.vid[i],{s[i],c[i]}); st.build(); ll q;cin >> q; while(q--){ ll t;cin >> t; if(!t){ ll v,w,z;cin >> v >> w >> z;v--;w--; hld.for_each(v,w,[&](int l,int r){st.update(l,r+1,z);}); } else{ mint ans=0; ll v,w;cin >> v >> w;v--;w--; hld.for_each(v,w,[&](int l,int r){ans+=st.query(l,r+1).fi;}); cout << ans <<endl; } } }