login
A302850
Permutation of natural numbers: a(1) = 1; a(2) = 2; for n > 2, a(n) = 2+A064989(A003961(n)-4).
3
1, 2, 3, 5, 4, 9, 7, 21, 12, 15, 6, 39, 13, 25, 31, 37, 8, 69, 19, 55, 28, 17, 11, 129, 14, 45, 51, 53, 10, 99, 16, 235, 61, 49, 73, 145, 33, 35, 18, 95, 24, 97, 43, 111, 70, 81, 27, 399, 46, 79, 57, 141, 23, 237, 48, 285, 64, 85, 36, 309, 22, 105, 271, 209, 59, 183, 63, 165, 88, 225, 40, 415, 20, 67, 241, 117, 139, 243, 75, 559, 154, 29, 41, 489, 84, 133
OFFSET
1,2
FORMULA
a(1) = 1; a(2) = 2; for n > 2, a(n) = 2+A064989(A003961(n)-4).
PROG
(PARI)
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
A064989(n) = {my(f); f = factor(n); if((n>1 && f[1, 1]==2), f[1, 2] = 0); for (i=1, #f~, f[i, 1] = precprime(f[i, 1]-1)); factorback(f)};
A302850(n) = if(n<=2, n, 2+A064989(A003961(n)-4));
CROSSREFS
Cf. A302849 (inverse).
Sequence in context: A259431 A085875 A330615 * A338253 A060030 A081025
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 26 2018
STATUS
approved