cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A235201 Self-inverse and multiplicative permutation of integers: a(0)=0, a(1)=1, a(2)=2, a(3)=4 and a(4)=3, a(p_i) = p_{a(i)} for primes with index i > 2, and for composites > 4, a(u * v) = a(u) * a(v) for u, v > 0.

Original entry on oeis.org

0, 1, 2, 4, 3, 7, 8, 5, 6, 16, 14, 17, 12, 19, 10, 28, 9, 11, 32, 13, 21, 20, 34, 53, 24, 49, 38, 64, 15, 43, 56, 59, 18, 68, 22, 35, 48, 37, 26, 76, 42, 67, 40, 29, 51, 112, 106, 107, 36, 25, 98, 44, 57, 23, 128, 119, 30, 52, 86, 31, 84, 131, 118, 80, 27, 133, 136, 41, 33
Offset: 0

Views

Author

Antti Karttunen, Jan 11 2014

Keywords

Comments

The permutation satisfies A000040(a(n)) = a(A000040(n)) for all positive n except n=2, and is self-inverse. It swaps 3 & 4, maps any prime p_i with index i > 2 to p_{a(i)}, and lets the multiplicativity take care of the rest.
This can be viewed also as a "signature-permutation" for a bijection on non-oriented rooted trees, mapped through the Matula-Goebel numbers (cf. A061773). This bijection will swap the subtrees encoded by numbers 3 and 4, wherever they occur as the terminal configurations anywhere in the tree:
....................
.o..................
.|..................
.o.............o...o
.|..............\./.
.x.....<--->.....x..
.3...............4..
That is, the last two edges of any branch which ends with at least in two edges long unbranched stem, will be changed to a V-branch (two single edges in parallel). Vice versa, any terminal configuration in the tree that consists of more than one single edges next to each other (in "parallel") will be transformed so that maximal even number (2k) of those single edges will be combined to k unbranching stems of two edges, and an extra odd edge, if present, will stay as it is.
This permutation commutes with A235199, i.e. a(A235199(n)) = A235199(a(n)) for all n. This can be easily seen, when comparing the above bijection to the one described in A235199. Composition A235199 o A235201 works as a "difference" of these two bijections, swapping the above subconfigurations only when they do not occur alone at the tips of singular edges. (Which cases are encoded by Matula-Goebel numbers 5 and 7, the third and fourth prime respectively).
Permutation fixes n! for n=0, 1, 2, 4, 7.
Note that a(5!) = a(120) = 168 = 120+(2*4!) and a(8!) = a(40320) = 30240 = 40320-(2*7!).

Crossrefs

Composition with A235487 gives A235485/A235486, composition with A235489 gives A235493/A235494.
List below gives similarly constructed permutations, which all force a swap of two small numbers, with (the rest of) primes permuted with the sequence itself and the new positions of composite numbers defined by the multiplicative property:
A234840 (swaps 2 & 3, conjugates A008578 back to itself).
A235200 (swaps 3 & 5, conjugates A065091 back to itself).
A235199 (swaps 5 & 7, conjugates A000040 back to itself).
A235487 (swaps 7 & 8, conjugates A000040 back to itself).
A235489 (swaps 8 & 9, conjugates A000040 back to itself).

Formula

Multiplicative with a(3^k) = 2^(2k), a(2^(2k)) = 3^k, a(2^(2k+1)) = 2*3^k, a(p_i) = p_{a(i)} for primes with index i > 2, and for composites > 4, a(u * v) = a(u) * a(v) for u, v > 0.

A235199 Self-inverse and multiplicative permutation of integers: For n < 4, a(n)=n, a(5)=7 and a(7)=5, a(p_i) = p_{a(i)} for primes with index i > 4, and a(u * v) = a(u) * a(v) for u, v > 0.

Original entry on oeis.org

0, 1, 2, 3, 4, 7, 6, 5, 8, 9, 14, 17, 12, 13, 10, 21, 16, 11, 18, 19, 28, 15, 34, 23, 24, 49, 26, 27, 20, 43, 42, 59, 32, 51, 22, 35, 36, 37, 38, 39, 56, 41, 30, 29, 68, 63, 46, 73, 48, 25, 98, 33, 52, 53, 54, 119, 40, 57, 86, 31, 84, 61, 118, 45, 64, 91, 102
Offset: 0

Views

Author

Antti Karttunen, Jan 04 2014

Keywords

Comments

The permutation satisfies A000040(a(n)) = a(A000040(n)) for all positive n except n=3 or 4, and is self-inverse. It swaps 5 & 7, maps all larger primes p_i (with index i > 4) to p_{a(i)}, and lets the multiplicativity take care of the rest.
It can be viewed also as a "signature-permutation" for a bijection of non-oriented rooted trees, mapped through Matula-Goebel numbers (cf. A061773). The bijection will swap the subtrees encoded by primes 5 and 7, wherever they occur as the terminal branches of the tree:
....................
.o..................
.|..................
.o.............o...o
.|..............\./.
.o.....<--->.....o..
.|...............|..
.x...............x..
.5...............7..
That is, any branch which ends at least in three edges long unbranched stem, will be changed so that its last two edges will become V-branch. Vice versa, any branch of the tree that ends with three edges in Y-formation, will be transformed so that those three edges will be straightened to an unbranching stem of three edges.
This permutation commutes with A235201, i.e. a(A235201(n)) = A235201(a(n)) for all n.
Permutation fixes n! for n=0, 1, 2, 3, 4, 7, 8 and 9.
Note also that a(5!) = a(120) = 168 = 120+(2*4!) and a(10!) = 5080320 = 3628800+(4*9!).

Crossrefs

Composition with A234840 gives A234743 & A234744.
List below gives similarly constructed permutations, which all force a swap of two small numbers, with (the rest of) primes permuted with the sequence itself and the new positions of composite numbers defined by the multiplicative property:
A234840 (swaps 2 & 3, conjugates A008578 back to itself).
A235200 (swaps 3 & 5, conjugates A065091 back to itself).
A235201 (swaps 3 & 4, conjugates A000040 back to itself).
A235487 (swaps 7 & 8, conjugates A000040 back to itself).
A235489 (swaps 8 & 9, conjugates A000040 back to itself).

Formula

For n < 4, a(n)=n, a(5)=7 and a(7)=5, a(p_i) = p_{a(i)} for primes with index i > 4, and a(u * v) = a(u) * a(v) for u, v > 0.
A000035(a(n)) = A000035(n) = (n mod 2) for all n. [Even terms occur only on even indices and odd terms only on odd indices, respectively]

A072029 Swap twin prime pairs of form (4*k+3,4*(k+1)+1) in prime factorization of n.

Original entry on oeis.org

1, 2, 5, 4, 3, 10, 7, 8, 25, 6, 13, 20, 11, 14, 15, 16, 17, 50, 19, 12, 35, 26, 23, 40, 9, 22, 125, 28, 29, 30, 31, 32, 65, 34, 21, 100, 37, 38, 55, 24, 41, 70, 43, 52, 75, 46, 47, 80, 49, 18, 85, 44, 53, 250, 39, 56, 95, 58, 61, 60, 59, 62, 175, 64, 33
Offset: 1

Views

Author

Reinhard Zumkeller, Jun 07 2002

Keywords

Examples

			a(42) = a(2*3*7) = a(2)*a(3)*a(7) = a(2)*a(4*0+3)*a(7) = 2*(4*1+1)*7 = 2*5*7 = 70.
		

Crossrefs

Programs

  • Mathematica
    a[n_] := Product[{p, e} = pe; Which[
         Mod[p, 4] == 3 && PrimeQ[p + 2], p + 2,
         Mod[p, 4] == 1 && PrimeQ[p - 2], p - 2,
         True, p]^e, {pe, FactorInteger[n]}];
    Array[a, 100] (* Jean-François Alcover, Nov 21 2021 *)
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, p = f[i,1]; if(p == 2, p, if(p%4 == 3 && isprime(p+2), p+2, if(p%4 == 1 && isprime(p-2), p-2, p)))^f[i,2]);} \\ Amiram Eldar, Feb 26 2024

Formula

Multiplicative with a(p) = (if p mod 4 = 3 and p+2 is prime then p+2 else (if p mod 4 = 1 and p-2 is prime then p-2 else p)), p prime.
a(a(n))=n, a self-inverse permutation of natural numbers.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{(p < q) swapped pair} ((p^2-p)*(q^2-q)/((p^2-q)*(q^2-p))) = 1.37140598833326962... . - Amiram Eldar, Feb 26 2024

A072027 Swap (2,3) and all twin prime pairs >(3,5) in prime factorization of n.

Original entry on oeis.org

1, 3, 2, 9, 7, 6, 5, 27, 4, 21, 13, 18, 11, 15, 14, 81, 19, 12, 17, 63, 10, 39, 23, 54, 49, 33, 8, 45, 31, 42, 29, 243, 26, 57, 35, 36, 37, 51, 22, 189, 43, 30, 41, 117, 28, 69, 47, 162, 25, 147, 38, 99, 53, 24, 91, 135, 34, 93, 61, 126, 59, 87, 20, 729, 77
Offset: 1

Views

Author

Reinhard Zumkeller, Jun 07 2002

Keywords

Examples

			a(143) = a(11*13) = a(11)*a(13) = 13*11 = 143.
a(77) = a(7*11) = a(7)*a(11) = 5*13 = 65.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[p < 5, 5 - p, If[PrimeQ[p + 2], p + 2, If[PrimeQ[p - 2], p - 2, p]]]^e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Feb 26 2024 *)
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, p = f[i,1]; if(p < 5, 5-p, if(isprime(p+2), p+2, if(isprime(p-2), p-2, p)))^f[i,2]);} \\ Amiram Eldar, Feb 26 2024

Formula

Multiplicative with a(p) = (if p<=3 then 5-p else (if p+2 is prime then p+2 else (if p-2 is prime then p-2 else p))), p prime.
a(a(n)) = n, self-inverse permutation of natural numbers.
a(n) = n for single primes (A007510) and products of twin prime pairs (A037074).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{(p < q) swapped pair} ((p^2-p)*(q^2-q)/((p^2-q)*(q^2-p))) = 1.832194438922717... . - Amiram Eldar, Feb 26 2024

A072028 Swap twin prime pairs of form (4*k+1,4*k+3) in prime factorization of n.

Original entry on oeis.org

1, 2, 3, 4, 7, 6, 5, 8, 9, 14, 11, 12, 13, 10, 21, 16, 19, 18, 17, 28, 15, 22, 23, 24, 49, 26, 27, 20, 31, 42, 29, 32, 33, 38, 35, 36, 37, 34, 39, 56, 43, 30, 41, 44, 63, 46, 47, 48, 25, 98, 57, 52, 53, 54, 77, 40, 51, 62, 59, 84, 61, 58, 45, 64, 91, 66
Offset: 1

Views

Author

Reinhard Zumkeller, Jun 07 2002

Keywords

Examples

			a(65) = a(5*13) = a(5)*a(13) = a(4*1+1)*a(13) = (4*1+3)*13 = 7*13 = 91.
		

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := If[p < 5, p, If[(m = Mod[p, 4]) == 1 && PrimeQ[p + 2], p + 2, If[m == 3 && PrimeQ[p - 2], p - 2, p]]]^e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Feb 26 2024 *)
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, p = f[i,1]; if(p < 5, p, if(p%4 == 1 && isprime(p+2), p+2, if(p%4 == 3 && isprime(p-2), p-2, p)))^f[i,2]);} \\ Amiram Eldar, Feb 26 2024

Formula

Multiplicative with a(p) = (if p mod 4 = 1 and p+2 is prime then p+2 else (if p mod 4 = 3 and p-2 is prime then p-2 else p)), p prime.
a(a(n)) = n, a self-inverse permutation of natural numbers.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = Product_{(p < q) swapped pair} ((p^2-p)*(q^2-q)/((p^2-q)*(q^2-p))) = 1.0627249749498993391... . - Amiram Eldar, Feb 26 2024
Showing 1-5 of 5 results.