atcoder 160 F

题目:https://atcoder.jp/contests/abc160/tasks/abc160_f

逆元 :https://blog.csdn.net/qq_35416331/article/details/81059747

换根dp:https://www.cnblogs.com/Kanoon/p/12644109.html 

#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
#include<math.h>
#include<string>
#include<string.h>
#include<map>
#include<unordered_map>
#include<unordered_set>
#include<set>
#include<stack>
#include<sstream>
//#include <Eigen/Dense>
//#include <Eigen/Geometry>
//using namespace Eigen;
using namespace std;
 
int n;
const int maxn = 2e5 + 100;
int mod = 1e9 + 7;
vector<int> g[maxn];
int vis[maxn], sz[maxn];
long long ans[maxn], inv[maxn];
 
void inv_n(int n)
{
	inv[1] = 1;
	for (long long i = 2; i <= n; i++)
	{
		inv[i] = inv[mod%i] * (mod - mod / i) % mod;
	}
}
 
 
int dfs1(int cur,int p)
{
	sz[cur] = 1;
	for (int i : g[cur])
	{
		if (i != p)	
			sz[cur] += dfs1(i, cur);
	}
	return sz[cur];
}
 
void dfs2(int cur, int p)
{
	ans[cur] = ans[p] * sz[cur] % mod * inv[(n - sz[cur])] % mod;//ans[p] * sz[cur] * inv[(n - sz[cur])] % mod溢出报错
	for (int i : g[cur])
	{
		if (i != p)
			dfs2(i,cur);
	}
}
int main()
{
	int u, v;
	cin >> n;
	inv_n(n);
	for (int i = 0; i < n - 1; i++)
	{
		cin >> u >> v;
		g[u].push_back(v);
		g[v].push_back(u);
	}
	ans[1] = 1;
	for (int i = 1; i <= n; i++)
		ans[1] = ans[1] * i%mod;
	dfs1(1, 0);
	for (int i = 1; i <= n; i++)
		ans[1] = ans[1] * inv[sz[i]] % mod;
	
	for (int i : g[1])
		dfs2(i, 1);
	for (int i = 1; i <= n; i++)
		cout << ans[i] << endl;
 
	return 0;
}

 


#include<iostream>
#include <algorithm>
#include<queue>
#include<cmath>
#include<math.h>
#include<string>
#include<string.h>
#include<map>
#include<set>
using LL = long long;
using namespace std;

const int M = 2e5 + 100;
const int mod = 1e9 + 7;

vector<int> e[M];
int n, sz[M], ans[M];
int fac[M], inv[M], inv2[M];

LL mul(LL a, LL b) {
	return a * b % mod;
}

void init() 
{
	fac[0] = 1;
	for (int i = 1; i <= n; i++) 
		fac[i] = mul(fac[i - 1], i);
	inv[1] = 1;
	for (int i = 2; i <= n; i++) 
		//常规公式为inv[i] = mul((- mod / i), inv[mod % i]),为使得逆元为正数,分子加上mod仍有i*inv[i] = 1(mod p)
		inv[i] = mul((mod - mod / i), inv[mod % i]);

}

void dfs(int u, int fa) 
{
	sz[u] = 1;
	for (int v : e[u]) 
	{
		if (v != fa) 
		{
			dfs(v, u);
			sz[u] += sz[v];
		}
	}
}

void reroot(int u, int fa) 
{
	ans[u] = mul(mul(ans[fa], sz[u]), inv[n - sz[u]]);
	for (int v : e[u]) 
	{
		if (v != fa) 
		{
			reroot(v, u);
		}
	}
}

int main() {
	cin >> n;
 	init();
	for (int i = 0; i < n - 1; i++) 
	{
		int u, v; 
		cin >> u >> v;
		e[u].push_back(v);
		e[v].push_back(u);
	}
	dfs(1, 0);
	ans[1] = fac[n];
	for (int i = 1; i <= n; i++) 
		ans[1] = mul(ans[1], inv[sz[i]]);
	for (int v : e[1]) 
		reroot(v, 1);
	for (int i = 1; i <= n; i++) 
		cout << ans[i] << "\n";
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值