#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
constexpr int dx[]={1,0,0,-1},dy[]={0,1,-1,0};
constexpr int Dx[]={1,1,0,-1,-1,-1,0,1},Dy[]={0,1,1,1,0,-1,-1,-1};
constexpr int mod=998244353,Mod=1e9+7,inf=Mod;
constexpr ll linf=(ll)inf*inf;
typedef pair<ll,int>P;
#define m_p make_pair
struct fastio{
	fastio(){
		cin.tie(0);
		ios::sync_with_stdio(0);
		cout<<fixed<<setprecision(20);
	};
}fio;
template<class T,class U>
bool chmax(T&a,const U&b){
	if(a<b){
		a=b;
		return 1;
	}
	return 0;
}
template<class T,class U>
bool chmin(T&a,const U&b){
	if(a>b){
		a=b;
		return 1;
	}
	return 0;
}
//*
#include<atcoder/all>
using namespace atcoder;
//*/
ll dp[2001][4001];
int main(){
	int n;
	cin>>n;
	dp[0][0]=1;
	for(int i=1;i<=n;i++){
		for(int j=0;j<=2*i;j++){
			if(j<2*i)
				(dp[i][j]+=dp[i-1][j]*(2*i-1-j)*(2*i-2-j)/2)%=mod;
			if(0<j&&j<2*i)
				(dp[i][j]+=dp[i-1][j-1]*((2*i-1-j)*(j-1)+1))%=mod;
			if(1<j)
				(dp[i][j]+=dp[i-1][j-2]*(j-1)*(j-2)/2)%=mod;
		}
	}
	for(int i=0;i<=2*n;i++)
		cout<<dp[n][i]<<endl;
}