OFFSET
0,2
COMMENTS
Equals A000217 (1, 3, 6, 10, 15, ...) convolved with A193356 (1, 0, 3, 0, 5, ...). - Gary W. Adamson, Feb 16 2009
F(1,4,n) is the number of bracelets with 1 blue, 4 red and n black beads. If F(1,4,1)=3 and F(1,4,2)=9 taken as a base;
F(1,4,n) = n(n+1)(n+2)/6+F(1,2,n) + F(1,4,n-2). [F(1,2,n) is the number of bracelets with 1 blue, 2 red and n black beads. If F(1,2,1)=2 and F(1,2,2)=4 taken as a base F(1,2,n)=n+1+F(1,2,n-2)]. - Ata Aydin Uslu and Hamdi G. Ozmenekse, Jan 11 2012
a(A254338(n)) = 6 for n > 0. - Reinhard Zumkeller, Feb 27 2015
REFERENCES
S. J. Cyvin et al., Polygonal systems including the corannulene and coronene homologs: novel applications of Pólya's theorem, Z. Naturforsch., 52a (1997), 867-873.
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber. 30 (1897), 1917-1926.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
T. D. Noe, Table of n, a(n) for n = 0..1000
Johann Cigler, Some remarks on Rogers-Szegö polynomials and Losanitsch's triangle, arXiv:1711.03340 [math.CO], 2017.
S. M. Losanitsch, Die Isomerie-Arten bei den Homologen der Paraffin-Reihe, Chem. Ber. 30 (1897), 1917-1926. (Annotated scanned copy)
N. J. A. Sloane, Classic Sequences
Index entries for linear recurrences with constant coefficients, signature (3, -1, -5, 5, 1, -3, 1).
"http://commons.wikimedia.org/wiki/File:Bracelet_Problem_(Bileklik_problemi).pdf" number of bracelets made with 1 blue, 4 red and n black beads [From Ata Aydin Uslu and Hamdi G. Ozmenekse, Jan 11 2012].
"http://commons.wikimedia.org/wiki/File:Bileklik_Problemi_(Bracelet_Problem).pdf" number of bracelets made with 1 blue, 2 red and n black beads [From Ata Aydin Uslu and Hamdi G. Ozmenekse, Jan 12 2012].
FORMULA
G.f.: (1+x^2)/((1-x)^3*(1-x^2)^2) = (1+x^2)/((1-x)^5*(1+x)^2).
l(c, r) = 1/2 C(c+r-3, r) + 1/2 d(c, r), where d(c, r) is C((c + r - 3)/2, r/2) if c is odd and r is even, 0 if c is even and r is odd, C((c + r - 4)/2, r/2) if c is even and r is even, C((c + r - 4)/2, (r - 1)/2) if c is odd and r is odd.
a(-5-n)=a(n). - Michael Somos, Mar 08 2007
Euler transform of length 4 sequence [3, 3, 0, -1]. - Michael Somos, Mar 08 2007
a(n) = 3a(n-1) - a(n-2) - 5a(n-3) + 5a(n-4) + a(n-5) - 3a(n-6) + a(n-7), with a(0)=1, a(1)=3, a(2)=9, a(4)=19, a(5)=38, a(6)=66, a(7)=110. - Harvey P. Dale, May 02 2011
a(n) = A006009(n)/2 - A000332(n+4) = ((1/2)*Sum_{i=1..n+1} (i+1)*floor((i+1)^2/2)) - binomial(n+4,4). - Enrique Pérez Herrero, May 11 2012
a(n) = (1/48)*(n+1)*(n+3)*((n+2)*(n+4)+3)+1/32*(2*n+5)*(1+(-1)^n). - Yosu Yurramendi, Jun 20 2013
MAPLE
a:= n -> (Matrix([[1, 0$4, 1, 3]]). Matrix(7, (i, j)-> if (i=j-1) then 1 elif j=1 then [3, -1, -5, 5, 1, -3, 1][i] else 0 fi)^n)[1, 1]: seq (a(n), n=0..40); # Alois P. Heinz, Jul 31 2008
MATHEMATICA
LinearRecurrence[{3, -1, -5, 5, 1, -3, 1}, {1, 3, 9, 19, 38, 66, 110}, 50] (* or *) CoefficientList[Series[(1+x^2)/((1-x)^3(1-x^2)^2), {x, 0, 50}], x] (* Harvey P. Dale, May 02 2011 *)
nn=45; With[{a=Accumulate[Range[nn]], b=Riffle[Range[1, nn, 2], 0]}, Flatten[ Table[ListConvolve[Take[a, n], Take[b, n]], {n, nn}]]] (* Harvey P. Dale, Nov 11 2011 *)
PROG
(PARI) {a(n)=if(n<-4, n=-5-n); polcoeff( (1+x^2)/((1-x)^3*(1-x^2)^2)+x*O(x^n), n)} /* Michael Somos, Mar 08 2007 */
(Haskell) Following Gary W. Adamson.
import Data.List (inits, intersperse)
a005994 n = a005994_list !! n
a005994_list = map (sum . zipWith (*) (intersperse 0 [1, 3 ..]) . reverse) $
tail $ inits $ tail a000217_list
-- Reinhard Zumkeller, Feb 27 2015
CROSSREFS
KEYWORD
nonn,easy,nice
AUTHOR
N. J. A. Sloane, Winston C. Yang (yang(AT)math.wisc.edu)
STATUS
approved