Juggler sequence: Difference between revisions

Line 482:

5812827 135 67 d(n) = 7996276 digits

7110201 205 119 d(n) = 89981517 digits

</pre>

=={{header|Phix}}==

with javascript_semantics

include mpfr.e

function juggler(integer n, bool bDigits=false)

assert(n>=1)

mpz a = mpz_init(n),

h = mpz_init(n)

integer l = 0,

hi = 0

while mpz_cmp_si(a,1)!=0 do

if mpz_odd(a) then

mpz_pow_ui(a,a,3)

end if

mpz_sqrt(a,a)

l += 1

if mpz_cmp(a,h)>0 then

mpz_set(h,a)

hi = l

end if

end while

atom hd = iff(bDigits?mpz_sizeinbase(h,10):mpz_get_atom(h))

return {n, l, hd, hi}

end function

procedure main()

printf(1," n l[n] h[n] i[n]\n")

printf(1,"-----------------------------------\n")

for n=20 to 39 do

printf(1,"%2d %2d %,18d %2d\n", juggler(n))

end for

sequence nums = {113, 173, 193, 2183, 11229, 15065, 15845, 30817}

-- alas mpz_sqrt() throws a wobbly over the rest:

-- 48443, 275485, 1267909, 2264915, 5812827, 7110201}

printf(1,"\n n l[n] d[n] i[n]\n")

printf(1,"-----------------------------\n")

-- ... and pwa/p2js(/mpfr.js) only copes with the first 3

for i=1 to iff(platform()=JS?3:length(nums)) do

printf(1,"%9d %3d %,6d %2d\n", juggler(nums[i],true))

end for

end procedure

main()

<pre>

n l[n] h[n] i[n]

-----------------------------------

20 3 20 0

21 9 140 4

22 3 22 0

23 9 110 1

24 3 24 0

25 11 52,214 3

26 6 36 3

27 6 140 1

28 6 36 3

29 9 156 1

30 6 36 3

31 6 172 1

32 6 36 3

33 8 2,598 2

34 6 36 3

35 8 2,978 2

36 3 36 0

37 17 24,906,114,455,136 8

38 3 38 0

39 14 233,046 3

n l[n] d[n] i[n]

-----------------------------

113 16 27 9

173 32 82 17

193 73 271 47

2183 72 5,929 32

11229 101 8,201 54

15065 66 11,723 25

15845 139 23,889 43

30817 93 45,391 39

</pre>