This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A093951 Sum of integers generated by n-1 substitutions, starting with 1, k -> k+1,k-1,..,1. 1
 1, 2, 4, 8, 17, 36, 80, 176, 403, 910, 2128, 4896, 11628, 27132, 65208, 153824, 373175, 888030, 2170740, 5202600, 12797265, 30853680, 76292736, 184863168, 459162452, 1117370696, 2786017120, 6804995008, 17024247304, 41717833740, 104673837384 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 COMMENTS Substitutions 1->{2}, 2->{3,1}, 3->{4,2}, 4->{5,3,1}, 5->{6,4,2}, 6->{7,5,3,1}, 7->{8,6,4,2}, etc. The function f[n] gives Det[IdentityMatrix[n]-x*A[n]] with A[n]=Table[If[j > i+1, 0, Mod[i+j,2]], {i,n}, {j,n}] and can be written in terms of Dickson polynomials as : g(w)= x D_(w-1)(1+x, x*(1+x)) +(1-2*x)*E_(w-1)(1+x, x*(1+x)) Francisco Salinas (franciscodesalinas(AT)hotmail.com), Apr 13 2004; Count of integers is A047749. Sum of integers with substitution starting from 0 is A084081. LINKS FORMULA GF up to n-th term = GF[n]= g[n]/f[n] with f[1]=1; f[2]=1-x^2; f[3]=1-2x^2; f[n]=f[n-1]-x^2 f[n-3] and g[1]=1; g[2]=1+2x; g[3]=1+2x+2x^2; g[n]=g[n-1]-x^2 g[n-3]+2 x^(n-1) a(2n) = 4*binomial(3n,n-1)/(n+1) = 2*A006629(n-1); a(2n+1) = 6*binomial(3n+2,n)/(2n+3) - binomial(3n+1,n)/(n+1) = A056096(n+3). - Paul D. Hanna, Apr 24 2006 EXAMPLE GF[12]=(1 +2*x -7*x^2 -14*x^3 +9*x^4 +20*x^5 +2*x^6 -2*x^7 +2*x^11)/(1 -11*x^2 +36*x^4 -35*x^6 +5*x^8) produces a[1] to a[12]. a(4)=8 since 4-1= 3 substitutions on 1 produce 1-> 2-> 3+1-> 4+2+2 =8. MATHEMATICA Plus@@@Flatten/@NestList[ #/.k_Integer:>Range[k+1, 1, -2]&, {1}, 8]; (*or for n>16 *); f[1]=1; f[2]=1-x^2; f[3]=1-2x^2; f[n_]:=f[n]=Expand[f[n-1]-x^2 f[n-3]]; g[1]=1; g[2]=1+2x; g[3]=1+2x+2x^2; g[n_]:=g[n]=Expand[g[n-1] -x^2 g[n-3]+2 x^(n-1)]; GF[n_]:=g[n]/f[n]; CoefficientList[Series[GF[36], {x, 0, 36}], x] PROG (PARI) {a(n)=if(n%2==0, 4*binomial(3*n/2, n/2-1)/(n/2+1), 6*binomial(3*(n\2)+2, n\2)/(2*(n\2)+3) - binomial(3*(n\2)+1, n\2)/(n\2+1))} - Paul D. Hanna, Apr 24 2006 CROSSREFS Cf. A084081, A047749. Cf. A006629, A056096. Sequence in context: A072925 A002955 A202844 * A137255 A247298 A325928 Adjacent sequences:  A093948 A093949 A093950 * A093952 A093953 A093954 KEYWORD nonn AUTHOR Wouter Meeussen, Apr 18 2004 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified August 22 02:44 EDT 2019. Contains 326169 sequences. (Running on oeis4.)