import math
ans=0
M=998244353
x,y=map(int,input().split())
n=min(x,y)
D=[0]+[1]*n
P=[0]+[1]*n
def s(i):
    return i*(i+1)//2
for i in range(2,n+1):
    if P[i]:
        for j in range(i,n+1,i):
            P[j]=0
            if j%(i**2)==0:
                D[j]=0
            else:
                D[j]*=-1
ans=s(x)*s(y)
for i in range(2,n+1):
    ans+=s(x//i)*s(y//i)*D[i]*(i**2)
print(ans%M)