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.

Previous Showing 51-60 of 93 results. Next

A060837 List the positive rationals in the canonical order A020652(n)/A020653(n) and apply the Sagher map to turn them into integers.

Original entry on oeis.org

1, 2, 4, 3, 9, 8, 12, 18, 16, 5, 25, 6, 20, 72, 48, 50, 36, 7, 45, 75, 49, 32, 28, 80, 200, 98, 64, 27, 63, 147, 81, 10, 108, 288, 112, 150, 180, 392, 192, 162, 100, 11, 175, 245, 121, 24, 44, 90, 432, 800, 252, 294, 320, 648, 300, 242, 144, 13, 99, 675, 405, 363, 169, 14
Offset: 1

Views

Author

N. J. A. Sloane, Jun 19 2002

Keywords

Comments

The Sagher map sends Product p_i^e_i / Product q_i^f_i (p_i and q_i being distinct primes) to Product p_i^(2e_i) * Product q_i^(2f_i-1). This map is multiplicative.

Examples

			The first few rationals and their images are 1/1 -> 1, 1/2 -> 2, 2/1 -> 4, 1/3 -> 3, 3/1 -> 9, 1/4 -> 8, ...
		

Crossrefs

Programs

  • Haskell
    a060837 n = (a020652 n ^ 2) *
       product (zipWith (^) (a027748_row m)
                            (map ((subtract 1) . (* 2)) (a124010_row m)))
       where m = a020653 n
    -- Reinhard Zumkeller, Feb 16 2014

Formula

a(n) = A020652(n)^2 * product(A027748(m,k)^(2*A124010(m,k)-1): m=a020653(n), k=1..A000005(m)). - Reinhard Zumkeller, Feb 16 2014

Extensions

More terms from Antonio G. Astudillo (afg_astudillo(AT)hotmail.com), Jan 12 2003
Corrected by Charles R Greathouse IV, Sep 02 2009
Definition slightly changed by Reinhard Zumkeller, Feb 16 2014

A060940 Triangle in which n-th row gives the phi(n) terms appearing as initial primes in arithmetic progressions with difference n, with initial term equal to the smallest positive residue coprimes to n.

Original entry on oeis.org

2, 3, 7, 5, 5, 7, 11, 7, 13, 19, 7, 11, 29, 23, 17, 11, 19, 13, 17, 11, 13, 23, 19, 11, 13, 23, 43, 17, 11, 13, 17, 19, 23, 13, 47, 37, 71, 17, 29, 19, 31, 43, 13, 17, 19, 23, 53, 41, 29, 17, 31, 19, 59, 47, 61, 23, 37, 103, 29, 17, 19, 23, 53, 41, 31, 17, 19, 37, 23, 41, 43, 29
Offset: 1

Views

Author

Labos Elemer, May 07 2001

Keywords

Examples

			For differences 1, 2, 3, 4, 5, 6, 7, .. the initial primes are 2; 3; 7, 5; 5, 7; 11, 7, 13, 19; 7, 11; 29, 23, 17, 11, 19, 13; ... etc. Suitable initial values (coprimes to difference) are in A038566. Position of end(start) of rows is given by values of A002088.
From _Seiichi Manyama_, Apr 02 2018: (Start)
   n | phi(n)|
  ---+-------+------------------------
   1 |   1   |  2;
   2 |   1   |  3;
   3 |   2   |  7,  5;
   4 |   2   |  5,  7;
   5 |   4   | 11,  7, 13, 19;
   6 |   2   |  7, 11;
   7 |   6   | 29, 23, 17, 11, 19, 13;
   8 |   4   | 17, 11, 13, 23;
   9 |   6   | 19, 11, 13, 23, 43, 17;
  10 |   4   | 11, 13, 17, 19;         (End)
		

Crossrefs

A069915 Sum of (1+phi)-divisors of n (cf. A061389).

Original entry on oeis.org

1, 3, 4, 3, 6, 12, 8, 7, 4, 18, 12, 12, 14, 24, 24, 11, 18, 12, 20, 18, 32, 36, 24, 28, 6, 42, 13, 24, 30, 72, 32, 31, 48, 54, 48, 12, 38, 60, 56, 42, 42, 96, 44, 36, 24, 72, 48, 44, 8, 18, 72, 42, 54, 39, 72, 56, 80, 90, 60, 72, 62, 96, 32, 35, 84, 144, 68, 54, 96, 144, 72
Offset: 1

Views

Author

Vladeta Jovovic, Apr 23 2002

Keywords

Crossrefs

Programs

  • Haskell
    a069915 n = product $ zipWith sum_1phi (a027748_row n) (a124010_row n)
       where sum_1phi p e = 1 + sum [p ^ k | k <- a038566_row e]
    -- Reinhard Zumkeller, Mar 13 2012
    
  • Mathematica
    a[1] = 1; a[p_?PrimeQ] = p + 1; a[n_] := Times @@ (1 + Sum[If[GCD[k, Last[#]] == 1, First[#]^k, 0], {k, 1, Last[#]}] & ) /@ FactorInteger[n]; Table[a[n], {n, 1, 71}] (* Jean-François Alcover, May 04 2012 *)
  • PARI
    a(n) = {my(f = factor(n)); prod(i = 1, #f~, 1 + sum(k = 1, f[i, 2], (gcd(k, f[i, 2]) == 1) * f[i, 1]^k));} \\ Amiram Eldar, Aug 15 2023

Formula

Multiplicative with a(p^e) = 1+Sum_{k=1..e, gcd(k, e)=1} p^k.

A071912 a(0) = 0, a(1) = 1; to get a(n+1) for n >= 1, let m = a(n) and consider in turn the numbers k = m-1, m-2, ..., 2, 1, m+1, m+2, m+3, ... until reach a k such that gcd(m,k) = 1 and m/k is different from all a(i)/a(i+1) for i = 0, ..., n-1.

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 2, 3, 1, 4, 3, 4, 1, 5, 4, 5, 3, 5, 2, 5, 1, 6, 5, 6, 1, 7, 6, 7, 5, 7, 4, 7, 3, 7, 2, 7, 1, 8, 7, 8, 5, 8, 3, 8, 1, 9, 8, 9, 7, 9, 5, 9, 4, 9, 2, 9, 1, 10, 9, 10, 7, 10, 3, 10, 1, 11, 10, 11, 9, 11, 8, 11, 7, 11, 6, 11, 5, 11, 4, 11, 3, 11, 2, 11, 1, 12, 11, 12, 7, 12
Offset: 0

Views

Author

N. J. A. Sloane, Jun 13 2002

Keywords

Comments

A version of a greedy construction of an integer-valued function a such that a(n)/a(n+1) runs through all nonnegative rationals exactly once.
After initial 0, odd-indexed terms are integers in order with m repeated phi(m) times; even-indexed terms are the corresponding numbers <= m and relatively prime to m, in descending order. - Franklin T. Adams-Watters, Dec 06 2006

Examples

			After [0 1 1 2 1 3 2] we have seen the fractions 0/1, 1/1, 1/2, 2/1, 1/3, 3/2; we consider k = 1, 3, 4, 5, ...; the first of these that gives a new ratio is k=3, giving 2/3, so the next term is 3.
		

Crossrefs

Cf. A002487.
Bisections: A038567 and essentially A020653.
Cf. A038566.

Programs

  • Haskell
    Following Franklin T. Adams-Watters's comment.
    import Data.List (transpose)
    a071912 n = a071912_list !! n
    a071912_list = 0 : concatMap f [1..] where
       f x = concat $ transpose [take (length tots) $ repeat x, reverse tots]
             where tots = a038566_row x
    -- Reinhard Zumkeller, Dec 16 2013
  • Mathematica
    a[0] = 0; a[1] = a[2] = 1; a[n_] := a[n] = Module[{m = a[n-1], ff = Table[ a[i]/a[i+1], {i, 0, n-2}], ok}, ok := GCD[m, k] == 1 && FreeQ[ff, m/k]; For[k = m-1, k >= 1, k--, If[ok, Return[k]]]; For[k = m+1, True, k++, If[ok, Return[k]]]]; Table[a[n], {n, 0, 89}] (* Jean-François Alcover, Oct 28 2017 *)

A116084 Number of partitions of 1 into distinct fractions i/j with 1<=i

Original entry on oeis.org

0, 0, 1, 2, 4, 6, 10, 15, 23, 36, 47, 70, 87, 132, 283, 434, 471, 772, 825, 1834, 4368, 5545, 5648, 9923, 16464, 19943, 32323, 75912, 76167, 140801, 141140, 238513, 537696, 598295, 2556064, 4674084, 4674843, 4985385, 9716586, 23983711, 23984970, 48523605, 48525214
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 04 2006

Keywords

Comments

Partial sums of A116085, which is more elementary to compute, cf. examples. Sequence A154888 has an equivalent definition except that i=j is allowed there, which yields the one-term sum 1/1 as an additional possibility, and thus A154888(n) = a(n)+1. Sequence A115855 is also about the same problem but does not require the fractions to be distinct. - M. F. Hasler, Jul 14 2016

Examples

			a(4) = # [1/3+2/3, 1/4+3/4] = 2;
a(5) = a(4) + # [1/5+4/5, 2/5+3/5] = 2 + 2 = 4;
a(6) = a(5) + # [1/6+5/6, 1/6+1/3+1/2] = 4 + 2 = 6.
		

Crossrefs

Equals A154888(n) - 1.

Programs

  • Mathematica
    Table[Length@ Select[Union /@ Flatten[Map[IntegerPartitions[1, {#}, Rest@ Union[Flatten@ TensorProduct[#, 1/#] &@ Range@ n /. {Integer -> 0, k /; k > 1 -> 0}]] &, Range@ n], 1], Total@# == 1 &], {n, 25}] (* Michael De Vlieger, Jul 14 2016, after Robert G. Wilson v at A154888 *)

Formula

A116085(n) = a(n+1) - a(n).
a(n) = Sum_{k=1..n-1} A116085(k), cf. examples. - M. F. Hasler, Jul 14 2016

Extensions

a(24)-a(34) from Don Reble, Jul 13 2016
a(35)-a(41) from Giovanni Resta, Jul 15 2016
a(42)-a(43) from Jinyuan Wang, Dec 12 2024

A189918 Sum of tetrahedral numbers A000292(k), with k in the reduced residue system modulo n.

Original entry on oeis.org

0, 1, 5, 11, 35, 36, 126, 130, 264, 260, 715, 406, 1365, 952, 1530, 1716, 3876, 1830, 5985, 3300, 5796, 5500, 12650, 5460, 15075, 10556, 16965, 12810, 31465, 9920, 40920, 24616, 34650, 30192, 49210, 26106, 82251, 46740, 67158, 47320
Offset: 1

Views

Author

Wolfdieter Lang, May 19 2011

Keywords

Comments

The reduced residue system modulo n used here is the set of numbers k from the set {0,1,...,n-1} which satisfy gcd(k,n)=1. There are phi(n) = A000010(n) such numbers k. Cf. A038566. See also the Apostol reference p. 133, and the Wikipedia link.
This is the m=3 member of a family of sequences, call them rmnS(m) (reduced mod n sum), with entries rmnS(m;n):=sum(binomial(k+m-1,m),0<=k<=n-1 with gcd(k,n)=1), m>=0, n>=1. Recall gcd(0,n)=n.
The members for m=0, 1, and 2 are A000010(n), A023896(n) and A127415(n), respectively, where in the last two the offset for n=1 should be taken as 0 (not 1).

Examples

			a(6) = A000292(1) + A000292(5)= 1 + 35 = 36.
a(6) = (6*8/4!)*(6*8 + 1*6)*((1/2)*(2/3)) = 36.
a(12) = A000292(1) + A000292(5) + A000292(7) + A000292(11) = 1 + 35 + 84 + +286 = 406.
a(12) = (12*14/4!)*(12*14 + 1*6)*((1/2)*(2/3)) = 406.
		

References

  • T. Apostol, Introduction to Analytic Number Theory, Springer, 1986.

Crossrefs

Programs

  • Maple
    A000292 := proc(n) binomial(n+2,3) ; end proc:
    A189918 := proc(n) local a; a := 0 ; for k from 0 to n-1 do if igcd(k,n) = 1 then a := a+A000292(k); end if; end do: a ; end proc:
    seq(A189918(n),n=1..40) ; # R. J. Mathar, Jun 13 2011
  • Mathematica
    a[n_] := Sum[ Boole[GCD[k, n] == 1]*k*(k+1)*(k+2)/6, {k, 0, n-1}]; Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Jul 12 2012 *)
  • PARI
    a(n) = sum(k=0, n-1, if (gcd(n,k)==1, k*(k+1)*(k+2)/6)); \\ Michel Marcus, Feb 01 2016

Formula

a(n) = Sum_{k=0..n-1, gcd(k,n)=1 } * A000292(k), n>=1.
a(n) = (n*(n+2)/4!) *{n*(n+2) + mu(rad(n))*rad(n)} *phi(n)/n, n>=2, with rad(n) = A007947(n) the squarefree kernel of n, mu(n)=A008683(n), and phi(n)= A000010(n).
Note that phi(n)/n = A076512(n)/A109395(n) = phi(rad(n))/rad(n).
Proof by principle of inclusion-exclusion.

A279399 Irregular triangle read by rows. Row n gives the primes of the smallest positive restricted residue system modulo A033949(n).

Original entry on oeis.org

3, 5, 7, 5, 7, 11, 2, 7, 11, 13, 3, 5, 7, 11, 13, 3, 7, 11, 13, 17, 19, 2, 5, 11, 13, 17, 19, 5, 7, 11, 13, 17, 19, 23, 3, 5, 11, 13, 17, 19, 23, 7, 11, 13, 17, 19, 23, 29, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 2, 5, 7, 13, 17, 19, 23, 29, 31, 2, 3, 11, 13, 17, 19, 23, 29, 31, 5, 7, 11, 13, 17, 19, 23, 29, 31, 2, 5, 7, 11, 17, 19, 23, 29, 31, 37, 3, 7, 11, 13, 17, 19, 23, 29, 31, 37
Offset: 1

Views

Author

Wolfdieter Lang, Jan 25 2017

Keywords

Comments

The length of row n is given by A279400(n)
For the restricted residue systems modulo n see A038566. For the primes of A038566 (for n >= 3) see A112484.
The primes of the restricted residue system modulo the (composite) positive numbers without a primitive root, given in A033949, are of interest for the determination of the Dirichlet characters modulo the A033949 numbers. For prime numbers (A000040) or for composite positive numbers that have prime primitive roots (A279398) the Dirichlet characters are determined from those of the prime primitive root.

Examples

			The triangle T(n, k) begins (here N = A033949(n)):
n,   N \ k 1  2  3  4  5  6  7  8  9 10 ...
1,   8:    3  5  7
2,  12:    5  7 11
3,  15:    2  7 11 13
4,  16:    3  5  7 11 13
5,  20:    3  7 11 13 17 19
6,  21:    2  5 11 13 17 19
7,  24:    5  7 11 13 17 19 23
8,  28:    3  5 11 13 17 19 23
9,  30:    7 11 13 17 19 23 29
10, 32:    3  5  7 11 13 17 19 23 29 31
11, 33:    2  5  7 13 17 19 23 29 31
12, 35:    2  3 11 13 17 19 23 29 31
13, 36:    5  7 11 13 17 19 23 29 31
14, 39:    2  5  7 11 17 19 23 29 31 37
15, 40:    3  7 11 13 17 19 23 29 31 37
...
		

Crossrefs

Formula

Row n of T is given by the primes of row A033949(n) of A038566, for n >= 1.
T(n, k) = A112484(A033949(n), k), n >= 1, k = 1..A279400(n).

A282600 a(n) = Sum_(k=1..phi(n)) floor(d_k/2) where d_k are the totatives of n.

Original entry on oeis.org

0, 0, 1, 1, 4, 2, 9, 6, 12, 8, 25, 10, 36, 18, 28, 28, 64, 24, 81, 36, 60, 50, 121, 44, 120, 72, 117, 78, 196, 56, 225, 120, 160, 128, 204, 102, 324, 162, 228, 152, 400, 120, 441, 210, 264, 242, 529, 184, 504, 240, 400, 300, 676, 234, 540, 324, 504, 392, 841, 232
Offset: 1

Views

Author

Michel Marcus, Feb 19 2017

Keywords

Comments

The totatives of n are the numbers k <= n with gcd(k,n) = 1.

Crossrefs

Programs

  • PARI
    a(n) = sum(k=1, n, (k\2)*(gcd(n, k)==1));

Formula

If n is odd, a(n) = A023896(n)/2 - A000010(n)/4;
If n is even, a(n) = A023896(n)/2 - A000010(n)/2.

A322937 Triangular array in which the n-th row lists the primes strongly prime to n (in ascending order). For the empty rows n = 1, 2, 3, 4 and 6 we set by convention 0.

Original entry on oeis.org

0, 0, 0, 0, 3, 0, 5, 3, 5, 5, 7, 7, 3, 7, 5, 7, 5, 7, 11, 3, 5, 11, 11, 13, 7, 11, 13, 3, 5, 7, 11, 13, 5, 7, 11, 13, 5, 7, 11, 13, 17, 3, 7, 11, 13, 17, 11, 13, 17, 19, 5, 13, 17, 19, 3, 5, 7, 13, 17, 19, 5, 7, 11, 13, 17, 19, 7, 11, 13, 17, 19, 23
Offset: 1

Views

Author

Peter Luschny, Apr 01 2019

Keywords

Comments

A number k is strongly prime to n if and only if k <= n is prime to n and k does not divide n-1. See the link to 'Strong Coprimality'. (Our terminology follows the plea of Knuth, Graham and Patashnik in Concrete Mathematics, p. 115.)

Examples

			The length of row n is A181834(n). The triangular array starts:
[1] {}
[2] {}
[3] {}
[4] {}
[5] {3}
[6] {}
[7] {5}
[8] {3, 5}
[9] {5, 7}
[10] {7}
[11] {3, 7}
[12] {5, 7}
[13] {5, 7, 11}
[14] {3, 5, 11}
[15] {11, 13}
[16] {7, 11, 13}
[17] {3, 5, 7, 11, 13}
[18] {5, 7, 11, 13}
[19] {5, 7, 11, 13, 17}
[20] {3, 7, 11, 13, 17}
		

Crossrefs

Programs

  • Maple
    Primes := n -> select(isprime, {$1..n}):
    StrongCoprimes := n -> select(k->igcd(k, n)=1, {$1..n}) minus numtheory:-divisors(n-1):
    StrongCoprimePrimes := n -> Primes(n) intersect StrongCoprimes(n):
    A322937row := proc(n) if n in {1, 2, 3, 4, 6} then return 0 else op(StrongCoprimePrimes(n)) fi end:
    seq(A322937row(n), n=1..25);
  • Mathematica
    Table[Select[Prime@ Range@ PrimePi@ n, And[GCD[#, n] == 1, Mod[n - 1, #] != 0] &] /. {} -> {0}, {n, 25}] // Flatten (* Michael De Vlieger, Apr 01 2019 *)
  • Sage
    def primes_primeto(n):
        return [p for p in prime_range(n) if gcd(p, n) == 1]
    def primes_strongly_primeto(n):
        return [p for p in set(primes_primeto(n)) - set((n-1).divisors())]
    def A322937row(n):
        if n in [1, 2, 3, 4, 6]: return [0]
        return sorted(primes_strongly_primeto(n))
    for n in (1..25): print(A322937row(n))

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).
Previous Showing 51-60 of 93 results. Next