結果

問題 No.235 めぐるはめぐる (5)
コンテスト
ユーザー chocorusk
提出日時 2019-10-04 14:09:15
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,436 ms / 10,000 ms
コード長 5,431 bytes
コンパイル時間 1,721 ms
コンパイル使用メモリ 133,152 KB
実行使用メモリ 37,284 KB
最終ジャッジ日時 2024-10-03 06:52:50
合計ジャッジ時間 7,977 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
template<typename Monoid, typename OperatorMonoid=Monoid>
struct LazySegmentTree{
	using F=function<Monoid(Monoid, Monoid)>;
	using G=function<Monoid(Monoid, OperatorMonoid, int, int)>;
	using H=function<OperatorMonoid(OperatorMonoid, OperatorMonoid)>;
	int sz;
	vector<Monoid> data;
	vector<OperatorMonoid> lazy;
	const F f;
	const G g;
	const H h;
	const Monoid e1;
	const OperatorMonoid e0;

	LazySegmentTree(int n, const F f, const G g, const H h, const Monoid &e1, const OperatorMonoid &e0): f(f), g(g), h(h), e1(e1), e0(e0){
		sz=1;
		while(sz<n) sz<<=1;
		data.resize(2*sz-1, e1);
		lazy.resize(2*sz-1, e0);
	}

	void build(vector<Monoid> v){
		for(int i=0; i<v.size(); i++) data[i+sz-1]=v[i];
		for(int i=sz-2; i>=0; i--) data[i]=f(data[2*i+1], data[2*i+2]);
	}

	void eval(int k, int l, int r){
		if(lazy[k]!=e0){
			data[k]=g(data[k], lazy[k], l, r);
			if(k<sz-1){
				lazy[2*k+1]=h(lazy[2*k+1], lazy[k]);
				lazy[2*k+2]=h(lazy[2*k+2], lazy[k]);
			}
		}
		lazy[k]=e0;
	}

	void update(int a, int b, const OperatorMonoid &x, int k, int l, int r){
		eval(k, l, r);
		if(r<=a || b<=l) return;
		if(a<=l && r<=b){
			lazy[k]=h(lazy[k], x);
			eval(k, l, r);
		}else{
			update(a, b, x, 2*k+1, l, (l+r)/2);
			update(a, b, x, 2*k+2, (l+r)/2, r);
			data[k]=f(data[2*k+1], data[2*k+2]);
		}
	}
	void update(int a, int b, const OperatorMonoid &x){
		return update(a, b, x, 0, 0, sz);
	}

	Monoid find(int a, int b, int k, int l, int r){
		eval(k, l, r);
		if(b<=l || r<=a) return e1;
		if(a<=l && r<=b) return data[k];
		else return f(find(a, b, 2*k+1, l, (l+r)/2), find(a, b, 2*k+2, (l+r)/2, r));
	}
	Monoid find(int a, int b){
		return find(a, b, 0, 0, sz);
	}
	Monoid operator[](const int &k){
		return find(k, k+1);
	}
};

template< typename G=vector<vector<int>> >
struct HeavyLightDecomposition {
  G &g;
  vector< int > sz, in, out, head, rev, par;

  HeavyLightDecomposition(G &g) :
      g(g), sz(g.size()), in(g.size()), out(g.size()), head(g.size()), rev(g.size()), par(g.size()) {}

  void dfs_sz(int idx, int p) {
    par[idx] = p;
    sz[idx] = 1;
    if(g[idx].size() && g[idx][0] == p) swap(g[idx][0], g[idx].back());
    for(auto &to : g[idx]) {
      if(to == p) continue;
      dfs_sz(to, idx);
      sz[idx] += sz[to];
      if(sz[g[idx][0]] < sz[to]) swap(g[idx][0], to);
    }
  }

  void dfs_hld(int idx, int par, int &times) {
    in[idx] = times++;
    rev[in[idx]] = idx;
    for(auto &to : g[idx]) {
      if(to == par) continue;
      head[to] = (g[idx][0] == to ? head[idx] : to);
      dfs_hld(to, idx, times);
    }
    out[idx] = times;
  }

  void build() {
    dfs_sz(0, -1);
    int t = 0;
    dfs_hld(0, -1, t);
  }

  int la(int v, int k) {
    while(1) {
      int u = head[v];
      if(in[v] - k >= in[u]) return rev[in[v] - k];
      k -= in[v] - in[u] + 1;
      v = par[u];
    }
  }

  int lca(int u, int v) {
    for(;; v = par[head[v]]) {
      if(in[u] > in[v]) swap(u, v);
      if(head[u] == head[v]) return u;
    }
  }

  template< typename T, typename Q, typename F >
  T query(int u, int v, const T &ti, const Q &q, const F &f, bool edge = false) {
    T l = ti;//, r = ti;
    for(;; v = par[head[v]]) {
      if(in[u] > in[v]) swap(u, v);//, swap(l, r);
      if(head[u] == head[v]) break;
      l = f(q(in[head[v]], in[v] + 1), l);
    }
    return f(q(in[u] + edge, in[v] + 1), l);
  }

  template< typename Q >
  void add(int u, int v, const Q &q, bool edge = false) {
    for(;; v = par[head[v]]) {
      if(in[u] > in[v]) swap(u, v);
      if(head[u] == head[v]) break;
      q(in[head[v]], in[v] + 1);
    }
    q(in[u] + edge, in[v] + 1);
  }
};

int main()
{
	int n;
	scanf("%d", &n);
	ll c1[200020], c[200020], cs[200020], s1[200020];
	vector<ll> s(n);
	for(int i=0; i<n; i++){
		scanf("%lld", &s1[i]);
	}
	const ll MOD=1e9+7;
	for(int i=0; i<n; i++){
		scanf("%lld", &c1[i]);
	}
	vector<vector<int>> G(n);
	for(int i=0; i<n-1; i++){
		int a, b;
		scanf("%d %d", &a, &b);
		a--; b--;
		G[a].push_back(b);
		G[b].push_back(a);
	}
	HeavyLightDecomposition<> hld(G);
	hld.build();
	for(int i=0; i<n; i++){
		s[hld.in[i]]=s1[i];
		c[hld.in[i]]=c1[i];
	}
	cs[0]=0;
	for(int i=0; i<n; i++){
		cs[i+1]=cs[i]+c[i];
		if(cs[i+1]>=MOD) cs[i+1]-=MOD;
	}
	auto f=[MOD](ll a, ll b){
		ll ret=a+b;
		if(ret>=MOD) ret-=MOD;
		return ret;
	};
	auto g=[&](ll a, ll x, int l, int r){
		return (a+x*(cs[r]-cs[l]+MOD))%MOD;
	};
	auto h=[MOD](ll x, ll y){
		ll ret=x+y;
		if(ret>=MOD) ret-=MOD;
		return ret;
	};
	LazySegmentTree<ll> seg(n, f, g, h, 0, 0);
	seg.build(s);
	int Q;
	scanf("%d", &Q);
	for(int z=0; z<Q; z++){
		int t;
		scanf("%d", &t);
		if(t==0){
			int x, y; ll z;
			scanf("%d %d %lld", &x, &y, &z);
			x--; y--;
			auto q=[&](int i, int j){ return seg.update(i, j, z);};
			hld.add(x, y, q);
		}else{
			int x, y;
			scanf("%d %d", &x, &y);
			x--; y--;
			auto q=[&](int i, int j){ return seg.find(i, j);};
			printf("%lld\n", hld.query(x, y, 0ll, q, f));
		}
	}
	return 0;
}
0