login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A268675
Self-inverse permutation of natural numbers: a(1) = 1; a(n) = A000079(A193231(A007814(n))) * A250469(a(A268674(n))).
2
1, 2, 3, 8, 5, 6, 7, 4, 21, 10, 11, 24, 13, 14, 15, 32, 17, 42, 19, 40, 9, 22, 23, 12, 55, 26, 27, 56, 29, 30, 31, 16, 69, 34, 35, 168, 37, 38, 39, 20, 41, 18, 43, 88, 93, 46, 47, 96, 91, 110, 123, 104, 53, 54, 25, 28, 117, 58, 59, 120, 61, 62, 63, 64, 65, 138, 67, 136, 33, 70, 71, 84, 73, 74, 75, 152, 77, 78, 79, 160
OFFSET
1,2
FORMULA
a(1) = 1, and for n > 1, a(n) = A000079(A193231(A007814(n))) * A250469(a(A268674(n))).
Other identities. For all n >= 1:
A000035(a(n)) = A000035(n). [This permutation preserves the parity of n.]
A020639(a(n)) = A020639(n). [More generally, it preserves the smallest prime dividing n.]
A055396(a(n)) = A055396(n).
PROG
(Scheme, with memoization-macro definec)
(definec (A268675 n) (if (= 1 n) 1 (* (A000079 (A193231 (A007814 n))) (A250469 (A268675 (A268674 n))))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 11 2016
STATUS
approved