#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=3000005,M=998244353;
int fac[N],inv[N],n,m,ans;
int ksm(int x,int y){
	int ans=1;
	for (;y;y>>=1,x=x*x%M)
		if (y&1)(ans*=x)%=M;
	return ans;	
}
int C(int x,int y){
	if (x<y)return 0;
	return fac[x]*inv[y]%M*inv[x-y]%M;
}
signed main(){
	scanf("%lld%lld",&n,&m);
	fac[0]=inv[0]=1;
	for (int i=1;i<N;i++)fac[i]=fac[i-1]*i%M,inv[i]=ksm(fac[i],M-2);
	for (int i=m&1;i<=m;i+=2)
		(ans+=C(n+(3*m-i)/2-1,n-1)*C(n,i))%=M;
	for (int i=m&1;i<=m;i+=2)
		(ans+=M-n*C(n+(m-i)/2-1,n-1)%M*C(n,i)%M)%=M;
	n--;
	for (int i=m&1;i<=m;i+=2)
		(ans+=(n+1)*C(n+(m-i)/2-1,n-1)%M*C(n,i)%M)%=M;		
	printf("%lld\n",ans);
	return 0;
} 