#include <bits/stdc++.h>
using namespace std;

const int MOD = 998244353;

long long modpow(long long x, long long p){
	if(p == 0)return 1;
	if(p % 2 == 1)return x * modpow(x, p - 1) % MOD;
	
	long long r = modpow(x, p / 2);
	return r * r % MOD;
}
long long modinv(long long x){
	return modpow(x, MOD - 2);
}

int main(){
	int n, m;
	scanf("%d %d", &n, &m);
	
	int k = 3 * m + n;
	vector <long long> fact(k + 1), invfact(k + 1);
	
	fact[0] = 1;
	for(int i = 1; i <= k; i++)
		fact[i] = i * fact[i - 1] % MOD;
		
	invfact[k] = modinv(fact[k]);
	for(int i = k - 1; i >= 0; i--)
		invfact[i] = invfact[i + 1] * (i + 1) % MOD;

	auto C = [&](int n, int r) -> long long{
				return fact[n] * invfact[r] % MOD * invfact[n - r] % MOD;
			};
	
	long long res = 0;
	
	for(int odd = 0; odd <= min(n, m); odd++){
		if(odd % 2 != m % 2)continue;
		res = (res + C(n, odd) * C((3 * m - odd) / 2 + n - 1, n - 1)) % MOD;
	}

	for(int mx = 2 * m + 1; mx <= 3 * m; mx++)
		res = (res + (MOD - n) * C(3 * m - mx + n - 2, n - 2)) % MOD;
	
	printf("%lld\n", res);
	
	return 0;
}
