#include<bits/stdc++.h>
using namespace std;
const int N=6060;
int f[N][N*2],mod,n,m,ans;
int main(){
	scanf("%d%d",&n,&mod);m=n;n*=3;f[n][N]=1;
	for(int i=n;i;i--){
		for(int j=-(n-i);j<=(n-i);j++){
			if(i>=3)f[i-3][j+N]=(f[i-3][j+N]+1ll*f[i][j+N]*(i-1)%mod*(i-2))%mod;
			if(i>=2)f[i-2][j+1+N]=(f[i-2][j+1+N]+1ll*f[i][j+N]*(i-1))%mod;
			f[i-1][j-1+N]=(f[i-1][j-1+N]+f[i][j+N])%mod;
		}
	}
	for(int i=-n;i<=0;i++)if(i%3==0)(ans+=f[0][i+N])%=mod;
	cout<<ans;
	return 0;
}