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-10 of 19 results. Next

A291755 Compound filter (multiplicative order of 2 mod 2n+1 & eulerphi(2n+1)): a(n) = P(A002326(n), A037225(n)), where P(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

1, 5, 25, 31, 61, 181, 265, 59, 261, 613, 142, 507, 761, 613, 1513, 566, 416, 607, 2521, 607, 1731, 1499, 607, 2301, 1912, 749, 5305, 1731, 1396, 6613, 7081, 826, 1723, 8581, 2102, 5391, 3169, 1731, 3946, 6709, 5725, 13285, 2493, 3431, 4764, 3415, 2356, 5707, 10201, 3946, 19801, 11527
Offset: 0

Views

Author

Antti Karttunen, Oct 02 2017

Keywords

Crossrefs

Cf. A000010, A000027, A002326, A037225, A291766 (rgs-version of this filter).
Cf. also A292249, A292268.

Programs

Formula

a(n) = (1/2)*(2 + ((A002326(n) + A000010(2n+1))^2) - A002326(n) - 3*A000010(2n+1)).

A291766 Restricted growth sequence transform of A291755; filter combining multiplicative order of 2 mod 2n+1 & eulerphi(2n+1) (A002326 & A037225).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 10, 14, 15, 16, 17, 18, 17, 19, 20, 17, 21, 22, 23, 24, 19, 25, 26, 27, 28, 29, 30, 31, 32, 33, 19, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 34, 45, 46, 29, 47, 48, 43, 49, 50, 41, 51, 52, 53, 45, 54, 55, 56, 57, 43, 58, 59, 60, 61, 49, 62, 63, 64, 51, 65, 66, 67, 68, 69, 53, 70, 71, 57, 72, 61, 73, 74, 75, 61
Offset: 0

Views

Author

Antti Karttunen, Oct 02 2017

Keywords

Crossrefs

Cf. A291769, A292267 for related filters.

Programs

  • PARI
    rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
    write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
    A002326(n) = if(n<0, 0, znorder(Mod(2, 2*n+1))); \\ This function from Michael Somos, Mar 31 2005
    A291755(n) = (1/2)*(2 + ((A002326(n)+eulerphi(n+n+1))^2) - A002326(n) - 3*eulerphi(n+n+1));
    write_to_bfile(0,rgs_transform(vector(32769,n,A291755(n-1))),"b291766_upto32768.txt");

A003972 Moebius transform of A003961; a(n) = phi(A003961(n)), where A003961 shifts the prime factorization of n one step towards the larger primes.

Original entry on oeis.org

1, 2, 4, 6, 6, 8, 10, 18, 20, 12, 12, 24, 16, 20, 24, 54, 18, 40, 22, 36, 40, 24, 28, 72, 42, 32, 100, 60, 30, 48, 36, 162, 48, 36, 60, 120, 40, 44, 64, 108, 42, 80, 46, 72, 120, 56, 52, 216, 110, 84, 72, 96, 58, 200, 72, 180, 88, 60, 60, 144, 66, 72, 200, 486, 96, 96, 70
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    b[1] = 1; b[p_?PrimeQ] := b[p] = Prime[ PrimePi[p] + 1]; b[n_] := b[n] = Times @@ (b[First[#]]^Last[#] &) /@ FactorInteger[n]; a[n_] := Sum[ MoebiusMu[n/d]*b[d], {d, Divisors[n]}]; Table[a[n], {n, 1, 70}]  (* Jean-François Alcover, Jul 18 2013 *)
  • PARI
    A003972(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); eulerphi(factorback(f)); }; \\ Antti Karttunen, Aug 20 2020
    
  • Python
    from math import prod
    from sympy import nextprime, factorint
    def A003972(n): return prod((q:=nextprime(p))**(e-1)*(q-1) for p, e in factorint(n).items()) # Chai Wah Wu, Jul 18 2022

Formula

Multiplicative with a(p^e) = (q-1)q^(e-1) where q = nextPrime(p). - David W. Wilson, Sep 01 2001
a(n) = A000010(A003961(n)) = A037225(A108228(n)) = A037225(A048673(n)-1). - Antti Karttunen, Aug 20 2020
Sum_{k=1..n} a(k) ~ c * n^2, where c = (3/Pi^2) * Product_{p prime} (p^2-p)/(p^2 - nextPrime(p)) = 1.2547593344... . - Amiram Eldar, Nov 18 2022

Extensions

More terms from David W. Wilson, Aug 29 2001
Secondary name added by Antti Karttunen, Aug 20 2020

A062570 a(n) = phi(2*n).

Original entry on oeis.org

1, 2, 2, 4, 4, 4, 6, 8, 6, 8, 10, 8, 12, 12, 8, 16, 16, 12, 18, 16, 12, 20, 22, 16, 20, 24, 18, 24, 28, 16, 30, 32, 20, 32, 24, 24, 36, 36, 24, 32, 40, 24, 42, 40, 24, 44, 46, 32, 42, 40, 32, 48, 52, 36, 40, 48, 36, 56, 58, 32, 60, 60, 36, 64, 48, 40, 66, 64, 44, 48, 70, 48, 72
Offset: 1

Views

Author

Jason Earls, Jul 03 2001

Keywords

Comments

a(n) is also the number of non-congruent solutions to x^2 - y^2 == 1 (mod n). - Yuval Dekel (dekelyuval(AT)hotmail.com), Sep 21 2003
a(n) is the size of a square companion matrix of the minimal cyclotomic polynomial of (-1)^(1/n). - Eric Desbiaux, Dec 08 2015
a(n) is the degree of the (2n)-th cyclotomic field Q(zeta_(2n)). Note that Q(zeta_n) = Q(zeta_(2n)) for odd n. - Jianing Song, May 17 2021
The number of integers k from 1 to n such that gcd(n,k) is a power of 2. - Amiram Eldar, May 18 2025

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, p. 28.

Crossrefs

Column 1 of A129559, column 2 of A372673.
Row 1 of A379010.
Row sums of A129558 and of A129564.

Programs

  • Maple
    [phi(2*n)$n=1..80]; # Muniru A Asiru, Mar 18 2019
  • Mathematica
    Table[EulerPhi[2 n], {n, 80}] (* Vincenzo Librandi, Aug 23 2013 *)
  • PARI
    a(n) = eulerphi(2*n)
    
  • Python
    from sympy import totient
    def A062570(n): return totient(n) if n&1 else totient(n)<<1 # Chai Wah Wu, Aug 04 2024
  • Sage
    [euler_phi(2*n) for n in range(1,74)] # Zerinvary Lajos, Jun 06 2009
    

Formula

a(n) = Sum_{d|n and d is odd} n/d*mu(d).
Multiplicative with a(2^e) = 2^e and a(p^e) = p^e-p^(e-1), p>2.
Dirichlet g.f.: zeta(s-1)/zeta(s)*2^s/(2^s-1). - Ralf Stephan, Jun 17 2007
a(n) = A000010(2*n).
a(n) = phi(n)*(1+((n+1) mod 2)). - Gary Detlefs, Jul 13 2011
a(n) = A173557(n)*b(n) where b(n) = 1, 2, 1, 4, 1, 2, 1, 8, 3, 2, 1, 4, 1, 2, ... is the multiplicative function defined by b(p^e) = p^(e-1) if p<>2 and b(2^e)=2^e. b(n) = n/A204455(n). - R. J. Mathar, Jul 02 2013
a(n) = -c_{2n}(n) where c_q(n) is Ramanujan's sum. - Michael Somos, Aug 23 2013
a(n) = A055034(2*n), for n >= 2. - Wolfdieter Lang, Nov 30 2013
O.g.f.: Sum_{n >= 1} mu(2*n-1)*x^(2*n-1)/(1 - x^(2*n-1))^2. - Peter Bala, Mar 17 2019
a(n) = A000010(4*n)/2, for n > = 1 (see Apostol, Theorem 2.5, (b), p. 28). - Wolfdieter Lang, Nov 17 2019
a(n) = n - Sum_{d|n, n/d odd, d < n} a(d). - Ilya Gutkovskiy, May 30 2020
Dirichlet convolution of A000010 and A209229. - Werner Schulte, Jan 17 2021
From Richard L. Ollerton, May 07 2021: (Start)
a(n) = Sum_{k=1..n} A209229(gcd(n,k)).
a(n) = Sum_{k=1..n} A209229(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). (End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = 4/Pi^2 = 0.405284... (A185199). - Amiram Eldar, Oct 22 2022
a(n) = A000034(n) * A000010(n). - Amiram Eldar, May 18 2025

Extensions

Corrected by Vladeta Jovovic, Dec 04 2002

A072451 Number of odd terms in the reduced residue system of 2*n-1.

Original entry on oeis.org

1, 1, 2, 3, 3, 5, 6, 4, 8, 9, 6, 11, 10, 9, 14, 15, 10, 12, 18, 12, 20, 21, 12, 23, 21, 16, 26, 20, 18, 29, 30, 18, 24, 33, 22, 35, 36, 20, 30, 39, 27, 41, 32, 28, 44, 36, 30, 36, 48, 30, 50, 51, 24, 53, 54, 36, 56, 44, 36, 48, 55, 40, 50, 63, 42, 65, 54, 36, 68, 69, 46, 60, 56
Offset: 1

Views

Author

Labos Elemer, Jun 19 2002

Keywords

Comments

Given the quasi-order and coach theorems of Hilton and Pedersen (cf. A003558): phi(b) = 2 * k * c with b odd. Dividing through by 2 we obtain phi(b)/2 = k * c which is the same as a(n) = A003558(n-1) * A135303(n-1). Here k refers to the "least exponent" (cf. A003558), while c is the number of coaches for odd b (cf. A135303). - Gary W. Adamson, Sep 25 2012
After the initial term, this is the totient function phi(2n-1)/2 (A037225(n-1)/2). Also a(n) is the number of partitions of the odd number (2n+1) into two ordered relatively prime parts. If p and q are such parts where p > q and p+q = 2n+1 then they can generate primitive Pythagorean triples of the form (p^2 - q^2, 2*p*q, p^2 + q^2). - Frank M Jackson, Oct 30 2012

Examples

			n=105: phi(105)=48 with 24 odd, 24 even; for even n=2k reduced residue system consists only of odd terms, so number of odd terms equals phi(n).
With odd integer 33, A072451(17) = 10 = A003558(16) * A135303(16); or 10 = 5 * 2.
		

References

  • Peter Hilton and Jean Pedersen, A Mathematical Tapestry, Demonstrating the Beautiful Unity of Mathematics, Cambridge University Press, 2010, p. 200.

Crossrefs

Programs

  • Maple
    A072451 := n -> ceil(numtheory:-phi(2*n-1)/2):
    seq(A072451(n), n=1..73); # Peter Luschny, Feb 24 2020
  • Mathematica
    gw[x_] := Table[GCD[x, w], {w, 1, x}] rrs[x_] := Flatten[Position[gw[x], 1]] Table[Count[OddQ[rrs[2*w-1]], True], {w, 1, 128}]
    (* Additional programs: *)
    Table[Count[Range[1, #, 2], k_ /; CoprimeQ[k, #]] &[2 n - 1], {n, 73}] (* or *) Array[If[# == 1, #, EulerPhi[2 # - 1]/2] &, 73] (* Michael De Vlieger, Jul 24 2017 *)
  • PARI
    A072451(n) = if(1==n,n,eulerphi(n+n-1)/2); \\ (After Benoit Cloitre's formula) - Antti Karttunen, Jul 24 2017

Formula

a(1) = 1, and for n > 1, a(n) = phi(2n-1)/2. - Benoit Cloitre, Oct 11 2002
It would appear that a(n) = Sum_{k=0..n} abs(Jacobi(k, 2n-2k+1)). - Paul Barry, Jul 20 2005
a(n) = A055034(2*n-1), n >= 1. - Wolfdieter Lang, Feb 07 2020
G.f.: (x + Sum_{n>=1} mu(2n-1) * x^n * (1 + x^(2n-1)) / (1 - x^(2n-1))^2) / 2. - Mamuka Jibladze, Dec 13 2022
Sum_{k=1..n} a(k) ~ c * n^2, where c = 4/Pi^2 = 0.405284... (A185199). - Amiram Eldar, Feb 11 2023

A037226 a(n) = phi(2n+1) / multiplicative order of 2 mod 2n+1.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 6, 2, 2, 1, 2, 2, 3, 2, 2, 2, 4, 1, 2, 2, 1, 1, 6, 4, 1, 2, 2, 8, 2, 2, 2, 1, 1, 8, 2, 8, 6, 6, 2, 2, 2, 1, 2, 4, 1, 3, 2, 4, 2, 6, 4, 1, 4, 1, 18, 6, 1, 6, 2, 2, 1, 2, 2, 4, 2, 1, 10, 4, 6, 3, 2, 4
Offset: 0

Views

Author

Keywords

Comments

Number of primitive irreducible factors of x^(2n+1) - 1 over integers mod 2. There are no primitive irreducible factors for x^(2n)-1 because it always has the same factors as x^n-1. Considering that A000374 also counts the cycles in the map f(x) = 2x mod n, a(n) is also the number of primitive cycles of that mapping. - T. D. Noe, Aug 01 2003
Equals number of irreducible factors of the cyclotomic polynomial Phi(2n+1,x) over Z/2Z. All factors have the same degree. - T. D. Noe, Mar 01 2008

Crossrefs

Cf. A000374 (number of irreducible factors of x^n - 1 over integers mod 2), A081844.
Cf. A006694 (cyclotomic cosets of 2 mod 2n+1).

Programs

Formula

a(n) = Sum{d|2n+1} mu((2n+1)/d) A000374(d), the inverse Mobius transform of A000374 - T. D. Noe, Aug 01 2003
a(n) = A037225(n)/A002326(n).

A099957 a(n) = Sum_{k=0..n-1} phi(2k+1).

Original entry on oeis.org

1, 3, 7, 13, 19, 29, 41, 49, 65, 83, 95, 117, 137, 155, 183, 213, 233, 257, 293, 317, 357, 399, 423, 469, 511, 543, 595, 635, 671, 729, 789, 825, 873, 939, 983, 1053, 1125, 1165, 1225, 1303, 1357, 1439, 1503, 1559, 1647, 1719, 1779, 1851, 1947
Offset: 1

Views

Author

Hugo Pfoertner, Nov 13 2004

Keywords

Comments

The n-th term is the number of notes of the (2n-1)-limit tonality diamond. This is a term from music theory and means the scale consisting of the rational numbers r, 1 <= r < 2, such that the odd part of both the numerator and the denominator of r, when reduced to lowest terms, is less than or equal to the fixed odd number 2n-1. - Gene Ward Smith, Mar 27 2006
(1/4)*Number of distinct angular positions under which an observer positioned at the center of a square of a square lattice can see the (2n) X (2n) points symmetrically surrounding his position.
(1/8)*number of distinct angular positions under which an observer positioned at a lattice point of a square lattice can see the (2n+1)X(2n+1) points symmetrically surrounding his position gives A002088.
(1/2)*number of distinct angular positions under which an observer positioned at the center of an edge of a square lattice can see the (2n)X(2n-1) points symmetrically surrounding his position gives A099958.

Crossrefs

Bisection of A274401.
Partial sums of A037225.

Programs

  • Mathematica
    Accumulate[EulerPhi[2*Range[0,50]+1]] (* Harvey P. Dale, Aug 20 2021 *)
  • PARI
    apply( {A099957(n)=sum(k=1,n, eulerphi(2*k-1))}, [1..55]) \\ M. F. Hasler, Apr 03 2023

Formula

a(n+1) - a(n) = phi(2n+1) (A037225).
a(n) = (8/Pi^2)*n^2 + O(n^(3/2+eps)) (Lemma 1 in Lv Chuan, 2004). - Amiram Eldar, Aug 02 2022, corrected by M. F. Hasler, Mar 26 2023
a(n) = A002088(2*n-1) - A049690(n-1). - Chai Wah Wu, Aug 04 2024

A325596 a(n) = Sum_{d|n} mu(n/d) * (-1)^(d + 1) * d.

Original entry on oeis.org

1, -3, 2, -2, 4, -6, 6, -4, 6, -12, 10, -4, 12, -18, 8, -8, 16, -18, 18, -8, 12, -30, 22, -8, 20, -36, 18, -12, 28, -24, 30, -16, 20, -48, 24, -12, 36, -54, 24, -16, 40, -36, 42, -20, 24, -66, 46, -16, 42, -60, 32, -24, 52, -54, 40, -24, 36, -84, 58, -16, 60, -90, 36, -32, 48
Offset: 1

Views

Author

Ilya Gutkovskiy, Sep 07 2019

Keywords

Comments

Moebius transform of A181983.

Crossrefs

Programs

  • Magma
    [&+[MoebiusMu(Floor(n/d))*(-1)^(d+1)*d:d in Divisors(n)]:n in [1..70]]; // Marius A. Burtea, Sep 07 2019
  • Mathematica
    a[n_] := Sum[MoebiusMu[n/d] (-1)^(d + 1) d, {d, Divisors[n]}]; Table[a[n], {n, 1, 65}]
    a[n_] := If[OddQ[n], EulerPhi[n], EulerPhi[n] - 4 EulerPhi[n/2]]; Table[a[n], {n, 1, 65}]
    nmax = 65; CoefficientList[Series[Sum[MoebiusMu[k] x^k/(1 + x^k)^2, {k, 1, nmax}], {x, 0, nmax}], x] // Rest
    f[p_, e_] := (p - 1)*p^(e - 1); f[2, 1] = -3; f[2, e_] := -2^(e - 1); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Nov 15 2022 *)
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*(-1)^(d+1)*d); \\ Michel Marcus, Sep 07 2019
    

Formula

G.f.: Sum_{k>=1} mu(k) * x^k / (1 + x^k)^2.
G.f. A(x) satisfies: A(x) = x / (1 + x)^2 - Sum_{k>=2} A(x^k).
a(n) = phi(n) if n odd, phi(n) - 4*phi(n/2) if n even, where phi = A000010.
a(n) = A319997(n) - A319998(n).
Multiplicative with a(2) = -3, a(2^e) = -2^(e-1) for e > 1, and a(p^e) = (p-1)*p^(e-1) for p > 2. - Amiram Eldar, Nov 15 2022

A072455 Number of totients in the reduced residue system of 2n-1.

Original entry on oeis.org

1, 2, 3, 4, 4, 6, 7, 4, 8, 9, 7, 11, 10, 8, 13, 14, 9, 11, 16, 10, 17, 18, 9, 20, 19, 13, 22, 17, 15, 25, 26, 14, 21, 28, 16, 29, 30, 14, 23, 31, 19, 33, 27, 19, 35, 28, 22, 29, 37, 19, 38, 39, 16, 41, 42, 26, 44, 33, 26, 38, 41, 27, 36, 47, 29, 49, 43, 22, 51, 52, 32, 43, 40, 27
Offset: 1

Views

Author

Labos Elemer, Jun 19 2002

Keywords

Examples

			For n=31: reduced residue system(31) = {1,...,30} with 15 odd and 15 even numbers. From the odd terms only the term 1 is totient, while from the 15 even terms, 2 terms, {14,26}, are nontotients, so 13 terms are totients. All totients count 1 + 13 = 14, thus a((31+1)/2) = a(16) = 14.
		

Crossrefs

Programs

  • PARI
    a(n) = {my(m = 2*n-1); sum(k = 1, m, gcd(m, k) == 1 && istotient(k));} \\ Amiram Eldar, Nov 07 2024

Formula

a(n) = phi(2*n-1) - A072454(n). [Corrected by Sean A. Irvine, Oct 04 2024]

A337712 Irregular triangle read by rows: row n gives the complete system of cycles of the doubling sequences modulo N = 2*n+1, for n >= 0.

Original entry on oeis.org

1, 2, 1, 2, 4, 3, 1, 2, 4, 3, 6, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 5, 10, 9, 7, 3, 6, 1, 2, 4, 8, 3, 6, 12, 11, 9, 5, 10, 7, 1, 2, 4, 8, 7, 14, 13, 11, 1, 2, 4, 8, 16, 15, 13, 9, 3, 6, 12, 7, 14, 11, 5, 10, 1, 2, 4, 8, 16, 13, 7, 14, 9, 18, 17, 15, 11, 3, 6, 12, 5, 10
Offset: 0

Views

Author

Gary W. Adamson and Wolfdieter Lang, Oct 14 2020

Keywords

Comments

The length of row n is A037225(n), for n >= 0.
The doubling sequence modulo N = 2*n+1, for n >= 0, has entries DS(N, s(N,i), j) = s(N,i)*2^j (mod N), with j >= 0, and certain positive odd integer seeds s(N, i), for i = 1, 2, ..., S(N) = A037226((N-1)/2), where gcd(s(N, i), N) = 1 (restricted seeds modulo N). These doubling sequences are periodic with period length P(N) = A002326((N-1)/2) (order of 2 modulo N). Only the periods (cycles) {DS(N, s(N, i), j)}_{j=0..P(N)-1}, for i = 1, 2, ..., S(N), are listed.
N = 1 (n=0) is special: one takes here the restricted residue system modulo N not as [0] but as [1]. The order of 2 modulo 1 is 1, because 2^1 == 1 (mod 1) (== 0 (mod 1)).
In order to obtain the complete system of doubling sequences one starts with seed s(N, 1) = 1, and if all numbers from the smallest positive reduced residue system modulo N (called RRS(N), given in row N of A038566) are obtained, i.e., if P(N) = #RRS(N) = phi(N) = A000010(N), then the system is complete. Otherwise the smallest missing number from RRS(N) is taken as new seed s(N, 2), etc. until the system is complete. This means that the number of seeds needed is S(N) = phi(N)/P(N) = A037226((N-1)/2)).
The irregular subtriangle where only seed s(N, 1) = 1 has been used is given in A201908. But there 0 (not 1) for N = 1 has been used.
From Gary W. Adamson and Wolfdieter Lang, Dec 15 2020: (Start)
The cycles in row n, for N = 2*n + 1, of period length P(N) = A002326((N-1)/2) give the periods of the iterated doubling function D(x) = frac(2*x) with seeds x = s(N, i)/N, for i = 1, 2, ..., S(N) = A037226((N-1)/2), after multiplication with N. This is the doubling function used in the Devaney reference, pp. 24-25, 27, 125. 132, 171,289.
Each cycle in row n can also be used to find from the base 2 version of its first entry (the seed s = s(N, i)) divided by N the other entries by repeated application of a cyclic left shift by one step (called sigma operation) to the period of the base 2 expression of s/N. E.g., n = 7, N = 15, P(N) = 4, s = 1: (1/15){10->2} = .repeat(0001), then (.repeat(0010)){2->10} = 2/10, (.repeat(0100)){2->10} = 4/10 and (.repeat(1000)){2->10} = 8/15. Similarly for s = 7: from (7/15)_{10->2} = .repeat(0111) one obtains by repeated sigma operations 14/15, 13/15 and 11/15. The proof uses the elementary formulas for the conversion from base 10 to base 2, and the reverse one, from base 2 to base 10. See also a comment on the period length P(N) given in A002326. (End)

Examples

			The irregular triangle T(n, k) begins (cycles are separated by a vertical bar)
n,  N \ k 1 2 3 4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 ...
0,  1:    1
1,  3:    1 2
2,  5:    1 2 4 3
3,  7:    1 2 4|3  6  5
4,  9:    1 2 4 8  7  5
5,  11:   1 2 4 8  5 10  9  7  3  6
6,  13:   1 2 4 8  3  6 12 11  9  5 10  7
7,  15:   1 2 4 8| 7 14 13 11
8,  17:   1 2 4 8 16 15 13  9| 7 14 11  5 10  3  6 12
9,  19:   1 2 4 8 16 13  7 14  9 18 17 15 11  3  6 12 5 10
10, 21:   1 2 4 8 16 11| 5 10 20 19 17 13
11, 23:   1 2 4 8 16  9 18 13  3  6 12| 5 10 20 17 11 22 21 19 15  7 14
12, 25:   1 2 4 8 16  7 14  3  6 12 24 23 21 17  9 18 11 22 19 13
13, 27:   1 2 4 8 16  5 10 20 13 26 25 23 19 11 22 17  7 14
...
n = 14, N = 29:  1 2 4 8 16  3  6 12 24 19  9 18  7 14 28 27 25 21 13 26 23 17  5 10 20 11 22 15,
n = 15, N = 31: 1 2 4 8 16|3 6 12 24 17|5 10 20 9 18|7 14 28 25 19|11 22 13 26 21|15 30 29 27 23.
		

References

  • Robert L. Devaney, A First Course in Chaotic Dynamical Systems, Addison-Wesley., 1992. pp. 24-25, 27, 125, 132, 171, 289. Second edition 2020.

Crossrefs

Cf. A000010, A002326, A037225, A037226, A201908, A038566, A334430 (modified doubling), A337936 (tripling), A339046 (quadrupling).

Programs

  • Mathematica
    Array[Block[{a = {}, k = 2, n = 2 # + 1, m}, m = EulerPhi[n]; While[Length@ Flatten@ a < m, AppendTo[a, Most@ NestWhileList[Mod[2 #, n] &, If[Length@ a == 0, 1, k], UnsameQ, All]]; Set[k, SelectFirst[Complement[Range[n], Union@ Flatten@ a], GCD[#, n] == 1 &] ]]; a] &, 9] // Flatten (* Michael De Vlieger, Nov 06 2020 *)

Formula

T(n, k) gives the k-th entry in the complete doubling system modulo N = 2*n+1, for n >= 0, with the S(N) = A037226((N-1)/2) cycles of length A002326((N-1)/2) written in row n. See the comment above for DS(N,s(N,i)), i = 1, 2, ..., S(N).
Showing 1-10 of 19 results. Next