A379649 a(0) = 1, a(n) = position of prime(n) in A379727, or a(n) = -1 if prime(n) is not in A379727.

1, -1, 2, 4, 3, 5, 9, 15, 8, 6, 13, 44, 10, 17, 40, 7, 38, 14, 134, 20, 16, 35, 37, 18, 46, 24, 12, 153, 39, 34, 27, 45, 71, 104, 43, 32, 11, 95, 70, 19, 82, 47, 22, 165, 98, 36, 179, 73, 178, 28, 159, 130, 554, 365, 67, 220, 72, 124, 21, 154, 58, 127, 83, 79

Offset

0,3

Comments

a(0) = 1 by definition to signify A379727(1) = 1.

Links

Michael De Vlieger, Table of n, a(n) for n = 0..10000

Mathematica

c[_] := True; j = 1; c[1] = False; r = 0; a[0] = 1;
Do[m = 2*j + 1;
  While[
    Set[k, SelectFirst[FactorInteger[m][[All, 1]], c]]; ! IntegerQ[k],
    m = 2*m + 1];
  c[k] = False; j = k; a[PrimePi[k]] = n;
  If[k > r, r = k], {n, 2, 600}];
{1, -1}~Join~TakeWhile[Array[a, PrimePi[r] + 1, 2], IntegerQ]

Crossrefs

Cf. A379727.

Sequence in context: A082329 A072799 A214928 * A240277 A091449 A100834

Adjacent sequences: A379646 A379647 A379648 * A379650 A379651 A379652

Keyword

sign,easy

Author

Michael De Vlieger, Jan 01 2025

Status

approved