%I #6 Jan 01 2025 01:48:38

%S 1,7,2,8,3,13,11,6,4,17,9,14,26,49,5,47,18,154,29,12,44,46,27,120,33,

%T 16,102,48,43,36,10,87,21,52,41,15,71,86,28,58,121,31,114,172,45,190,

%U 89,189,37,108,150,19,230,83,247,88,146,30,103,74,149,59,55,119

%N a(n) = position of prime(n) in A379652, or a(n) = -1 if prime(n) is not in A379652.

%H Michael De Vlieger, <a href="/A379648/b379648.txt">Table of n, a(n) for n = 1..10000</a>

%t c[_] := True; j = 2; c[2] = False; r = 0; a[1] = PrimePi[j];

%t Do[m = 2*j + 1;

%t While[

%t Set[k, SelectFirst[FactorInteger[m][[All, 1]], c]]; ! IntegerQ[k],

%t m = 2*m + 1];

%t c[k] = False; j = k; a[PrimePi[k]] = n;

%t If[k > r, r = k], {n, 2, 256}];

%t TakeWhile[Array[a, PrimePi[r]], IntegerQ]

%Y Cf. A379652.

%K nonn,easy

%O 1,2

%A _Michael De Vlieger_, Jan 01 2025