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.

A224809 Number of permutations (p(1), p(2), ..., p(n)) satisfying -k <= p(i)-i <= r and p(i)-i not in the set I, i=1..n, with k=2, r=4, I={-1,1,2,3}.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 4, 6, 9, 12, 16, 24, 36, 54, 81, 117, 169, 247, 361, 532, 784, 1148, 1681, 2460, 3600, 5280, 7744, 11352, 16641, 24381, 35721, 52353, 76729, 112462, 164836, 241570, 354025, 518840, 760384, 1114416, 1633284
Offset: 0

Views

Author

Vladimir Baltic, May 16 2013

Keywords

Comments

Number of subsets of {1,2,...,n-4} without differences equal to 2 or 4.

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[-(x-1)*(1+x+x^2)/((x^3+x-1)*(x^6-x^4-1)), {x, 0, 50}], x] (* G. C. Greubel, Apr 28 2017 *)
  • PARI
    N = 42; x = 'x + O('x^N);
    Vec(Ser(-(x-1)*(1+x+x^2)/((x^3+x-1)*(x^6-x^4-1))))  \\ Gheorghe Coserea, Nov 11 2016

Formula

a(n) = a(n-1) + a(n-3) - a(n-4) + a(n-5) + a(n-6) - a(n-9).
G.f.: -(x-1)*(1+x+x^2) / ( (x^3+x-1)*(x^6-x^4-1) ).
a(2*k) = (A000930(k))^2, a(2*k+1) = A000930(k) * A000930(k+1).