#include <fstream>
 
using namespace std;
 
ifstream in("kperm.in");
ofstream out("kperm.out");
 
typedef long long i64;
 
const i64 mod= 666013;
const i64 nmax= 5000;
 
i64 f[nmax+1];
 
i64 lgpow( i64 x, i64 y ) {
    i64 sol= 1;
    while ( y>0 ) {
        if ( y%2==1 ) {
            sol= (sol*x)%mod;
            --y;
        }
        x= (x*x)%mod;
        y/= 2;
    }
 
    return sol;
}
 
int main(  ) {
    f[0]= 1;
    for ( i64 i= 1; i<=nmax; ++i ) {
        f[i]= (f[i-1]*i)%mod;
    }
 
    i64 n, k;
    in>>n>>k;
 
    i64 sol= 0;
    if ( k%2==1 ) {
        sol= (f[k-n%k]*f[n%k]%mod)*(lgpow(f[n/k+1], n%k)*lgpow(f[n/k], k-n%k)%mod)%mod;
    }
    out<<sol<<"\n";
 
    return 0;
}