python3で

def gcd(m,n):
if n==0:
return m
return gcd(n,m%n)

def lcm(m,n):
g=gcd(m,n)
return (m/g)*n

a=1000000000000000123600000
b=100000000000000012367890
int(lcm(a,b))


で、casioのサイトと比較しても違う（多分桁が足りない）。gcdは合ってる。
何故かnumpyのlog10もデカすぎるのか不機嫌になってる（floatにすればOK)。