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.

Showing 1-2 of 2 results.

A334882 Numbers k such that k and k+2 are both primitive practical numbers (A267124).

Original entry on oeis.org

28, 304, 306, 340, 460, 462, 858, 868, 1482, 1768, 1974, 2440, 2728, 2838, 2860, 3318, 3738, 4134, 4264, 4288, 4420, 4422, 5236, 5694, 6100, 6102, 7590, 8814, 9040, 9042, 10218, 11128, 11620, 11778, 12558, 12978, 13110, 14320, 14382, 14670, 15568, 16048, 16110
Offset: 1

Views

Author

Amiram Eldar, May 14 2020

Keywords

Examples

			28 is a term since 28 and 28 + 2 = 30 are both primitive practical numbers.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := (p^(e + 1) - 1)/(p - 1); pracQ[fct_] := (ind = Position[fct[[;; , 1]]/(1 + FoldList[Times, 1, f @@@ Most@fct]), ?(# > 1 &)]) == {}; pracTestQ[fct, k_] := Module[{f = fct}, f[[k, 2]] -= 1; pracQ[f]]; primPracQ[n_] := Module[{fct = FactorInteger[n]}, pracQ[fct] && AllTrue[Range@Length[fct], fct[[#, 2]] == 1 || ! pracTestQ[fct, #] &]]; Select[Range[2, 16200, 2], primPracQ[#] && primPracQ[# + 2] &]

A364975 Admirable numbers (A111592) with a record gap to the next admirable number.

Original entry on oeis.org

12, 30, 42, 88, 120, 140, 186, 534, 678, 6774, 7962, 77118, 94108, 152826, 478194, 662154, 935564, 1128174, 2028198, 6934398, 7750146, 8330924, 9984738, 10030804, 22956114, 62062566, 151040622, 284791602, 732988732, 804394974, 1151476732, 9040886574, 31302713634
Offset: 1

Views

Author

Amiram Eldar, Aug 15 2023

Keywords

Comments

The corresponding record gaps are 8, 10, 12, 14, 18, 34, 36, 48, 84, 132, 204, 216, 254, 312, 348, 360, 392, 468, 516, 528, 552, 598, 624, 638, 828, 852, 936, 1056, 1082, 1128, 1454, 1692, 1752, ... .

Examples

			The first 5 admirable numbers are 12, 20, 24, 30 and 40. The differences between these terms are 8, 4, 6 and 10. The record gaps, 8 and 10, occur after the terms 12 and 30, which are the first two terms of this sequence.
		

Crossrefs

Similar sequences: A306953, A330870, A334418, A334419, A334883, A363296.

Programs

  • Mathematica
    admQ[n_] := (ab = DivisorSigma[1, n] - 2 n) > 0 && EvenQ[ab] && ab/2 < n && Divisible[n, ab/2];
    seq[kmax_] := Module[{s = {}, m = 12, dm = 0}, Do[If[admQ[k], d = k - m; If[d > dm, dm = d; AppendTo[s, m]]; m = k], {k, m + 1, kmax}]; s]; seq[10^6]
  • PARI
    isadm(n) = {my(ab=sigma(n)-2*n); ab>0 && ab%2 == 0 && ab/2 < n && n%(ab/2) == 0; }
    lista(kmax) = {my(m = 12, dm = 0); for(k = m+1, kmax, if(isadm(k), d = k - m; if(d > dm, dm = d; print1(m, ", ")); m = k));}
Showing 1-2 of 2 results.