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.

User: Benoit Jubin

Benoit Jubin's wiki page.

Benoit Jubin has authored 48 sequences. Here are the ten most recent ones:

A306954 Triangle read by rows: T(n, k), for 1 <= k <= n, is the maximum integer q such that k non-attacking armies of q queens can be placed on an n X n chessboard.

Original entry on oeis.org

1, 4, 0, 9, 1, 0, 16, 2, 1, 1, 25, 4, 1, 1, 1, 36, 5, 2, 1, 1, 1, 49, 7
Offset: 1

Author

Benoit Jubin, Mar 17 2019

Keywords

Comments

One has T(n, k) = 0 exactly when (n, k) = (2, 2) or (3, 3).
One has T(n, n) = 1 except when n = 2 or 3 (that is, when A000170(n) = 0).
For a fixed k, the sequence T(-, k) is nondecreasing.
For a fixed n, the sequence T(n, -) is nonincreasing.
For a fixed nonzero p, the sequence (T(k + p, k))_k is nonincreasing. Indeed, given a configuration with k+1 armies on a k+p+1 chessboard, remove the row and column containing a given queen; this row and column can contain only queens of one army, so this yields a configuration of k armies on a k+p chessboard.
One has T(n, 1) = n^2 and T(n, 2) = A250000(n).
For a fixed k, one has, asymptotically in n, that 1/2.(n/k)^2 <= T(n, k) <= (n/k)^2, which can be proved as follows.
For the upper bound, let R(n, k) be defined as T(n, k) with rooks instead of queens. Then, T(n, k) <= R(n, k) ~ (n/k)^2. Indeed, for 1 <= i <= k, say the i^th army controls a_i rows and b_i columns. The sum of the a_i's, as well as the sum of the b_i's, is at most n. The size of the i^th army is at most a_i b_i. Therefore, one wants to maximize min(a_i b_i, i = 1..k). Ignoring rounding, the maximum is reached when all the a_i's and b_i's are equal.
For the lower bound, consider the n X n chessboard as a k X k grid with cells of size (n/k) X (n/k). Consider a configuration of k non-attacking queens on the k X k chessboard as in A000170(k). Place each of the k armies inside one cell occupied in that configuration. The precise placement is as follows: the army occupies the square whose vertices are the midpoints of the sides of the cell, hence has size 1/2.(n/k)^2. These armies are non-attacking.

Examples

			Triangle begins:
   1
   4  0
   9  1  0
  16  2  1  1
  25  4  1  1  1
  36  5  2  1  1  1
  49  7  .  .  1  1  1
  64  9  .  .  .  1  1  1
  81 12  .  .  .  .  1  1  1
		

Crossrefs

Column 1 gives A000290, n >= 1.
Cf. A000170 (non-attacking queens).
Column 2 gives A250000 (case of two armies).

A261070 Irregular triangle read by rows: T(n,k) is the number of arrangements of n circles with 2k intersections (using the same rules as A250001).

Original entry on oeis.org

1, 1, 2, 1, 4, 4, 2, 4, 9, 15, 15, 31, 24, 35, 44, 20, 50
Offset: 0

Author

Benoit Jubin, Aug 08 2015

Keywords

Comments

Length of n-th row: 1 + (n-1)n/2 (for a configuration for T(n,(n-1)n/2), consider n circles of radius 1 and centers at (k/n,0) for 1<=k<=n).
The generating function down the column k=1 is 1+z^2 *C^2(z) *[C^2(z)+C(z^2)]/ (2*[1-z*C(z)]) = 1+ z^2 +4*z^3 +15*z^4+ 50*z^5+...where C(z) = 1+z+2*z^2+4*z^3+... is the g.f. of A000081 divided by z; eq. (78) in arXiv:1603.00077. - R. J. Mathar, Mar 05 2016

Examples

			n\k 0  1  2  3  4  5  6
0   1
1   1
2   2  1
3   4  4  2  4
4   9 15 15 31 24 35 44
5  20 50  .  .  .  .  .  .  .  .  .
		

Crossrefs

Row sums give A250001.
Cf. A000081, A152947, A249752, A252158, A280786 (column k=1)

Formula

A250001(n) = Sum_{k>=0} T(n,k).
A000081(n+1) = T(n,0).

Extensions

T(4,2)..T(5,0) (6 terms) from Travis Vasquez, Nov 28 2024

A246878 a(0) = 1, then a(n) = sum(a(k), k = floor(log_2(n)) .. n - 1).

Original entry on oeis.org

1, 1, 1, 2, 3, 6, 12, 24, 47, 94, 188, 376, 752, 1504, 3008, 6016, 12030, 24060, 48120, 96240, 192480, 384960, 769920, 1539840, 3079680, 6159360, 12318720, 24637440, 49274880, 98549760, 197099520, 394199040, 788398077
Offset: 0

Author

Benoit Jubin, Sep 06 2014

Keywords

Comments

a(n) = Theta(2^n), and more precisely, for n >= 4, (13/16)*(3/16)2^n <= a(n) <= (3/16)*2^n.
Indeed, from the formula, one gets a(n) <= (3/16)*2^n, and injecting this in the formula, one gets a(n) >= 2*a(n - 1) - (3/32)*n. Then by induction, and using the formula sum(k*2^k, k = 1 .. n) = (n - 1)*2^(n + 1) + 2, one obtains a(n) >= (13/16)*(3/16)2^n + (3/32)*n + (3/8).

Examples

			a(2) = a(1) = a(0) = 1.
a(3) = a(2) + a(1) = 2.
a(4) = a(3) + a(2) = 3.
a(5) = a(4) + a(3) + a(2) = 6.
		

Crossrefs

Cf. A000523.

Programs

  • Haskell
    import Data.List (genericDrop)
    a246878 n = a246878_list !! n
    a246878_list = 1 : f [1] a000523_list where
       f xs (k:ks) = y : f (xs ++ [y]) ks where y = sum $ genericDrop k xs
    -- Reinhard Zumkeller, Sep 16 2014

Formula

If n >= 1 is not a power of 2, then a(n) = 2*a(n - 1), and if k >= 1, then a(2^k) = 2*a(2^k - 1) - a(k - 1).

A179101 Numbers which are not the sum of exactly one positive square and a sum of squares-minus-1.

Original entry on oeis.org

2, 3, 5, 6, 8, 11, 14
Offset: 1

Author

Benoit Jubin, Jun 29 2010

Keywords

Comments

Sequence motivated by the study of certain replicate tilings, where each tile can be replaced by a square number of tiles.
Adding multiples of 3=2^2-1 to the numbers 1=1^2, 9=3^2 and 17=3^2+(3^2-1), one obtains all the integers not in the sequence.

Crossrefs

Cf. A000290 (squares), A005563 (squares-minus-1), A123120, A078135.

A162187 a(0)=4 and a(n+1)=a(n)!

Original entry on oeis.org

4, 24, 620448401733239439360000
Offset: 0

Author

Benoit Jubin, Jun 27 2009

Keywords

Comments

Initializing at 0, 1 or 2 gives constant sequences (except for the first term when initializing at 0).
Initializing at 3 gives 3,6,720,720!,...

Crossrefs

Cf. A000142 (factorials).

A165732 Decimal expansion of sum(e^(1/n!)-1,n>=0).

Original entry on oeis.org

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

Author

Benoit Jubin, Sep 25 2009

Keywords

Examples

			4.319176544623038682474939196864331221480491851277749028168381742824336001505931423481751288467967...
		

Crossrefs

Cf. A165733.

A165733 Decimal expansion of sum(e^(1/n!)-1,n>=1).

Original entry on oeis.org

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

Author

Benoit Jubin, Sep 25 2009

Keywords

Examples

			2.600894716163993447114651725511668723723244757577789453201414115100259371152383828910369109942801...
		

Crossrefs

Cf. A165732.

A162414 a(0)=3 and a(n+1)=a(n)^a(n).

Original entry on oeis.org

3, 27, 443426488243037769948249630619149892803
Offset: 0

Author

Benoit Jubin, Jul 02 2009

Keywords

Comments

If the initial term is 2 the sequence is 2, 4, 256, 256^256, ...

Crossrefs

Cf. A000312.

Programs

A168341 The lexicographically earliest injective sequence of nonnegative integers such that a(a(n)) is a square for all n>=0.

Original entry on oeis.org

0, 1, 3, 4, 9, 6, 16, 8, 25, 36, 11, 49, 13, 64, 15, 81, 100, 18, 121, 20, 144, 22, 169, 24, 196, 225, 27, 256, 29, 289, 31, 324, 33, 361, 35, 400, 441, 38, 484, 40, 529, 42, 576, 44, 625, 46, 676, 48, 729, 784, 51, 841, 53, 900, 55, 961, 57, 1024, 59, 1089, 61, 1156, 63
Offset: 0

Author

Eric Angelini and Benoit Jubin, Nov 23 2009

Keywords

Comments

The term a(n) is either n+1 or a square. All the squares appear and they appear in increasing order. Every other term is a square, except when the index is a square, in which case, the corresponding term is also a square (which shifts the pattern). See FORMULA for a more precise statement.

Examples

			For n=6, we have k=floor(sqrt(6))=2; since 6-2=4 is even, a(6)=((6+2)/2)^2=16.
		

Formula

To define a(n), let k = floor(sqrt(n)). Then a(n) = n+1 if n-k^2 is odd and ((n+k)/2)^2 if n-k^2 is even.
Note that k^2 is the largest square which is at most n.

A162185 a(0)=3 and a(n+1)=2^a(n).

Original entry on oeis.org

3, 8, 256, 115792089237316195423570985008687907853269984665640564039457584007913129639936
Offset: 0

Author

Benoit Jubin, Jun 27 2009

Keywords

Comments

Initializing at 0, 1, 2 or 4 gives A014221.

Crossrefs

A000079 (powers of 2).

Programs