login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A371191
Nonsquare terms of A349062.
1
1372, 465125, 4879688, 6272006419, 3533294646441, 405211279147678088
OFFSET
1,1
COMMENTS
There are no more terms below 10^18.
There are 276407671 terms of A349062 below 10^18, of them only 6 are nonsquare numbers.
EXAMPLE
1372 = 2^2 * 7^3 is a term since it is a term of A349062 (the gap to the next powerful number, 1444, is 72, which is a record) and it is not a square.
MATHEMATICA
seq[max_] := Module[{pows = Union[Flatten[Table[i^2*j^3, {j, 1, Surd[max, 3]}, {i, 1, Sqrt[max/j^3]}]]], gapmax = 0, gap, s = {}}, Do[gap = pows[[k+1]] - pows[[k]]; If[gap > gapmax, gapmax = gap; If[!IntegerQ[Sqrt[pows[[k]]]], AppendTo[s, pows[[k]]]]], {k, 1, Length[pows] - 1}]; s]; seq[10^10]
PROG
(PARI) lista(mx) = {my(s = List(), gap, gapmax = 0); for(j = 1, sqrtnint(mx, 3), for(i = 1, sqrtint(mx\j^3), listput(s, i^2 * j^3))); s = Set(s); for(k = 1, #s - 1, gap = s[k+1] - s[k]; if(gap > gapmax, gapmax = gap; if(!issquare(s[k]), print1(s[k], ", ")))); }
CROSSREFS
Intersection of A000037 and A349062.
Cf. A227297.
Sequence in context: A031535 A031715 A135819 * A181969 A139414 A155925
KEYWORD
nonn,hard,more
AUTHOR
Amiram Eldar, Mar 14 2024
STATUS
approved