Submission #1679455


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<=3*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 0
Code Size 2312 Byte
Status RE
Exec Time 249 ms
Memory 37504 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 0 / 1900
Status
AC × 3
AC × 13
RE × 35
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 RE 249 ms 34040 KB
doublestar1 RE 194 ms 33528 KB
doublestar2 RE 193 ms 33912 KB
doublestar3 RE 193 ms 33528 KB
doublestar4 RE 193 ms 33528 KB
example0 AC 6 ms 12544 KB
example1 AC 6 ms 12544 KB
example2 AC 6 ms 12544 KB
line0 RE 207 ms 37504 KB
line1 RE 207 ms 37504 KB
line2 RE 208 ms 37504 KB
line3 RE 206 ms 37504 KB
line4 RE 209 ms 37504 KB
maxrand0 RE 231 ms 33408 KB
maxrand1 RE 220 ms 33792 KB
maxrand10 RE 214 ms 33792 KB
maxrand11 RE 215 ms 34304 KB
maxrand12 RE 214 ms 34048 KB
maxrand13 RE 216 ms 33408 KB
maxrand14 RE 235 ms 35072 KB
maxrand15 RE 222 ms 33408 KB
maxrand16 RE 237 ms 33408 KB
maxrand17 RE 215 ms 34176 KB
maxrand18 RE 214 ms 33536 KB
maxrand19 RE 214 ms 33408 KB
maxrand2 RE 213 ms 33408 KB
maxrand3 RE 221 ms 33408 KB
maxrand4 RE 220 ms 34944 KB
maxrand5 RE 227 ms 33408 KB
maxrand6 RE 214 ms 33408 KB
maxrand7 RE 213 ms 33408 KB
maxrand8 RE 213 ms 33408 KB
maxrand9 RE 218 ms 33792 KB
rand0 AC 11 ms 14976 KB
rand1 AC 6 ms 12544 KB
rand2 AC 9 ms 14720 KB
rand3 AC 10 ms 14848 KB
rand4 AC 7 ms 12672 KB
rand5 AC 9 ms 14848 KB
rand6 AC 9 ms 14720 KB
rand7 AC 9 ms 14848 KB
rand8 AC 6 ms 12544 KB
rand9 AC 7 ms 12672 KB
star0 RE 227 ms 34804 KB
star1 RE 186 ms 33908 KB
star2 RE 187 ms 34036 KB
star3 RE 188 ms 34036 KB
star4 RE 187 ms 34036 KB