#include using namespace __gnu_pbds; #pragma GCC optimize("Ofast,unroll-loops") #include using namespace std; #include using namespace atcoder; using ll = long long int; using ull = unsigned long long int; using ld = long double; constexpr ll MAX = 2000000000000000000; constexpr ld PI = 3.14159265358979; constexpr ll MOD = 0;//2024948111; ld dotorad(ld K){ return PI * K / 180.0; } ld radtodo(ld K){ return K * 180.0 / PI; } mt19937 mt; void randinit(){ srand((unsigned)time(NULL));mt = mt19937(rand()); } __int128_t modinv(__int128_t A,__int128_t M){ //A*r+B*n=gcd(A,B) A %= M; if(A == 0 || __gcd(A,M) != 1){ //cout << "Error modinv(" << A << "," << M << ")" << endl; return -1; } __int128_t B = M,U = 1,V = 0; while(B){ __int128_t T = A / B; A -= T * B; swap(A,B); U -= T * V; swap(U,V); } U %= M; if(U < 0) U += M; return U; } int main(){ ll N,M; cin >> N >> M; vector pow3(60,1),pow4(60,1),pow6(60,1),pow2(60,1); for(ll i = 1;i < 60;i++){ pow2[i] = pow2[i - 1] * 2; pow3[i] = pow3[i - 1] * 3; pow4[i] = pow4[i - 1] * 4; pow6[i] = pow6[i - 1] * 6; } ll l = pow2[N * 2] * pow3[N]; if(N >= 13){ __int128_t t1 = (modinv(pow3[N],pow2[N * 2]) * pow3[N]) % l; __int128_t t2 = (modinv(pow2[N * 2],pow3[N]) * pow2[N * 2]) % l; ll ID = 0; vector memo2(pow3[N - 1],0); for(ll bb = 0;bb < pow3[N];bb += 3){ ll c = bb; ll p4 = 0,p6 = 0; for(ll i = 0;i < N;i++){ ll x = c % 3; c /= 3; p4 += x * pow4[i]; p6 += x * pow6[i]; } ll x = (bb * t2 + p4 * t1) % l; ll zurasi = p6 - x; zurasi %= pow6[N]; if(zurasi < 0) zurasi += pow6[N]; memo2[bb / 3] = (zurasi); } gp_hash_table X; ll ans = 0; for(ll id = 0;id < pow3[N - 1];id++){ ll bb = 3 * id; ll z = memo2[id]; if(X.find(z) != X.end()){ ll z1 = X[z]; ll z2 = bb; if(z2 < z1) swap(z1,z2); ll d1 = (z2 - (z1 + 3)) % l; if(d1 >= M) ans++; ll d2 = (z1 - (z2 + 3)) % l; if(d2 < 0) d2 += l; if(d2 >= M) ans++; X.erase(z); continue; } else{ X[z] = bb; } } for(auto [z,bb] : X){ if(M <= l - 3) ans++; } cout << ans << endl; } else{ map> A,B; vector Z; for(ll bb = 0;bb < pow3[N];bb++){ ll c = bb; ll p3 = 0,p4 = 0,p6 = 0; for(ll i = 0;i < N;i++){ ll x = c % 3; c /= 3; p3 += x * pow3[i]; p4 += x * pow4[i]; p6 += x * pow6[i]; } tuple P = crt({p3,p4},{pow3[N],pow4[N]}); ll x = get<0>(P),m = get<1>(P); ll zurasi = p6 - x; zurasi %= pow6[N]; if(zurasi < 0) zurasi += pow6[N]; Z.emplace_back(zurasi); //assert(m == l); A[zurasi].emplace_back(x); A[zurasi].emplace_back(x + l); A[zurasi].emplace_back(x + l + l); x -= M; x %= l; if(x < 0) x += l; B[zurasi].emplace_back(x); } sort(Z.begin(),Z.end()); Z.erase(unique(Z.begin(),Z.end()),Z.end()); ll ans = 0; for(ll z : Z){ vector a = A[z],b = B[z]; sort(a.begin(),a.end()); sort(b.begin(),b.end()); for(ll i = 0;i < b.size();i++){ ll x = b[i]; auto it = lower_bound(a.begin(),a.end(),x); if(*it == x + M) ans++; } } cout << ans << endl; } }