This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A370071 #64 Mar 28 2024 00:01:03 %S A370071 0,1,8,144,610,46368,1346269,14930352,63245986,165580141,267914296, %T A370071 2971215073,4807526976,7778742049,12586269025,139583862445, %U A370071 365435296162,591286729879,956722026041,1548008755920,2504730781961,4052739537881,6557470319842,17167680177565,27777890035288,190392490709135 %N A370071 Fibonacci numbers such that some permutation of their digits is a perfect power. %C A370071 0, 1, 8 and 144 are the only Fibonacci numbers that are themselves perfect powers (see A227875). %C A370071 A term might have multiple permutations which are perfect powers. %C A370071 Leading 0 digits are allowed in the perfect power, so that 610 is a term since 016 = 2^4. %C A370071 From _David A. Corneth_, Feb 17 2024: (Start) %C A370071 Four ideas to find terms: %C A370071 1. For each Fibonacci number, check each anagram to determine whether it is a perfect power. A Fibonacci number is a term if and only if such an anagram exists. %C A370071 2. Let q be the number of digits of some Fibonacci number F. Then check all perfect powers m < 10^q if the frequency of positive digits matches the corresponding positive digits of F and F has at least as many 0's as m. E.g., 610 is a term as the perfect power 16 has the same number of 1's, and the same number of 6's, and 610 has at least as many 0's as 16. %C A370071 3. Match the last digits of perfect powers and whittle down the number of candidates. So for 46368, a perfect square must end in 4 or 6 for the last digit, 36, 64 or 84 for the last two digits and so on. %C A370071 4. If some Fibonacci number F has q digits then we could list the possible strings the first floor((q + 1)/2) such numbers could form and see what squares start with those digits. (End) %C A370071 a(46) = 1500520536206896083277. Some other terms: 3928413764606871165730, 6356306993006846248183, 10284720757613717413913, 16641027750620563662096, 26925748508234281076009, 43566776258854844738105, 3311648143516982017180081, 5358359254990966640871840, 8670007398507948658051921, 14028366653498915298923761, 155576970220531065681649693, 659034621587630041982498215, 30960598847965113057878492344. - _Chai Wah Wu_, Mar 27 2024 %H A370071 Michael S. Branicky, <a href="/A370071/b370071.txt">Table of n, a(n) for n = 1..45</a> %e A370071 46368 is a term because it is a Fibonacci number whose digits can be permuted to 36864 = 192^2 (and also to 86436 = 294^2). %o A370071 (PARI) isok(f) = my(d=digits(f), n=#d); for (i=1, n!, my(p=numtoperm(n, i), dd=vector(#d, i, d[p[i]])); if (ispower(fromdigits(dd)), return(1));); %o A370071 lista(nn) = my(list = List([0, 1])); for (n=3, nn, my(f=fibonacci(n)); if (isok(f), listput(list, f));); Vec(list); \\ _Michel Marcus_, Feb 17 2024 %o A370071 (Python) %o A370071 from itertools import count, islice %o A370071 from sympy import integer_log %o A370071 def A370071_gen(): # generator of terms %o A370071 a, b = 1, 2 %o A370071 yield from (0,1) %o A370071 while True: %o A370071 s = sorted(str(b)) %o A370071 l = len(s) %o A370071 m = int(''.join(s[::-1])) %o A370071 u = int(''.join(s)) %o A370071 for i in count(2): %o A370071 if i**2 > m: %o A370071 break %o A370071 for j in count(max(2,integer_log(u,i)[0])): %o A370071 if (k:=i**j) > m: %o A370071 break %o A370071 t = sorted(str(k)) %o A370071 if ['0']*(l-len(t))+t == s: %o A370071 yield b %o A370071 break %o A370071 else: %o A370071 continue %o A370071 if k<=m: %o A370071 break %o A370071 a, b = b, a+b %o A370071 A370071_list = list(islice(A370071_gen(),20)) # _Chai Wah Wu_, Mar 27 2024 %Y A370071 Cf. A000045, A227875, A001597, A118715. %K A370071 nonn,base %O A370071 1,3 %A A370071 _Gonzalo MartÃnez_, Feb 08 2024 %E A370071 a(19) inserted and a(21)-a(23) by _Michael S. Branicky_, Feb 17 2024 %E A370071 More terms from _David A. Corneth_, Feb 17 2024