login
A277911
Self-inverse permutation of natural numbers induced when A118306 is restricted to A007310.
6
1, 3, 2, 5, 4, 7, 6, 10, 17, 8, 13, 26, 11, 15, 14, 18, 9, 16, 31, 21, 20, 40, 24, 23, 27, 12, 25, 30, 45, 28, 19, 54, 34, 33, 36, 35, 38, 37, 68, 22, 57, 115, 44, 43, 29, 47, 46, 74, 73, 51, 50, 87, 55, 32, 53, 58, 41, 56, 180, 61, 60, 96, 83, 65, 64, 67, 66, 39, 101, 100, 75, 110, 49, 48, 71, 77, 76, 80, 124, 78, 84, 283, 63, 81, 126, 88, 52
OFFSET
1,2
FORMULA
a(n) = A126760(A118306(A007310(n))).
PROG
(PARI)
A118306(n) = { if(1==n, 1, my(f = factor(n)); my(d = (-1)^primepi(f[1, 1])); for(i=1, #f~, f[i, 1] = prime(primepi(f[i, 1])-d)); factorback(f)); };
A007310(n) = n\2*6-(-1)^n; \\ From M. F. Hasler, Oct 31 2014
A126760(n) = {n&&n\=3^valuation(n, 3)<<valuation(n, 2); n%3+n\6*2}; \\ From M. F. Hasler, Jan 19 2016
for(n=1, 10001, write("b277911.txt", n, " ", A277911(n)));
(Scheme) (define (A277911 n) (A126760 (A118306 (A007310 n))))
CROSSREFS
Cf. A277908, A277909, A277910 (for "almost fixed points"), also A277907.
Sequence in context: A282664 A085230 A093715 * A085229 A093714 A353904
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 05 2016
STATUS
approved