#include<bits/stdc++.h>
#define LL long long
#define DB double
int MOD;
#define ls(x) x << 1
#define rs(x) x << 1 | 1
#define lowbit(x) x & (-x)
#define PII pair<int, int>
#define MP make_pair
#define VI vector<int>
#define VII vector<int>::iterator
#define all(x) x.begin(), x.end()
#define EB emplace_back
#define SI set<int>
#define SII set<int>::iterator
#define QI queue<int>
using namespace std;
template<typename T> void chkmn(T &a, const T &b) { (a > b) && (a = b); }
template<typename T> void chkmx(T &a, const T &b) { (a < b) && (a = b); }
int inc(const int &a, const int &b) { return a + b >= MOD ? a + b - MOD : a + b; }
int dec(const int &a, const int &b) { return a - b < 0 ? a - b + MOD : a - b; }
int mul(const int &a, const int &b) { return 1LL * a * b % MOD; }
int sqr(const int &a) { return 1LL * a * a % MOD; }
void Inc(int &a, const int &b) { ((a += b) >= MOD) && (a -= MOD); }
void Dec(int &a, const int &b) { ((a -= b) < 0) && (a += MOD); }
void Mul(int &a, const int &b) { a = 1LL * a * b % MOD; }
void Sqr(int &a) { a = 1LL * a * a % MOD; }
int qwqmi(int x, int k = MOD - 2)
{
int res = 1;
while(k)
{
if(k & 1) Mul(res, x);
Sqr(x), k >>= 1;
}
return res;
}
const int N = 6e3 + 5;
int n;
int f[N][N << 1], ans;
int main()
{
scanf("%d %d", &n, &MOD);
n = n * 3;
f[0][N] = 1;
for(int i = 0; i < n; ++i)
for(int j = N - i; j <= N + i; ++j)
{
Inc(f[i + 1][j + 1], f[i][j]);
Inc(f[i + 2][j - 1], mul(f[i][j], i + 1));
Inc(f[i + 3][j], mul(f[i][j], mul(i + 1, i + 2)));
}
for(int i = N; i < (N << 1); ++i)
Inc(ans, f[n][i]);
printf("%d\n", ans);
return 0;
}