#include <algorithm>
#include <iostream>
constexpr int N(2e3+10);
using ull=unsigned long long;
using L=__uint128_t;
struct FM
{
	ull b,m;
	FM(ull p):b(p),m((L(1)<<64)/p){}
	ull R(ull t)
	{
		ull a((L(t)*m)>>64);
		ull c=t-a*b;
		return c>=b?c-b:c;
	}
}F(2);
ull f[N][N],h[N][N][2],g[N][N][2][2][2][2];
int main()
{
	std::ios::sync_with_stdio(false);int n,m,p;
	std::cin.tie(nullptr);std::cin>>n>>m>>p;F=FM(p);
	f[0][0]=1;
	for(int i(0);i<=n;++i) for(int j(0);j<=m;++j)
	{
		f[i][j]=F.R(f[i][j]);for(int k(0);k<2;++k){h[i][j][k]=F.R(h[i][j][k]);}
		ull v1(1ull*f[i][j]*(n-i-1)),v2(1ull*f[i][j]*(m-j));
		f[i+1][j]+=1ull*f[i][j]*(i+j+1)+h[i][j][0]+h[i][j][1];
		h[i+1][j][0]+=v1*(i+j+1)+1ull*(h[i][j][0]+h[i][j][1])*(n-i-1);
		h[i][j+1][1]+=v2*i+1ull*h[i][j][0]*(m-j);
		g[i+1][j][0][0][1][0]+=v1;
		g[i][j+1][1][1][0][1]+=v2;
		for(int k(0);k<2;++k) for(int t(0);t<2;++t) for(int l(0);l<2;++l) for(int r(0);r<2;++r)
		{
			ull a(F.R(g[i][j][k][t][l][r])),b(1ull*a*(n-i-1)),c(1ull*a*(m-j));
			if(!((l^1)&r)) f[i+1][j]+=a,h[i+1][j][0]+=b;
			if(!k && !t && !(l&(r^1))) h[i][j+1][1]+=c;
			g[i+1][j][k][0][l^1][r]+=b;
			if(!t) g[i][j+1][k][1][l][r^1]+=c;
		}
	}
	int ans(0);
	for(int i(m),t(1);i>=0;--i,t=F.R(1ull*t*n)) ans=F.R(ans+1ull*f[n][i]*t);
	std::cout<<ans<<std::endl;
	return 0;
}