C,X=map(int,input().split()) P=998244353 Y=X+1 D=[0]*Y I=range(Y) for i in I: if i:D[i]=[(D[i-1][x-1]*(C+i+1-x)+D[i-1][x]*x)*pow(C+i,-1,P)%P for x in I] else:D[i]=[1]+[0]*X print([sum(D[max(0,X-C)][x]*(X-x)*pow(Y-x,-1,P)for x in I)%P,pow(C+1,-1,P)][X<1])