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.

A002714 Number of different keys with n cuts, depths between 1 and 7 and depth difference at most 1 between adjacent cut depths.

Original entry on oeis.org

1, 7, 19, 53, 149, 421, 1193, 3387, 9627, 27383, 77923, 221805, 631469, 1797957, 5119593, 14578387, 41514003, 118218823, 336653331, 958698053, 2730124261, 7774706437, 22140438345, 63050541515, 179552587883, 511322221559, 1456121982755, 4146683677885
Offset: 0

Views

Author

Keywords

Comments

Also number of base 7 n-digit numbers with adjacent digits differing by one or less.
[Empirical] a(base,n)=a(base-1,n)+3^(n-1) for base>=n; a(base,n)=a(base-1,n)+3^(n-1)-2 when base=n-1. - R. H. Hardin, Dec 26 2006

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Programs

  • Maple
    A002714:=-(7-9*z-9*z**2+3*z**3)/(-1+4*z-2*z**2-4*z**3+z**4); # conjectured by Simon Plouffe in his 1992 dissertation; correct up to offset
    T := proc(d,n) option remember ; if n = 1 then 1; else if d = 7 then T(d,n-1)+T(d-1,n-1) ; elif d = 1 then T(d,n-1)+T(d+1,n-1) ; else T(d-1,n-1)+T(d,n-1)+T(d+1,n-1) ; fi ; fi ; end: A002714 := proc(n) local d ; add( T(d,n),d=1..7) ; end: seq(A002714(n),n=1..35) ; # R. J. Mathar, Jun 15 2008
  • Mathematica
    CoefficientList[Series[(2*x^4-5*x^3-7*x^2+3*x+1)/(-x^4+4*x^3+2*x^2-4*x+1),{x,0,200}],x] (* Vincenzo Librandi, Aug 13 2012 *)
    Join[{1}, LinearRecurrence[{4, -2, -4, 1}, {7, 19, 53, 149}, 30]] (* Jean-François Alcover, Jan 07 2019 *)
  • PARI
    /* from the Knopfmacher et al. reference */
    default(realprecision,99); /* using floats */
    sn(n,k)=1/n*sum(i=1,k,sumdiv(n,j,eulerphi(j)*(1+2*cos(i*Pi/(k+1)))^(n/j)));
    vector(66,n, if (n==1,1,round(sn(n-1,7))) )
    /* Joerg Arndt, Aug 13 2012 */

Formula

G.f.: (2*x^4 - 5*x^3 - 7*x^2 + 3*x + 1)/(-x^4 + 4*x^3 + 2*x^2 - 4*x + 1); (from the Knopfmacher et al. reference). - Joerg Arndt, Aug 10 2012

Extensions

Information added from A126361, offset changed to 0 by Joerg Arndt, Aug 13 2012