#pragma GCC optimize("O3") #include using namespace std; using ll=long long; using P=pair; template using V=vector; #define fi first #define se second #define all(v) (v).begin(),(v).end() const ll inf=(1e18); //const ll mod=998244353; // const ll mod=1000000007; const vector dy={-1,0,1,0},dx={0,-1,0,1}; ll GCD(ll a,ll b) {return b ? GCD(b,a%b):a;} ll LCM(ll c,ll d){return c/GCD(c,d)*d;} struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout< bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a=b; return 1; } return 0; } templatevoid debag(const vector &a){cerr<<"debag :";for(auto v:a)cerr<void print(const vector &a){for(auto v:a)cout<>=1ll; } return res; } mint inv()const{return pow(mod-2);} //拡張ユークリッドの互除法 /* mint inv()const{ int x,y; int g=extgcd(v,mod,x,y); assert(g==1); if(x<0)x+=mod; return mint(x); }*/ friend ostream& operator<<(ostream&os,const mint&val){ return os<(const mint&val)const{return v>val.v;} }; const ll MAX = 2000010;//設定 mint fac[MAX], finv[MAX], inv[MAX]; // テーブルを作る前処理 void init(){ fac[0] = fac[1] = 1; for(int i=1;i=0;i--)finv[i]=finv[i+1]*(i+1); for(int i=MAX-2;i>=1;i--)inv[i]=finv[i]+fac[i-1]; } //階乗 mint factor(ll n,ll k){ if (n>n>>mod; init(); V> dp(n+5,V(n+5,mint(0))); dp[0][0]=mint(1); cout<<1<<"\n"; V inv(n+5,mint(1)); for(int i=1;i<=n;i++)inv[i]/=mint(i); for(int i=2;i<=n;i++){ mint res=mint(i).pow(i-2); V> ndp=dp;//すでにとった個数,分けた個数,辺の本数 for(int j=0;jn)break; ans+=dp[j][i+j-1]; } cout<