#include using namespace std; typedef long long ll; templatebool chmax(T &a, const T &b) { if (abool chmin(T &a, const T &b) { if (b #define vl vector #define vii vector> #define vll vector> #define vvi vector> #define vvl vector> #define vvii vector>> #define vvll vector>> #define vst vector #define pii pair #define pll pair #define pb push_back #define all(x) (x).begin(),(x).end() #define mkunique(x) sort(all(x));(x).erase(unique(all(x)),(x).end()) #define fi first #define se second #define mp make_pair #define si(x) int(x.size()) const int mod=998244353,MAX=300005,INF=15<<26; // BIT セグ木 遅延セグ木 のみ // from: https://gist.github.com/yosupo06/ddd51afb727600fd95d9d8ad6c3c80c9 // (based on AtCoder STL) #include #include #ifdef _MSC_VER #include #endif namespace atcoder { namespace internal { int ceil_pow2(int n) { int x = 0; while ((1U << x) < (unsigned int)(n)) x++; return x; } int bsf(unsigned int n) { #ifdef _MSC_VER unsigned long index; _BitScanForward(&index, n); return index; #else return __builtin_ctz(n); #endif } } // namespace internal } // namespace atcoder #include #include #include namespace atcoder { namespace internal { #ifndef _MSC_VER template using is_signed_int128 = typename std::conditional::value || std::is_same::value, std::true_type, std::false_type>::type; template using is_unsigned_int128 = typename std::conditional::value || std::is_same::value, std::true_type, std::false_type>::type; template using make_unsigned_int128 = typename std::conditional::value, __uint128_t, unsigned __int128>; template using is_integral = typename std::conditional::value || is_signed_int128::value || is_unsigned_int128::value, std::true_type, std::false_type>::type; template using is_signed_int = typename std::conditional<(is_integral::value && std::is_signed::value) || is_signed_int128::value, std::true_type, std::false_type>::type; template using is_unsigned_int = typename std::conditional<(is_integral::value && std::is_unsigned::value) || is_unsigned_int128::value, std::true_type, std::false_type>::type; template using to_unsigned = typename std::conditional< is_signed_int128::value, make_unsigned_int128, typename std::conditional::value, std::make_unsigned, std::common_type>::type>::type; #else template using is_integral = typename std::is_integral; template using is_signed_int = typename std::conditional::value && std::is_signed::value, std::true_type, std::false_type>::type; template using is_unsigned_int = typename std::conditional::value && std::is_unsigned::value, std::true_type, std::false_type>::type; template using to_unsigned = typename std::conditional::value, std::make_unsigned, std::common_type>::type; #endif template using is_signed_int_t = std::enable_if_t::value>; template using is_unsigned_int_t = std::enable_if_t::value>; template using to_unsigned_t = typename to_unsigned::type; } // namespace internal } // namespace atcoder #include #include namespace atcoder { template struct fenwick_tree { using U = internal::to_unsigned_t; public: fenwick_tree() : _n(0) {} fenwick_tree(int n) : _n(n), data(n) {} void add(int p, T x) { assert(0 <= p && p < _n); p++; while (p <= _n) { data[p - 1] += U(x); p += p & -p; } } T sum(int l, int r) { assert(0 <= l && l <= r && r <= _n); return sum(r) - sum(l); } private: int _n; std::vector data; U sum(int r) { U s = 0; while (r > 0) { s += data[r - 1]; r -= r & -r; } return s; } }; } // namespace atcoder #include #include #include #include namespace atcoder { template struct lazy_segtree { public: lazy_segtree() : lazy_segtree(0) {} lazy_segtree(int n) : lazy_segtree(std::vector(n, e())) {} lazy_segtree(const std::vector& v) : _n(int(v.size())) { log = internal::ceil_pow2(_n); size = 1 << log; d = std::vector(2 * size, e()); lz = std::vector(size, id()); for (int i = 0; i < _n; i++) d[size + i] = v[i]; for (int i = size - 1; i >= 1; i--) { update(i); } } void set(int p, S x) { assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); d[p] = x; for (int i = 1; i <= log; i++) update(p >> i); } S get(int p) { assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); return d[p]; } S prod(int l, int r) { assert(0 <= l && l <= r && r <= _n); if (l == r) return e(); l += size; r += size; for (int i = log; i >= 1; i--) { if (((l >> i) << i) != l) push(l >> i); if (((r >> i) << i) != r) push(r >> i); } S sml = e(), smr = e(); while (l < r) { if (l & 1) sml = op(sml, d[l++]); if (r & 1) smr = op(d[--r], smr); l >>= 1; r >>= 1; } return op(sml, smr); } S all_prod() { return d[1]; } void apply(int p, F f) { assert(0 <= p && p < _n); p += size; for (int i = log; i >= 1; i--) push(p >> i); d[p] = mapping(f, d[p]); for (int i = 1; i <= log; i++) update(p >> i); } void apply(int l, int r, F f) { assert(0 <= l && l <= r && r <= _n); if (l == r) return; l += size; r += size; for (int i = log; i >= 1; i--) { if (((l >> i) << i) != l) push(l >> i); if (((r >> i) << i) != r) push((r - 1) >> i); } { int l2 = l, r2 = r; while (l < r) { if (l & 1) all_apply(l++, f); if (r & 1) all_apply(--r, f); l >>= 1; r >>= 1; } l = l2; r = r2; } for (int i = 1; i <= log; i++) { if (((l >> i) << i) != l) update(l >> i); if (((r >> i) << i) != r) update((r - 1) >> i); } } template int max_right(int l) { return max_right(l, [](S x) { return g(x); }); } template int max_right(int l, G g) { assert(0 <= l && l <= _n); assert(g(e())); if (l == _n) return _n; l += size; for (int i = log; i >= 1; i--) push(l >> i); S sm = e(); do { while (l % 2 == 0) l >>= 1; if (!g(op(sm, d[l]))) { while (l < size) { push(l); l = (2 * l); if (g(op(sm, d[l]))) { sm = op(sm, d[l]); l++; } } return l - size; } sm = op(sm, d[l]); l++; } while ((l & -l) != l); return _n; } template int min_left(int r) { return min_left(r, [](S x) { return g(x); }); } template int min_left(int r, G g) { assert(0 <= r && r <= _n); assert(g(e())); if (r == 0) return 0; r += size; for (int i = log; i >= 1; i--) push((r - 1) >> i); S sm = e(); do { r--; while (r > 1 && (r % 2)) r >>= 1; if (!g(op(d[r], sm))) { while (r < size) { push(r); r = (2 * r + 1); if (g(op(d[r], sm))) { sm = op(d[r], sm); r--; } } return r + 1 - size; } sm = op(d[r], sm); } while ((r & -r) != r); return 0; } private: int _n, size, log; std::vector d; std::vector lz; void update(int k) { d[k] = op(d[2 * k], d[2 * k + 1]); } void all_apply(int k, F f) { d[k] = mapping(f, d[k]); if (k < size) lz[k] = composition(f, lz[k]); } void push(int k) { all_apply(2 * k, lz[k]); all_apply(2 * k + 1, lz[k]); lz[k] = id(); } }; } // namespace atcoder #include #include #include namespace atcoder { template struct segtree { public: segtree() : segtree(0) {} segtree(int n) : segtree(std::vector(n, e())) {} segtree(const std::vector& v) : _n(int(v.size())) { log = internal::ceil_pow2(_n); size = 1 << log; d = std::vector(2 * size, e()); for (int i = 0; i < _n; i++) d[size + i] = v[i]; for (int i = size - 1; i >= 1; i--) { update(i); } } void set(int p, S x) { assert(0 <= p && p < _n); p += size; d[p] = x; for (int i = 1; i <= log; i++) update(p >> i); } S get(int p) { assert(0 <= p && p < _n); return d[p + size]; } S prod(int l, int r) { assert(0 <= l && l <= r && r <= _n); S sml = e(), smr = e(); l += size; r += size; while (l < r) { if (l & 1) sml = op(sml, d[l++]); if (r & 1) smr = op(d[--r], smr); l >>= 1; r >>= 1; } return op(sml, smr); } S all_prod() { return d[1]; } template int max_right(int l) { return max_right(l, [](S x) { return f(x); }); } template int max_right(int l, F f) { assert(0 <= l && l <= _n); assert(f(e())); if (l == _n) return _n; l += size; S sm = e(); do { while (l % 2 == 0) l >>= 1; if (!f(op(sm, d[l]))) { while (l < size) { l = (2 * l); if (f(op(sm, d[l]))) { sm = op(sm, d[l]); l++; } } return l - size; } sm = op(sm, d[l]); l++; } while ((l & -l) != l); return _n; } template int min_left(int r) { return min_left(r, [](S x) { return f(x); }); } template int min_left(int r, F f) { assert(0 <= r && r <= _n); assert(f(e())); if (r == 0) return 0; r += size; S sm = e(); do { r--; while (r > 1 && (r % 2)) r >>= 1; if (!f(op(d[r], sm))) { while (r < size) { r = (2 * r + 1); if (f(op(d[r], sm))) { sm = op(d[r], sm); r--; } } return r + 1 - size; } sm = op(d[r], sm); } while ((r & -r) != r); return 0; } private: int _n, size, log; std::vector d; void update(int k) { d[k] = op(d[2 * k], d[2 * k + 1]); } }; } // namespace atcoder using T=int; T f(T a,T b){ return min(a,b); } T ti(){ return INF; } //HL分解(構造体) struct HeavyLightDecomposition{ int n; vector sz,in,out,nxt,par,depth; vector> G; HeavyLightDecomposition(){} HeavyLightDecomposition(int n_){ n=n_; sz.assign(n,0); in.assign(n,0); out.assign(n,0); nxt.assign(n,0); par.assign(n,0); depth.assign(n,0); G.assign(n,vector()); } void add_edge(int u,int v){ G[u].push_back(v); G[v].push_back(u); } void dfs_sz(int u,int p){ par[u]=p; sz[u]=1; if(G[u].size()&&G[u][0]==p) swap(G[u][0],G[u].back()); for(auto &a:G[u]){ if(a==p) continue; depth[a]=depth[u]+1; dfs_sz(a,u); sz[u]+=sz[a]; if(sz[a]>sz[G[u][0]]){ swap(a,G[u][0]); } } } void dfs_hld(int u,int p,int &t){ in[u]=t++; for(auto a:G[u]){ if(a==p) continue; nxt[a]=(a==G[u][0] ? nxt[u] : a); dfs_hld(a,u,t); } out[u]=t; } void build(int u){ int t=0; dfs_sz(u,-1); dfs_hld(u,-1,t); } int lca(int u,int v){ if(in[u]>in[v]) swap(u,v); if(nxt[u]==nxt[v]) return u; return lca(u,par[nxt[v]]); } int mov1(int a,int b){ if(a==b) return a; int c=lca(a,b); if(c==a){ int l=0,r=si(G[a]); while(r-l>1){ int m=(l+r)/2; if(par[a]==G[a][m]){ if(m+1 seg; HeavyLightDecomposition hld; //hld=HeavyLightDecomposition(N)、のようにmainに書く; //buildするときは根が0じゃないとバグるらしい、適宜番号振りを頑張ってくれ ///ここまで貼る void UPD(int l,int r,ll x){ if(l+1==r){ seg.set(l,x); } } void QUE(int l,int r,int &x,bool type){ chmin(x,seg.prod(l,r)); } void updateE(int a,int b,int s){ int c=hld.lca(a,b); while(a!=c&&a!=-1){ if(hld.nxt[a]==hld.nxt[c]){ UPD(max(1,hld.in[c]+1),hld.in[a]+1,s); break; } UPD(max(1,hld.in[hld.nxt[a]]),hld.in[a]+1,s); a=hld.par[hld.nxt[a]]; } while(b!=c&&b!=-1){ if(hld.nxt[b]==hld.nxt[c]){ UPD(max(1,hld.in[c]+1),hld.in[b]+1,s); break; } UPD(max(1,hld.in[hld.nxt[b]]),hld.in[b]+1,s); b=hld.par[hld.nxt[b]]; } } ll queryE(int a,int b){ int c=hld.lca(a,b); int res=INF; while(a!=c&&a!=-1){ if(hld.nxt[a]==hld.nxt[c]){ QUE(max(1,hld.in[c]+1),hld.in[a]+1,res,0); break; } QUE(max(1,hld.in[hld.nxt[a]]),hld.in[a]+1,res,0); a=hld.par[hld.nxt[a]]; } while(b!=c&&b!=-1){ if(hld.nxt[b]==hld.nxt[c]){ QUE(max(1,hld.in[c]+1),hld.in[b]+1,res,0); break; } QUE(max(1,hld.in[hld.nxt[b]]),hld.in[b]+1,res,0); b=hld.par[hld.nxt[b]]; } return res; } //辺属性 int main(){ std::ifstream in("text.txt"); std::cin.rdbuf(in.rdbuf()); cin.tie(0); ios::sync_with_stdio(false); int N,M;cin>>N>>M; hld=HeavyLightDecomposition(N); vii E(N-1); for(int i=0;i>a>>b;a--;b--; E[i]=mp(a,b); hld.add_edge(a,b); } hld.build(0); seg=atcoder::segtree(N); for(int i=0;i>a>>b;a--;b--; int res=queryE(a,b); if(res) imo[res-1]++; } for(int i=N-3;i>=0;i--) imo[i]+=imo[i+1]; for(int i=0;i