tps/AlgoFermatFacto.py~

   1 import math as m
   2
   3 def estUnCarre(k):
   4     return (int(m.sqrt(k)))**2 == k
   5
   6 def facto(n):
   7     t= int(m.sqrt(n))+1
   8     sp = t**2 - n
   9     while not estUnCarre(sp):
  10         t += 1
  11         sp = t**2 -n
  12         print "+",
  13     return(int(t+m.sqrt(sp)),int(t-m.sqrt(sp)))
  14
  15
  16
  17
  18
  19 print facto(5959)
  20 print facto(2953*3037)
  21 print facto(6197*6299)
  22 print facto(2953*3037)
  23 print facto(17729*17939)
  24
  25
  26 print facto(9623827)
  27 print facto(343570291)
  28