from math import lcm
from atcoder.convolution import convolution

N,M=map(int,input().split())
ans=0
mod=998244353
for a in range(1,N+1):
    for b in range(1,M+1):
        ans+=lcm(a,b)
        ans%=mod
print(ans)