﻿import math as M
import numpy as N

end=200

for n in range(end):
    m=N.sqrt(N.array(M.factorial(n)+1))
    if int(m)==m and M.factorial(n)+1==m**2:
        print(n,'&',m)