A387072 Index of n in A386482, or -1 if n does not appear in A386482.

1, 2, 5, 3, 17, 4, 11, 9, 6, 8, 24, 7, 44, 10, 16, 14, 52, 13, 31, 15, 12, 23, 73, 22, 18, 21, 26, 20, 126, 19, 57, 29, 25, 28, 41, 27, 115, 30, 43, 40, 809, 39, 821, 38, 42, 37, 76, 36, 49, 35, 51, 34, 652, 33, 70, 48, 32, 47, 144, 46, 189, 56, 50, 55, 45, 54

Offset

1,2

Links

Michael De Vlieger, Table of n, a(n) for n = 1..20000

Mathematica

Block[{a, j, k, p, m, s, nn},
  nn = 2^10; a[_] := 0; m[_] := 1; a[1] = 1; j = a[2] = 2;
  Monitor[Do[
    If[PrimePowerQ[j],
      Set[{p, k, m}, {#1, #1^(#2 - 1), #1^(#2 - 1)}] & @@
        FactorInteger[j][[1]]; While[And[a[k*p] != 0, k != 0], k--];
        If[k == 0, k = m; While[a[k*p] != 0, k++] ]; k *= p,
      k = j - 1; While[And[Or[a[k] != 0, CoprimeQ[j, k]], k != 1], k--];
      If[k == 1, k += j; While[Or[a[k] != 0, CoprimeQ[j, k] ], k++] ] ];
    Set[{a[k], j}, {n, k}], {n, 3, nn}], n];
  k = 1; Reap[While[a[k] > 0, Sow[a[k]]; k++] ][[-1, 1]] ]

Crossrefs

Cf. A386482.

Sequence in context: A063703 A109619 A248576 * A087228 A285743 A077216

Adjacent sequences: A387069 A387070 A387071 * A387073 A387074 A387075

Keyword

nonn

Author

Michael De Vlieger, Aug 15 2025

Status

approved