#include<bits/stdc++.h>
using namespace std;
using i64=long long;
using u32=unsigned;
using poly=array<u32,32>;
poly mul(poly f,poly g){
   poly h{};
   for(int i=0;i<32;i++)for(int j=0;i+j<32;j++)
   h[i+j]+=f[i]*g[j];
   return h;
}
inline i64 calc(poly f,u32 x){
   u32 res=0;
   for(int i=31;i>=0;i--)res=res*x+f[i];
   return res;
}
u32 binom[32][32];
poly comp(poly f,u32 c){
   poly g{};
   u32 pw[32]{1};for(int i=1;i<32;i++)pw[i]=pw[i-1]*c;
   for(int i=0;i<32;i++)for(int j=0;j<=i;j++)
   g[j]+=f[i]*binom[i][j]*pw[i-j];
   return g;
}
poly f[64];
inline i64 calct(i64 n){
   i64 res=0;for(;n;res+=n/=2);
   return res;
}
inline u32 inv(u32 x){
   u32 r=1;for(int t=5;t--;r*=2-r*x);
   return r;
}
u32 fac(i64 n){
   u32 res=1;
   for(;n;n/=2){
      u32 w=0;
      for(int i=63;i>=0;i--)if(n>>i&1)res*=calc(f[i],w),w+=1ull<<i;
   }return res;
}
int main(){
   for(int i=0;i<32;i++)binom[i][0]=1;
   for(int i=1;i<32;i++)for(int j=1;j<32;j++)binom[i][j]=binom[i-1][j]+binom[i-1][j-1];
   f[0]=f[1]={1,1};for(int i=2;i<64;i++){
      f[i]=mul(f[i-1],comp(f[i-1],1ull<<i-1));
   }
   i64 n,m;cin>>n>>m;
   i64 t=calct(n)-calct(m)-calct(n-m);u32 val=fac(n)*inv(fac(m))*inv(fac(n-m));
   cout<<(t>=32?0:val<<t);
   return 0;
}