|
|
A267099
|
|
Fully multiplicative involution swapping the positions of 4k+1 and 4k+3 primes: a(1) = 1; a(prime(k)) = A267101(k), a(x*y) = a(x)*a(y) for x, y > 1.
|
|
32
|
|
|
1, 2, 5, 4, 3, 10, 13, 8, 25, 6, 17, 20, 7, 26, 15, 16, 11, 50, 29, 12, 65, 34, 37, 40, 9, 14, 125, 52, 19, 30, 41, 32, 85, 22, 39, 100, 23, 58, 35, 24, 31, 130, 53, 68, 75, 74, 61, 80, 169, 18, 55, 28, 43, 250, 51, 104, 145, 38, 73, 60, 47, 82, 325, 64, 21, 170, 89, 44, 185, 78, 97, 200, 59, 46, 45, 116, 221, 70, 101
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Lexicographically earliest self-inverse permutation of natural numbers where each prime of the form 4k+1 is replaced by a prime of the form 4k+3 and vice versa, with the composite numbers determined by multiplicativity.
Sequences A072202 and A078613 are closed with respect to this permutation.
|
|
LINKS
|
|
|
FORMULA
|
Other identities. For all n >= 1:
a(2*n) = 2*a(n).
a(3*n) = 5*a(n).
a(5*n) = 3*a(n).
a(7*n) = 13*a(n).
a(11*n) = 17*a(n).
|
|
PROG
|
(PARI)
up_to = 2^16;
A267097list(up_to) = { my(v=vector(up_to), i=0, c=0); forprime(p=2, prime(up_to), if(1==(p%4), c++); i++; v[i] = c); (v); };
v267097 = A267097list(up_to);
list_primes_of_the_form(up_to, m, k) = { my(v=vector(up_to), i=0); forprime(p=2, , if(k==(p%m), i++; v[i] = p; if(i==up_to, return(v)))); };
v002144 = list_primes_of_the_form(2*up_to, 4, 1);
v002145 = list_primes_of_the_form(2*up_to, 4, 3);
(Scheme, with memoization-macro definec)
|
|
CROSSREFS
|
Cf. A000035, A000040, A000720, A010051, A020639, A032742, A267100, A267101, A354102 (Möbius transform), A354103 (inverse Möbius transform), A354192 (fixed points).
|
|
KEYWORD
|
nonn,mult
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|