Codeforces Round 891 (Div. 3) G. Counting Graphs

Problem Description

给定一棵 n n n 个节点的树,树的每条边权为 w i w_i wi

对一个图有以下 4 4 4 种限制:

  1. 图没有自环和重边。
  2. 图中边权为整数,且不超过 S S S
  3. 该图只有一棵最小生成树。
  4. 该图的最小生成树为给定的树。

求可行图的个数。

Input

第一行输入一个整数 t ( 1 ≤ t ≤ 1 0 4 ) t(1 \le t \le 10^4) t(1t104) ,表示测试组数。

对于每组测试,

第一行输入两个整数 n , S ( 2 ≤ n ≤ 2 × 1 0 5 , 1 ≤ S ≤ 1 0 9 ) n,S(2 \le n \le 2 \times 10^5,1 \le S \le 10^9) n,S(2n2×105,1S109) ,表示节点个数和边权上界。

接下来 n − 1 n-1 n1 ,每行输入三个整数 u i , v i , w i ( 1 ≤ u i , v i ≤ n , u i ≠ v i , 1 ≤ w i ≤ S ) u_i,v_i,w_i(1 \le u_i,v_i \le n,u_i \not= v_i,1 \le w_i \le S) ui,vi,wi(1ui,vin,ui=vi,1wiS) ,表示 u i u_i ui v i v_i vi 之间有一条边权为 w i w_i wi 的边。

题目保证 ∑ n ≤ 2 × 1 0 5 \sum n \le 2 \times 10^5 n2×105

Output

输出模 998244353 998244353 998244353 后的可行图的个数。

Solution

我们考虑 K r u s k a l Kruskal Kruskal 生成的过程中,会发现当添加了一条权值为 w i w_i wi 的边后,会合并两个连通块得到一个新的连通块 b b b ,而在 b b b 中添加任意边权 w j ∈ ( w i , S ] w_j \in (w_i,S] wj(wi,S] 后,并不会妨碍最小生成树的唯一性。

那么这题的解法便变成了在做 K r u s k a l Kruskal Kruskal 每次合并两个连通块 x , y x,y x,y 时,计算能够加边对于 a n s ans ans 的贡献,贡献即为加边方案数 c n t cnt cnt ,然后令 a n s = a n s × c n t ans=ans \times cnt ans=ans×cnt 即可。但是我们会发现由于之前生成 x , y x,y x,y 连通块时已经计算了他们各自能够提供的方案数贡献,所以我们在每次计算方案数时将能加的边不考虑 x , y x,y x,y 内部生成的,那么新生成的边数就变成了 s i z x × s i z y − 1 siz_x \times siz_y - 1 sizx×sizy1 c n t = ( S − w i + 1 ) s i z x × s i z y − 1 cnt=(S-w_i+1)^{siz_x \times siz_y -1} cnt=(Swi+1)sizx×sizy1 ,这样我们就能不重不漏地计算出答案。

最终时间复杂度为 O ( n l o g n ) O(nlogn) O(nlogn)

Code

#include <bits/stdc++.h>
#define endl '\n'
using namespace std;
typedef long long ll;
constexpr int mod = 998244353;
ll qmi(ll a, ll k) {
	ll res = 1 % mod, t = a % mod;
	while (k) {
		if (k & 1)res = res * t % mod;
		t = t * t % mod;
		k >>= 1;
	}
	return res;
}
void solve() {
	int n, S;
	cin >> n >> S;
	vector<ll>p(n + 1), siz(n + 1, 1); iota(p.begin(), p.end(), 0);
	vector<tuple<int, int, int>>edge;
	auto leader = [&](int x) {
		while (x != p[x])x = p[x] = p[p[x]];
		return x;
	};
	for (int i = 1; i <= n - 1; i++) {
		int u, v, w;
		cin >> u >> v >> w;
		edge.emplace_back(w, u, v);
	}
	sort(edge.begin(), edge.end());
	ll ans = 1;
	for (auto [w, u, v] : edge) {
		u = leader(u), v = leader(v);
		ans = ans * qmi(S - w + 1, siz[u] * siz[v] - 1) % mod;
		siz[u] += siz[v];
		p[v] = u;
	}
	cout << ans << endl;
}
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int t;
	cin >> t;
	while (t--)solve();
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值