#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]; __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; }