from sys import stdin input=lambda :stdin.readline()[:-1] mod=998244353 n,m=map(int,input().split()) def divisors(n): lower,upper=[],[] i=1 while i*i<=n: if n%i==0: lower.append(i) if i!=n//i: upper.append(n//i) i+=1 return lower+upper[::-1] xy=[] if m==0: for i in range(-n,n+1): xy.append((i,i)) else: for d in divisors(abs(m)): if (d+m//d)%2==0: # x + y = d, x - y = m // d x=(d+(m//d))//2 y=(d-(m//d))//2 if abs(x)<=n and abs(y)<=n: xy.append((x,y)) xy.append((-x,-y)) mod=998244353 table_size=2*10**5 fac=[1]*(table_size+1) finv=[1]*(table_size+1) for i in range(2,table_size+1): fac[i]=fac[i-1]*i%mod finv[table_size]=pow(fac[table_size],mod-2,mod) for i in range(table_size-1,-1,-1): finv[i]=finv[i+1]*(i+1)%mod def rebuild(n): global table_size,fac,finv fac+=[0]*(n-table_size) fac+=[0]*(n-table_size) finv+=[0]*(n-table_size) for i in range(table_size+1,n+1): fac[i]=fac[i-1]*i%mod finv[n]=inv(fac[n]) for i in range(n-1,table_size,-1): finv[i]=finv[i+1]*(i+1)%mod table_size=n def binom(n,k): if n<0 or k<0: return 0 if k>n: return 0 if n>table_size: rebuild(n+10**4) return (fac[n]*finv[k]%mod)*finv[n-k]%mod def fpow(x,k): res=1 while k: if k&1: res=res*x%mod x=x*x%mod k>>=1 return res def inv(a): if a