結果

問題 No.399 動的な領主
ユーザー SnowBeenDidingSnowBeenDiding
提出日時 2022-02-26 18:34:59
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 430 ms / 2,000 ms
コード長 10,447 bytes
コンパイル時間 5,843 ms
コンパイル使用メモリ 319,960 KB
実行使用メモリ 22,204 KB
最終ジャッジ日時 2023-09-17 20:57:05
合計ジャッジ時間 11,058 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,384 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 3 ms
4,380 KB
testcase_05 AC 25 ms
4,692 KB
testcase_06 AC 398 ms
17,352 KB
testcase_07 AC 399 ms
17,224 KB
testcase_08 AC 405 ms
17,280 KB
testcase_09 AC 400 ms
17,336 KB
testcase_10 AC 4 ms
4,380 KB
testcase_11 AC 20 ms
4,704 KB
testcase_12 AC 305 ms
17,680 KB
testcase_13 AC 292 ms
17,964 KB
testcase_14 AC 79 ms
22,204 KB
testcase_15 AC 106 ms
22,128 KB
testcase_16 AC 159 ms
19,632 KB
testcase_17 AC 430 ms
17,228 KB
testcase_18 AC 415 ms
17,544 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <atcoder/all>
using namespace atcoder;
using mint = modint;

#include <bits/stdc++.h>
#define rep(i,a,b) for(ll i=(ll)(a);i<(ll)(b);i++)
#define repeq(i,a,b) for(ll i=(ll)(a);i<=(ll)(b);i++)
#define repreq(i,a,b) for(ll i=(ll)(a);i>=(ll)(b);i--)
#define each(a,b) for(auto&(a):(b))
#define endl '\n'
#define cYes cout<<"Yes"<<endl
#define cNo cout<<"No"<<endl
#define sortr(v) sort(v,greater<>())
#define pb push_back
#define mp make_pair
#define mt make_tuple
#define tget(a,b) get<b>(a)
#define FI first
#define SE second
#define ALL(v) (v).begin(),(v).end()
#define INFLL 3000000000000000000LL
#define INF 1000000010
#define PI acos(-1.0L)

using namespace std;

typedef long long ll;
typedef pair<ll,ll> Pll;
typedef tuple<ll,ll,ll> Tlll;
typedef vector<int> Vi;
typedef vector<Vi> VVi;
typedef vector<ll> Vl;
typedef vector<Vl> VVl;
typedef vector<VVl> VVVl;
typedef vector<mint> Vm;
typedef vector<Vm> VVm;
typedef vector<string> Vs;
typedef vector<double> Vd;
typedef vector<char> Vc;
typedef vector<bool> Vb;
typedef vector<Pll> VPll;
typedef priority_queue<ll> PQl;
typedef priority_queue<ll,vector<ll>,greater<ll>> PQlr;

/* print */
template <typename T>
ostream& operator<<(ostream& os,const vector<T> &V){
int N=V.size();if(N==0){os<<endl;return os;}
rep(i,0,N-1){os<<V[i]<<' ';}os<<V[N-1]<<endl;return os;}
template <typename T>
ostream& operator<<(ostream& os,const vector<vector<T>> &V){
int N=V.size();rep(i,0,N)os<<V[i];return os;}
template <typename T,typename S>
ostream& operator<<(ostream&os,pair<T,S>const&P){os<<P.FI<<' '<<P.SE;return os;}
ostream& operator<<(ostream&os,mint const&M){os<<M.val();return os;}

/* useful */
template<typename T>void Vin(vector<T> &v){int n=v.size();rep(i,0,n)cin>>v[i];}
template<typename T>int SMALLER(vector<T> &a,T x){
return lower_bound(a.begin(),a.end(),x)-a.begin();}
template<typename T>int orSMALLER(vector<T> &a,T x){
return upper_bound(a.begin(),a.end(),x)-a.begin();}
template<typename T>int BIGGER(vector<T>&a,T x){return a.size()-orSMALLER(a,x);}
template<typename T>int orBIGGER(vector<T>&a,T x){return a.size()-SMALLER(a,x);}
template<typename T>int COUNT(vector<T> &a,T x){
return upper_bound(ALL(a),x)-lower_bound(ALL(a),x);}
template<typename T>bool chmax(T &a,T b) {if(a<b){a=b;return 1;}return 0;}
template<typename T>bool chmin(T &a,T b) {if(a>b){a=b;return 1;}return 0;}
template<typename T>void press(T &v){v.erase(unique(ALL(v)),v.end());}
template<typename T>vector<int> zip(vector<T> b){pair<T,int> p[b.size()+10];
int a=b.size();vector<int> l(a);for(int i=0;i<a;i++) p[i]=mp(b[i],i);
sort(p,p+a);int w=0;for(int i=0;i<a;i++){if(i&&p[i].first!=p[i-1].first)w++;
l[p[i].second]=w;}return l;}
template<typename T>vector<T> vis(vector<T> &v){
vector<T> S(v.size()+1);rep(i,1,S.size())S[i]+=v[i-1]+S[i-1];return S;}

ll dem(ll a,ll b){return((a+b-1)/(b));}
ll dtoll(double d,int g){return round(d*pow(10,g));}

//const long long MOD = 1000000007;
const long long MOD = 998244353;
const double EPS = 1e-10;

void init(){mint::set_mod(MOD);
cin.tie(0);cout.tie(0);ios::sync_with_stdio(0);cout<<fixed<<setprecision(12);}

/********************************** START **********************************/


class HLdecomposition {
public:
    int V;
    vector<vector<int> > G;
    vector<int> stsize, parent, pathtop, in, out;
    int root;
    void BuildStsize(int u, int p){
        stsize[u] = 1, parent[u] = p;
        for(int& v : G[u]){
            if(v == p){
                if(v == G[u].back()) break;
                else swap(v, G[u].back());
            }
            BuildStsize(v, u);
            stsize[u] += stsize[v];
            if(stsize[v] > stsize[G[u][0]]){
                swap(v, G[u][0]);
            }
        }
    }
    void BuildPath(int u, int p, int& tm){
        in[u] = tm++;
        for(int v : G[u]){
            if(v == p) continue;
            pathtop[v] = (v == G[u][0] ? pathtop[u] : v);
            BuildPath(v, u, tm);
        }
        out[u] = tm;
    }
 
public:
    void add_edge(int u, int v){
        G[u].push_back(v), G[v].push_back(u);
    }
    void build(int _root){
        root = _root;
        int tm = 0;
        BuildStsize(root, -1);
        pathtop[root] = root;
        BuildPath(root, -1, tm);
    }
    //元の頂点のインデックスの配列上でのidを返す
    inline int get(int a){
        return in[a];
    }
    int lca(int a, int b){
        int pa = pathtop[a], pb = pathtop[b];
        while(pathtop[a] != pathtop[b]){
            if(in[pa] > in[pb]){
                a = parent[pa], pa = pathtop[a];
            }else{
                b = parent[pb], pb = pathtop[b];
            }
        }
        if(in[a] > in[b]) swap(a, b);
        return a;
    }
    void subtree_query(int a, const function< void(int, int) > &func){
        func(in[a], out[a]);
    }
    void path_query(int a, int b, const function< void(int, int) > &func){
        int pa = pathtop[a], pb = pathtop[b];
        while(pathtop[a] != pathtop[b]){
            if(in[pa] > in[pb]){
                func(in[pa], in[a] + 1);
                a = parent[pa], pa = pathtop[a];
            }else{
                func(in[pb], in[b] + 1);
                b = parent[pb], pb = pathtop[b];
            }
        }
        if(in[a] > in[b]) swap(a, b);
        func(in[a], in[b] + 1);
    }
    HLdecomposition() {}
    HLdecomposition(int node_size) : V(node_size), G(V), stsize(V, 0), parent(V, -1),
        pathtop(V, -1), in(V, -1), out(V, -1){}
};


struct Tree {

  int n;
  vector<vector<int>> gr;
  vector<vector<int>> chi;
  vector<vector<int>> par_exp;
  vector<int> al;
  vector<int> euler;
  vector<int> par;
  vector<int> dep;
  vector<int> subsize;
  HLdecomposition hl;
  
  Tree() {}
  Tree(int input) {
    n = input;
    gr = vector<vector<int>>(n);
    par = vector<int>(n);
    dep = vector<int>(n);
    chi = vector<vector<int>>(n);
  }
  Tree(vector<vector<int>> gr) {
    par = vector<int>(n);
    dep = vector<int>(n);
    chi = vector<vector<int>>(n);
    int st = 0;
    int now = st;
    vector<bool> al(n);
    al[now] = 1;
    queue<int> q;
    q.push(now);
    while(!q.empty()) {
      now = q.front();
      q.pop();
      rep(i,0,(ll)gr[now].size()) {
        int ne = gr[now][i];
        if(!al[ne]) {
          connect(now, ne);
          q.push(ne);
          al[ne] = 1;
        }
      }
    }
  }
  void connect(int s,int t) {
    gr[s].push_back(t);
    gr[t].push_back(s);
  }
  void init_dfs (int d) {
    int len = gr[d].size();
    int sum_subsize = 0;
    rep(i,0,len) {
      if(al[gr[d][i]]) continue;
      al[gr[d][i]] = 1;
      par[gr[d][i]] = d;
      dep[gr[d][i]] = dep[d] + 1;
      chi[d].push_back(gr[d][i]);
      euler.push_back(d);
      init_dfs(gr[d][i]);
      sum_subsize += subsize[gr[d][i]];
    }
    subsize[d] = sum_subsize + 1;
    euler.push_back(d);
  }
  void make(int root) {
    rep(i,0,n) sort(ALL(gr[i]));
    al = vector<int>(n,0);
    par = vector<int>(n);
    dep = vector<int>(n);
    subsize = vector<int>(n);
    chi = vector<vector<int>>(n);
    euler = {};
    al[root] = 1;
    par[root] = -1;
    dep[root] = 0;
    init_dfs(root);
  }
  int rad() {
    int max_ind = -1, max_dep = -1;
    int s,t;
    rep(i,0,n) if(chmax(max_dep, dep[i])) max_ind = i;
    make(max_ind);
    s = max_ind;
    max_ind = -1,max_dep = -1;
    rep(i,0,n) if(chmax(max_dep, dep[i])) max_ind = i;
    t = max_ind;
    swap(s,t);
    return max_dep;
    //return make_pair(s,t);
  }
  void init_lca() {
    int ct = 0;
    while(n >= (1<<ct)) ct++;
    ct++;
    par_exp = vector<vector<int>>(ct,vector<int>(n,-1));
    rep(i,0,n) par_exp[0][i] = par[i];
    rep(i,0,ct-1){
      rep(j,0,n){
        if(par_exp[i][j] == -1) par_exp[i+1][j] = -1;
        else par_exp[i+1][j] = par_exp[i][par_exp[i][j]];
      }
    }
  }
  int lca(int u,int v) {
    int h = par_exp.size();
    if(dep[u] > dep[v]) swap(u,v);
    for(int i=h-1; i>=0; i--) {
      if(((dep[v]-dep[u])>>i)&1) v = par_exp[i][v];
    }
    if(u == v) return u;
    for(int i=h-1;i>=0;i--) {
      if(par_exp[i][u] != par_exp[i][v]){
        u = par_exp[i][u];
        v = par_exp[i][v];
      }
    }
    return par_exp[0][u];
  }
  void dfs_(int d) { // dfs template
    int len = gr[d].size();
    rep(i,0,len) {
      if(al[gr[d][i]]) continue;
      al[gr[d][i]] = 1;
      dfs_(gr[d][i]);
    }
  }
  void dfs(int root) {
    al = vector<int>(n);
    al[root] = 1;
    dfs_(root);
  }
  void bfs(int root) { // bfs template
    al = vector<int>(n);
    al[root] = 1;
    queue<int> q;
    stack<int> st;
    q.push(root);
    st.push(root);
    vector<int> order;
    while(!q.empty()) {
      int d = q.front();
      order.pb(d);
      // 今見ている頂点
      q.pop();
      int len = gr[d].size();
      rep(i,0,len) {
        if(al[gr[d][i]]) continue;
        al[gr[d][i]] = 1;
        q.push(gr[d][i]);
        st.push(gr[d][i]);
        // ここで次の頂点を見ている
      }
    }
    while(!st.empty()) {
      int d = st.top();
      st.pop();
      order.pb(d);
      // ここでbfsの逆順に頂点をたどっている
    }
  }

};
/* 使い方 */
/*
  Tree tr(n) 	      // 宣言: n は頂点数
  tr.connect(a,b)   // 頂点 a と b をつなぐ
  tr.make(k) 		    // 頂点 k を根としてpar,dep,chiを作り上げる
  tr.init_lca()     // lcaを使う前に必ずこれをする
  tr.lca(a,b)       // 頂点aとbの共通の"先祖"で一番近い"頂点"を返す
  tr.rad()          // 木の直径を返す
*/

struct S {ll value, size;};
using F = ll;

S fx(S l, S r) { return {l.value + r.value, l.size + r.size}; }
S ex() { return {0, 0}; }
S fa(F f, S x) { return {x.value + x.size * f, x.size}; }
F fm(F f, F g) { return (f + g); } // gにfを加えるイメージ
F em() { return 0; }

void sol() {
  int n;
  cin >> n;
  
  HLdecomposition hl(n);
  rep(i,0,n-1) {
    int qw,er;
    cin >> qw >> er;
    qw--,er--;
    hl.add_edge(qw,er);
  }
  
  hl.build(0);
  vector<S> v(n, {0, 1});
  lazy_segtree<S, fx, ex, F, fa, fm, em> st(v);
  auto f = [&](int l, int r){
    st.apply(l,r,1);
  };
  int q;cin >> q;
  while(q--) {
    int qw,er;
    cin >> qw >> er;
    qw--,er--;
    hl.path_query(qw,er,f);
  }
  ll ans = 0;
  rep(i,0,n) {
    ll qw = st.get(i).value;
    //cerr << qw << endl;
    ans += qw * (qw + 1) / 2;
  }
  cout << ans << endl;
  
}

int main() {
  init();
  int q = 1;
  //cin >> q;
  while(q--) sol();
  return 0;
}
0