結果

問題 No.235 めぐるはめぐる (5)
コンテスト
ユーザー 沙耶花
提出日時 2021-11-18 22:12:20
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2,104 ms / 10,000 ms
コード長 3,031 bytes
コンパイル時間 5,477 ms
コンパイル使用メモリ 271,548 KB
最終ジャッジ日時 2025-01-25 19:23:13
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint1000000007;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000001

struct HLD{
	vector<int> sz,parent,depth,root,pos;
	vector<int> arr;
	HLD(vector<vector<int>> &E){
		sz.resize(E.size(),1);
		parent.resize(E.size(),0);
		depth.resize(E.size(),0);
		root.resize(E.size(),0);
		pos.resize(E.size(),0);
		
		dfs(0,-1,E);
		dfs2(0,-1,E,0);
	}
	
	void dfs(int now,int p,vector<vector<int>> &E){
		parent[now] = p;
		if(p==-1){
			depth[now] = 0;
		}
		else{
			depth[now] = depth[p]+1;
		}
		for(int i=0;i<E[now].size();i++){
			int to = E[now][i];
			if(to==p)continue;
			dfs(to,now,E);
			sz[now] += sz[to];
		}
	}
	
	void dfs2(int now,int p,vector<vector<int>> &E,int r){
		pos[now] = arr.size();
		arr.push_back(now);
		root[now] = r;
		int maxi = 0;
		int ind = -1;
		for(int i=0;i<E[now].size();i++){
			int to = E[now][i];
			if(to==p)continue;
			if(maxi<sz[to]){
				maxi = sz[to];
				ind = to;
			}
		}
		if(ind==-1)return;
		dfs2(ind,now,E,r);
		for(int i=0;i<E[now].size();i++){
			int to = E[now][i];
			if(to==p||to==ind)continue;
			dfs2(to,now,E,to);
		}
	}
	
	vector<pair<int,int>> query(int u,int v){
		vector<pair<int,int>> ret;
		int t = 0;
		while(root[u]!=root[v]){
			if(depth[root[u]] <= depth[root[v]]){
				ret.insert(ret.begin()+t,{pos[root[v]], pos[v]});
				v = parent[root[v]];
			}
			else{
				ret.insert(ret.begin()+t,{pos[u],pos[root[u]]});
				u = parent[root[u]];
				t++;
			}
		}
		ret.insert(ret.begin()+t,{pos[u],pos[v]});
		return ret;
	}
	
	int lca(int u,int v){
		for(;;v=parent[root[v]]){
			if(pos[u]>pos[v])swap(u,v);
			if(root[u]==root[v])return u;
		}
	}
	
	int get_distance(int u,int v){
		return depth[u] + depth[v] - 2 * depth[lca(u,v)];
	}
	
};

using P = pair<mint,mint>;

P op(P a,P b){
	a.first += b.first;
	a.second += b.second;
	return a;
}

P e(){
	return make_pair(0,0);
}

P mapping(mint a,P b){
	b.first += a * b.second;
	return b;
}

mint composition(mint a,mint b){
	return a+b;
}

mint id(){
	return 0;
}

int main(){
	
	int N;
	cin>>N;
	
	vector<int> S(N),C(N);
	rep(i,N)cin>>S[i];
	rep(i,N)cin>>C[i];
	
	vector<vector<int>> E(N);
	rep(i,N-1){
		int u,v;
		cin>>u>>v;
		
		u--;v--;
		E[u].push_back(v);
		E[v].push_back(u);
	}
	
	HLD H(E);
	
	lazy_segtree<P,op,e,mint,mapping,composition,id> seg(N);
	
	rep(i,N){
		seg.set(H.pos[i],make_pair(S[i],C[i]));
	}
	

	int Q;
	cin>>Q;
	
	rep(_,Q){
		int t;
		cin>>t;
		if(t==0){
			int X,Y,Z;
			cin>>X>>Y>>Z;
			
			X--;Y--;
			auto ret = H.query(X,Y);
			rep(i,ret.size()){
				int l = ret[i].first,r = ret[i].second;
				if(l>r)swap(l,r);
				seg.apply(l,r+1,Z);
			}
		}
		else{
			int X,Y;
			cin>>X>>Y;
			
			X--;Y--;
			auto ret = H.query(X,Y);
			mint ans = 0;
			rep(i,ret.size()){
				int l = ret[i].first,r = ret[i].second;
				if(l>r)swap(l,r);
				ans += seg.prod(l,r+1).first;
				//cout<<l<<','<<r<<endl;
			}
			cout<<ans.val()<<endl;
		}
	}
		
	
	return 0;
}
0