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 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