Pisano period: Difference between revisions

=={{header|Phix}}==

<lang Phix>function pisano_period(integer m)

-- Calculates the Pisano period of 'm' from first principles. (copied from Go)

integer p = 0, c = 1

for i=0 to m*m-1 do

{p, c} = {c, mod(p+c,m)}

if p == 0 and c == 1 then

return i + 1

end if

end for

return 1

end function

function pisanoPrime(integer p, k)

if not is_prime(p) or k=0 then ?9/0 end if

return power(p,k-1)*pisano_period(p)

end function

function pisano(integer m)

-- Calculates the Pisano period of 'm' using pisanoPrime.

sequence s = prime_factors(m, true, get_maxprime(m))&0,

pps = {}

integer k = 1, p = s[1]

for i=2 to length(s) do

integer n = s[i]

if n!=p then

pps = append(pps,pisanoPrime(p,k))

{k,p} = {1,n}

else

k += 1

end if

end for

return lcm(pps)

end function

procedure p(integer k, lim)

printf(1,"pisanoPrimes")

integer pdx = 1, c = 0

while true do

integer p = get_prime(pdx)

if p>=lim then exit end if

c += 1

if c=7 then puts(1,"\n ") c = 1

elsif pdx>1 then puts(1,", ") end if

printf(1,"(%3d,%d)=%3d",{p,k,pisanoPrime(p,k)})

pdx += 1

end while

printf(1,"\n")

end procedure

p(2,15)

p(1,180)

sequence p180 = {}

for n=2 to 180 do p180 &= pisano(n) end for

printf(1,"pisano(2..180):\n")

pp(p180,{pp_IntFmt,"%4d",pp_IntCh,false})</lang>

{{out}}

<pre>

pisanoPrimes( 2,2)= 6, ( 3,2)= 24, ( 5,2)=100, ( 7,2)=112, ( 11,2)=110, ( 13,2)=364

pisanoPrimes( 2,1)= 3, ( 3,1)= 8, ( 5,1)= 20, ( 7,1)= 16, ( 11,1)= 10, ( 13,1)= 28

( 17,1)= 36, ( 19,1)= 18, ( 23,1)= 48, ( 29,1)= 14, ( 31,1)= 30, ( 37,1)= 76

( 41,1)= 40, ( 43,1)= 88, ( 47,1)= 32, ( 53,1)=108, ( 59,1)= 58, ( 61,1)= 60

( 67,1)=136, ( 71,1)= 70, ( 73,1)=148, ( 79,1)= 78, ( 83,1)=168, ( 89,1)= 44

( 97,1)=196, (101,1)= 50, (103,1)=208, (107,1)= 72, (109,1)=108, (113,1)= 76

(127,1)=256, (131,1)=130, (137,1)=276, (139,1)= 46, (149,1)=148, (151,1)= 50

(157,1)=316, (163,1)=328, (167,1)=336, (173,1)=348, (179,1)=178

pisano(2..180):

{ 3, 8, 6, 20, 24, 16, 12, 24, 60, 10, 24, 28, 48, 40, 24,

36, 24, 18, 60, 16, 30, 48, 24, 100, 84, 72, 48, 14, 120, 30,

48, 40, 36, 80, 24, 76, 18, 56, 60, 40, 48, 88, 30, 120, 48,

32, 24, 112, 300, 72, 84, 108, 72, 20, 48, 72, 42, 58, 120, 60,

30, 48, 96, 140, 120, 136, 36, 48, 240, 70, 24, 148, 228, 200, 18,

80, 168, 78, 120, 216, 120, 168, 48, 180, 264, 56, 60, 44, 120, 112,

48, 120, 96, 180, 48, 196, 336, 120, 300, 50, 72, 208, 84, 80, 108,

72, 72, 108, 60, 152, 48, 76, 72, 240, 42, 168, 174, 144, 120, 110,

60, 40, 30, 500, 48, 256, 192, 88, 420, 130, 120, 144, 408, 360, 36,

276, 48, 46, 240, 32, 210, 140, 24, 140, 444, 112, 228, 148, 600, 50,

36, 72, 240, 60, 168, 316, 78, 216, 240, 48, 216, 328, 120, 40, 168,

336, 48, 364, 180, 72, 264, 348, 168, 400, 120, 232, 132, 178, 120}

</pre>