from math import factorial
from math import sqrt
import time
start,end=0,1001
Time_start=time.time()

def is_square(m):
  x = m // 2
  e = set([x])
  while x * x != m:
    x = (x + (m // x)) // 2
    if x in e: return False
    e.add(x)
  return True


print("N & M")
for n in range(start,end):
    m=factorial(n)+1
    if is_square(m):
        print(n,'&',sqrt(m))

Time_end=time.time()
temps=Time_end-Time_start
minute=temps//60
H=minute//60
minute-=H*60
print('===========================')
print('Temps : %s heure %s mins '% (H,minute))