cp's OEIS Frontend

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.

A147802 Least m coprime to 3 minimizing A007947(m*(3^n-m)).

This page as a plain text file.
%I A147802 #19 Aug 11 2020 09:47:07
%S A147802 1,1,2,1,1,25,139,289,31,1,16096,49,7424,588665,619115,83521,8000000,
%T A147802 1515625,505620842,17643776,244140625,5443635008
%N A147802 Least m coprime to 3 minimizing A007947(m*(3^n-m)).
%C A147802 Related to the abc conjecture: Since m is coprime to 3, it is also coprime to 3^n and thus to 3^n-m. Thus A007947(m*(3^n-m)*3^n) = 3*A007947(m(3^n-m)).
%t A147802 rad[n_] := Times @@ FactorInteger[n][[All, 1]];
%t A147802 a[n_] := MinimalBy[Select[Range[3^n - 2], CoprimeQ[#, 3] &], rad[# (3^n - #)] &][[1]];
%t A147802 Reap[Do[Print[n, " ", an = a[n]]; Sow[an], {n, 1, 16}]][[2, 1]] (* _Jean-François Alcover_, Mar 27 2020 *)
%o A147802 (PARI) A147802(n, p=3) = {local(b, m=n=p^n); for(a=1, (n-1)\2, a%p || next; A007947(n-a)*A007947(a)<m || next; m=A007947((n-a)*b=a)); b; }
%Y A147802 Cf. A007947, A143700 (analog for 2^n), A147300 (general case).
%K A147802 nonn,hard,more
%O A147802 1,3
%A A147802 _M. F. Hasler_, Nov 13 2008
%E A147802 a(17) from _Jean-François Alcover_, Mar 28 2020
%E A147802 a(18)-a(22) from _Giovanni Resta_, Mar 29 2020