Submission #1679457


Source Code Expand

#include<bits/stdc++.h>
#define mo 924844033ll
using namespace std;
long long a[800010],b[800010],c[800010],ym[800010],fac[200010],nifac[200010];
vector<int> des[200010];
int n;
int getmi(int a,int x)
{
	int ans=1;
	while (x)
	{
		if (x&1) ans=(long long)ans*a%mo;
		a=(long long)a*a%mo;
		x>>=1;
	}
	return ans;
}
long long C(int n,int m)
{
	if ((n<0)||(m<0)||(n<m)) return 0;
	return (long long)fac[n]*nifac[m]%mo*nifac[n-m]%mo;
}
int dfs(int s,int pre=0)
{
	int ans=1;
	for (int k=0;k<des[s].size();k++) if (des[s][k]!=pre) ans+=dfs(des[s][k],s);
	if (pre) {a[ans]++;a[n-ans]++;}
	return ans;
}
void fft(long long a[],int n)
{
	if (n==1) return;
	long long l[n>>1],r[n>>1];
	for (int i=0;i<(n>>1);i++)
	{
		l[i]=a[i<<1];
		r[i]=a[(i<<1)+1];
	}
	fft(l,n>>1);fft(r,n>>1);
	long long wn=ym[n],w=1;
	for (int i=0;i<(n>>1);i++)
	{
		a[i]=(l[i]+w*r[i])%mo;
		a[i+(n>>1)]=(l[i]-w*r[i])%mo;
		w=w*wn%mo;
	}
}
int main()
{
	scanf("%d",&n);
	fac[0]=nifac[0]=1;
	for (int i=1;i<=n;i++)
	{
		fac[i]=(long long)fac[i-1]*i%mo;
		nifac[i]=getmi(fac[i],mo-2);
	}
	if (fac[n]*nifac[n]%mo!=1) assert(0);
	for (int i=1;i<n;i++)
	{
		int u,v;scanf("%d%d",&u,&v);
		des[u].push_back(v);des[v].push_back(u);
	}
	dfs(1);
	for (int i=0;i<=n;i++)
	{
		a[i]=(long long)a[i]*fac[i]%mo;
		c[i]=nifac[n-i];
	}
	//for (int i=0;i<=2*n;i++) cout<<i<<' '<<a[i]<<' '<<c[i]<<endl;
	/*for (int i=0;i<=n;i++)
	for (int j=0;j<=2*n;j++)
		b[i+j]=(b[i+j]+a[i]*c[j])%mo;*/
	//for (int i=0;i<=3*n;i++) cout<<i<<' '<<b[i]<<endl;
	//for (int i=1;i<=n;i++) printf("%lld\n",(mo-b[i+n]*nifac[i]%mo+n*C(n-1,i-1)+n*C(n-1,i))%mo);
	
	int mi=1;while (mi<=2*n+10) mi<<=1;
	for (int i=1;i<=mi;i<<=1) ym[i]=getmi(5,(mo-1)/i);
	fft(a,mi);fft(c,mi);
	for (int i=0;i<mi;i++) a[i]=a[i]*c[i]%mo;
	fft(a,mi);
	for (int i=1;i<(mi>>1);i++) swap(a[i],a[mi-i]);
	int ni=getmi(mi,mo-2);
	for (int i=0;i<mi;i++) a[i]=a[i]*ni%mo;
	for (int i=0;i<mi;i++) a[i]+=mo;
	//for (int i=0;i<=n;i++) cout<<a[i+n]*nifac[i]%mo<<endl;
	for (int i=1;i<=n;i++) printf("%lld\n",(mo-a[i+n]*nifac[i]%mo+n*C(n-1,i-1)+n*C(n-1,i))%mo);
	/*for (int i=1;i<=n;i++)
	{
		int ans=(n*C(n-1,i-1)+n*C(n-1,i))%mo;
		for (int j=0;j<=n;j++)
		{
			ans-=(a[j]+b[j])*C(j,i);
		}
		cout<<(ans%mo+mo)%mo<<endl;
	}*/
}

Submission Info

Submission Time
Task F - Many Easy Problems
User yx_lu
Language C++14 (GCC 5.4.1)
Score 1900
Code Size 2312 Byte
Status AC
Exec Time 358 ms
Memory 43648 KB

Compile Error

./Main.cpp: In function ‘int main()’:
./Main.cpp:50:16: warning: ignoring return value of ‘int scanf(const char*, ...)’, declared with attribute warn_unused_result [-Wunused-result]
  scanf("%d",&n);
                ^
./Main.cpp:60:30: warning: ignoring return value of ‘int scanf(const char*, ...)’, declared with attribute warn_unused_result [-Wunused-result]
   int u,v;scanf("%d%d",&u,&v);
                              ^

Judge Result

Set Name Sample All
Score / Max Score 0 / 0 1900 / 1900
Status
AC × 3
AC × 48
Set Name Test Cases
Sample example0, example1, example2
All doublestar0, doublestar1, doublestar2, doublestar3, doublestar4, example0, example1, example2, line0, line1, line2, line3, line4, maxrand0, maxrand1, maxrand10, maxrand11, maxrand12, maxrand13, maxrand14, maxrand15, maxrand16, maxrand17, maxrand18, maxrand19, maxrand2, maxrand3, maxrand4, maxrand5, maxrand6, maxrand7, maxrand8, maxrand9, rand0, rand1, rand2, rand3, rand4, rand5, rand6, rand7, rand8, rand9, star0, star1, star2, star3, star4
Case Name Status Exec Time Memory
doublestar0 AC 261 ms 41592 KB
doublestar1 AC 261 ms 41592 KB
doublestar2 AC 261 ms 41592 KB
doublestar3 AC 261 ms 41592 KB
doublestar4 AC 260 ms 41592 KB
example0 AC 6 ms 12544 KB
example1 AC 6 ms 12544 KB
example2 AC 6 ms 12544 KB
line0 AC 276 ms 43520 KB
line1 AC 276 ms 43520 KB
line2 AC 277 ms 43648 KB
line3 AC 277 ms 43648 KB
line4 AC 276 ms 43520 KB
maxrand0 AC 289 ms 41344 KB
maxrand1 AC 290 ms 41344 KB
maxrand10 AC 288 ms 41344 KB
maxrand11 AC 289 ms 41344 KB
maxrand12 AC 289 ms 41344 KB
maxrand13 AC 289 ms 41344 KB
maxrand14 AC 289 ms 41344 KB
maxrand15 AC 290 ms 41344 KB
maxrand16 AC 289 ms 41344 KB
maxrand17 AC 289 ms 41344 KB
maxrand18 AC 288 ms 41344 KB
maxrand19 AC 292 ms 41344 KB
maxrand2 AC 291 ms 41344 KB
maxrand3 AC 293 ms 41344 KB
maxrand4 AC 293 ms 41344 KB
maxrand5 AC 295 ms 41344 KB
maxrand6 AC 293 ms 41344 KB
maxrand7 AC 289 ms 41344 KB
maxrand8 AC 303 ms 41344 KB
maxrand9 AC 292 ms 41344 KB
rand0 AC 10 ms 14848 KB
rand1 AC 6 ms 12544 KB
rand2 AC 8 ms 12672 KB
rand3 AC 10 ms 14848 KB
rand4 AC 6 ms 12544 KB
rand5 AC 10 ms 14848 KB
rand6 AC 8 ms 12672 KB
rand7 AC 9 ms 14848 KB
rand8 AC 6 ms 12544 KB
rand9 AC 7 ms 12672 KB
star0 AC 358 ms 41972 KB
star1 AC 255 ms 41972 KB
star2 AC 256 ms 41972 KB
star3 AC 254 ms 41972 KB
star4 AC 255 ms 41972 KB