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-9 of 9 results.

A135141 a(1)=1, a(p_n)=2*a(n), a(c_n)=2*a(n)+1, where p_n = n-th prime, c_n = n-th composite number.

Original entry on oeis.org

1, 2, 4, 3, 8, 5, 6, 9, 7, 17, 16, 11, 10, 13, 19, 15, 12, 35, 18, 33, 23, 21, 14, 27, 39, 31, 25, 71, 34, 37, 32, 67, 47, 43, 29, 55, 22, 79, 63, 51, 20, 143, 26, 69, 75, 65, 38, 135, 95, 87, 59, 111, 30, 45, 159, 127, 103, 41, 24, 287, 70, 53, 139, 151, 131, 77, 36, 271, 191
Offset: 1

Views

Author

Katarzyna Matylla, Feb 13 2008

Keywords

Comments

A permutation of the positive integers, related to A078442.
a(p) is even when p is prime and is divisible by 2^(prime order of p).
From Robert G. Wilson v, Feb 16 2008: (Start)
What is the length of the cycle containing 10? Is it infinite? The cycle begins 10, 17, 12, 11, 16, 15, 19, 18, 35, 29, 34, 43, 26, 31, 32, 67, 36, 55, 159, 1055, 441, 563, 100, 447, 7935, 274726911, 1013992070762272391167, ... Implementation in Mmca: NestList[a(AT)# &, 10, 26] Furthermore, it appears that any non-single-digit number has an infinite cycle.
Records: 1, 2, 4, 8, 9, 17, 19, 35, 39, 71, 79, 143, 159, 287, 319, 575, 639, 1151, 1279, 2303, 2559, 4607, 5119, 9215, 10239, 18431, 20479, 36863, 40959, 73727, 81919, 147455, 163839, 294911, 327679, 589823, 655359, ..., . (End)

Examples

			a(20) = 33 = 2*16 + 1 because 20 is 11th composite and a(11)=16. Or, a(20)=33=100001(bin). In other words it is a composite number, its index is a prime number, whose index is a prime....
		

Crossrefs

Cf. A246346, A246347 (record positions and values).
Cf. A227413 (inverse).
Cf. A071574, A245701, A245702, A245703, A245704, A246377, A236854, A237427 for related and similar permutations.

Programs

  • Haskell
    import Data.List (genericIndex)
    a135141 n = genericIndex a135141_list (n-1)
    a135141_list = 1 : map f [2..] where
       f x | iprime == 0 = 2 * (a135141 $ a066246 x) + 1
           | otherwise   = 2 * (a135141 iprime)
           where iprime = a049084 x
    -- Reinhard Zumkeller, Jan 29 2014
    
  • Mathematica
    a[1] = 1; a[n_] := If[PrimeQ@n, 2*a[PrimePi[n]], 2*a[n - 1 - PrimePi@n] + 1]; Array[a, 69] (* Robert G. Wilson v, Feb 16 2008 *)
  • Maxima
    /* Let pc = prime count (which prime it is), cc = composite count: */
    pc[1]:0;
    cc[1]:0;
    pc[2]:1;
    cc[4]:1;
    pc[n]:=if primep(n) then 1+pc[prev_prime(n)] else 0;
    cc[n]:=if primep(n) then 0 else if primep(n-1) then 1+cc[n-2] else 1+cc[n-1];
    a[1]:1;
    a[n]:=if primep(n) then 2*a[pc[n]] else 1+2*a[cc[n]];
    
  • PARI
    A135141(n) = if(1==n, 1, if(isprime(n), 2*A135141(primepi(n)), 1+(2*A135141(n-primepi(n)-1)))); \\ Antti Karttunen, Dec 09 2019
  • Python
    from sympy import isprime, primepi
    def a(n): return 1 if n==1 else 2*a(primepi(n)) if isprime(n) else 2*a(n - 1 - primepi(n)) + 1 # Indranil Ghosh, Jun 11 2017, after Mathematica code
    

Formula

a(n) = 2*A135141((A049084(n))*chip + A066246(n)*(1-chip)) + 1 - chip, where chip = A010051(n). - Reinhard Zumkeller, Jan 29 2014
From Antti Karttunen, Dec 09 2019: (Start)
A007814(a(n)) = A078442(n).
A070939(a(n)) = A246348(n).
A080791(a(n)) = A246370(n).
A054429(a(n)) = A246377(n).
A245702(a(n)) = A245703(n).
a(A245704(n)) = A245701(n). (End)

A245703 Permutation of natural numbers: a(1) = 1, a(p_n) = A014580(a(n)), a(c_n) = A091242(a(n)), where p_n = n-th prime, c_n = n-th composite number and A014580(n) and A091242(n) are binary codes for n-th irreducible and n-th reducible polynomials over GF(2), respectively.

Original entry on oeis.org

1, 2, 3, 4, 7, 5, 11, 6, 8, 12, 25, 9, 13, 17, 10, 14, 47, 18, 19, 34, 15, 20, 31, 24, 16, 21, 62, 26, 55, 27, 137, 45, 22, 28, 42, 33, 37, 23, 29, 79, 59, 35, 87, 71, 36, 166, 41, 58, 30, 38, 54, 44, 61, 49, 32, 39, 99, 76, 319, 46, 91, 108, 89, 48, 200, 53, 97, 75, 40, 50, 203, 70, 67, 57, 78, 64, 43, 51
Offset: 1

Views

Author

Antti Karttunen, Aug 02 2014

Keywords

Comments

All the permutations A091202, A091204, A106442, A106444, A106446, A235041 share the same property that primes (A000040) are mapped bijectively to the binary representations of irreducible GF(2) polynomials (A014580) but while they determine the mapping of composites (A002808) to the corresponding binary codes of reducible polynomials (A091242) by a simple multiplicative rule, this permutation employs index-recursion also in that case.

Crossrefs

Programs

  • PARI
    allocatemem(123456789);
    a014580 = vector(2^18);
    a091242 = vector(2^22);
    isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
    i=0; j=0; n=2; while((n < 2^22), if(isA014580(n), i++; a014580[i] = n, j++; a091242[j] = n); n++)
    A245703(n) = if(1==n, 1, if(isprime(n), a014580[A245703(primepi(n))], a091242[A245703(n-primepi(n)-1)]));
    for(n=1, 10001, write("b245703.txt", n, " ", A245703(n)));
    
  • Scheme
    ;; With memoization-macro definec.
    (definec (A245703 n) (cond ((= 1 n) n) ((= 1 (A010051 n)) (A014580 (A245703 (A000720 n)))) (else (A091242 (A245703 (A065855 n))))))

Formula

a(1) = 1, a(p_n) = A014580(a(n)) and a(c_n) = A091242(a(n)), where p_n is the n-th prime, A000040(n) and c_n is the n-th composite, A002808(n).
a(1) = 1, after which, if A010051(n) is 1 [i.e. n is prime], then a(n) = A014580(a(A000720(n))), otherwise a(n) = A091242(a(A065855(n))).
As a composition of related permutations:
a(n) = A245702(A135141(n)).
a(n) = A091204(A245821(n)).
Other identities. For all n >= 1, the following holds:
a(A007097(n)) = A091230(n). [Maps iterates of primes to the iterates of A014580. Permutation A091204 has the same property]
A091225(a(n)) = A010051(n). [Maps primes to binary representations of irreducible GF(2) polynomials, A014580, and nonprimes to union of {1} and the binary representations of corresponding reducible polynomials, A091242. The permutations A091202, A091204, A106442, A106444, A106446 and A235041 have the same property.]

A245704 Permutation of natural numbers: a(1) = 1, a(A014580(n)) = A000040(a(n)), a(A091242(n)) = A002808(a(n)), where A000040(n) = n-th prime, A002808(n) = n-th composite number, and A014580(n) and A091242(n) are binary codes for n-th irreducible and n-th reducible polynomial over GF(2), respectively.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 5, 9, 12, 15, 7, 10, 13, 16, 21, 25, 14, 18, 19, 22, 26, 33, 38, 24, 11, 28, 30, 34, 39, 49, 23, 55, 36, 20, 42, 45, 37, 50, 56, 69, 47, 35, 77, 52, 32, 60, 17, 64, 54, 70, 78, 94, 66, 51, 29, 105, 74, 48, 41, 84, 53, 27, 88, 76, 95, 106, 73, 125, 91, 72, 44, 140, 97, 100, 68, 58, 115, 75, 40
Offset: 1

Views

Author

Antti Karttunen, Aug 02 2014

Keywords

Comments

All the permutations A091203, A091205, A106443, A106445, A106447, A235042 share the same property that the binary representations of irreducible GF(2) polynomials (A014580) are mapped bijectively to the primes (A000040) but while they determine the mapping of corresponding reducible polynomials (A091242) to the composite numbers (A002808) by a simple multiplicative rule, this permutation employs index-recursion also in that case.

Crossrefs

Programs

Formula

a(1) = 1, after which, if A091225(n) is 1 [i.e. n is in A014580], then a(n) = A000040(a(A091226(n))), otherwise a(n) = A002808(a(A091245(n))).
As a composition of related permutations:
a(n) = A227413(A245701(n)).
a(n) = A245822(A091205(n)).
Other identities. For all n >= 1, the following holds:
a(A091230(n)) = A007097(n). [Maps iterates of A014580 to the iterates of primes. Permutation A091205 has the same property].
A010051(a(n)) = A091225(n). [After a(1)=1, maps binary representations of irreducible GF(2) polynomials (= A014580) to primes and the corresponding representations of reducible polynomials to composites].

A245702 Permutation of natural numbers: a(1) = 1, a(2n) = A014580(a(n)), a(2n+1) = A091242(a(n)), where A014580(n) = binary code for n-th irreducible polynomial over GF(2) and A091242(n) = binary code for n-th reducible polynomial over GF(2).

Original entry on oeis.org

1, 2, 4, 3, 5, 11, 8, 7, 6, 13, 9, 47, 17, 31, 14, 25, 12, 19, 10, 59, 20, 37, 15, 319, 62, 87, 24, 185, 42, 61, 21, 137, 34, 55, 18, 97, 27, 41, 16, 415, 76, 103, 28, 229, 49, 67, 22, 3053, 373, 433, 79, 647, 108, 131, 33, 1627, 222, 283, 54, 425, 78, 109, 29, 1123, 166, 203, 45, 379, 71, 91, 26, 731, 121, 145, 36, 253, 53, 73, 23
Offset: 1

Views

Author

Antti Karttunen, Aug 02 2014

Keywords

Crossrefs

Inverse: A245701.
Similar entanglement permutations: A193231, A227413, A237126, A243288, A245703, A245704.

Programs

  • PARI
    allocatemem(123456789);
    a014580 = vector(2^18);
    a091242 = vector(2^22);
    isA014580(n)=polisirreducible(Pol(binary(n))*Mod(1, 2)); \\ This function from Charles R Greathouse IV
    i=0; j=0; n=2; while((n < 2^22), if(isA014580(n), i++; a014580[i] = n, j++; a091242[j] = n); n++)
    A245702(n) = if(1==n, 1, if(0==(n%2), a014580[A245702(n/2)], a091242[A245702((n-1)/2)]));
    for(n=1, 383, write("b245702.txt", n, " ", A245702(n)));
    
  • Scheme
    ;; With memoizing definec-macro.
    (definec (A245702 n) (cond ((< n 2) n) ((even? n) (A014580 (A245702 (/ n 2)))) (else (A091242 (A245702 (/ (- n 1) 2))))))

Formula

a(1) = 1, a(2n) = A014580(a(n)), a(2n+1) = A091242(a(n)).
As a composition of related permutations:
a(n) = A245703(A227413(n)).
Other identities:
For all n >= 1, 1 - A091225(a(n)) = A000035(n). [Maps even numbers to binary representations of irreducible GF(2) polynomials (= A014580) and odd numbers to the corresponding representations of reducible polynomials].

A246201 Permutation of natural numbers: a(1) = 1, a(A014580(n)) = (2*a(n))+1, a(A091242(n)) = 2*a(n), where A014580(n) = binary code for n-th irreducible polynomial over GF(2), A091242(n) = binary code for n-th reducible polynomial over GF(2).

Original entry on oeis.org

1, 3, 7, 2, 6, 14, 15, 4, 12, 28, 5, 30, 13, 8, 24, 56, 10, 60, 29, 26, 16, 48, 112, 20, 31, 120, 58, 52, 32, 96, 9, 224, 40, 62, 240, 116, 25, 104, 64, 192, 57, 18, 448, 80, 124, 480, 11, 232, 50, 208, 128, 384, 114, 36, 61, 896, 160, 248, 27, 960, 17, 22, 464, 100, 416, 256, 49, 768, 228, 72, 122, 1792, 113, 320, 496, 54, 1920, 34, 44
Offset: 1

Views

Author

Antti Karttunen, Aug 19 2014

Keywords

Comments

Because 2 is the only even term in A014580, it implies that, apart from a(2)=3, odd numbers occur in odd positions only (along with many even numbers that also occur in odd positions).
Note that for any value k in A246156, "Odd reducible polynomials over GF(2)": 5, 9, 15, 17, 21, 23, ..., a(k) will be even, and apart from 2, all other even numbers are mapped to some even number, so all those terms reside in infinite cycles. Furthermore, apart from 5 and 15, all of them reside in separate cycles. The infinite cycle containing 5 and 15 goes as: ..., 47, 11, 5, 6, 14, 8, 4, 2, 3, 7, 15, 24, 20, 26, 120, 7680, ... and it is only because a(2) = 3, that it can temporarily switch back from even terms to odd terms, until after a(15) = 24 it is finally doomed to the eternal evenness.
(Compare also to the comments given at A246161).

Crossrefs

Inverse: A246202.
Similar or related permutations: A245701, A246161, A006068, A054429, A193231, A246163, A246203, A237427.

Formula

a(1) = 1, and for n > 1, if A091225(n) = 1 [i.e. when n is in A014580], a(n) = 1 + (2*a(A091226(n))), otherwise a(n) = 2*a(A091245(n)).
As a composition of related permutations:
a(n) = A054429(A245701(n)).
a(n) = A006068(A246161(n)).
a(n) = A193231(A246163(n)).
a(n) = A246203(A193231(n)).
Other identities:
For all n > 1, A000035(a(n)) = A091225(n). [After 1 maps binary representations of reducible GF(2) polynomials to even numbers and the corresponding representations of irreducible polynomials to odd numbers, in some order. A246203 has the same property].

A246163 Permutation of natural numbers: a(1) = 1, a(A014580(n)) = A065621(1+a(n)), a(A091242(n)) = A048724(a(n)), where A065621(n) and A048724(n) give the reversing binary representation of n and -n, respectively, and A014580 resp. A091242 are the binary coded irreducible resp. reducible polynomials over GF(2).

Original entry on oeis.org

1, 2, 7, 3, 6, 9, 8, 5, 10, 27, 4, 24, 11, 15, 30, 45, 12, 40, 26, 29, 17, 34, 119, 20, 25, 120, 46, 39, 51, 102, 14, 153, 60, 43, 136, 114, 31, 105, 85, 170, 44, 18, 427, 68, 125, 408, 13, 150, 33, 187, 255, 510, 116, 54, 41, 765, 204, 135, 28, 680, 16, 23, 442, 99, 461, 257, 35, 514, 156, 90, 123, 1799, 118, 340, 393, 36
Offset: 1

Views

Author

Antti Karttunen, Aug 19 2014

Keywords

Comments

This is an instance of entanglement permutation, where the two complementary pairs to be entangled with each other are A014580/A091242 (binary codes for irreducible and reducible polynomials over GF(2)) and A065621/A048724, the latter which themselves are permutations of A000069/A001969 (odious and evil numbers), which means that this permutation shares many properties with A246161.
Because 3 is the only evil number in A014580, it implies that, apart from a(3)=7, odious numbers occur in odious positions only (along with many evil numbers that also occur in odious positions).
Furthermore, all terms of A246158 reside in infinite cycles, and apart from 4 and 8, all of them reside in separate cycles. The infinite cycle containing 4 and 8 goes as: ..., 2091, 97, 47, 13, 11, 4, 3, 7, 8, 5, 6, 9, 10, 27, 46, 408, 2535, ... and it is only because a(3) = 7, that it can temporarily switch back from evil terms to odious terms, until right after a(8) = 5 it is finally doomed to the eternal evilness.
Please see also the comments at A246201 and A246161.

Crossrefs

Inverse: A246164.
Similar or related permutations: A246205, A193231, A246201, A234026, A245701, A234612, A246161, A246203.

Formula

a(1) = 1, and for n > 1, if n is in A014580, a(n) = A065621(1+a(A091226(n))), otherwise a(n) = A048724(a(A091245(n))).
As a composition of related permutations:
a(n) = A193231(A246201(n)).
a(n) = A234026(A245701(n)).
a(n) = A234612(A246161(n)).
a(n) = A193231(A246203(A193231(n))).
Other identities:
For all n > 1, A010060(a(n)) = A091225(n). [Maps binary representations of irreducible GF(2) polynomials (A014580) to odious numbers and the corresponding representations of reducible polynomials (A091242) to evil numbers, in some order].

A091230 Iterates of A014580, starting with a(0) = 1, a(n) = A014580^(n)(1). [Here A014580^(n) means the n-th fold application of A014580].

Original entry on oeis.org

1, 2, 3, 7, 25, 137, 1123, 13103, 204045, 4050293, 99440273
Offset: 0

Views

Author

Antti Karttunen, Jan 03 2004

Keywords

Crossrefs

Programs

Formula

a(0)=1, a(n) = A014580(a(n-1)). [The defining recurrence].
From Antti Karttunen, Aug 03 2014: (Start)
Other identities. For all n >= 0, the following holds:
A091238(a(n)) = n+1.
a(n) = A091204(A007097(n)) and A091205(a(n)) = A007097(n).
a(n) = A245703(A007097(n)) and A245704(a(n)) = A007097(n).
a(n) = A245702(A000079(n)) and A245701(a(n)) = A000079(n).
(End)

Extensions

Terms a(8)-a(10) computed by Antti Karttunen, Aug 02 2014

A246161 Permutation of positive integers: a(1) = 1, a(A014580(n)) = A000069(1+a(n)), a(A091242(n)) = A001969(1+a(n)), where A000069 and A001969 are the odious and evil numbers, and A014580 resp. A091242 are the binary coded irreducible resp. reducible polynomials over GF(2).

Original entry on oeis.org

1, 2, 4, 3, 5, 9, 8, 6, 10, 18, 7, 17, 11, 12, 20, 36, 15, 34, 19, 23, 24, 40, 72, 30, 16, 68, 39, 46, 48, 80, 13, 144, 60, 33, 136, 78, 21, 92, 96, 160, 37, 27, 288, 120, 66, 272, 14, 156, 43, 184, 192, 320, 75, 54, 35, 576, 240, 132, 22, 544, 25, 29, 312, 86, 368, 384, 41
Offset: 1

Views

Author

Antti Karttunen, Aug 17 2014

Keywords

Comments

This is an instance of entanglement permutation, where the two complementary pairs to be entangled with each other are A014580/A091242 (binary codes for irreducible and reducible polynomials over GF(2)) and A000069/A001969 (odious and evil numbers).
Because 3 is the only evil number in A014580, it implies that, apart from a(3)=4, odious numbers occur in odious positions only (along with many evil numbers that also occur in odious positions).
Note that the two values n=21 and n=35 given in the Example section both encode polynomials reducible over GF(2) and have an odd number of 1-bits in their binary representation (that is, they are both terms of A246158). As this permutation maps all terms of A091242 to the terms of A001969, and apart from a single exception 3 (which here is in a closed cycle: a(3) = 4, a(4) = 3), no term of A001969 is a member of A014580, so they must be members of A091242, thus successive iterations a(21), a(a(21)), a(a(a(21))), etc. always yield some evil number (A001969), so the cycle can never come back to 21 as it is an odious number, so that cycle must be infinite.
On the other hand, when we iterate with the inverse of this permutation, A246162, starting from 21, we see that its successive pre-images 37, 41, 67, 203, 5079 [e.g., 21 = a(a(a(a(a(5079)))))] are all irreducible and thus also odious.
In each such infinite cycle, there can be at most one term which is both reducible (in A091242) and odious (in A000069), i.e. in A246158, thus 21 and 35 must reside in different infinite cycles.
The sequence of fixed points begin as: 1, 2, 5, 19, 54, 71, 73, 865.
Question: apart from them and transposition (3 4) are there any more instances of finite cycles?

Examples

			Consider n=21. In binary it is 10101, encoding for polynomial x^4 + x^2 + 1, which factorizes as (x^2 + x + 1)(x^2 + x + 1) over GF(2), in other words, 21 = A048720(7,7). As such, it occurs as the 14th term in A091242, reducible polynomials over GF(2), coded in binary.
By definition of this permutation, a(21) is thus obtained as A001969(1+a(14)). 14 in turn is 8th term in A091242, thus a(14) = A001969(1+a(8)). In turn, 8 = A091242(4), thus a(8) = A001969(1+a(4)), and 4 = A091242(1).
By working the recursion back towards the toplevel, the result is a(21) = A001969(1+A001969(1+A001969(1+A001969(1+1)))) = 24.
Consider n=35. In binary it is 100011, encoding for polynomial x^5 + x + 1, which factorizes as (x^2 + x + 1)(x^3 + x^2 + 1) over GF(2), in other words, 35 = A048720(7,13). As such, it occurs as the 26th term in A091242, thus a(35) = A001969(1+a(26)), and as 26 = A091242(18) and 18 = A091242(12) and 12 = A091242(7) and 7 = A014580(3) [the polynomial x^2 + x + 1 is irreducible over GF(2)], and 3 = A014580(2) and 2 = A014580(1), we obtain the result as a(35) = A001969(1+A001969(1+A001969(1+A001969(1+A000069(1+A000069(1+A000069(2))))))) = 136.
		

Crossrefs

Formula

a(1) = 1, and for n > 1, if n is in A014580, a(n) = A000069(1+a(A091226(n))), otherwise a(n) = A001969(1+a(A091245(n))).
As a composition of related permutations:
a(n) = A233280(A245701(n)).
a(n) = A003188(A246201(n)).
a(n) = A234612(A246163(n)).
Other identities:
For all n > 1, A010060(a(n)) = A091225(n). [Maps binary representations of irreducible GF(2) polynomials (A014580) to odious numbers and the corresponding representations of reducible polynomials (A091242) to evil numbers, in some order].

A246205 Permutation of natural numbers: a(1) = 1, a(A014580(n)) = A117968(a(n)), a(A091242(n)) = A117967(1+a(n)), where A117967 and A117968 give positive and negative parts of inverse of balanced ternary enumeration of integers, and A014580 resp. A091242 are the binary coded irreducible resp. reducible polynomials over GF(2).

Original entry on oeis.org

1, 2, 7, 5, 3, 11, 23, 15, 4, 12, 22, 33, 6, 52, 17, 13, 35, 43, 25, 16, 137, 45, 53, 36, 58, 155, 29, 47, 462, 154, 66, 135, 37, 152, 426, 30, 8, 156, 1273, 428, 24, 148, 460, 41, 423, 1426, 71, 31, 9, 427, 4283, 1410, 34, 431, 75, 1274, 159, 1423, 21, 3707, 194, 99, 44, 10, 1412, 11115, 64, 3850, 38, 1404, 103, 4281, 26, 412, 3722, 49
Offset: 1

Views

Author

Antti Karttunen, Aug 19 2014

Keywords

Crossrefs

Inverse: A246206.
Similar or related entanglement permutations: A246163, A245701, A246201, A246207, A246209.

Formula

a(1) = 1, and for n > 1, if A091225(n) = 1 [i.e. n is in A014580], a(n) = A117968(a(A091226(n))), otherwise a(n) = A117967(1+a(A091245(n))).
As a composition of related permutations:
a(n) = A246207(A245701(n)).
a(n) = A246209(A246201(n)).
Showing 1-9 of 9 results.