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.

A375622 Numbers k such that k gives the maximum quality for abc triples of the form (1, k^n-1, k^n) where n is the sequence index.

This page as a plain text file.
%I A375622 #11 Sep 05 2024 03:27:53
%S A375622 4375,49,361,7,9,19,129,7,28,3,243,19,625,26,9,3
%N A375622 Numbers k such that k gives the maximum quality for abc triples of the form (1, k^n-1, k^n) where n is the sequence index.
%C A375622 An abc triple is defined as (a, b, c) with a + b = c, gcd(a, b) = 1 and radical of a*b*c, rad(a*b*c) < c. The quality of an abc triple is q = log(c)/log(rad(a*b*c)). For each n, a sample of the first 100 abc triples of the form (1, k^n-1, k^n) is compared to find the value of k that gives the abc triple maximum quality. The sample size s = 100 of abc triples appears adequate to identify the maximum quality because the quality term tends rapidly towards the lim sup(q) = 1 as s -> oo.
%H A375622 Elise Alvarez-Salazar, Alexander J. Barrios, Calvin Henaku, and Summer Soller, <a href="https://arxiv.org/abs/2301.01376">On abc triples of the form (1,c-1,c)</a>, arXiv:2301.01376 [math.NT], 2023.
%e A375622 a(2) = 49 because from the sample of 100 abc triples of the form (1, k^2-1, k^2) (see A375019) where k takes values 3, 7,..., 49,..., 3362, 3375, when k = 49 = A375019(12), we get maximum quality q = 1.45567... with triple (1, 2400, 2401).
%t A375622 Rad[n_] := Module[{lst=FactorInteger[n]}, Times@@(First/@lst)]; Table[(lst={}; k=2; While[Length@lst<100, If[Rad[(k^n-1)*k]<k^n, AppendTo[lst, k]]; k++]; lst1=Table[k=lst[[j]]; r=Rad[k(k^n-1)]; {n*Log[k]/Log[r], k}, {j, 1, Length@lst}]; (Flatten@Select[lst1, #[[1]]==Max@(First /@ lst1) &])[[2]]), {n, 1, 16}]
%Y A375622 Cf. A007947, A216323, A375019.
%K A375622 nonn,more
%O A375622 1,1
%A A375622 _Frank M Jackson_, Aug 21 2024