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

A297845 Encoded multiplication table for polynomials in one indeterminate with nonnegative integer coefficients. Symmetric square array T(n, k) read by antidiagonals, n > 0 and k > 0. See comment for details.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 5, 4, 1, 1, 5, 9, 9, 5, 1, 1, 6, 7, 16, 7, 6, 1, 1, 7, 15, 25, 25, 15, 7, 1, 1, 8, 11, 36, 11, 36, 11, 8, 1, 1, 9, 27, 49, 35, 35, 49, 27, 9, 1, 1, 10, 25, 64, 13, 90, 13, 64, 25, 10, 1, 1, 11, 21, 81, 125, 77, 77, 125, 81
Offset: 1

Views

Author

Rémy Sigrist, Jan 10 2018

Keywords

Comments

For any number n > 0, let f(n) be the polynomial in a single indeterminate x where the coefficient of x^e is the prime(1+e)-adic valuation of n (where prime(k) denotes the k-th prime); f establishes a bijection between the positive numbers and the polynomials in a single indeterminate x with nonnegative integer coefficients; let g be the inverse of f; T(n, k) = g(f(n) * f(k)).
This table has many similarities with A248601.
For any n > 0 and m > 0, f(n * m) = f(n) + f(m).
Also, f(1) = 0 and f(2) = 1.
The function f can be naturally extended to the set of positive rational numbers: if r = u/v (not necessarily in reduced form), then f(r) = f(u) - f(v); as such, f is a homomorphism from the multiplicative group of positive rational numbers to the additive group of polynomials of a single indeterminate x with integer coefficients.
See A297473 for the main diagonal of T.
As a binary operation, T(.,.) is related to A306697(.,.) and A329329(.,.). When their operands are terms of A050376 (sometimes called Fermi-Dirac primes) the three operations give the same result. However the rest of the multiplication table for T(.,.) can be derived from these results because T(.,.) distributes over integer multiplication (A003991), whereas for A306697 and A329329, the equivalent derivation uses distribution over A059896(.,.) and A059897(.,.) respectively. - Peter Munn, Mar 25 2020
From Peter Munn, Jun 16 2021: (Start)
The operation defined by this sequence can be extended to be the multiplicative operator of a ring over the positive rationals that is isomorphic to the polynomial ring Z[x]. The extended function f (described in the author's original comments) is the isomorphism we use, and it has the same relationship with the extended operation that exists between their unextended equivalents.
Denoting this extension of T(.,.) as t_Q(.,.), we get t_Q(n, 1/k) = t_Q(1/n, k) = 1/T(n, k) and t_Q(1/n, 1/k) = T(n, k) for positive integers n and k. The result for other rationals is derived from the distributive property: t_Q(q, r*s) = t_Q(q, r) * t_Q(q, s); t_Q(q*r, s) = t_Q(q, s) * t_Q(r, s). This may look unusual because standard multiplication of rational numbers takes on the role of the ring's additive group.
There are many OEIS sequences that can be shown to be a list of the integers in an ideal of this ring. See the cross-references.
There are some completely additive sequences that similarly define by extension completely additive functions on the positive rationals that can be shown to be homomorphisms from this ring onto the integer ring Z, and these functions relate to some of the ideals. For example, the extended function of A048675, denoted A048675_Q, maps i/j to A048675(i) - A048675(j) for positive integers i and j. For any positive integer k, the set {r rational > 0 : k divides A048675_Q(r)} forms an ideal of the ring; for k=2 and k=3 the integers in this ideal are listed in A003159 and A332820 respectively.
(End)

Examples

			Array T(n, k) begins:
  n\k|  1   2   3    4    5    6    7     8    9    10
  ---+------------------------------------------------
    1|  1   1   1    1    1    1    1     1    1     1  -> A000012
    2|  1   2   3    4    5    6    7     8    9    10  -> A000027
    3|  1   3   5    9    7   15   11    27   25    21  -> A003961
    4|  1   4   9   16   25   36   49    64   81   100  -> A000290
    5|  1   5   7   25   11   35   13   125   49    55  -> A357852
    6|  1   6  15   36   35   90   77   216  225   210  -> A191002
    7|  1   7  11   49   13   77   17   343  121    91
    8|  1   8  27   64  125  216  343   512  729  1000  -> A000578
    9|  1   9  25   81   49  225  121   729  625   441
   10|  1  10  21  100   55  210   91  1000  441   550
From _Peter Munn_, Jun 24 2021: (Start)
The encoding, n, of polynomials, f(n), that is used for the table is further described in A206284. Examples of encoded polynomials:
   n      f(n)        n           f(n)
   1         0       16              4
   2         1       17            x^6
   3         x       21        x^3 + x
   4         2       25           2x^2
   5       x^2       27             3x
   6     x + 1       35      x^3 + x^2
   7       x^3       36         2x + 2
   8         3       49           2x^3
   9        2x       55      x^4 + x^2
  10   x^2 + 1       64              6
  11       x^4       77      x^4 + x^3
  12     x + 2       81             4x
  13       x^5       90   x^2 + 2x + 1
  15   x^2 + x       91      x^5 + x^3
(End)
		

Crossrefs

Row n: n=1: A000012, n=2: A000027, n=3: A003961, n=4: A000290, n=5: A357852, n=6: A191002, n=8: A000578.
Main diagonal: A297473.
Functions f satisfying f(T(n,k)) = f(n) * f(k): A001222, A048675 (and similarly, other rows of A104244), A195017.
Powers of k: k=3: A000040, k=4: A001146, k=5: A031368, k=6: A007188 (see also A066117), k=7: A031377, k=8: A023365, k=9: main diagonal of A329050.
Integers in the ideal of the related ring (see Jun 2021 comment) generated by S: S={3}: A005408, S={4}: A000290\{0}, S={4,3}: A003159, S={5}: A007310, S={5,4}: A339690, S={6}: A325698, S={6,4}: A028260, S={7}: A007775, S={8}: A000578\{0}, S={8,3}: A191257, S={8,6}: A332820, S={9}: A016754, S={10,4}: A340784, S={11}: A008364, S={12,8}: A145784, S={13}: A008365, S={15,4}: A345452, S={15,9}: A046337, S={16}: A000583\{0}, S={17}: A008366.
Equivalent sequence for polynomial composition: A326376.
Related binary operations: A003991, A306697/A059896, A329329/A059897.

Programs

  • PARI
    T(n,k) = my (f=factor(n), p=apply(primepi, f[, 1]~), g=factor(k), q=apply(primepi, g[, 1]~)); prod (i=1, #p, prod(j=1, #q, prime(p[i]+q[j]-1)^(f[i, 2]*g[j, 2])))

Formula

T is completely multiplicative in both parameters:
- for any n > 0
- and k > 0 with prime factorization Prod_{i > 0} prime(i)^e_i:
- T(prime(n), k) = T(k, prime(n)) = Prod_{i > 0} prime(n + i - 1)^e_i.
For any m > 0, n > 0 and k > 0:
- T(n, k) = T(k, n) (T is commutative),
- T(m, T(n, k)) = T(T(m, n), k) (T is associative),
- T(n, 1) = 1 (1 is an absorbing element for T),
- T(n, 2) = n (2 is an identity element for T),
- T(n, 2^i) = n^i for any i >= 0,
- T(n, 4) = n^2 (A000290),
- T(n, 8) = n^3 (A000578),
- T(n, 3) = A003961(n),
- T(n, 3^i) = A003961(n)^i for any i >= 0,
- T(n, 6) = A191002(n),
- A001221(T(n, k)) <= A001221(n) * A001221(k),
- A001222(T(n, k)) = A001222(n) * A001222(k),
- A055396(T(n, k)) = A055396(n) + A055396(k) - 1 when n > 1 and k > 1,
- A061395(T(n, k)) = A061395(n) + A061395(k) - 1 when n > 1 and k > 1,
- T(A000040(n), A000040(k)) = A000040(n + k - 1),
- T(A000040(n)^i, A000040(k)^j) = A000040(n + k - 1)^(i * j) for any i >= 0 and j >= 0.
From Peter Munn, Mar 13 2020 and Apr 20 2021: (Start)
T(A329050(i_1, j_1), A329050(i_2, j_2)) = A329050(i_1+i_2, j_1+j_2).
T(n, m*k) = T(n, m) * T(n, k); T(n*m, k) = T(n, k) * T(m, k) (T distributes over multiplication).
A104244(m, T(n, k)) = A104244(m, n) * A104244(m, k).
For example, for m = 2, the above formula is equivalent to A048675(T(n, k)) = A048675(n) * A048675(k).
A195017(T(n, k)) = A195017(n) * A195017(k).
A248663(T(n, k)) = A048720(A248663(n), A248663(k)).
T(n, k) = A306697(n, k) if and only if T(n, k) = A329329(n, k).
A007913(T(n, k)) = A007913(T(A007913(n), A007913(k))) = A007913(A329329(n, k)).
(End)

Extensions

New name from Peter Munn, Jul 17 2021

A055777 a(n) = 3^(3^n).

Original entry on oeis.org

3, 27, 19683, 7625597484987, 443426488243037769948249630619149892803
Offset: 0

Views

Author

Henry Bottomley, Jul 12 2000

Keywords

Comments

Next term is too big to include.
a(n+1) = a(n) written in base 3 and read as if in base 27 (and recorded in base 10).
Number of distinct n-ary operators in a ternary logic. - Ross Drewe, Feb 13 2008
The next term has 116 digits. - Harvey P. Dale, Mar 28 2019

Crossrefs

Programs

Formula

a(n) = a(n-1)^3.
Sum_{n>=0} 1/a(n) = A383817. - Amiram Eldar, May 16 2025

A211669 Number of iterations f(f(f(...(n)...))) such that the result is < 2, where f(x) = cube root of x.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2
Offset: 1

Views

Author

Hieronymus Fischer, Apr 30 2012

Keywords

Comments

For the general case of "Number of iterations f(f(f(...(n)...))) such that the result is < q, where f(x) = x^(1/p)", with p > 1, q > 1, the resulting g.f. is g(x) = 1/(1 - x)*Sum_{k>=0} x^(q^(p^k))
= (x^q + x^(q^p) + x^(q^(p^2)) + x^(q^(p^3)) + ...)/(1 - x).
The first term that equals 3 is a(512). - Harvey P. Dale, Jan 02 2015

Examples

			a(n) = 1, 2, 3, 4, 5, ... for n = 2^1, 2^3, 2^9, 2^27, 2^81, ..., i.e., n = 2, 8, 512, 134217728, 2417851639229258349412352, ... = A023365.
		

Crossrefs

Programs

  • Mathematica
    Table[Length[NestWhileList[Surd[#,3]&,n,#>=2&]],{n,90}]-1 (* Harvey P. Dale, Jan 02 2015 *)
  • PARI
    a(n,c=0)={while(n>=2, n=sqrtnint(n,3); c++);c} \\ M. F. Hasler, Dec 07 2018

Formula

a(2^(3^n)) = a(2^(3^(n-1))) + 1, for n >= 1.
G.f.: 1/(1-x)*Sum_{k>=0} x^(2^(3^k))
= (x^2 + x^8 + x^512 + x^134217728 + ...)/(1 - x).

Extensions

Edited by M. F. Hasler, Dec 07 2018

A376867 Reduced numerators of Newton's iteration for 1/sqrt(2), starting with 1/2.

Original entry on oeis.org

1, 5, 355, 94852805, 1709678476417571835487555, 9994796326591347130392203807311551183419838794447313956622219314498503205
Offset: 0

Views

Author

Steven Finch, Oct 07 2024

Keywords

Comments

An explicit formula for a(n) is not known, although it arises from a recurrence and the corresponding denominators are simply 2^(3^n) = A023365(n+1).
Next term is too large to include.

Examples

			a(1) = 5 because b(1) = (1/2)*(3/2 - 1/4) = 5/8.
1/2, 5/8, 355/512, 94852805/134217728, ... = a(n)/A023365(n+1).
		

Crossrefs

Programs

  • Maple
    b:= proc(n) b(n):= `if`(n=0, 1/2, b(n-1)*(3/2-b(n-1)^2)) end:
    a:= n-> numer(b(n)):
    seq(a(n), n=0..5);  # Alois P. Heinz, Oct 07 2024
  • Mathematica
    a[0]=1/2; a[n_]:=a[n-1](3/2-a[n-1]^2); Numerator[Array[a,6,0]] (* Stefano Spezia, Oct 15 2024 *)
  • Python
    from itertools import count, islice
    def A376867_gen(): # generator of terms
        p = 1
        for k in count(0):
            yield p
            p *= ((3<<((3**k<<1)-1))-p**2)
    A376867_list = list(islice(A376867_gen(),6)) # Chai Wah Wu, Oct 11 2024

Formula

a(n) is the reduced numerator of b(n) = b(n-1)*(3/2 - b(n-1)^2); b(0) = 1/2.
Limit_{n -> oo} a(n)/A023365(n+1) = 1/sqrt(2) = A010503.
a(n+1) = a(n)*(3*2^(2*3^n-1)-a(n)^2). - Chai Wah Wu, Oct 11 2024

A343786 Decimal expansion of Sum_{k>=0} 2^(3^k)/(2^(2*(3^k)) + 2^(3^k) + 1).

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Apr 29 2021

Keywords

Comments

Proven to be a transcendental number by Borwein and Coons (2008).

Examples

			0.39725264457801453528444061025296764648419265353335...
		

Crossrefs

Programs

  • Mathematica
    f[x_] := x/(x^2 + x + 1); RealDigits[Sum[f[2^(3^n)], {n, 0, 10}], 10, 100][[1]]

Formula

Equals Sum_{k>=1} A343785(k)/2^k.

A381705 Length of iteration sequence of shortest unimodal Collatz (3x+1)/2 sequence that begins with exactly n increases and ends with continuous decreases until reaching 1.

Original entry on oeis.org

3, 6, 13, 32, 87, 250, 737, 2196, 6571, 19694, 59061, 177160, 531455, 1594338, 4782985, 14348924, 43046739, 129140182, 387420509, 1162261488, 3486784423, 10460353226, 31381059633, 94143178852, 282429536507, 847288609470, 2541865828357, 7625597485016, 22876792454991
Offset: 1

Views

Author

David Dewan, Mar 04 2025

Keywords

Comments

A unimodal Collatz sequence has one peak because it starts with only odd numbers (which increase) followed by only even numbers (which decrease). It uses the rule odd x -> (3x+1)/2.
A sequence of length a(n) starts with exactly n odd numbers and ends with 3^(n-1) even numbers and the final 1 for a total length of n + 3^(n-1) + 1.
The peak of a given sequence is 2^(3^(n-1)). See A023365.

Examples

			For n=2, the shortest unimodal sequence has length a(2) = 6 terms and one such sequence is
  3 -> 5 ->  8  -> 4 -> 2 -> 1
    \-----/     \----------/
  2 increases, then decreases
		

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{5,-7,3},{3,6,13},29] (* James C. McMahon, Apr 02 2025 *)

Formula

a(n) = n + 3^(n-1) + 1.
From Stefano Spezia, Mar 07 2025: (Start)
G.f.: x*(3 - 9*x + 4*x^2)/((1 - x)^2*(1 - 3*x)).
E.g.f.: (exp(3*x) + 3*exp(x)*(1 + x) - 4)/3. (End)

A383420 Maximum (equal) number of red and blue tiles on an n X n matrix, where opposite colors cannot be adjacent diagonally or edgewise, and where a cluster of the same color can be no greater than n.

Original entry on oeis.org

0, 0, 6, 8, 16, 24, 30, 38
Offset: 1

Views

Author

Jakub Buczak, Apr 26 2025

Keywords

Comments

a(n) refers to the total number of both red and blue tiles covering the n X n matrix, and thus all the terms are even.

Examples

			a(2) = 0, since red and blue tiles, regardless of the arrangement will always either border sideways or diagonally.
a(3) = 6, a simple example could be:
  [R R R]
  [     ]
  [B B B]
		

Crossrefs

Formula

The lower and upper bounds apply to all n > 3.
k is the root of the smallest square greater than n, b = ceiling((n+1)/(k+1)). The variable x is defined as 5*((n+1)/(k+1) - b), if (n+1)/(k+1) - b > 0, otherwise x=0.
a(n) <= (b*k)^2 + n^2 - (n - x)^2 for all n.
a(n) => n^2/2, for n == 0 mod 4.
a(n) => (n^2 + n - 2)/2, for n == 1 mod 4.
a(n) => n^2/2 + n - 4, for n == 2 mod 4.
a(n) => (n^2 + n - 4)/2, for n == 3 mod 4.

A185981 a(n) = 2^(2*(5^(n-1) - 1)).

Original entry on oeis.org

1, 256, 281474976710656, 452312848583266388373324160190187140051835877600158453279131187530910662656
Offset: 1

Views

Author

Wolfdieter Lang, Feb 24 2011

Keywords

Comments

The number of digits of a(n) is 1, 3, 15, 75, 376, 1881, 9407, 47036, 235180, 1175898, ....
-1/(4*a(n)) is the coefficient of x^0 of the minimal polynomial Psi(5^n,x) of cos(2*Pi/5^n). Hence 4*a(n)*Psi(5^n,x) is the integer polynomial with coefficient -1 of x^0. E.g., Psi(5,1)= x^2 + (1/2)*x -1/4, Psi(25,x)= x^10 + ... -1/1024. See A181875/A181876, A181877 and the W. Lang link under A181875.

Crossrefs

Cf. A023365.

Programs

Formula

a(n) = 2^(2*(5^(n-1) - 1)).
Showing 1-8 of 8 results.