rt。

样例过不去，调了好久了，哪个大佬可以帮忙看看啊！

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#define ll long long
//#define int ll
using namespace std;
const int N=6e5+10;
const int mod=1004535809,gen=3,invg=334845270;
inline int read(){
	register int d=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){d=(d<<1)+(d<<3)+(ch^48);ch=getchar();}
	return d*f;
}
inline int qpow(int a,int b,int p){
	int res=1;
	while(b){
		if(b&1) res=1ll*res*a%p;
		a=1ll*a*a%p;
		b>>=1;
	}
	return res;
}
int n,tr[N],W[N],f[N],g[N],h[N];
inline int init(int n){
	int Limit=1,len=0,up=(n<<1);
	while(Limit<up) Limit<<=1,++len;
	for(int i=0;i<Limit;i++) tr[i]=(tr[i>>1]>>1|((i&1)<<len-1));
	return Limit;
}
inline void NTT(int *f,int n,bool opt){
	for(int i=0;i<n;i++) if(i<tr[i])swap(f[i],f[tr[i]]);
	for(int mid=1;mid<n;mid<<=1){
		int wn=qpow(opt?gen:invg,(mod-1)/(mid<<1),mod);
		W[0]=1;for(int i=1;i<mid;i++) W[i]=1ll*W[i-1]*wn%mod;
		for(int len=mid<<1,j=0;j<n;j+=len)
			for(int k=0;k<mid;k++){
				int x=f[j+k],y=1ll*f[j+mid+k]*W[k]%mod;
				f[j+k]=(1ll*x+y)%mod,f[j+mid+k]=(1ll*x-y+mod)%mod;
			}
	}
	if(!opt){
		int invn=qpow(n,mod-2,mod);
		for(int i=0;i<n;i++) f[i]=1ll*f[i]*invn%mod;
	}
}
inline void Mul(int *f,int *g,int *h,int n,int m){
	static int A[N],B[N];
	int len=n+m-1;
	int limit=init(len);
	for(int i=0;i<limit;i++) A[i]=B[i]=0;
	for(int i=0;i<n;i++) A[i]=f[i];
	for(int i=0;i<m;i++) B[i]=g[i];
	NTT(A,limit,1);NTT(B,limit,1);
	for(int i=0;i<limit;i++) A[i]=1ll*A[i]*B[i]%mod;
	NTT(A,limit,0);
	for(int i=0;i<len;i++) h[i]=A[i];
}
int d[N];
void Inv(int *f,int *g,int len){
	if(len==1){g[0]=qpow(f[0],mod-2,mod);return;}
	Inv(f,g,(len+1)>>1);
	int limit=init(len);
	for(int i=0;i<len;i++) d[i]=f[i];
	for(int i=len;i<limit;i++) d[i]=0;
	NTT(d,limit,1);NTT(g,limit,1);
	for(int i=0;i<limit;i++) g[i]=((2ll-1ll*g[i]*d[i]%mod)+mod)%mod*g[i]%mod;
	NTT(g,limit,0);
	for(int i=len;i<limit;i++) g[i]=0;
//	memset(d,0,sizeof(d));
}
int fac[N],ifac[N],ig[N];
signed main(){
	n=read()+1;
	fac[0]=1;for(int i=1;i<=n;i++) fac[i]=1ll*i*fac[i-1]%mod;
	ifac[n]=qpow(fac[n],mod-2,mod);ifac[0]=1;ifac[1]=1;
	for(int i=n-1;i>=2;i--) ifac[i]=1ll*(i+1)*ifac[i+1]%mod;
	g[0]=1;h[0]=0;
	for(int i=1;i<=n;i++){
		int tmp=qpow(2,(1ll*i*(i-1)/2)%(mod-1),mod);
		g[i]=1ll*tmp*ifac[i]%mod;h[i]=1ll*tmp*ifac[i-1]%mod;
//		printf("%d %d\n",g[i],h[i]);
	}
	Inv(g,ig,n);
	int limit=init(n);
//	for(int i=0;i<=n;i++){
//		int res=0;
//		for(int k=0;k<=i;k++){
//			res=(1ll*res+1ll*g[k]*ig[i-k]%mod)%mod;
//		}
//		 printf("%d\n",res);
//	}
//	for(int i=0;i<n;i++) printf("%d\n",ig[i]);
	NTT(ig,limit,1);NTT(h,limit,1);
	for(int i=0;i<=limit;i++) f[i]=1ll*ig[i]*h[i]%mod;
	NTT(f,limit,0);
	printf("%d",1ll*f[n-1]*(ifac[n-2])%mod);
	return 0;
}