edofrhasa - 経済統計練習帳

import numpy as np
import sympy as sym
 
pr=[]
 
for i in range(100):
    m=sym.prime(i+1)
    pr.append(m)
 
p=np.random.choice(pr)
 
pr.remove(p)
 
q=np.random.choice(pr)
 
print(p)
print(q)
 
n=p*q
 
print(n)
 
lx=sym.lcm(p-1, q-1)
 
print(lx)