CF1554E You

在到校两个小时二十分钟后，我终于开始决定写题。

我们考虑实际上这个 \(a_i\) 数组，必然与我们直接定向每一条边之后每个点的出度数组一一对应。

然后就去吃饭了好耶！

在到校四个小时四十分钟后，我终于开始继续写题。

我们再来考虑如何求出这个 \(\gcd(a_1,a_2,\cdots,a_n)=k\) 。

感觉可以考虑莫反，我们令 \(g(i)\) 表示恰好 \(\gcd\) 为 \(i\) 的方案数，\(f(i)\) 表示 \(\gcd\) 为 \(i\) 的倍数的方案数。

\[f(k)=\sum_{k|d}g(d)\\ g(k)=\sum_{k|d}\mu(\frac{d}{k})f(d)\\ \]

我们考虑 \(f(d)\) 怎么算。我们发现，实际上我们只需要满足每一棵子树中的 \(siz_u\) 加上该子树连向父亲的边的对该子树根的出度贡献是 \(d\) 的倍数，就是一种方案。

我们就直接考虑算出所有的可能的 \(d\) 的 \(f(d)\) ，然后用调和级数来计算 \(g(k)\) 就行了。

我们再来考虑可能的 \(d\) 的范围，实际上就是 \(n-1\) 的因数。我们考虑暴力枚举因数然后判断就行了吧。

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+5;
const int MOD=998244353;
int ADD(int x,int y){return x+y>=MOD?x+y-MOD:x+y;}
int TIME(int x,int y){return (int)(1ll*x*y%MOD);}
int n,siz[N],f[N],g[N];
struct Edge{int nxt,to;}e[N<<1];int fir[N];
void add(int u,int v,int i){e[i]=(Edge){fir[u],v},fir[u]=i;}
void dfs(int u,int fa){
	siz[u]=1;
	for(int i=fir[u];i;i=e[i].nxt){
		int v=e[i].to;if(v==fa) continue;
		dfs(v,u),siz[u]+=siz[v];
	}
}
int check(int k){
	int res=1;
	for(int i=2;i<=n;++i)
	res=TIME(res,(siz[i]%k==0)+((siz[i]-1)%k==0));
	return res;
}
int mu[N];bool tag[N];
vector<int> pri;
void init(){
	mu[1]=1;
	for(int i=2;i<N;++i){
		if(!tag[i]) pri.push_back(i),mu[i]=MOD-1;
		for(int j=0;j<(int)pri.size();++j){
			if(i*pri[j]>=N) break;
			tag[i*pri[j]]=true;
			if(i%pri[j]==0){mu[i*pri[j]]=0;break;}
			else mu[i*pri[j]]=TIME(mu[i],mu[pri[j]]);
		}
	}
}
void solve(){
	cin>>n;
	for(int i=1;i<=n;++i) fir[i]=0;
	for(int i=1;i<n;++i){
		int u,v;scanf("%d%d",&u,&v);
		add(u,v,i<<1),add(v,u,i<<1|1);
	}
	dfs(1,0);
	for(int i=1;i<=n;++i) f[i]=0,g[i]=0;
	for(int i=1;i*i<=n-1;++i){
		if((n-1)%i) continue;
		f[(n-1)/i]=ADD(f[(n-1)/i],check((n-1)/i));
		if(i*i!=n-1) f[i]=ADD(f[i],check(i));
	}
	for(int i=1;i<=n;++i){
		for(int j=i;j<=n;j+=i)
		g[i]=ADD(g[i],TIME(mu[j/i],f[j]));
	}
	for(int i=1;i<=n;++i) printf("%d ",g[i]);
	return printf("\n"),void();
}
int main(){
	init();int T;cin>>T;
	while(T--) solve();
	return 0;
}