题解 P3412 【仓鼠找sugar II】

\(\huge\texttt{P3412}\)

题意

给定一棵树,求任意一条路径从起点随机游走到终点的期望距离的期望。

思路

讨论求出每条边的贡献,当且仅当每条路径两个端点分别在这条边分成两个连通块中。

其中向上的期望可以从 \(u->fa[u]\) 或者 \(u->v->u->fa[u]\)。

向下的期望可以从 \(u->v\) 或者 \(u->v'->u->v\) 或者 \(u->fa[u]->u->v\)(\(v'\in u\&v'\ne v\))。

然后发现这些都可以通过已知项求得,对其讨论即可。

然后列式子求解每条边向上走向下走的期望步数即可,为了方便标记每条边,\(up[u]\) 和 \(down[u]\) 都表示连 \(u\) 的父边,而 \(up\) 表示向上,\(down\) 对应向下。

则会有下面的式子(感觉我推烦了,但还是比较清晰的吧):

令 \(sum[u]=\sum_{v\in u} up[v]\),\(v\in u\) 表示v是u的子节点,\(p[u]\) 表示 \(u\) 的连边数量。

\[up[u]=\frac{1}{p[u]}+\frac{1}{p[u]}(\sum_{v\in u}(1+up[v]+up[u]))\\ up[u]=\frac{1}{p[u]}+\frac{(p[u]-1)up[u]}{p[u]}+\frac{1}{p[u]}(\sum_{v\in u}(1+up[v]))\\ up[u]=1+p[u]-1+\sum_{v\in u}up[v]\\ up[u]=p[u]+\sum_{v\in u}up[v] \]

\[down[u]=\frac{1}{p[fa[u]]}+\frac{1}{p[fa[u]]}(\sum_{v\in fa[u]\&v\ne u}1+up[v]+down[u])+\frac{1}{p[fa[u]]}(1+down[fa[u]]+down[u])\\ donw[u]=1+\frac{(p-1)(down[u])}{p[fa[u]]}+\frac{1}{p[fa[u]]}(sum[fa[u]]-up[u]+down[fa[u]])\\ down[u]=p[fa[u]]+sum[fa[u]]-up[u]+down[fa[u]] \]

然后就很容易了。

代码

int a, b, up[N + 5], sum[N + 5], down[N + 5], ans, siz[N + 5];
vector<int> st[N + 5];

void dfs(int n, int fa) {
	siz[n] = 1;
	rep(i, 0, siz(st[n]) - 1) {
		int v = st[n][i];
		if (v == fa) continue;
		dfs(v, n);
		siz[n] += siz[v];
		(sum[n] += up[v]) %= mod;
	}
	up[n] = (siz(st[n]) + sum[n]) % mod;
}

void Dfs(int n, int fa) {
	if (n != 1) down[n] = (siz(st[fa]) + sum[fa] - up[n] + down[fa] + mod) % mod;
	(ans += siz[n] * (a - siz[n]) % mod * (down[n] + up[n]) % mod) %= mod;
	rep(i, 0, siz(st[n]) - 1) {
		int v = st[n][i];
		if (v == fa) continue;
		Dfs(v, n);
	}
}

int qpow(int n, int m = mod - 2) {
	int res = 1;
	for (; m; m >>= 1) {
		if (m & 1) res = res * n % mod;
		n = n * n % mod;
	}
	return res;
}

signed main() {
	// freopen("in1.in", "r", stdin);
	// freopen("out.out", "w", stdout);
	a = read();
	int x, y;
	rep(i, 1, a - 1) {
		x = read();
		y = read();
		st[x].PB(y);
		st[y].PB(x);
	}
	dfs(1, 0);
	Dfs(1, 0);
	printf("%lld", ans * qpow(a * a % mod) % mod);
	return 0;
}
上一篇:格式化namenode时:SHUTDOWN_MSG: Shutting down NameNode at java.net.UnknownHostException: xxx


下一篇:顺时针打印矩阵(剑指offer)