#380816 (C++14) No.901 K-ary εxtrεεmε

提出ソース

結果

問題	No.901 K-ary εxtrεεmε
ユーザー	polylogK
提出日時	2019-09-19 23:27:33
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	884 ms / 3,000 ms
コード長	6,975 bytes
コンパイル時間	2,545 ms
コンパイル使用メモリ	198,100 KB
実行使用メモリ	30,588 KB
最終ジャッジ日時	2024-04-15 00:03:38
合計ジャッジ時間	20,237 ms
ジャッジサーバーID （参考情報）	judge4 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	270 ms 30,588 KB
testcase_01	AC	2 ms 6,940 KB
testcase_02	AC	6 ms 6,940 KB
testcase_03	AC	7 ms 6,940 KB
testcase_04	AC	7 ms 6,944 KB
testcase_05	AC	7 ms 6,940 KB
testcase_06	AC	7 ms 6,944 KB
testcase_07	AC	798 ms 28,208 KB
testcase_08	AC	794 ms 28,200 KB
testcase_09	AC	768 ms 28,068 KB
testcase_10	AC	796 ms 28,260 KB
testcase_11	AC	783 ms 28,076 KB
testcase_12	AC	641 ms 28,144 KB
testcase_13	AC	675 ms 28,288 KB
testcase_14	AC	679 ms 28,084 KB
testcase_15	AC	640 ms 28,096 KB
testcase_16	AC	674 ms 28,196 KB
testcase_17	AC	848 ms 28,060 KB
testcase_18	AC	860 ms 28,084 KB
testcase_19	AC	854 ms 28,172 KB
testcase_20	AC	884 ms 28,128 KB
testcase_21	AC	816 ms 28,060 KB
testcase_22	AC	454 ms 28,200 KB
testcase_23	AC	443 ms 28,144 KB
testcase_24	AC	448 ms 28,272 KB
testcase_25	AC	440 ms 28,204 KB
testcase_26	AC	452 ms 28,088 KB
testcase_27	AC	525 ms 28,128 KB
testcase_28	AC	523 ms 28,184 KB
testcase_29	AC	516 ms 28,296 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std::literals::string_literals;
using i64 = long long;
using std::cout;
using std::endl;
using std::cin;

template<typename T>
std::vector<T> make_v(size_t a){return std::vector<T>(a);}

template<typename T,typename... Ts>
auto make_v(size_t a,Ts... ts){
  return std::vector<decltype(make_v<T>(ts...))>(a,make_v<T>(ts...));
}

// build: you have to run it before you use
// for_each: process [l, r]
// for_each_edge: process [l, r]
// distance: returns the dist (l, r)
class heavy_light_decomposition {
	using size_type = size_t;

	using F = std::function<void (int, int)>;

	public:
	std::vector<std::vector<int>> g;
	std::vector<int> vid, inv;

	private:
	std::vector<int> head, sz, heavy, par, dep, type;

	void dfs(int root) {
		using P = std::pair<int, int>;
		
		par[root] = -1;
		dep[root] = 0;
		std::stack<P> st;
		st.push({root, 0});
		while(!st.empty()) {
			int v = st.top().first;
			int & i = st.top().second;
		
			if(i < g[v].size()) {
				int u = g[v][i++];
				if(u == par[v]) continue;
				par[u] = v;
				dep[u] = dep[v] + 1;
				st.push({u, 0});
			} else {
				st.pop();
				int tmp = 0;
				for(auto e: g[v]) {
					if(e == par[v]) continue;
					sz[v] += sz[e];
					if(tmp < sz[e]) {
						tmp = sz[e];
						heavy[v] = e;
					}
				}
			}
		}
	}
	void bfs(int root, int c, int & k) {
		std::queue<int> qu({root});
		while(!qu.empty()) {
			int h = qu.front(); qu.pop();

			for(int v = h; v != -1; v = heavy[v]) {
				type[v] = c;
				vid[v] = k++;
				inv[vid[v]] = v;
				head[v] = h;
				for(int e: g[v])
					if(e != par[v] and e != heavy[v]) qu.push(e);
			}
		}
	}

	public:
	heavy_light_decomposition() {}
	heavy_light_decomposition(int n) :
		g(n), vid(n, -1), head(n), sz(n, 1), heavy(n, -1),
		par(n), dep(n), inv(n), type(n) {}
	void add_edge(int a, int b) {
		g[a].push_back(b);
		g[b].push_back(a);
	}
	void build(std::vector<int> rs = std::vector<int>(1, 0)) {
		int c = 0, k = 0;
		for(int r: rs) {
			dfs(r);
			bfs(r, c++, k);
		}
	}
	int lca(int u, int v) {
		while(true) {
			if(vid[u] > vid[v]) std::swap(u, v);
			if(head[u] == head[v]) return u;
			v = par[head[v]];
		}
	}
	void for_each(int u, int v, const F & f) {
		while(true) {
			if(vid[u] > vid[v]) std::swap(u, v);
			f(std::max(vid[head[v]], vid[u]), vid[v]);
			if(head[u] != head[v]) v = par[head[v]];
			else break;
		}
	}
	void for_each_edge(int u, int v, const F & f) {
		while(true) {
			if(vid[u] > vid[v]) std::swap(u, v);
			if(head[u] != head[v]) {
				f(vid[head[v]], vid[v]);
				v = par[head[v]];
			} else {
				if(u != v) f(vid[u] + 1, vid[v]);
				break;
			}
		}
	}

	int distance(int u, int v) {
		return dep[u] + dep[v] - 2 * dep[lca(u, v)];
	}
};

using namespace std;
template <typename T,typename E>
struct SegmentTree{
  using F = function<T(T,T)>;
  using G = function<T(T,E)>;
  using H = function<E(E,E)>;
  int n,height;
  F f;
  G g;
  H h;
  T ti;
  E ei;
  vector<T> dat;
  vector<E> laz;
  SegmentTree(F f,G g,H h,T ti,E ei):
    f(f),g(g),h(h),ti(ti),ei(ei){}

  void init(int n_){
    n=1;height=0;
    while(n<n_) n<<=1,height++;
    dat.assign(2*n,ti);
    laz.assign(2*n,ei);
  }
  void build(const vector<T> &v){
    int n_=v.size();
    init(n_);
    for(int i=0;i<n_;i++) dat[n+i]=v[i];
    for(int i=n-1;i;i--)
      dat[i]=f(dat[(i<<1)|0],dat[(i<<1)|1]);
  }
  inline T reflect(int k){
    return laz[k]==ei?dat[k]:g(dat[k],laz[k]);
  }
  inline void eval(int k){
    if(laz[k]==ei) return;
    laz[(k<<1)|0]=h(laz[(k<<1)|0],laz[k]);
    laz[(k<<1)|1]=h(laz[(k<<1)|1],laz[k]);
    dat[k]=reflect(k);
    laz[k]=ei;
  }
  inline void thrust(int k){
    for(int i=height;i;i--) eval(k>>i);
  }
  inline void recalc(int k){
    while(k>>=1)
      dat[k]=f(reflect((k<<1)|0),reflect((k<<1)|1));
  }
  void update(int a,int b,E x){
    thrust(a+=n);
    thrust(b+=n-1);
    for(int l=a,r=b+1;l<r;l>>=1,r>>=1){
      if(l&1) laz[l]=h(laz[l],x),l++;
      if(r&1) --r,laz[r]=h(laz[r],x);
    }
    recalc(a);
    recalc(b);
  }
  void set_val(int a,T x){
    thrust(a+=n);
    dat[a]=x;laz[a]=ei;
    recalc(a);
  }
  T query(int a,int b){
    thrust(a+=n);
    thrust(b+=n-1);
    T vl=ti,vr=ti;
    for(int l=a,r=b+1;l<r;l>>=1,r>>=1) {
      if(l&1) vl=f(vl,reflect(l++));
      if(r&1) vr=f(reflect(--r),vr);
    }
    return f(vl,vr);
  }
};

struct node {
	int used;
	int depth;
	int label;
	int depth_all;
	int label_all;
	
	node() : used(1), depth(-1), label(-1), depth_all(-1), label_all(-1) {}
	node(int used, int depth, int label, int depth_all, int label_all) : used(used), depth(depth), label(label), depth_all(depth_all), label_all(label_all) {}
};

int main() {
	int n; scanf("%d", &n);

	std::vector<std::vector<std::pair<int, i64>>> g(n);
	heavy_light_decomposition hld(n);
	for(int i = 0; i < n - 1; i++) {
		int a, b, c; scanf("%d%d%d", &a, &b, &c);

		g[a].push_back({b, c});
		g[b].push_back({a, c});
		hld.add_edge(a, b);
	}
	hld.build();

	std::vector<i64> dist(n), depth(n);
	auto dfs = [&](auto && dfs, int v, int par) -> void {
		for(auto e: g[v]) {
			if(e.first == par) continue;
			dist[e.first] = dist[v] + e.second;
			depth[e.first] = depth[v] + 1;
			dfs(dfs, e.first, v);
		}
	};
	dfs(dfs, 0, -1);
	
	auto calc = [dist, &hld](int u, int v) {
		return dist[u] + dist[v] - 2LL * dist[hld.lca(u, v)];
	};

	auto f = [](node a, node b) {
		node ret;
		ret.used = 0;
		if(a.depth_all > b.depth_all) {
			ret.label_all = a.label_all;
			ret.depth_all = a.depth_all;
		} else {
			ret.label_all = b.label_all;
			ret.depth_all = b.depth_all;
		}

		int A = a.label, B = a.depth, C = b.label, D = b.depth;
		if(a.used) {
			A = a.label_all;
			B = a.depth_all;
		}
		if(b.used) {
			C = b.label_all;
			D = b.depth_all;
		}

		if(B > D) {
			ret.label = A;
			ret.depth = B;
		} else {
			ret.label = C;
			ret.depth = D;
		}
		return ret;
	};
	auto gg = [](node a, int b) {
		if(b) {
			a.used = 1;
			a.depth = a.depth_all;
			a.label = a.label_all;
		} else {
			a.used = 0;
			a.depth = -1;
			a.label = -1;
		}
		return a;
	};
	auto h = [](int a, int b) { return b; };
	
	SegmentTree<node, int> A(f, gg, h,  node(), -1);
	A.init(n);
	for(int i = 0; i < n; i++) {
		int v = hld.inv[i];

		A.set_val(i, node(0, -1, -1, depth[v], v));
	}
	
	auto kiri = [&A](int l, int r) {
		A.update(l, r + 1, 1);
	};
	auto irik = [&A](int l, int r) {
		A.update(l, r + 1, 0);
	};

	int q; scanf("%d", &q);
	while(q--) {
		int k; scanf("%d", &k);
		std::vector<int> vec(k);
		for(int i = 0; i < k; i++) scanf("%d", &vec[i]);

		i64 ans = 0;
		int c = vec.front();
		for(auto v: vec) {
			node tmp = node();
			auto tanpo = [&](int l, int r) {
				tmp = f(tmp, A.query(l, r + 1));
			};
			hld.for_each(0, v, tanpo);

			if(tmp.depth == -1) {
				ans += calc(v, c);
				hld.for_each(v, c, kiri);
				c = hld.lca(v, c);
			} else {
				int u = tmp.label;
				ans += calc(u, v);
				hld.for_each(u, v, kiri);
			}
		}

		A.update(0, n, 0);
		printf("%lld\n", ans);
	}
	return 0;
}

yukicoder

結果

テストケース

ソースコード