#include <fstream>
#define NMax 3005
#define modulo 666013
using namespace std;
ifstream f("perm.in");
ofstream g("perm.out");
int n, k, mat[3][NMax], i, j;
int main()
{
    f>>n>>k;
    mat[0][0]=1;
    
    for (i=2; i<=n; i++) {
        mat[1-i%2][0] = (mat[i%2][1] + mat[i%2][0]*(i-1)) % modulo;
        for (j=1; j<i; j++)
            mat[1-i%2][j] = (mat[i%2][j+1]*(j+1) + mat[i%2][j-1] + mat[i%2][j]*(i-j-1)) % modulo;
    }
    g<<mat[1][k]%modulo;
    return 0;
}