0, 1, 2, 6, 15, 40, 104, 273, 714, 1870, 4895, 12816, 33552, 87841, 229970, 602070, 1576239, 4126648, 10803704, 28284465, 74049690, 193864606, 507544127, 1328767776, 3478759200, 9107509825, 23843770274, 62423800998, 163427632719, 427859097160, 1120149658760
Offset: 0
A192744
Constant term of the reduction by x^2->x+1 of the polynomial p(n,x) defined below in Comments.
Original entry on oeis.org
1, 1, 3, 8, 29, 133, 762, 5215, 41257, 369032, 3676209, 40333241, 483094250, 6271446691, 87705811341, 1314473334832, 21017294666173, 357096406209005, 6424799978507178, 122024623087820183, 2439706330834135361, 51219771117454755544
Offset: 0
The first five polynomials and their reductions:
1 -> 1
1+x -> 1+x
2+x+x^2 -> 3+2x
6+2x+x^2+x^3 -> 8+5x
24+6x+2x^2+x^3+x^4 -> 29+13x, so that
A192744=(1,1,3,8,29,...) and A192745=(0,1,2,5,13,...).
-
A192744p := proc(n,x)
option remember;
if n = 0 then
1;
else
x*procname(n-1,x)+n! ;
expand(%) ;
end if;
end proc:
A192744 := proc(n)
local p;
p := A192744p(n,x) ;
while degree(p,x) > 1 do
p := algsubs(x^2=x+1,p) ;
p := expand(p) ;
end do:
coeftayl(p,x=0,0) ;
end proc: # R. J. Mathar, Dec 16 2015
-
q = x^2; s = x + 1; z = 40;
p[0, n_] := 1; p[n_, x_] := x*p[n - 1, x] + n!;
Table[Expand[p[n, x]], {n, 0, 7}]
reduce[{p1_, q_, s_, x_}] :=
FixedPoint[(s PolynomialQuotient @@ #1 +
PolynomialRemainder @@ #1 &)[{#1, q, x}] &, p1]
t = Table[reduce[{p[n, x], q, s, x}], {n, 0, z}];
u1 = Table[Coefficient[Part[t, n], x, 0], {n, 1, z}]
(* A192744 *)
u2 = Table[Coefficient[Part[t, n], x, 1], {n, 1, z}]
(* A192745 *)
A035317
Pascal-like triangle associated with A000670.
Original entry on oeis.org
1, 1, 1, 1, 2, 2, 1, 3, 4, 2, 1, 4, 7, 6, 3, 1, 5, 11, 13, 9, 3, 1, 6, 16, 24, 22, 12, 4, 1, 7, 22, 40, 46, 34, 16, 4, 1, 8, 29, 62, 86, 80, 50, 20, 5, 1, 9, 37, 91, 148, 166, 130, 70, 25, 5, 1, 10, 46, 128, 239, 314, 296, 200, 95, 30, 6, 1, 11, 56, 174, 367, 553, 610, 496, 295, 125
Offset: 0
Triangle begins:
1;
1, 1;
1, 2, 2;
1, 3, 4, 2;
1, 4, 7, 6, 3;
1, 5, 11, 13, 9, 3;
1, 6, 16, 24, 22, 12, 4;
1, 7, 22, 40, 46, 34, 16, 4;
1, 8, 29, 62, 86, 80, 50, 20, 5;
1, 9, 37, 91, 148, 166, 130, 70, 25, 5;
1, 10, 46, 128, 239, 314, 296, 200, 95, 30, 6;
...
- Vincenzo Librandi, Rows n = 0..100, flattened
- Joseph Briggs, Alex Parker, Coy Schwieder, and Chris Wells, Frogs, hats and common subsequences, arXiv preprint arXiv:2404.07285 [math.CO], 2024. See p. 28.
- A. Hlavác, M. Marvan, Nonlocal conservation laws of the constant astigmatism equation, arXiv preprint arXiv:1602.06861 [nlin.SI], 2016.
- E. Mendelson, Races with Ties, Math. Mag. 55 (1982), 170-175.
- Index entries for triangles and arrays related to Pascal's triangle
Triangle sums (see the comments):
A000975 (Row1),
A059841 (Row2),
A080239 (Kn11),
A052952 (Kn21),
A129696 (Kn22),
A001906 (Kn3),
A001654 (Kn4),
A001045 (Fi1, Fi2),
A023435 (Ca2), Gi2 (
A193146),
A190525 (Ze2),
A193147 (Ze3),
A181532 (Ze4). -
Johannes W. Meijer, Jul 20 2011
-
a035317 n k = a035317_tabl !! n !! k
a035317_row n = a035317_tabl !! n
a035317_tabl = map snd $ iterate f (0, [1]) where
f (i, row) = (1 - i, zipWith (+) ([0] ++ row) (row ++ [i]))
-- Reinhard Zumkeller, Jul 09 2012
-
A035317 := proc(n,k): add((-1)^(i+k) * binomial(i+n-k+1, i), i=0..k) end: seq(seq(A035317(n,k), k=0..n), n=0..10); # Johannes W. Meijer, Jul 20 2011
A035317 := proc(n,k): coeff(coeftayl(1/((y*x-1)*(y*x+1)*((y+1)*x-1)), x=0, n), y, k) end: seq(seq(A035317(n,k), k=0..n), n=0..10); # Johannes W. Meijer, Jul 20 2011
-
t[n_, k_] := (-1)^k*(((-1)^k*(n+2)!*Hypergeometric2F1[1, n+3, k+2, -1])/((k+1)!*(n-k+1)!) + 2^(k-n-2)); Flatten[ Table[ t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, Dec 14 2011, after Johannes W. Meijer *)
-
{T(n,k)=if(n==k,(n+2)\2,if(k==0,1,if(n>k,T(n-1,k-1)+T(n-1,k))))}
for(n=0,12,for(k=0,n,print1(T(n,k),","));print("")) \\ Paul D. Hanna, Jul 18 2012
-
def A035317_row(n):
@cached_function
def prec(n, k):
if k==n: return 1
if k==0: return 0
return -prec(n-1,k-1)-sum(prec(n,k+i-1) for i in (2..n-k+1))
return [(-1)^k*prec(n+2, k) for k in (1..n)]
for n in (1..11): print(A035317_row(n)) # Peter Luschny, Mar 16 2016
A104763
Triangle read by rows: Fibonacci(1), Fibonacci(2), ..., Fibonacci(n) in row n.
Original entry on oeis.org
1, 1, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 2, 3, 5, 1, 1, 2, 3, 5, 8, 1, 1, 2, 3, 5, 8, 13, 1, 1, 2, 3, 5, 8, 13, 21, 1, 1, 2, 3, 5, 8, 13, 21, 34, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233
Offset: 1
First few rows of the triangle are:
1;
1, 1;
1, 1, 2;
1, 1, 2, 3;
1, 1, 2, 3, 5;
1, 1, 2, 3, 5, 8;
1, 1, 2, 3, 5, 8, 13; ...
-
Flat(List([1..15], n -> List([1..n], k -> Fibonacci(k)))); # G. C. Greubel, Jul 13 2019
-
a104763 n k = a104763_tabl !! (n-1) !! (k-1)
a104763_row n = a104763_tabl !! (n-1)
a104763_tabl = map (flip take $ tail a000045_list) [1..]
-- Reinhard Zumkeller, Aug 15 2013
-
[Fibonacci(k): k in [1..n], n in [1..15]]; // G. C. Greubel, Jul 13 2019
-
Table[Fibonacci[k], {n,15}, {k,n}]//Flatten (* G. C. Greubel, Jul 13 2019 *)
-
for(n=1,15, for(k=1,n, print1(fibonacci(k), ", "))) \\ G. C. Greubel, Jul 13 2019
-
[[fibonacci(k) for k in (1..n)] for n in (1..15)] # G. C. Greubel, Jul 13 2019
A124502
a(1)=a(2)=1; thereafter, a(n+1) = a(n) + a(n-1) + 1 if n is a multiple of 5, otherwise a(n+1) = a(n) + a(n-1).
Original entry on oeis.org
1, 1, 2, 3, 5, 9, 14, 23, 37, 60, 98, 158, 256, 414, 670, 1085, 1755, 2840, 4595, 7435, 12031, 19466, 31497, 50963, 82460, 133424, 215884, 349308, 565192, 914500, 1479693, 2394193, 3873886, 6268079, 10141965, 16410045, 26552010, 42962055, 69514065, 112476120
Offset: 1
a(6) = a(5) + a(4) + 1 = 5 + 3 + 1 = 9 because n=5 is a multiple of 5.
a(7) = a(6) + a(5) = 9 + 5 = 14 because n=6 is not a multiple of 5.
-
A124502:=proc(n) option remember; local t1; if n <= 2 then return 1; fi: if n mod 5 = 1 then t1:=1 else t1:=0; fi: procname(n-1)+procname(n-2)+t1; end proc; [seq(A124502(n), n=1..100)]; # N. J. A. Sloane, May 25 2008
-
a=0; b=0; lst={a,b}; Do[z=a+b+1; AppendTo[lst,z]; a=b; b=z; z=a+b; AppendTo[lst,z]; a=b; b=z; z=a+b; AppendTo[lst,z]; a=b; b=z; z=a+b; AppendTo[lst,z]; a=b; b=z; z=a+b; AppendTo[lst,z]; a=b; b=z,{n,4!}]; lst (* Vladimir Joseph Stephan Orlovsky, Feb 16 2010 *)
nxt[{n_,a_,b_}]:={n+1,b,If[Divisible[n,5],a+b+1,a+b]}; NestList[nxt,{2,1,1},40][[All,2]] (* or *) LinearRecurrence[{1,1,0,0,1,-1,-1},{1,1,2,3,5,9,14},40] (* Harvey P. Dale, Jun 15 2017 *)
A131132
a(n) = a(n-1) + a(n-2) + 1 if n is a multiple of 6, otherwise a(n) = a(n-1) + a(n-2).
Original entry on oeis.org
1, 1, 2, 3, 5, 8, 14, 22, 36, 58, 94, 152, 247, 399, 646, 1045, 1691, 2736, 4428, 7164, 11592, 18756, 30348, 49104, 79453, 128557, 208010, 336567, 544577, 881144, 1425722, 2306866, 3732588, 6039454, 9772042, 15811496, 25583539, 41395035, 66978574, 108373609
Offset: 0
Since 5 is not a multiple of 6, a(5) = a(4) + a(3) = 5 + 3 = 8. Since 6 is a multiple of 6, a(6) = a(5) + a(4) + 1 = 8 + 5 + 1 = 14. - _Michael B. Porter_, Jul 26 2016
- H. Matsui et al., Problem B-1035, Fibonacci Quarterly, Vol. 45, Number 2; 2007; p. 182.
- Index entries for linear recurrences with constant coefficients, signature (1,1,0,0,0,1,-1,-1).
-
A131132:=proc(n) option remember; local t1; if n <= 2 then RETURN(1); fi: if n mod 6 = 1 then t1:=1 else t1:=0; fi: procname(n-1)+procname(n-2)+t1; end; [seq(A131132(n), n=1..100)]; # N. J. A. Sloane, May 25 2008; Typo corrected by R. J. Mathar, May 30 2008
-
Print[Table[Round[(1 + Sqrt[5])/8 Fibonacci[n + 3]], {n, 0, 50}]] ;
Print[RecurrenceTable[{c[n] == c[n - 1] + c[n - 2] + c[n - 6] - c[n - 7] - c[n - 8], c[0] == 1, c[1] == 1, c[2] == 2, c[3] == 3, c[4] == 5, c[5] == 8, c[6] == 14, c[7] == 22}, c, {n, 0, 50}]] ; (* John M. Campbell, Jul 07 2016 *)
LinearRecurrence[{1, 1, 0, 0, 0, 1, -1, -1}, {1, 1, 2, 3, 5, 8, 14, 22}, 40] (* Vincenzo Librandi, Jul 07 2016 *)
A192758
Coefficient of x in the reduction by x^2->x+1 of the polynomial p(n,x) defined below in Comments.
Original entry on oeis.org
0, 1, 2, 4, 7, 13, 22, 37, 61, 101, 165, 269, 437, 710, 1151, 1865, 3020, 4890, 7915, 12810, 20730, 33546, 54282, 87834, 142122, 229963, 372092, 602062, 974161, 1576231, 2550400, 4126639, 6677047, 10803695, 17480751, 28284455, 45765215
Offset: 1
-
q = x^2; s = x + 1; z = 40;
p[0, n_] := 1; p[n_, x_] := x*p[n - 1, x] + Floor[(n + 4)/4] /; n > 0;
Table[Expand[p[n, x]], {n, 0, 7}]
reduce[{p1_, q_, s_, x_}] :=
FixedPoint[(s PolynomialQuotient @@ #1 +
PolynomialRemainder @@ #1 &)[{#1, q, x}] &, p1]
t = Table[reduce[{p[n, x], q, s, x}], {n, 0, z}];
u1 = Table[Coefficient[Part[t, n], x, 0], {n, 1, z}]
(* A080239 *)
u2 = Table[Coefficient[Part[t, n], x, 1], {n, 1, z}]
(* A192758 *)
A230448
T(n, k) = T(n-1, k-1) + T(n-1, k) with T(n, 0) = 1 and T(n, n) = A226205(n+1), n >= 0 and 0 <= k <= n.
Original entry on oeis.org
1, 1, 0, 1, 1, 3, 1, 2, 4, 5, 1, 3, 6, 9, 16, 1, 4, 9, 15, 25, 39, 1, 5, 13, 24, 40, 64, 105, 1, 6, 18, 37, 64, 104, 169, 272, 1, 7, 24, 55, 101, 168, 273, 441, 715, 1, 8, 31, 79, 156, 269, 441, 714, 1156, 1869, 1, 9, 39, 110, 235, 425, 710, 1155, 1870, 3025, 4896
Offset: 0
The first few rows of triangle T(n, k), n >= 0 and 0 <= k <= n.
n/k 0 1 2 3 4 5 6 7
------------------------------------------------
0| 1
1| 1, 0
2| 1, 1, 3
3| 1, 2, 4, 5
4| 1, 3, 6, 9, 16
5| 1, 4, 9, 15, 25, 39
6| 1, 5, 13, 24, 40, 64, 105
7| 1, 6, 18, 37, 64, 104, 169, 272
The triangle as a square array Tsq(n, k) = T(n+k, k), n >= 0 and k >= 0.
n/k 0 1 2 3 4 5 6 7
------------------------------------------------
0| 1, 0, 3, 5, 16, 39, 105, 272
1| 1, 1, 4, 9, 25, 64, 169, 441
2| 1, 2, 6, 15, 40, 104, 273, 714
3| 1, 3, 9, 24, 64, 168, 441, 1155
4| 1, 4, 13, 37, 101, 269, 710, 1865
5| 1, 5, 18, 55, 156, 425, 1135, 3000
6| 1, 6, 24, 79, 235, 660, 1795, 4795
7| 1, 7, 31, 110, 345, 1005, 2800, 7595
-
T := proc(n, k) option remember: if k=0 then return(1) elif k=n then return(combinat[fibonacci](n+2)*combinat[fibonacci](n-1)) else procname(n-1, k-1) + procname(n-1, k) fi: end: seq(seq(T(n, k), k=0..n), n=0..10); # End first program.
T := proc(n, k): add(A035317(n+k-p-2, p), p=0..k) end: A035317 := proc(n, k): add((-1)^(i+k) * binomial(i+n-k+1, i), i=0..k) end: seq(seq(T(n, k), k=0..n), n=0..10); # End second program.
A260710
Expansion of 1/(1 - x - x^2 - x^4 + x^5 + x^7).
Original entry on oeis.org
1, 1, 2, 3, 6, 9, 16, 25, 43, 69, 116, 188, 313, 511, 846, 1386, 2288, 3756, 6191, 10174, 16756, 27552, 45357, 74604, 122787, 201996, 332414, 546901, 899946, 1480699, 2436459, 4008858, 6596366, 10853563, 17858788, 29384804, 48350401, 79555943, 130902711
Offset: 0
There are 25 partially ordered partitions of 7, i.e., a(7) = 25. These are (43=34),(421=412),(124=214),(241),(142),(4111),(1411),(1141),(1114),(331),(313),(133),(1132=1123),(2131=1231),(1312=1321),(2311=3211),(31111),(13111),(11311),(11131),(11113),(2221=four),(22111=ten),(211111=six),(1111111).
-
I:=[1,1,2,3,6,9,16]; [n le 7 select I[n] else Self(n-1)+Self(n-2)+Self(n-4)-Self(n-5)-Self(n-7): n in [1..40]]; // Vincenzo Librandi, Aug 04 2015
-
LinearRecurrence[{1, 1, 0, 1, -1, 0, -1}, {1, 1, 2, 3, 6, 9, 16}, 50] (* Vincenzo Librandi, Aug 04 2015 *)
-
Vec(1/(1 - x - x^2 - x^4 + x^5 + x^7) + O(x^50)) \\ Michel Marcus, Aug 06 2015
Comments