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

A244930 Indices of primes in A214831.

Original entry on oeis.org

3, 4, 7, 8, 16, 26, 34, 42, 78, 94, 101, 107, 216, 255, 543, 562, 851, 981, 1099, 1528, 1824, 1955, 2122, 2488, 2500, 15331, 15961, 24107, 24938, 26051, 58504, 61617, 81034, 85119, 89768, 90597, 97191, 116899, 195346
Offset: 1

Views

Author

Robert Price, Jul 08 2014

Keywords

Comments

a(40) > 2*10^5.

Crossrefs

Programs

A244931 Prime terms in A214831.

Original entry on oeis.org

19, 37, 223, 409, 53617, 23757289, 3111662089, 407556643177, 1372675688565303822697, 23548271681390871672120649, 1676892190264006259992141409, 64923481849284379431377700019
Offset: 1

Views

Author

Robert Price, Jul 08 2014

Keywords

Comments

a(13) has 58 digits and thus is too large to display here. It corresponds to A214831(216).

Crossrefs

Programs

  • Mathematica
    a={1,9,9}; For[n=3, n<=1000, n++, sum=Plus@@a; If[PrimeQ[sum], Print[sum]]; a=RotateLeft[a]; a[[3]]=sum]

A214727 a(n) = a(n-1) + a(n-2) + a(n-3) with a(0) = 1, a(1) = a(2) = 2.

Original entry on oeis.org

1, 2, 2, 5, 9, 16, 30, 55, 101, 186, 342, 629, 1157, 2128, 3914, 7199, 13241, 24354, 44794, 82389, 151537, 278720, 512646, 942903, 1734269, 3189818, 5866990, 10791077, 19847885, 36505952, 67144914, 123498751, 227149617, 417793282
Offset: 0

Views

Author

Abel Amene, Jul 27 2012

Keywords

Comments

Part of a group of sequences defined by a(0), a(1)=a(2), a(n) = a(n-1) + a(n-2) + a(n-3) which is a subgroup of sequences with linear recurrences and constant coefficients listed in the index.
Note: A000073 (with offset=1), 1 followed by A000073, A000213, A141523, A214727, A214825 to A214831 completely define possible sequences with a(0)=0,1,2...9 and a(1)=a(2)=0,1,2...9 excluding any multiples of these sequences and the trivial case of a(0)=a(1)=a(2)=0.
Note: allowing a(0)=0 and a(1)=a(2)=1,2,3....9 leads to A000073 (with offset=1) and its multiples.
Note: allowing a(0)=1,2,3....9 a(1)=a(2)=0 leads to 1 followed by A000073 and its multiples.
With offset of 6 this sequence is the 8th row of tribonacci array A136175.

Examples

			G.f. = 1 + 2*x + 2*x^2 + 5 x^3 + 9*x^4 + 16*x^5 + 30*x^6 + 55*x^7 + ...
		

Crossrefs

Programs

  • GAP
    a:=[1,2,2];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 23 2019
  • Haskell
    a214727 n = a214727_list !! n
    a214727_list = 1 : 2 : 2 : zipWith3 (\x y z -> x + y + z)
       a214727_list (tail a214727_list) (drop 2 a214727_list)
    -- Reinhard Zumkeller, Jul 31 2012
    
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 23 2019
    
  • Mathematica
    LinearRecurrence[{1,1,1},{1,2,2},40] (* Ray Chandler, Dec 08 2013 *)
  • PARI
    a(n)=([0,1,0; 0,0,1; 1,1,1]^n*[1;2;2])[1,1] \\ Charles R Greathouse IV, Mar 22 2016
    
  • PARI
    my(x='x+O('x^40)); Vec((1+x-x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 23 2019
    
  • SageMath
    ((1+x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 23 2019
    

Formula

G.f.: (1+x-x^2)/(1-x-x^2-x^3).
a(n) = K(n) -2*T(n+1) + 3*T(n), where K(n) = A001644(n), T(n) = A000073(n+1). - G. C. Greubel, Apr 23 2019
a(n) = Sum_{r root of x^3-x^2-x-1} r^n/(-r^2+2*r+1). - Fabian Pereyra, Nov 20 2024

A214825 a(n) = a(n-1) + a(n-2) + a(n-3), with a(0) = 1, a(1) = a(2) = 3.

Original entry on oeis.org

1, 3, 3, 7, 13, 23, 43, 79, 145, 267, 491, 903, 1661, 3055, 5619, 10335, 19009, 34963, 64307, 118279, 217549, 400135, 735963, 1353647, 2489745, 4579355, 8422747, 15491847, 28493949, 52408543, 96394339, 177296831, 326099713, 599790883, 1103187427
Offset: 0

Views

Author

Abel Amene, Jul 28 2012

Keywords

Comments

Part of a group of sequences defined by a(0), a(1)=a(2), a(n) = a(n-1) + a(n-2) + a(n-3) which is a subgroup of sequences with linear recurrences and constant coefficients listed in the index. See Comments in A214727.

Crossrefs

Programs

  • GAP
    a:=[1,3,3];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 23 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+2*x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 23 2019
    
  • Mathematica
    LinearRecurrence[{1,1,1},{1,3,3},40] (* Harvey P. Dale, Oct 05 2013 *)
  • PARI
    a(n)=([0,1,0; 0,0,1; 1,1,1]^n*[1;3;3])[1,1] \\ Charles R Greathouse IV, Mar 22 2016
    
  • PARI
    my(x='x+O('x^40)); Vec((1+2*x-x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 23 2019
    
  • SageMath
    ((1+2*x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 23 2019
    

Formula

G.f.: (1+2*x-x^2)/(1-x-x^2-x^3).
a(n) = K(n) - 2*T(n+1) + 4*T(n), where K(n) = A001644(n), and T(n) = A000073(n+1). - G. C. Greubel, Apr 23 2019

A214827 a(n) = a(n-1) + a(n-2) + a(n-3), with a(0) = 1, a(1) = a(2) = 5.

Original entry on oeis.org

1, 5, 5, 11, 21, 37, 69, 127, 233, 429, 789, 1451, 2669, 4909, 9029, 16607, 30545, 56181, 103333, 190059, 349573, 642965, 1182597, 2175135, 4000697, 7358429, 13534261, 24893387, 45786077, 84213725, 154893189, 284892991, 523999905
Offset: 0

Views

Author

Abel Amene, Jul 29 2012

Keywords

Comments

See comments in A214727.

Crossrefs

Programs

  • GAP
    a:=[1,5,5];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 24 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+4*x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 24 2019
    
  • Mathematica
    LinearRecurrence[{1,1,1},{1,5,5},40] (* Ray Chandler, Dec 08 2013 *)
  • PARI
    my(x='x+O('x^40)); Vec((1+4*x-x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 24 2019
    
  • Sage
    ((1+4*x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 24 2019
    

Formula

G.f.: (x^2-4*x-1)/(x^3+x^2+x-1).
a(n) = -A000073(n) + 4*A000073(n+1) + A000073(n+2). - R. J. Mathar, Jul 29 2012

A214828 a(n) = a(n-1) + a(n-2) + a(n-3), with a(0) = 1, a(1) = a(2) = 6.

Original entry on oeis.org

1, 6, 6, 13, 25, 44, 82, 151, 277, 510, 938, 1725, 3173, 5836, 10734, 19743, 36313, 66790, 122846, 225949, 415585, 764380, 1405914, 2585879, 4756173, 8747966, 16090018, 29594157, 54432141, 100116316, 184142614, 338691071, 622950001
Offset: 0

Views

Author

Abel Amene, Jul 30 2012

Keywords

Comments

See comments in A214727.

Crossrefs

Programs

  • GAP
    a:=[1,6,6];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 24 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+5*x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 24 2019
    
  • Mathematica
    LinearRecurrence[{1,1,1},{1,6,6},33] (* Ray Chandler, Dec 08 2013 *)
  • PARI
    my(x='x+O('x^40)); Vec((1+5*x-x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 24 2019
    
  • Sage
    ((1+5*x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 24 2019
    

Formula

G.f.: (1+5*x-x^2)/(1-x-x^2-x^3).
a(n) = -A000073(n) + 5*A000073(n+1) + A000073(n+2). - G. C. Greubel, Apr 24 2019

A214829 a(n) = a(n-1) + a(n-2) + a(n-3), with a(0) = 1, a(1) = a(2) = 7.

Original entry on oeis.org

1, 7, 7, 15, 29, 51, 95, 175, 321, 591, 1087, 1999, 3677, 6763, 12439, 22879, 42081, 77399, 142359, 261839, 481597, 885795, 1629231, 2996623, 5511649, 10137503, 18645775, 34294927, 63078205, 116018907, 213392039, 392489151, 721900097, 1327781287, 2442170535
Offset: 0

Views

Author

Abel Amene, Aug 07 2012

Keywords

Comments

See comments in A214727.

Crossrefs

Programs

  • GAP
    a:=[1,7,7];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 24 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+6*x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 24 2019
    
  • Mathematica
    LinearRecurrence[{1,1,1}, {1,7,7}, 40] (* G. C. Greubel, Apr 24 2019 *)
  • PARI
    Vec((x^2-6*x-1)/(x^3+x^2+x-1) + O(x^40)) \\ Michel Marcus, Jun 04 2017
    
  • Sage
    ((1+6*x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 24 2019
    

Formula

G.f.: (1+6*x-x^2)/(1-x-x^2-x^3).
a(n) = -A000073(n) + 6*A000073(n+1) + A000073(n+2). - G. C. Greubel, Apr 24 2019

A214826 a(n) = a(n-1) + a(n-2) + a(n-3), with a(0) = 1, a(1) = a(2) = 4.

Original entry on oeis.org

1, 4, 4, 9, 17, 30, 56, 103, 189, 348, 640, 1177, 2165, 3982, 7324, 13471, 24777, 45572, 83820, 154169, 283561, 521550, 959280, 1764391, 3245221, 5968892, 10978504, 20192617, 37140013, 68311134, 125643764, 231094911, 425049809
Offset: 0

Views

Author

Abel Amene, Jul 29 2012

Keywords

Comments

See Comments in A214727.

Crossrefs

Programs

  • GAP
    a:=[1,4,4];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 23 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+3*x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 23 2019
    
  • Mathematica
    LinearRecurrence[{1,1,1},{1,4,4},33] (* Ray Chandler, Dec 08 2013 *)
  • PARI
    my(x='x+O('x^40)); Vec((1+3*x-x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 23 2019
    
  • Sage
    ((1+3*x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 23 2019
    

Formula

G.f.: (1+3*x-x^2)/(1-x-x^2-x^3).
a(n) = K(n) - 2*T(n+1) + 5*T(n), where K(n) = A001644(n) and T(n) = A000073(n+1). - G. C. Greubel, Apr 23 2019

A214830 a(n) = a(n-1) + a(n-2) + a(n-3), with a(0) = 1, a(1) = a(2) = 8.

Original entry on oeis.org

1, 8, 8, 17, 33, 58, 108, 199, 365, 672, 1236, 2273, 4181, 7690, 14144, 26015, 47849, 88008, 161872, 297729, 547609, 1007210, 1852548, 3407367, 6267125, 11527040, 21201532, 38995697, 71724269, 131921498, 242641464, 446287231, 820850193, 1509778888
Offset: 0

Views

Author

Abel Amene, Aug 07 2012

Keywords

Comments

See comments in A214727.

Crossrefs

Programs

  • GAP
    a:=[1,8,8];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 24 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 40); Coefficients(R!( (1+7*x-x^2)/(1-x-x^2-x^3) )); // G. C. Greubel, Apr 24 2019
    
  • Mathematica
    CoefficientList[Series[(x^2-7*x-1)/(x^3+x^2+x-1), {x, 0, 40}], x] (* Wesley Ivan Hurt, Jun 18 2014 *)
    LinearRecurrence[{1,1,1}, {1,8,8}, 40] (* G. C. Greubel, Apr 24 2019 *)
  • PARI
    my(x='x+O('x^40)); Vec((1+7*x-x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 24 2019
    
  • Sage
    ((1+7*x-x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 24 2019
    

Formula

G.f.: (1+7*x-x^2)/(1-x-x^2-x^3).
a(n) = -A000073(n) + 7*A000073(n+1) + A000073(n+2). - G. C. Greubel, Apr 24 2019

A246517 Indices of primes in A141036.

Original entry on oeis.org

0, 5, 14, 26, 33, 222, 234, 482, 937, 1170, 1290, 1877, 1897, 3413, 6017, 9365, 47470, 48254, 97421, 102057, 119689, 132418, 192517, 194442
Offset: 1

Views

Author

Robert Price, Aug 28 2014

Keywords

Comments

a(25) > 2*10^5.
A141036(a(n)) = A246518(n).

Crossrefs

Programs

  • Haskell
    a246517 n = a246517_list !! (n-1)
    a246517_list = filter ((== 1) . a010051'' . a141036) [0..]
    -- Reinhard Zumkeller, Sep 15 2014
  • Mathematica
    a={2,1,1}; Print[0]; For[n=3, n<=1000, n++, sum=Plus@@a; If[PrimeQ[sum], Print[n]]; a=RotateLeft[a]; a[[3]]=sum]
Showing 1-10 of 12 results. Next