login
A332807
a(n) = A000720(A108546(n)).
3
1, 2, 3, 4, 6, 5, 7, 8, 10, 9, 12, 11, 13, 14, 16, 15, 18, 17, 21, 19, 24, 20, 25, 22, 26, 23, 29, 27, 30, 28, 33, 31, 35, 32, 37, 34, 40, 36, 42, 38, 44, 39, 45, 41, 50, 43, 51, 46, 53, 47, 55, 48, 57, 49, 59, 52, 60, 54, 62, 56, 65, 58, 66, 61, 68, 63, 70, 64, 71, 67, 74, 69, 77, 72, 78, 73, 79, 75, 80, 76, 82, 81, 84, 83, 87, 85, 88, 86
OFFSET
1,2
FORMULA
a(n) = A000720(A108546(n)).
PROG
(PARI)
up_to = 50507;
A332807list(up_to) = { my(v=vector(up_to), p, q); v[1] = 2; v[2] = 3; v[3] = 5; for(n=4, up_to, p = v[n-2]; q = nextprime(1+p); while(q%4 != p%4, q=nextprime(1+q)); v[n] = q); apply(primepi, v); };
v332807 = A332807list(up_to);
A332807(n) = v332807[n];
CROSSREFS
Cf. A332805 (inverse permutation).
Fixed points are given by A000720(A108547(n)), n>=1.
Cf. also A267100.
Sequence in context: A361940 A361926 A332805 * A373038 A267106 A275119
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 27 2020
STATUS
approved