题意
给一棵树,要求对于任意的k,输出在树上任选k个点构成的虚树大小的和。
n<=200000
分析
考虑每个点对ans[k]的贡献。
容斥一下可以得到一个点的贡献就是
Ckn−∑x∈sonCksize[x]
C
n
k
−
∑
x
∈
s
o
n
C
s
i
z
e
[
x
]
k
其中son表示以当前点为根时的儿子集合。
那么有
ans[k]=∑i=1n(Ckn−∑x∈son[i]Cksize[x])
a
n
s
[
k
]
=
∑
i
=
1
n
(
C
n
k
−
∑
x
∈
s
o
n
[
i
]
C
s
i
z
e
[
x
]
k
)
设a[r]表示 Ckr C r k 在ans[k]中的系数,不难发现对于每一个k,他们的a数组都是一样的。
那么
ans[k]=∑i=kna[i]∗Cki
a
n
s
[
k
]
=
∑
i
=
k
n
a
[
i
]
∗
C
i
k
ans[k]=1k!∗∑i=kna[i]∗i!(i−k)!
a
n
s
[
k
]
=
1
k
!
∗
∑
i
=
k
n
a
[
i
]
∗
i
!
(
i
−
k
)
!
一开始我还以为这东西不能卷积,但实际上是可以的。
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int N=200005;
const int MOD=924844033;
int n,cnt,last[N],size[N],a[N*4],b[N*4],rev[N*4],L,jc[N],ny[N];
struct edge{int to,next;}e[N*2];
int read()
{
int x=0,f=1;char ch=getchar();
while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void addedge(int u,int v)
{
e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt;
}
void dfs(int x,int fa)
{
size[x]=1;a[n]++;
for (int i=last[x];i;i=e[i].next)
{
if (e[i].to==fa) continue;
dfs(e[i].to,x);
size[x]+=size[e[i].to];
a[size[e[i].to]]--;
}
if (n-size[x]) a[n-size[x]]--;
}
int C(int n,int m)
{
return (LL)jc[n]*ny[m]%MOD*ny[n-m]%MOD;
}
int ksm(int x,int y)
{
int ans=1;
while (y)
{
if (y&1) ans=(LL)ans*x%MOD;
x=(LL)x*x%MOD;y>>=1;
}
return ans;
}
void NTT(int *a,int f)
{
for (int i=0;i<L;i++) if (i<rev[i]) swap(a[i],a[rev[i]]);
for (int i=1;i<L;i<<=1)
{
int wn=ksm(5,f==1?(MOD-1)/i/2:MOD-1-(MOD-1)/i/2);
for (int j=0;j<L;j+=(i<<1))
{
int w=1;
for (int k=0;k<i;k++)
{
int u=a[j+k],v=(LL)a[j+k+i]*w%MOD;
a[j+k]=u+v;a[j+k]-=a[j+k]>=MOD?MOD:0;
a[j+k+i]=u+MOD-v;a[j+k+i]-=a[j+k+i]>=MOD?MOD:0;
w=(LL)w*wn%MOD;
}
}
}
int ny=ksm(L,MOD-2);
if (f==-1) for (int i=0;i<L;i++) a[i]=(LL)a[i]*ny%MOD;
}
int main()
{
n=read();
for (int i=1;i<n;i++)
{
int x=read(),y=read();
addedge(x,y);
}
jc[0]=jc[1]=ny[0]=ny[1]=1;
for (int i=2;i<=n;i++) jc[i]=(LL)jc[i-1]*i%MOD,ny[i]=(LL)(MOD-MOD/i)*ny[MOD%i]%MOD;
for (int i=2;i<=n;i++) ny[i]=(LL)ny[i-1]*ny[i]%MOD;
dfs(1,0);
for (int i=0;i<=n;i++) b[i]=ny[i],a[i]+=a[i]<0?MOD:0,a[i]=(LL)a[i]*jc[i]%MOD;
reverse(b,b+n+1);
int lg=0;
for (L=1;L<=n*2;L<<=1,lg++);
for (int i=0;i<L;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<(lg-1));
NTT(a,1);NTT(b,1);
for (int i=0;i<L;i++) a[i]=(LL)a[i]*b[i]%MOD;
NTT(a,-1);
for (int i=1;i<=n;i++) printf("%d\n",(LL)ny[i]*a[n+i]%MOD);
return 0;
}