#include<bits/stdc++.h>
#define MO 998244353
typedef long long ll;
const int maxn = 3000000;

ll n,m,ans;
ll inv[maxn],fac[maxn];

ll C(ll n, ll m)
{
	if (n < m) return 0;
	return fac[n]*inv[m]%MO*inv[n-m]%MO;
}
ll f(ll n, ll m, ll k)
{
	ll ret = 0;
	for (ll p=0; p<=k; p++)
		if (((m-p)&1)==0){
			ret = (ret+C((m-p)/2ll+n-1, n-1)*C(n, p)%MO)%MO;
		}
	return ret;
}
int main()
{
	scanf("%lld%lld",&n,&m);
	fac[0] = inv[0] = inv[1] = 1;
	for (int i=2; i<=n+3*m; i++)
		inv[i] = (MO-1ll*MO/i*inv[MO%i]%MO);
	for (int i=1; i<=n+3*m; i++)
		inv[i] = inv[i-1]*inv[i]%MO, fac[i] = fac[i-1]*i%MO;
	ans = (f(n, 3*m, m)-n*(f(n, m-1, m))%MO+MO)%MO;
	printf("%lld\n",ans);
	return 0;
}