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 18 results. Next

A167198 Fractal sequence of the interspersion A083047.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 4, 2, 3, 5, 1, 4, 6, 2, 7, 3, 5, 8, 1, 9, 4, 6, 10, 2, 7, 11, 3, 12, 5, 8, 13, 1, 9, 14, 4, 15, 6, 10, 16, 2, 17, 7, 11, 18, 3, 12, 19, 5, 20, 8, 13, 21, 1, 22, 9, 14, 23, 4, 15, 24, 6, 25, 10, 16, 26, 2, 17, 27, 7, 28, 11, 18, 29, 3, 30, 12, 19, 31, 5, 20, 32, 8, 33, 13
Offset: 1

Views

Author

Clark Kimberling, Oct 30 2009

Keywords

Comments

As a fractal sequence, if the first occurrence of each term is deleted, the remaining sequence is the original. In general, the interspersion of a fractal sequence is constructed by rows: row r consists of all n, such that a(n)=r; in particular, A083047 is constructed in this way from A167198.
a(n-1) gives the row number which contains n in the dual Wythoff array A126714 (beginning the row count at 1), see also A223025 and A019586. - Casey Mongoven, Mar 11 2013

Examples

			To produce row 5, first write row 4: 2,3,1, then place 4 right before 2, and then place 5 right before 1, getting 4,2,3,5,1.
		

References

  • Clark Kimberling, Stolarsky interspersions, Ars Combinatoria 39 (1995), 129-138.

Crossrefs

Formula

Following is a construction that avoids reference to A083047.
Write initial rows:
Row 1: .... 1
Row 2: .... 1
Row 3: .... 2..1
Row 4: .... 2..3..1
For n>=4, to form row n+1, let k be the least positive integer not yet used; write row n, and right before the first number that is also in row n-1, place k; right before the next number that is also in row n-1, place k+1, and continue. A167198 is the concatenation of the rows. (If "before" is replaced by "after", the resulting fractal sequence is A003603, and the associated interspersion is the Wythoff array, A035513.)

A083048 Main diagonal of table A083047.

Original entry on oeis.org

1, 5, 17, 36, 80, 164, 300, 575, 1020, 1884, 3426, 5921, 10568, 18697, 31850, 55716, 94332, 163579, 282388, 474625, 814328, 1363979, 2328358, 3963781, 6609951, 11209355, 18969158, 31524783, 53186480, 88235842, 148471479, 249459365
Offset: 0

Views

Author

Paul D. Hanna, Apr 18 2003

Keywords

Crossrefs

Formula

a(n) = T(n, n), where T(n, 0) = floor(n*x/(x-1)) + 1, T(n, k+1) = ceiling(x*T(n, k)), for n >= 0, k >= 0, with x = (sqrt(5)+1)/2.

A083049 Antidiagonal sums of table A083047.

Original entry on oeis.org

1, 5, 15, 34, 68, 127, 225, 387, 652, 1084, 1787, 2927, 4775, 7769, 12616, 20462, 33160, 53709, 86962, 140769, 227834, 368711, 596658, 965488, 1562270, 2527887, 4090292, 6618319, 10708756, 17327225, 28036136, 45363522, 73399824, 118763517
Offset: 0

Views

Author

Paul D. Hanna, Apr 18 2003

Keywords

Crossrefs

Formula

a(n) = Sum_{k=0..n} T(k, n-k), where T(n, 0) = floor(n*x/(x-1))+1, T(n, k+1) = ceiling(x*T(n, k)), for n>=0, k>=0, with x = (sqrt(5)+1)/2.

A048739 Expansion of 1/((1 - x)*(1 - 2*x - x^2)).

Original entry on oeis.org

1, 3, 8, 20, 49, 119, 288, 696, 1681, 4059, 9800, 23660, 57121, 137903, 332928, 803760, 1940449, 4684659, 11309768, 27304196, 65918161, 159140519, 384199200, 927538920, 2239277041, 5406093003, 13051463048, 31509019100, 76069501249
Offset: 0

Views

Author

Keywords

Comments

Partial sums of Pell numbers A000129.
W(n){1,3;2,-1,1} = Sum_{i=1..n} W(i){1,2;2,-1,0}, where W(n){a,b; p,q,r} implies x(n) = p*x(n-1) - q*x(n-2) + r; x(0)=a, x(1)=b.
Number of 2 X (n+1) binary arrays with path of adjacent 1's from upper left to lower right corner. - R. H. Hardin, Mar 16 2002
Binomial transform of A029744. - Paul Barry, Apr 23 2004
Number of (s(0), s(1), ..., s(n+2)) such that 0 < s(i) < 4 and |s(i) - s(i-1)| <= 1 for i = 1,2,...,n+2, s(0) = 1, s(n+2) = 3. - Herbert Kociemba, Jun 16 2004
Equals row sums of triangle A153346. - Gary W. Adamson, Dec 24 2008
Equals the sum of the terms of the antidiagonals of A142978. - J. M. Bergot, Nov 13 2012
a(p-2) == 0 mod p where p is an odd prime, see A270342. - Altug Alkan, Mar 15 2016
Also, the lexicographically earliest sequence of positive integers such that for n > 3, {sqrt(2)*a(n)} is located strictly between {sqrt(2)*a(n-1)} and {sqrt(2)*a(n-2)} where {} denotes the fractional part. - Ivan Neretin, May 02 2017
a(n+1) is the number of weak orderings on {1,...,n} that are weakly single-peaked w.r.t. the total ordering 1 < ... < n. - J. Devillet, Oct 06 2017

References

  • Allombert, Bill, Nicolas Brisebarre, and Alain Lasjaunias. "On a two-valued sequence and related continued fractions in power series fields." The Ramanujan Journal 45.3 (2018): 859-871. See Theorem 3, d_{4n+3}.

Crossrefs

First row of table A083087.
With a different offset, a(4n)=A008843(n), a(4n-2)=8*A001110(n), a(2n-1)=A001652(n).

Programs

  • Maple
    a:=n->sum(fibonacci(i,2), i=0..n): seq(a(n), n=1..29); # Zerinvary Lajos, Mar 20 2008
  • Mathematica
    Join[{a=1,b=3},Table[c=2*b+a+1;a=b;b=c,{n,60}]] (* Vladimir Joseph Stephan Orlovsky, Feb 01 2011 *)
    CoefficientList[Series[1/(1-3x+x^2+x^3),{x,0,30}],x] (* or *) LinearRecurrence[{3,-1,-1},{1,3,8},30] (* Harvey P. Dale, Jun 13 2011 *)
  • PARI
    a(n)=local(w=quadgen(8));-1/2+(3/4+1/2*w)*(1+w)^n+(3/4-1/2*w)*(1-w)^n
    
  • PARI
    vector(100, n, n--; floor((1+sqrt(2))^(n+2)/4)) \\ Altug Alkan, Oct 07 2015
    
  • PARI
    Vec(1/((1-x)*(1-2*x-x^2)) + O(x^40)) \\ Michel Marcus, May 06 2017

Formula

a(n) = 2*a(n-1) + a(n-2) + 1 with n > 1, a(0)=1, a(1)=3.
a(n) = ((2 + (3*sqrt(2))/2)*(1 + sqrt(2))^n - (2 - (3*sqrt(2))/2)*(1 - sqrt(2))^n )/(2*sqrt(2)) - 1/2.
a(0)=1, a(n+1) = ceiling(x*a(n)) for n > 0, where x = 1+sqrt(2). - Paul D. Hanna, Apr 22 2003
a(n) = 3*a(n-1) - a(n-2) - a(n-3). With two leading zeros, e.g.f. is exp(x)(cosh(sqrt(2)x)-1)/2. a(n) = Sum_{k=0..floor((n+2)/2)} binomial(n+2, 2k+2)2^k. - Paul Barry, Aug 16 2003
-a(-3-n) = A077921(n). - N. J. A. Sloane, Sep 13 2003
E.g.f.: exp(x)(cosh(x/sqrt(2)) + sqrt(2)sinh(x/sqrt(2)))^2. - N. J. A. Sloane, Sep 13 2003
a(n) = floor((1+sqrt(2))^(n+2)/4). - Bruno Berselli, Feb 06 2013
a(n) = (((1-sqrt(2))^(n+2) + (1+sqrt(2))^(n+2) - 2) / 4). - Altug Alkan, Mar 16 2016
2*a(n) = A001333(n+2)-1. - R. J. Mathar, Oct 11 2017
a(n) = Sum_{k=0..n} binomial(n+1,k+1)*2^floor(k/2). - Tony Foster III, Oct 12 2017

Extensions

Corrected and extended by Larry Reeves (larryr(AT)acm.org), Jun 11 2002

A083044 Square table read by antidiagonals forms a permutation of the natural numbers: T(n,0) = floor(n*x/(x-1))+1, T(n,k+1) = ceiling(x*T(n,k)), where x=3/2, n >= 0, k >= 0.

Original entry on oeis.org

1, 2, 4, 3, 6, 7, 5, 9, 11, 10, 8, 14, 17, 15, 13, 12, 21, 26, 23, 20, 16, 18, 32, 39, 35, 30, 24, 19, 27, 48, 59, 53, 45, 36, 29, 22, 41, 72, 89, 80, 68, 54, 44, 33, 25, 62, 108, 134, 120, 102, 81, 66, 50, 38, 28, 93, 162, 201, 180, 153, 122, 99, 75, 57, 42, 31, 140, 243
Offset: 0

Views

Author

Paul D. Hanna, Apr 18 2003

Keywords

Comments

First row is A061419, first column is T(n,0) = A016777(n) = 3n+1 (namely the numbers not of the form ceiling(3*k/2) for any natural number k, in increasing order), main diagonal is A083045, antidiagonal sums give A083046. [further detail on first column added by Glen Whitney, Aug 03 2018]
A083044 is the dispersion of the sequence A007494 of positive integers congruent to (0 or 2) mod 3; see A191655. - Clark Kimberling, Jun 10 2011
If T(n+1,k) - T(n,k) = 2m, then T(n+1,k+1) - T(n,k+1) = ceiling(3T(n+1,k)/2) - ceiling(3T(n,k)/2) = ceiling(3T(n,k)/2 + 3m) - ceiling(3T(n,k)/2) = 3m. Similarly, if T(n+1,k) - T(n,k) = 2m+1, then T(n+1,k+1) - T(n,k+1) = ceiling(3T(n,k)/2 + 3m + 3/2) - ceiling(3T(n,k)/2) = {3m+1 or 3m+2, according to whether T(n,k) is even or odd}. The first differences of the first column T(n,0) are periodic: (3)*. The parities of the first column T(n,0) are periodic: (odd,even)*. Hence by induction using the prior two observations, the first differences and parities of every column will be periodic; e.g., for the second column T(n,2): the first differences are (4,5)* and the parities are (even,even,odd,odd)*; for the third column T(n,3): (6,8,6,7)* and (odd,odd,odd,odd,even,even,even,even)*; for the fourth column T(n,4): (9,12,9,10,9,12,9,11)* and (o,e,e,o,o,e,e,o,e,o,o,e,e,o,o,e)*. Is the period length of the first differences of column k always 2^{k-1}? And is the period length of parities always 2^k? Does every integer > 2 occur as T(n+1,k) - T(n,k) for some n and k? Is the smallest first difference in column k always A061418(k+1)? And is the largest first difference in column k always A061419(k+2)? - Glen Whitney, Aug 03 2018
Consider the following two-player game: Start with two nonempty piles of counters. Players alternate taking turns consisting of first discarding one of the piles and then dividing the remaining pile into two nonempty piles. The smaller pile may always be discarded; the larger pile may only be discarded if the smaller pile is at least half as large. (Either pile may be discarded if they are equal.) The player who cannot move (because the configuration has reached two piles of one counter each) loses. Then the numbers c for which two piles of size c is a losing configuration (for the player whose turn it is) are exactly T(4,k) for k > 1, together with 1,3,5, and 9. - Glen Whitney, Aug 03 2018

Examples

			Table begins:
   1  2  3   5   8  12  18  27  41   62   93  140 ...
   4  6  9  14  21  32  48  72 108  162  243  365 ...
   7 11 17  26  39  59  89 134 201  302  453  680 ...
  10 15 23  35  53  80 120 180 270  405  608  912 ...
  13 20 30  45  68 102 153 230 345  518  777 1166 ...
  16 24 36  54  81 122 183 275 413  620  930 1395 ...
  19 29 44  66  99 149 224 336 504  756 1134 1701 ...
  22 33 50  75 113 170 255 383 575  863 1295 1943 ...
  25 38 57  86 129 194 291 437 656  984 1476 2214 ...
  28 42 63  95 143 215 323 485 728 1092 1638 2457 ...
  31 47 71 107 161 242 363 545 818 1227 1841 2762 ...
		

Crossrefs

Row in which a number occurs: A163491.
Column in which a number occurs: A087088.

Formula

T(A163491(n)-1, A087088(n)-1) = n. - Peter Munn, Jul 16 2020 [corrected Peter Munn, Aug 02 2020]

A083050 Square table read by antidiagonals forms a permutation of the natural numbers: T(n,0) = floor(n*x/(x-1))+1, T(n,k+1) = ceiling(x*T(n,k)), where x = sqrt(2), n>=0, k>=0.

Original entry on oeis.org

1, 2, 4, 3, 6, 7, 5, 9, 10, 11, 8, 13, 15, 16, 14, 12, 19, 22, 23, 20, 18, 17, 27, 32, 33, 29, 26, 21, 25, 39, 46, 47, 42, 37, 30, 24, 36, 56, 66, 67, 60, 53, 43, 34, 28, 51, 80, 94, 95, 85, 75, 61, 49, 40, 31, 73, 114, 133, 135, 121, 107, 87, 70, 57, 44, 35, 104, 162, 189
Offset: 0

Views

Author

Paul D. Hanna, Apr 18 2003

Keywords

Comments

First column is A083051, main diagonal is A083052, antidiagonal sums give A083053.
A083050 is the dispersion of the sequence given by floor(1+n*sqrt(2)); for a discussion of dispersions, see A191429.

Examples

			Table begins:
   1  2  3   5   8  12  17  25  36  51   73  104 ...
   4  6  9  13  19  27  39  56  80 114  162  230 ...
   7 10 15  22  32  46  66  94 133 189  268  380 ...
  11 16 23  33  47  67  95 135 191 271  384  544 ...
  14 20 29  42  60  85 121 172 244 346  490  693 ...
  18 26 37  53  75 107 152 215 305 432  611  865 ...
  21 30 43  61  87 124 176 249 353 500  708 1002 ...
  24 34 49  70  99 141 200 283 401 568  804 1138 ...
  28 40 57  81 115 163 231 327 463 655  927 1311 ...
  31 44 63  90 128 182 258 365 517 732 1036 1466 ...
  35 50 71 101 143 203 288 408 577 817 1156 1635 ...
		

Crossrefs

Programs

  • Mathematica
    (* program generates the dispersion array T of the complement of increasing sequence f[n] *)
    r = 40; r1 = 12; (* r=# rows of T, r1=# rows to show *)
    c = 40; c1 = 12; (* c=# cols of T, c1=# cols to show *)
    x = Sqrt[2]; f[n_] := Floor[n*x + 1]
    (* f(n) is complement of column 1 *)
    mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
    rows = {NestList[f, 1, c]};
    Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
    t[i_, j_] := rows[[i, j]];
    TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]]
    (* A083050 as an array *)
    Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* array as a sequence *)
    (* Program by Peter J. C. Moses, Jun 01 2011, added here Jun 03 2011 by Clark Kimberling *)

A083087 Square table read by antidiagonals which forms a permutation of the natural numbers: T(n,0) = floor(n*x/(x-1))+1, T(n,k+1) = ceiling(x*T(n,k)), for n>=0, k>=0, where x = 1 + sqrt(2).

Original entry on oeis.org

1, 3, 2, 8, 5, 4, 20, 13, 10, 6, 49, 32, 25, 15, 7, 119, 78, 61, 37, 17, 9, 288, 189, 148, 90, 42, 22, 11, 696, 457, 358, 218, 102, 54, 27, 12, 1681, 1104, 865, 527, 247, 131, 66, 29, 14, 4059, 2666, 2089, 1273, 597, 317, 160, 71, 34, 16, 9800, 6437, 5044, 3074
Offset: 0

Views

Author

Paul D. Hanna, Apr 21 2003

Keywords

Comments

The array in A083087 is the dispersion of the sequence given floor(n+1+n*sqrt(2)). The Mathematica program at A191438 generates A083087 using f[n_]:=Floor[n*x+n+1] instead of f[n_]:=Floor[n*x+n]. - Clark Kimberling, Jun 04 2011

Examples

			Table begins:
   1  3   8  20  49  119  288 ...
   2  5  13  32  78  189  457 ...
   4 10  25  61 148  358  865 ...
   6 15  37  90 218  527 1273 ...
   7 17  42 102 247  597 1442 ...
   9 22  54 131 317  766 1850 ...
  11 27  66 160 387  935 2258 ...
  12 29  71 172 416 1005 2427 ...
  14 34  83 201 486 1174 2835 ...
  16 39  95 230 556 1343 3243 ...
  18 44 107 259 626 1512 3651 ...
  19 46 112 271 655 1582 3820 ...
  21 51 124 300 725 1751 4228 ...
  23 56 136 329 795 1920 4636 ...
  24 58 141 341 824 1990 4805 ...
		

Crossrefs

Cf. A083088 (first column), A048739 (first row), A083090 (diagonal), A083091 (antidiagonal sums), A083044, A083047, A083050.

Programs

  • Magma
    z:=10; x:=1+Sqrt(2); S:=[]; for n in [0..z] do for k in [0..n] do if n-k eq 0 then Append(~S, Floor(n*x/(x-1))+1); else Append(~S, Ceiling(x*S[k+1+(n*(n-1) div 2)])); end if; end for; end for; S; // Klaus Brockhaus, Jan 04 2011
  • Mathematica
    (See Comments.)

Formula

T(n,k+1) = 2*T(n,k) + T(n,k-1) + 1 for n>=0, k>=1.

A083051 First column of table A083050.

Original entry on oeis.org

1, 4, 7, 11, 14, 18, 21, 24, 28, 31, 35, 38, 41, 45, 48, 52, 55, 59, 62, 65, 69, 72, 76, 79, 82, 86, 89, 93, 96, 100, 103, 106, 110, 113, 117, 120, 123, 127, 130, 134, 137, 140, 144, 147, 151, 154, 158, 161, 164, 168, 171, 175, 178, 181, 185, 188, 192, 195, 199, 202
Offset: 0

Views

Author

Paul D. Hanna, Apr 18 2003

Keywords

Crossrefs

Formula

a(n)=floor(n*x/(x-1))+1, x=sqrt(2).

A343152 Reverse the order of all but the most significant bits in the maximal Fibonacci expansion of n.

Original entry on oeis.org

1, 2, 3, 4, 6, 5, 7, 8, 11, 10, 9, 12, 16, 14, 19, 13, 18, 17, 15, 20, 21, 29, 27, 24, 32, 26, 23, 31, 22, 30, 28, 25, 33, 42, 37, 50, 35, 48, 45, 40, 53, 34, 47, 44, 39, 52, 43, 38, 51, 36, 49, 46, 41, 54, 55, 76, 71, 63, 84, 69, 61, 82, 58, 79, 74, 66, 87
Offset: 1

Views

Author

J. Parker Shectman, Apr 07 2021

Keywords

Comments

A self-inverse permutation of the natural numbers.
Analogous to A059893 with binary expansion replaced by maximal Fibonacci expansion.
Analogous to A343150 with minimal Fibonacci expansion replaced by maximal Fibonacci expansion.
For n=1, the expansion equals 1. For n>=2, the expansion equals A104326(n-1) with a 1 appended. The 1 corresponds to a digit (always equal to 1) for F(1)=1, in addition to the digit for F(2)=1. (This expansion is NOT a representation, see reference in link, pp. 106 and 137.)
Write the sequence as a (right-justified) "tetrangle" or "irregular triangle" tableau with F(t) (Fibonacci number) entries on each row, for t=1,2,3,.... Then, columns of the tableau equal rows of the array A083047 (see reference in link, p. 131):
1
2
3, 4
6, 5, 7
8, 11, 10, 9, 12
16, 14, 19, 13, 18, 17, 15, 20
...

Examples

			For an example of calculation by reversing Fibonacci binary digits, see reference in link, p. 144:
On the basis (1,1,2,3,5,8) n=13 is written as 110101, Reversing all but the most AND least significant digits gives 101011, which evaluates to 16, so a(13)=16.
On the basis (1,1,2,3,5,8) n=14 is written as 101101, Reversing all but the most AND least significant digits gives 101101, which evaluates to 14, so a(14)=14.
		

Crossrefs

Programs

  • Mathematica
    (*Produce indices of maximal Fibonacci expansion (recursively)*)
    MaxFibInd[n_] := Module[{t = Floor[Log[GoldenRatio, Sqrt[5]*n + 1]] - 1}, Piecewise[{{{1}, n == 1}, {Append[MaxFibInd[n - Fibonacci[t]], t], n > 1}},]];
    (*Define a(n)*)
    a[n_] := Module[{MFI = MaxFibInd[n]}, Apply[Plus, Fibonacci[Last[MFI] - MFI + 1]]];
    (*Generate DATA*)
    Array[a, 67]

A083088 First column of table A083087.

Original entry on oeis.org

1, 2, 4, 6, 7, 9, 11, 12, 14, 16, 18, 19, 21, 23, 24, 26, 28, 30, 31, 33, 35, 36, 38, 40, 41, 43, 45, 47, 48, 50, 52, 53, 55, 57, 59, 60, 62, 64, 65, 67, 69, 70, 72, 74, 76, 77, 79, 81, 82, 84, 86, 88, 89, 91, 93, 94, 96, 98, 100, 101, 103, 105, 106, 108, 110, 111, 113, 115
Offset: 0

Views

Author

Paul D. Hanna, Apr 21 2003

Keywords

Comments

It appears that A188937 gives the positions of 0 in the zero-one sequence A188037; complement of A080754. - Clark Kimberling, Mar 19 2011
Is this (apart from the prefixed a(0)) the same as A080755? - R. J. Mathar, Jul 31 2025

Crossrefs

Programs

  • Magma
    z:=70; x:=1+Sqrt(2); [ Floor(n*x/(x-1))+1: n in [0..z] ]; // Klaus Brockhaus, Jan 04 2011
  • Mathematica
    f[n_] := Floor[n/Sqrt@2 + n + 1]; Array[f, 68, 0]

Formula

a(n) = floor(n*x/(x-1)) + 1, n>=0, where x=1+sqrt(2).
a(n) = floor(n/sqrt(2)) + n + 1 = 1+n+A049472(n).

Extensions

This entry formerly contained an erroneous comment, which was deleted by N. J. A. Sloane, Jan 30 2008
Showing 1-10 of 18 results. Next