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.

Previous Showing 21-30 of 38 results. Next

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

A061451 Array T(n,k) of k-th order Fibonacci numbers read by antidiagonals in up-direction.

Original entry on oeis.org

2, 3, 3, 4, 5, 5, 5, 7, 9, 8, 6, 9, 13, 17, 13, 7, 11, 17, 25, 31, 21, 8, 13, 21, 33, 49, 57, 34, 9, 15, 25, 41, 65, 94, 105, 55, 10, 17, 29, 49, 81, 129, 181, 193, 89, 11, 19, 33, 57, 97, 161, 253, 349, 355, 144, 12, 21, 37, 65, 113, 193, 321, 497, 673, 653, 233
Offset: 1

Views

Author

Frank Ellermann, Jun 11 2001

Keywords

Examples

			2, 3, 5, 8 ... first order, a(1)=2, a(3)=3, a(6)=5, a(10)=8, ...
3, 5, 9,17 ... 2nd order, a(2)=3, a(5)=5, a(9)=9, ...
4, 7,13,25 ... 3rd order, a(4)=4, a(8)=7, ...
5, 9,17,33 ... 4th order, a(7)=5, ...
		

References

  • N. Wirth, Algorithmen und Datenstrukturen, 1975, table 2.15 (ch. 2.3.4)

Crossrefs

Programs

  • Mathematica
    max = 12; Clear[f]; f[k_, n_] /; n > k := f[k, n] = Sum[f[k, n - j], {j, 1, k + 1}]; f[k_, n_] = 1; t = Table[ Table[ f[k, n], {n, k + 1, max}], {k, 1, max}]; Table[ t[[k - n + 1, n]], {k, 1, max - 1}, {n, 1, k}] // Flatten (* Jean-François Alcover, Apr 10 2013 *)

A253705 Indices of primes in the 8th-order Fibonacci number sequence, A079262.

Original entry on oeis.org

9, 17, 25, 125, 350, 1322, 108935, 199528
Offset: 1

Views

Author

Robert Price, Jan 09 2015

Keywords

Comments

a(9) > 2*10^5.

Crossrefs

Programs

  • Mathematica
    a={0,0,0,0,0,0,0,1}; step=8; lst={}; For[n=step,n<=1000,n++, sum=Plus@@a; If[PrimeQ[sum], AppendTo[lst,n]]; a=RotateLeft[a]; a[[step]]=sum]; lst
    Flatten[Position[LinearRecurrence[{1,1,1,1,1,1,1,1},{0,0,0,0,0,0,0,1},200000],?PrimeQ]]-1 (* The program takes a long time to run *) (* _Harvey P. Dale, Apr 26 2018 *)
  • PARI
    lista(nn) = {gf = x^7/(1-x-x^2-x^3-x^4-x^5-x^6-x^7-x^8); for (n=0, nn, if (isprime(polcoeff(gf+O(x^(n+1)), n)), print1(n, ", ")););} \\ Michel Marcus, Jan 12 2015

A163551 13th-order Fibonacci numbers: a(n) = a(n-1) + ... + a(n-13) with a(1)=...=a(13)=1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 13, 25, 49, 97, 193, 385, 769, 1537, 3073, 6145, 12289, 24577, 49153, 98305, 196597, 393169, 786289, 1572481, 3144769, 6289153, 12577537, 25153537, 50304001, 100601857, 201191425, 402358273, 804667393
Offset: 1

Views

Author

Jainit Purohit (mjainit(AT)gmail.com), Jul 30 2009

Keywords

Crossrefs

Cf. A000045 (Fibonacci numbers), A000213 (tribonacci), A000288 (tetranacci), A000322 (pentanacci), A000383 (hexanacci), A060455 (heptanacci), A123526 (octanacci), A127193 (nonanacci), A127194 (decanacci), A127624 (undecanacci), A207539 (dodecanacci).

Programs

  • Mathematica
    With[{c=Table[1,{13}]},LinearRecurrence[c,c,40]] (* Harvey P. Dale, Aug 09 2013 *)
  • PARI
    x='x+O('x^50); Vec((1-x^2 -2*x^3-3*x^4 -4*x^5-5*x^6 -6*x^7-7*x^8 -8*x^9 -9*x^10 -10*x^11 -11*x^12) / (1-x-x^2-x^3-x^4-x^5-x^6-x^7-x^8-x^9-x^10-x^11-x^12-x^13)) \\ G. C. Greubel, Jul 28 2017

Formula

a(n) = a(n-1)+a(n-2)+...+a(n-13) for n > 12, a(0)=a(1)=...=a(12)=1.
G.f.: (-1)*(-1+x^2+2*x^3+3*x^4+4*x^5+5*x^6+6*x^7+7*x^8+8*x^9+9*x^10 +10*x^11 +11*x^12) / (1-x-x^2-x^3-x^4-x^5-x^6-x^7-x^8-x^9-x^10-x^11-x^12-x^13). - Michael Burkhart, Feb 18 2012

Extensions

Values adapted to the definition by R. J. Mathar, Aug 01 2009

A254412 Indices of primes in the 8th-order Fibonacci number sequence, A123526.

Original entry on oeis.org

11, 13, 15, 24, 30, 33, 57, 104, 121, 132, 149, 158, 178, 220, 295, 389, 1070, 1101, 1373, 1761, 1778, 2333, 2731, 4541, 5189, 5237, 5738, 8025, 8787, 10561, 11783, 13435, 14638, 15337, 20985, 21722, 24770, 31009, 57367, 65877, 129773, 134630, 167020
Offset: 1

Views

Author

Robert Price, Jan 30 2015

Keywords

Comments

a(44) > 2*10^5.

Crossrefs

Programs

  • Mathematica
    a={1,1,1,1,1,1,1,1}; step=8; lst={}; For[n=step+1,n<=1000,n++, sum=Plus@@a; If[PrimeQ[sum], AppendTo[lst,n]]; a=RotateLeft[a]; a[[step]]=sum]; lst

A254413 Primes in the 8th-order Fibonacci numbers A123526.

Original entry on oeis.org

29, 113, 449, 226241, 14307889, 113783041, 1820091580429249, 233322881089059894782836851617, 29566627412209231076314948970028097, 59243719929958343565697204780596496129, 7507351981539044730893385057192143660843521
Offset: 1

Views

Author

Robert Price, Jan 30 2015

Keywords

Comments

a(12) is too large to display here. It has 46 digits and is the 158th term in A123526.

Crossrefs

Programs

  • Mathematica
    a={1,1,1,1,1,1,1,1}; step=8; lst={}; For[n=step+1,n<=1000,n++, sum=Plus@@a; If[PrimeQ[sum], AppendTo[lst,sum]]; a=RotateLeft[a]; a[[step]]=sum]; lst
    Select[With[{lr=PadRight[{},8,1]},LinearRecurrence[lr,lr,200]],PrimeQ] (* Harvey P. Dale, Dec 03 2022 *)

A174323 Numbers n such that omega(Fibonacci(n)) is a square.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 11, 13, 17, 20, 23, 24, 27, 28, 29, 32, 43, 47, 52, 55, 74, 77, 80, 83, 84, 85, 87, 88, 91, 93, 96, 97, 100, 108, 115, 123, 131, 132, 137, 138, 143, 146, 149, 156, 157, 161, 163, 178, 184, 187, 189, 196, 197, 209, 211, 214, 215, 221, 222, 223, 232
Offset: 1

Views

Author

Michel Lagneau, Mar 15 2010

Keywords

Comments

Numbers n such that omega(A000045(n)) is a square, where omega(p) is the number of distinct prime factors of p (A001221). Remark: for the larger Fibonacci numbers F(n) (n > 300), the Maple program (below) is very slow. So we use a two-step process: factoring F(n) with the elliptic curve method, and then calculate the distinct prime factors.

Examples

			omega(Fibonacci(1)) = omega(Fibonacci(2)) = omega(1) = 0,
omega(Fibonacci(3)) = omega(2) = 1,
omega(Fibonacci(20)) = omega(6765) = 4,
omega(Fibonacci(80)) = omega(23416728348467685) = 9.
		

References

  • Majorie Bicknell and Verner E Hoggatt, Fibonacci's Problem Book, Fibonacci Association, San Jose, Calif., 1974.
  • Alfred Brousseau, Fibonacci and Related Number Theoretic Tables, The Fibonacci Association, 1972, pages 1-8.

Crossrefs

Cf. A038575 (number of prime factors of n-th Fibonacci number, with multiplicity).
Cf. A022307 (number of distinct prime factors of n-th Fibonacci number), A086597 (number of primitive prime factors).

Programs

  • Magma
    [k:k in [1..240]| IsSquare(#PrimeDivisors(Fibonacci(k)))]; // Marius A. Burtea, Oct 15 2019
  • Maple
    with(numtheory):u0:=0:u1:=1:for p from 2 to 400 do :s:=u0+u1:u0:=u1:u1:=s: s1:=nops( ifactors(s)[2]): w1:=sqrt(s1):w2:=floor(w1):if w1=w2 then print (p): else fi:od:
    # alternative:
    P[1]:= {}: count:= 1: res:= 1:
    for i from 2 to 300 do
      pn:= map(t -> i/t, numtheory:-factorset(i));
      Cprimes:= `union`(seq(P[t],t=pn));
      f:= combinat:-fibonacci(i);
      for p in Cprimes do f:= f/p^padic:-ordp(f,p) od;
      P[i]:= Cprimes union numtheory:-factorset(f);
      if issqr(nops(P[i])) then
         count:= count+1;
         res:= res, i;
      fi;
    od:
    res; # Robert Israel, Oct 13 2016
  • Mathematica
    Select[Range[200], IntegerQ[Sqrt[PrimeNu[Fibonacci[#]]]] &] (* G. C. Greubel, May 16 2017 *)
  • PARI
    is(n)=issquare(omega(fibonacci(n))) \\ Charles R Greathouse IV, Oct 13 2016
    

A253706 Primes in the 8th-order Fibonacci numbers A079262.

Original entry on oeis.org

2, 509, 128257, 133294824621464999938178340471931877, 4596852049500861351052672455121859744010232939954169259264638023409631672658340253083284317818242062413
Offset: 1

Views

Author

Robert Price, Jan 09 2015

Keywords

Comments

a(6) is too large to display here. It has 395 digits and is the 1322nd term in A079262.

Crossrefs

Programs

  • Mathematica
    a={0,0,0,0,0,0,0,1}; step=8; lst={}; For[n=step,n<=1000,n++, sum=Plus@@a; If[PrimeQ[sum], AppendTo[lst,sum]]; a=RotateLeft[a]; a[[step]]=sum]; lst
  • PARI
    lista(nn) = {gf = x^7/(1-x-x^2-x^3-x^4-x^5-x^6-x^7-x^8); for (n=0, nn, if (isprime(p=polcoeff(gf+O(x^(n+1)), n)), print1(p, ", ")););} \\ Michel Marcus, Jan 12 2015

A255529 Indices of primes in the 9th-order Fibonacci number sequence, A104144.

Original entry on oeis.org

10, 19, 878
Offset: 1

Views

Author

Robert Price, Feb 24 2015

Keywords

Comments

a(4) > 2*10^5.

Crossrefs

Programs

  • Mathematica
    a={0,0,0,0,0,0,0,0,1}; step=9; lst={}; For[n=step,n<=1000,n++, sum=Plus@@a; If[PrimeQ[sum], AppendTo[lst,n]]; a=RotateLeft[a]; a[[step]]=sum]; lst
  • PARI
    a104144(n) = polcoeff(x^8/(1-x-x^2-x^3-x^4-x^5-x^6-x^7-x^8-x^9) + O(x^(n+1)), n);
    lista(nn) = {for (n=1, nn, if (isprime(a104144(n)), print1(n, ", ")););} \\ Michel Marcus, Feb 27 2015

A255530 Indices of primes in the 9th-order Fibonacci number sequence, A251746.

Original entry on oeis.org

10, 19, 59, 79, 12487
Offset: 1

Views

Author

Robert Price, Feb 24 2015

Keywords

Comments

a(6) > 2*10^5.

Crossrefs

Programs

  • Mathematica
    a={0,0,0,0,0,0,0,1,0}; step=9; lst={}; For[n=step,n<=1000,n++, sum=Plus@@a; If[PrimeQ[sum], AppendTo[lst,n]]; a=RotateLeft[a]; a[[step]]=sum]; lst
Previous Showing 21-30 of 38 results. Next