#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=2e7+10,mod=998244353;
int n,m,fac[maxn],ifac[maxn],F[maxn];
int qpow(int x,int y){
    int res=1;
    for(;y;y>>=1){
        if(y&1)res=res*x%mod;
        x=x*x%mod;
    }
    return res;
}
void init(){
    fac[0]=1;
    for(int i=1;i<maxn;i++)fac[i]=fac[i-1]*i%mod;
    ifac[maxn-1]=qpow(fac[maxn-1],mod-2);
    for(int i=maxn-2;i>=0;i--)ifac[i]=ifac[i+1]*(i+1)%mod;
}
int C(int x,int y){
    if(x<y)return 0;
    return fac[x]*ifac[y]%mod*ifac[x-y]%mod;
}
signed main(){
    cin>>n>>m;
    init();
    int ans=C(n+m-1,n);
    for(int k=1;k*(3*k-1)/2<=m-1;k++){
        if(k&1){
            int s=k*(3*k-1)/2;
            ans=(ans-C(n+m-1-s,n)+mod)%mod;
            s=k*(3*k+1)/2;
            ans=(ans-C(n+m-1-s,n)+mod)%mod;
        } else {
            int s=k*(3*k-1)/2;
            ans=(ans+C(n+m-1-s,n))%mod;
            s=k*(3*k+1)/2;
            ans=(ans+C(n+m-1-s,n))%mod;
        }
    }
    cout<<ans*ifac[n]%mod;
    return 0;
}