#include #include #include using namespace atcoder; using mint = modint998244353; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf32 1000000001 #define Inf64 1000000000000000001 long long N,K; long long En = 1; vector la,lb; long long L; map mp[25]; map mp2[25]; void dfs(int ci,int ck,long long cx,__int128 cv){ if(cx > L)return; if(ck>=la.size()){ __int128 t = 1; rep(i,K){ t *= i+1; } mp[ci][cx] += t/cv; return; } __int128 t = 1LL; rep(i,K-ci+1){ if(i==0){ dfs(ci,ck+1,cx,cv); } else{ long long ncx = cx; ncx += (L / la[ck])*i; dfs(ci+i,ck+1,ncx,cv*t); } t *= i+1; } } int main(){ cin>>K>>N; rep(i,K)En *= i+1; L = 1; rep(i,N)L = lcm(L,i+1); rep(i,N){ if(i%2)la.push_back(i+1); else lb.push_back(i+1); } dfs(0,0,0,1); swap(la,lb); swap(mp,mp2); dfs(0,0,0,1); __int128 ans =0; rep(i,K+1){ for(auto x:mp[i]){ __int128 A = x.first,B = x.second; __int128 X = B; __int128 Y = mp2[K-i][L-A]; __int128 rem = 1LL; rep(j,K+1){ if(j==0)continue; int cc = j; if(X%j==0){ X /= j; continue; } if(Y%j==0){ Y /= j; continue; } rem *= j; } X *= Y; X /= rem; ans += X; } } long long aa = ans; cout<