A103377 -id:A103377 - cp's OEIS frontend

A103387 Primes in A103377.

2, 3, 5, 7, 17, 31, 71, 127, 157, 227, 257, 293, 349, 419, 503, 16007, 32783, 33329, 80429, 220919, 786433, 878831, 85790399, 10666238497, 394760836573, 22265159019759079, 617347651958903360669, 98459899570803393815965301

Offset: 1

Jonathan Vos Post, Feb 15 2005

Cf. A103377, A103397.

SemiprimeQ[n_]:=Plus@@FactorInteger[n][[All, 2]]?2; Clear[a]; k=9; Do[a[n]=1, {n, k+1}]; a[n_]:=a[n]=a[n-k]+a[n-k-1]; A103377=Array[a, 100] A103387=Union[Select[Array[a, 1000], PrimeQ]] A103397=Union[Select[Array[a, 300], SemiprimeQ]] N[Solve[x^10 - x - 1 == 0, x], 111][[2]]

Intersection of A103377 with A000040.

A103397 Semiprimes in A103377.

Original entry on oeis.org

4, 9, 15, 21, 33, 38, 58, 65, 86, 106, 121, 129, 265, 511, 8114, 8193, 16307, 16853, 17855, 19857, 31298, 68037, 104739, 124205, 131209, 134149, 140457, 152849, 252914, 259918, 265358, 274606, 417527, 2498871, 5291863, 8424051, 8743821

Offset: 1

Jonathan Vos Post, Feb 15 2005

			2071468241 is an element of A103377 and 2071468241= 17 * 121851073 which shows that it is a semiprime.

Cf. A001358, A000931, A079398, A103372-103381, A103377, A103397.

SemiprimeQ[n_]:=Plus@@FactorInteger[n][[All, 2]]?2; Clear[a]; k=9; Do[a[n]=1, {n, k+1}]; a[n_]:=a[n]=a[n-k]+a[n-k-1]; A103377=Array[a, 100] A103387=Union[Select[Array[a, 1000], PrimeQ]] A103397=Union[Select[Array[a, 300], SemiprimeQ]] N[Solve[x^10 - x - 1 == 0, x], 111][[2]]

Intersection of A103377 with A001358.

A103373 a(1) = a(2) = a(3) = a(4) = a(5) = a(6) = 1 and for n>6: a(n) = a(n-5) + a(n-6).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 4, 4, 4, 4, 5, 7, 8, 8, 8, 9, 12, 15, 16, 16, 17, 21, 27, 31, 32, 33, 38, 48, 58, 63, 65, 71, 86, 106, 121, 128, 136, 157, 192, 227, 249, 264, 293, 349, 419, 476, 513, 557, 642, 768, 895, 989, 1070, 1199, 1410, 1663, 1884, 2059, 2269

Offset: 1

Jonathan Vos Post, Feb 03 2005

k=5 case of the family of sequences whose k=1 case is the Fibonacci sequence A000045, k=2 case is the Padovan sequence A000931 (offset so as to begin 1,1,1), k=3 case is A079398 (offset so as to begin 1,1,1,1) and k=4 case is A103372.
The general case for integer k>1 is defined: a(1) = a(2) = ... = a(k+1) and for n>(k+1) a(n) = a(n-k) + a(n-[k+1]).
For this k=5 case, the ratio of successive terms a(n)/a(n-1) approaches the unique positive root of the characteristic polynomial: x^6 - x - 1 = 0. This is the real constant 1.1347241384015194926054460545064728402796672263828014859251495516682....
The sequence of prime values in this k=5 case is A103383; the sequence of semiprime values in this k=5 case is A103393.

			a(22) = 9 because a(22) = a(22-5) + a(22-6) = a(17) + a(16) = 5 + 4 = 9.

Zanten, A. J. van, "The golden ratio in the arts of painting, building and mathematics", Nieuw Archief voor Wiskunde, 4 (17) (1999) 229-245.

G. C. Greubel, Table of n, a(n) for n = 1..1000
J.-P. Allouche and T. Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [The hal link does not always work. - _N. J. A. Sloane_, Feb 19 2025]
J.-P. Allouche and T. Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [Local copy with annotations and a correction from _N. J. A. Sloane_, Feb 19 2025]
Richard Padovan, Dom Hans van der Laan and the Plastic Number.
E. S. Selmer, On the irreducibility of certain trinomials, Math. Scand., 4 (1956) 287-302.
J. Shallit, A generalization of automatic sequences, Theoretical Computer Science, 61 (1988) 1-16.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,1,1).

Cf. A000045, A000931, A079398, A103383, A103393.

Cf. A103372, A103374, A103375, A103376, A103377, A103378, A103379, A103380

k = 5; Do[a[n] = 1, {n, k + 1}]; a[n_] := a[n] = a[n - k] + a[n - k - 1]; Array[a, 65]
RecurrenceTable[{a[n] == a[n - 5] + a[n - 6], a[1] == a[2] == a[3] == a[4] == a[5] == a[6] == 1}, a, {n, 65}] (* or *)
Rest@ CoefficientList[Series[-x (1 + x + x^2 + x^3 + x^4)/(-1 + x^5 + x^6), {x, 0, 65}], x] (* Michael De Vlieger, Oct 03 2016 *)
LinearRecurrence[{0,0,0,0,1,1},{1,1,1,1,1,1},70] (* Harvey P. Dale, Jul 20 2019 *)

a(n)=([0,1,0,0,0,0; 0,0,1,0,0,0; 0,0,0,1,0,0; 0,0,0,0,1,0; 0,0,0,0,0,1; 1,1,0,0,0,0]^(n-1)*[1;1;1;1;1;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

x='x+O('x^50); Vec(x*(1+x+x^2+x^3+x^4)/(1-x^5-x^6 )) \\ G. C. Greubel, May 01 2017

G.f.: x*(1+x+x^2+x^3+x^4) / (1-x^5-x^6 ). - R. J. Mathar, Aug 26 2011

Edited by Ray Chandler and Robert G. Wilson v, Feb 06 2005

A103374 a(1) = a(2) = a(3) = a(4) = a(5) = a(6) = a(7) = 1 and for n>7: a(n) = a(n-6) + a(n-7).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 4, 4, 4, 4, 4, 5, 7, 8, 8, 8, 8, 9, 12, 15, 16, 16, 16, 17, 21, 27, 31, 32, 32, 33, 38, 48, 58, 63, 64, 65, 71, 86, 106, 121, 127, 129, 136, 157, 192, 227, 248, 256, 265, 293, 349, 419, 475, 504, 521, 558, 642, 768, 894, 979, 1025, 1079

Offset: 1

Jonathan Vos Post, Feb 03 2005

k=6 case of the family of sequences whose k=1 case is the Fibonacci sequence A000045, k=2 case is the Padovan sequence A000931 (offset so as to begin 1,1,1), k=3 case is A079398 (offset so as to begin 1,1,1,1), k=4 case is A103372 and k=5 case is A103373.
The general case for integer k>1 is defined: a(1) = a(2) = ... = a(k+1) and for n>(k+1) a(n) = a(n-k) + a(n-[k+1]).
For this k=6 case, the ratio of successive terms a(n)/a(n-1) approaches the unique positive root of the characteristic polynomial: x^7 - x - 1 = 0. This is the real constant 1.1127756842787... (see A230160).
The sequence of prime values in this k=6 case is A103384; the sequence of semiprime values in this k=6 case is A103394.

			a(32) = 17 because a(32) = a(32-6) + a(32-7) = a(26) + a(25) = 9 + 8 = 17.

Zanten, A. J. van, "The golden ratio in the arts of painting, building and mathematics", Nieuw Archief voor Wiskunde, 4 (17) (1999) 229-245.

G. C. Greubel, Table of n, a(n) for n = 1..1000
J.-P. Allouche and T. Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [The hal link does not always work. - _N. J. A. Sloane_, Feb 19 2025]
J.-P. Allouche and T. Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [Local copy with annotations and a correction from _N. J. A. Sloane_, Feb 19 2025]
Richard Padovan, Dom Hans van der Laan and the Plastic Number.
E. S. Selmer, On the irreducibility of certain trinomials, Math. Scand., 4 (1956) 287-302.
J. Shallit, A generalization of automatic sequences, Theoretical Computer Science, 61 (1988) 1-16.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,1,1).

Cf. A000045, A000931, A079398, A103384, A103394, A230160.

Cf. A103372, A103373, A103375, A103376, A103377, A103378, A103379, A103380

k = 6; Do[a[n] = 1, {n, k + 1}]; a[n_] := a[n] = a[n - k] + a[n - k - 1]; Array[a, 70]
RecurrenceTable[{a[n] == a[n - 6] + a[n - 7], a[1] == a[2] == a[3] == a[4] == a[5] == a[6] == a[7] == 1}, a, {n, 70}] (* or *)
Rest@ CoefficientList[Series[-x (1 + x) (1 + x + x^2) (x^2 - x + 1)/(-1 + x^6 + x^7), {x, 0, 70}], x] (* Michael De Vlieger, Oct 03 2016 *)
LinearRecurrence[{0,0,0,0,0,1,1},{1,1,1,1,1,1,1},80] (* Harvey P. Dale, Sep 02 2024 *)

a(n)=([0,1,0,0,0,0,0; 0,0,1,0,0,0,0; 0,0,0,1,0,0,0; 0,0,0,0,1,0,0; 0,0,0,0,0,1,0; 0,0,0,0,0,0,1; 1,1,0,0,0,0,0]^(n-1)*[1;1;1;1;1;1;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

x='x+O('x^50); Vec(x*(1+x)*(1+x+x^2)*(x^2-x+1)/(1-x^6-x^7)) \\ G. C. Greubel, May 01 2017

G.f.: x*(1+x)*(1+x+x^2)*(x^2-x+1) / ( 1-x^6-x^7 ). - R. J. Mathar, Aug 26 2011

Edited by Ray Chandler and Robert G. Wilson v, Feb 06 2005

A103378 a(n) = a(n-10) + a(n-11) for n > 11, and a(n) = 1 for 1 <= n <= 11.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 7, 8, 8, 8, 8, 8, 8, 8, 8, 9, 12, 15, 16, 16, 16, 16, 16, 16, 16, 17, 21, 27, 31, 32, 32, 32, 32, 32, 32, 33, 38, 48, 58, 63, 64, 64, 64, 64, 64, 65, 71, 86, 106, 121, 127, 128

Offset: 1

Jonathan Vos Post, Feb 15 2005

			a(52)=17 because a(52)=a(52-10)+a(52-11) = a(42)+a(41) = 9 + 8.

Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,1,1).

Cf. A079398, A103372-A103377, A103379-A103380, A103388, A103398.

A103378 := proc(n) option remember; if n <= 11 then 1 ; else A103378(n-10)+A103378(n-11) ; fi ; end: seq(A103378(n),n=1..78) ; # R. J. Mathar, Nov 22 2007

k=10; Do[a[n]=1, {n, k+1}]; a[n_]:=a[n]=a[n-k]+a[n-k-1]; A103377=Array[a, 100] N[Solve[x^10 - x - 1 == 0, x], 111][[2]]
LinearRecurrence[Join[Table[0,{9}],{1,1}],Table[1,{11}],80] (* Harvey P. Dale, Aug 14 2013 *)

Vec((x^10-1)/(x-1)/(1-x^10-x^11)+O(x^80)) \\ M. F. Hasler, Sep 19 2015

G.f.: x*(1+x+x^2+x^3+x^4+x^5+x^6+x^7+x^8+x^9)/(1-x^10-x^11). - R. J. Mathar, Nov 22 2007

Corrected and extended by R. J. Mathar, Nov 22 2007

Edited by M. F. Hasler, Sep 19 2015

A103376 a(1) = a(2) = a(3) = a(4) = a(5) = a(6) = a(7) = a(8) = a(9) = 1 and for n>9: a(n) = a(n-8) + a(n-9).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 4, 4, 4, 4, 4, 4, 4, 5, 7, 8, 8, 8, 8, 8, 8, 9, 12, 15, 16, 16, 16, 16, 16, 17, 21, 27, 31, 32, 32, 32, 32, 33, 38, 48, 58, 63, 64, 64, 64, 65, 71, 86, 106, 121, 127, 128, 128, 129, 136, 157, 192, 227, 248, 255, 256, 257, 265, 293

Offset: 1

Jonathan Vos Post, Feb 05 2005

k=8 case of the family of sequences whose k=1 case is the Fibonacci sequence A000045, k=2 case is the Padovan sequence A000931 (offset so as to begin 1,1,1), k=3 case is A079398 (offset so as to begin 1,1,1,1), k=4 case is A103372, k=5 case is A103373, k=6 case is A103374 and k=7 case is A103375.
The general case for integer k>1 is defined: a(1) = a(2) = ... = a(k+1) and for n>(k+1) a(n) = a(n-k) + a(n-[k+1]).
For this k=8 case, the ratio of successive terms a(n)/a(n-1) approaches the unique positive root of the characteristic polynomial: x^9 - x - 1 = 0. This is the real constant (to 50 digits accuracy): 1.0850702454914508283368958640973142340506536310308 = A230162. Note that x = (1 + x)^(1/9) = (1 + (1 + (1 + ...)^(1/9))^(1/9))^(1/9).
The sequence of prime values in this k=8 case is A103386; The sequence of semiprime values in this k=8 case is A103396.

			a(93) = 1200 because a(93) = a(93-8) + a(93-9) = a(85) + a(84) = 642 + 558.

Zanten, A. J. van, "The golden ratio in the arts of painting, building and mathematics", Nieuw Archief voor Wiskunde, 4 (17) (1999) 229-245.

J.-P. Allouche and T. Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [The hal link does not always work. - _N. J. A. Sloane_, Feb 19 2025]
J.-P. Allouche and T. Johnson, Narayana's cows and delayed morphisms, in G. Assayag, M. Chemillier, and C. Eloy, Troisièmes Journées d'Informatique Musicale, JIM '96, Île de Tatihou, France, 1996, pp. 2-7. [Local copy with annotations and a correction from _N. J. A. Sloane_, Feb 19 2025]
Richard Padovan, Dom Hans van der Laan and the Plastic Number.
E. S. Selmer, On the irreducibility of certain trinomials, Math. Scand., 4 (1956) 287-302.
J. Shallit, A generalization of automatic sequences, Theoretical Computer Science, 61 (1988), 1-16.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,1,1).

Cf. A000045, A000931, A079398, A103372-A103375, A103377-A103380, A103386, A103396.

k = 8; Do[a[n] = 1, {n, k + 1}]; a[n_] := a[n] = a[n - k] + a[n - k - 1]; Array[a, 76]
LinearRecurrence[{0,0,0,0,0,0,0,1,1},{1,1,1,1,1,1,1,1,1},80] (* Harvey P. Dale, May 07 2015 *)

a(n)=([0,1,0,0,0,0,0,0,0; 0,0,1,0,0,0,0,0,0; 0,0,0,1,0,0,0,0,0; 0,0,0,0,1,0,0,0,0; 0,0,0,0,0,1,0,0,0; 0,0,0,0,0,0,1,0,0; 0,0,0,0,0,0,0,1,0; 0,0,0,0,0,0,0,0,1; 1,1,0,0,0,0,0,0,0]^(n-1)*[1;1;1;1;1;1;1;1;1])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

G.f.: x*(1+x)*(1+x^2)*(1+x^4)/(1-x^8-x^9). - R. J. Mathar, Dec 14 2009

a(1)=1, a(2)=1, a(3)=1, a(4)=1, a(5)=1, a(6)=1, a(7)=1, a(8)=1, a(9)=1, a(n)=a(n-8)+a(n-9). - Harvey P. Dale, May 07 2015

Edited by Ray Chandler, Feb 10 2005

A103388 Primes in A103378.

Original entry on oeis.org

2, 3, 5, 7, 17, 31, 71, 127, 157, 227, 257, 293, 349, 419, 503, 32299, 33343, 72421, 80429, 134269, 140473, 252761, 2499061, 201329923, 607488611, 1005428989, 2920552289, 8185638173, 10676478541, 14058719281, 15985335181, 34020175663, 159315910211, 1448256661853

Offset: 1

Jonathan Vos Post, Feb 15 2005

Cf. A103378, A103398.

A103378 := proc(n) option remember ; if n <= 11 then 1; else procname(n-10)+procname(n-11) ; fi; end: isA103378 := proc(n) option remember ; local i ; for i from 1 do if A103378(i) = n then RETURN(true) ; elif A103378(i) > n then RETURN(false) ; fi; od: end: A103388 := proc(n) option remember ; local a; if n = 1 then 2; else a := nextprime(procname(n-1)) ; while true do if isA103378(a) then RETURN(a) ; fi; a := nextprime(a) ; od: fi; end: for n from 1 to 37 do printf("%d, ",A103388(n)) ; od: # R. J. Mathar, Aug 30 2008

Clear[a]; k=10; Do[a[n]=1, {n, k+1}]; a[n_]:=a[n]=a[n-k]+a[n-k-1]; A103387=Union[Select[Array[a, 1000], PrimeQ]] (* See A103377 and A103397 for code related to those. - M. F. Hasler, Sep 19 2015, . *)

{a=vector(m=10, n, 1); L=0; for(n=m, 10^5, isprime(a[i=n%m+1]+=a[(n+1)%m+1]) && LM. F. Hasler, Sep 19 2015

Intersection of A103378 with A000040.

Corrected from a(16) on by R. J. Mathar, Aug 30 2008

Edited and more terms added by M. F. Hasler, Sep 19 2015

A152847 Triangle read by rows, A007318 rows repeated nine times .

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 3, 3, 1, 1, 4, 6, 4, 1, 1, 4, 6, 4, 1, 1, 4, 6, 4, 1

Offset: 0

Philippe Deléham, Dec 14 2008

Diagonal sums : A103377 .

			Triangle begins : 1 ; 1 ; 1 ; 1 ; 1 ; 1 ; 1 ; 1 ; 1 ; 1,1 ; 1,1 ; 1,1 ; 1,1 ; 1,1 ; 1,1 ; 1,1 ; 1,1 ; 1,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,2,1 ; 1,3,3,1 ; ...

Cf. A007318, A152198, A152815, A152828, A152830, A152831, A152844, A152845, A152846.

Flatten[Table[#, {9}] & /@ Table[Binomial[n, k], {n, 0, 10}, {k, 0, n}]] (* G. C. Greubel, May 03 2017 *)

A103398 Semiprimes in A103378.

Original entry on oeis.org

4, 9, 15, 21, 33, 38, 58, 65, 86, 106, 121, 129, 265, 511, 2047, 2049, 4109, 16293, 16489, 17855, 19857, 32678, 34709, 66217, 104739, 220918, 240367, 262298, 293323, 954413, 2082999, 3145729, 3498467, 4296813, 16442015, 18037939, 21317326

Offset: 1

Jonathan Vos Post, Feb 15 2005

Cf. A001358, A103378, A103388.

A103378 := proc(n) option remember; if n <= 11 then 1 ; else procname(n-10)+procname(n-11) ; fi ; end proc:
a78prev := -1 ; for n from 1 to 400 do a78 := A103378(n) ; if numtheory[bigomega](a78) = 2 and a78 <> a78prev then printf("%d,",a78) ; end if; a78prev := a78 ; end do: # R. J. Mathar, Jun 11 2010

SemiprimeQ[n_]:=Plus@@FactorInteger[n][[All, 2]]?2; Clear[a]; k=10; Do[a[n]=1, {n, k+1}]; a[n_]:=a[n]=a[n-k]+a[n-k-1]; A103377=Array[a, 100] A103387=Union[Select[Array[a, 1000], PrimeQ]] A103397=Union[Select[Array[a, 300], SemiprimeQ]] N[Solve[x^11 - x - 1 == 0, x], 111][[2]] (* Ray Chandler and Robert G. Wilson v *)

Intersection of A103378 with A001358.

Edited and extended by Ray Chandler and Robert G. Wilson v

Entries >511 corrected by R. J. Mathar, Jun 11 2010

cp's OEIS Frontend

A103387 Primes in A103377.

Views

Author

Keywords

Crossrefs

Programs

Mathematica

Formula

A103397 Semiprimes in A103377.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Formula

A103373 a(1) = a(2) = a(3) = a(4) = a(5) = a(6) = 1 and for n>6: a(n) = a(n-5) + a(n-6).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A103374 a(1) = a(2) = a(3) = a(4) = a(5) = a(6) = a(7) = 1 and for n>7: a(n) = a(n-6) + a(n-7).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A103378 a(n) = a(n-10) + a(n-11) for n > 11, and a(n) = 1 for 1 <= n <= 11.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A103376 a(1) = a(2) = a(3) = a(4) = a(5) = a(6) = a(7) = a(8) = a(9) = 1 and for n>9: a(n) = a(n-8) + a(n-9).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A103388 Primes in A103378.

Views

Author

Keywords

Crossrefs