結果

問題 No.1030 だんしんぐぱーりない
ユーザー 👑 tute7627
提出日時 2020-04-17 22:26:44
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 188 ms / 2,000 ms
コード長 6,890 bytes
コンパイル時間 2,130 ms
コンパイル使用メモリ 186,784 KB
実行使用メモリ 46,784 KB
最終ジャッジ日時 2024-10-03 13:58:30
合計ジャッジ時間 8,120 ms
ジャッジサーバーID
(参考情報)
judge5 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 40
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

//#define _GLIBCXX_DEBUG
#include<bits/stdc++.h>
using namespace std;
#define endl '\n'
#define lfs cout<<fixed<<setprecision(10)
#define ALL(a) (a).begin(),(a).end()
#define ALLR(a) (a).rbegin(),(a).rend()
#define spa << " " <<
#define fi first
#define se second
#define MP make_pair
#define MT make_tuple
#define PB push_back
#define EB emplace_back
#define rep(i,n,m) for(ll i = (n); i < (ll)(m); i++)
#define rrep(i,n,m) for(ll i = (ll)(m) - 1; i >= (ll)(n); i--)
using ll = long long;
using ld = long double;
const ll MOD1 = 1e9+7;
const ll MOD9 = 998244353;
const ll INF = 1e18;
using P = pair<ll, ll>;
template<typename T1, typename T2>
bool chmin(T1 &a,T2 b){if(a>b){a=b;return true;}else return false;}
template<typename T1, typename T2>
bool chmax(T1 &a,T2 b){if(a<b){a=b;return true;}else return false;}
ll median(ll a,ll b, ll c){return a+b+c-max({a,b,c})-min({a,b,c});}
void ans1(bool x){if(x) cout<<"Yes"<<endl;else cout<<"No"<<endl;}
void ans2(bool x){if(x) cout<<"YES"<<endl;else cout<<"NO"<<endl;}
void ans3(bool x){if(x) cout<<"Yay!"<<endl;else cout<<":("<<endl;}
template<typename T1,typename T2>
void ans(bool x,T1 y,T2 z){if(x)cout<<y<<endl;else cout<<z<<endl;}
template<typename T>
void debug(vector<vector<T>>&v,ll h,ll w){for(ll i=0;i<h;i++)
{cout<<v[i][0];for(ll j=1;j<w;j++)cout spa v[i][j];cout<<endl;}};
void debug(vector<string>&v,ll h,ll w){for(ll i=0;i<h;i++)
{for(ll j=0;j<w;j++)cout<<v[i][j];cout<<endl;}};
template<typename T>
void debug(vector<T>&v,ll n){if(n!=0)cout<<v[0];
for(ll i=1;i<n;i++)cout spa v[i];cout<<endl;};
template<typename T>
vector<vector<T>>vec(ll x, ll y, T w){
vector<vector<T>>v(x,vector<T>(y,w));return v;}
ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;}
vector<ll>dx={1,0,-1,0,1,1,-1,-1};
vector<ll>dy={0,1,0,-1,1,-1,1,-1};
template<typename T>
vector<T> make_v(size_t a,T b){return vector<T>(a,b);}
template<typename... Ts>
auto make_v(size_t a,Ts... ts){
return vector<decltype(make_v(ts...))>(a,make_v(ts...));
}
template<typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2>&p){
return os << p.first << " " << p.second;
}
template< typename Monoid >
struct SegmentTree {
using F = function< Monoid(Monoid, Monoid) >;
int sz;
vector< Monoid > seg;
const F f;
const Monoid M1;
SegmentTree(int n, const F f, const Monoid &M1) : f(f), M1(M1) {
sz = 1;
while(sz < n) sz <<= 1;
seg.assign(2 * sz, M1);
}
void set(int k, const Monoid &x) {
seg[k + sz] = x;
}
void build() {
for(int k = sz - 1; k > 0; k--) {
seg[k] = f(seg[2 * k + 0], seg[2 * k + 1]);
}
}
void update(int k, const Monoid &x) {
k += sz;
seg[k] = x;
while(k >>= 1) {
seg[k] = f(seg[2 * k + 0], seg[2 * k + 1]);
}
}
Monoid query(int a, int b) {
Monoid L = M1, R = M1;
for(a += sz, b += sz; a < b; a >>= 1, b >>= 1) {
if(a & 1) L = f(L, seg[a++]);
if(b & 1) R = f(seg[--b], R);
}
return f(L, R);
}
Monoid operator[](const int &k) const {
return seg[k + sz];
}
template< typename C >
int find_subtree(int a, const C &check, Monoid &M, bool type) {
while(a < sz) {
Monoid nxt = type ? f(seg[2 * a + type], M) : f(M, seg[2 * a + type]);
if(check(nxt)) a = 2 * a + type;
else M = nxt, a = 2 * a + 1 - type;
}
return a - sz;
}
template< typename C >
int find_first(int a, const C &check) {
Monoid L = M1;
if(a <= 0) {
if(check(f(L, seg[1]))) return find_subtree(1, check, L, false);
return -1;
}
int b = sz;
for(a += sz, b += sz; a < b; a >>= 1, b >>= 1) {
if(a & 1) {
Monoid nxt = f(L, seg[a]);
if(check(nxt)) return find_subtree(a, check, L, false);
L = nxt;
++a;
}
}
return -1;
}
template< typename C >
int find_last(int b, const C &check) {
Monoid R = M1;
if(b >= sz) {
if(check(f(seg[1], R))) return find_subtree(1, check, R, true);
return -1;
}
int a = sz;
for(b += sz; a < b; a >>= 1, b >>= 1) {
if(b & 1) {
Monoid nxt = f(seg[--b], R);
if(check(nxt)) return find_subtree(b, check, R, true);
R = nxt;
}
}
return -1;
}
};
struct LCA{
ll datasize = 18;
ll n,root;//rootn
vector<ll>depth;//0
vector<vector<ll>>G,data;
LCA(vector<vector<ll>>g,ll r):G(g),n(g.size()),root(r){
data.assign(n,vector<ll>(datasize,-1));
depth.assign(n,-1LL);
build();
}
void build(){
dfs(root,0);
for(ll i=1;i<datasize;i++){
for(ll j=0;j<n;j++){
if(data[j][i-1]==-1)data[j][i]=-1;
else data[j][i]=data[data[j][i-1]][i-1];
}
}
}
void dfs(ll k,ll d){
depth[k]=d;
for(ll i=0;i<G[k].size();i++){
if(depth[G[k][i]]==-1){
data[G[k][i]][0]=k;
dfs(G[k][i],d+1);
}
}
}
ll deep(ll x){
return x>=0?depth[x]:-1;
}
ll query(ll x,ll y){
if(depth[x]<depth[y])swap(x,y);
ll tmp = datasize - 1;
while(depth[x]!=depth[y]){
while(tmp>=1&&depth[y]>=deep(data[x][tmp]))tmp--;
x = data[x][tmp];
}
tmp = datasize - 1;
while(x!=y){
while(tmp>=1&&data[x][tmp]==data[y][tmp])tmp--;
x = data[x][tmp];
y = data[y][tmp];
}
return x;
}
ll path(ll x,ll y){
return deep(x)+deep(y)-2*deep(query(x,y));
}
ll par(ll x){
return data[x][0];
}
};
vector<ll>c(100005);
void dfs(ll k,vector<vector<ll>>&g,ll par){
for(auto to:g[k]){
if(to != par){
//cout<<to spa k<<endl;
chmax(c[to],c[k]);
dfs(to,g,k);
}
}
}
void dfs_tree(ll k,vector<vector<ll>>&g,ll par,vector<ll>&tree){
tree.PB(k);
for(auto to:g[k]){
if(to != par){
dfs_tree(to,g,k,tree);
}
}
}
int main(){
cin.tie(nullptr);
ios_base::sync_with_stdio(false);
ll res=0,buf=0;
bool judge = true;
ll n,k,q;cin>>n>>k>>q;
rep(i,0,n)cin>>c[i];
vector<ll>a(k);
rep(i,0,k)cin>>a[i],a[i]--;
vector<ll>dt;
vector<vector<ll>>g(n);
rep(i,0,n-1){
ll e,f;cin>>e>>f;e--;f--;
g[e].PB(f);
g[f].PB(e);
}
dfs(0,g,-1);
dfs_tree(0,g,-1,dt);
vector<ll>idx(n);
rep(i,0,n)idx[dt[i]]=i;
auto f1=[&](ll x,ll y){return min(x,y);};
auto f2=[&](ll x,ll y){return max(x,y);};
SegmentTree<ll>segmi(k,f1,INT_MAX);
SegmentTree<ll>segmx(k,f2,-INT_MAX);
rep(i,0,k){
segmi.set(i,idx[a[i]]);
segmx.set(i,idx[a[i]]);
}
segmi.build();
segmx.build();
LCA lca(g,0);
while(q--){
ll t;cin>>t;
if(t==1){
ll x,y;cin>>x>>y;x--;y--;
segmi.update(x,idx[y]);
segmx.update(x,idx[y]);
}
else{
ll l,r;cin>>l>>r;l--;
ll mi=dt[segmi.query(l,r)];
ll mx=dt[segmx.query(l,r)];
ll p=lca.query(mi,mx);
//cout<<mi spa mx spa p<<endl;
//debug(c,n);
cout<<c[p]<<endl;
}
}
return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0