|
%I #70 Sep 08 2022 08:46:16
%S 0,1,2,6,17,50,147,434,1282,3789,11200,33109,97878,289354,855413,
%T 2528850,7476023,22101326,65338038,193158521,571033600,1688143881,
%U 4990651642,14753839486,43616704857,128943855250,381196100507,1126928202714,3331532438042,9848993360069
%N Self-composition of the Fibonacci sequence.
%C This sequence has the same relation to the Fibonacci numbers A000045 as A030267 has to the natural numbers A000027.
%C From _Oboifeng Dira_, Jun 28 2020: (Start)
%C This sequence can be generated from a family of composition pairs of generating functions g(f(x)), where k is an integer and where
%C f(x) = x/(1-k*x-x^2) and g(x) = (x+(k-1)*x^2)/(1-(3-2*k)*x-(3*k-k^2-1)*x^2).
%C Some cases of k values are:
%C k=-5, f(x) g.f. 0,A052918(-1)^n and g(x) g.f. 0,A081571
%C k=-4, f(x) g.f. A001076(-1)^(n+1) and g(x) g.f. 0,A081570
%C k=-3, f(x) g.f. A006190(-1)^(n+1) and g(x) g.f. 0,A081569
%C k=-2, f(x) g.f. A215936(n+2) and g(x) g.f. 0,A081568
%C k=-1, f(x) g.f. A039834(n+2) and g(x) g.f. 0,A081567
%C k=0, f(x) g.f. A000035 and g(x) g.f. 0,A001519(n+1)
%C k=1, f(x) g.f. A000045 and g(x) g.f. A000045
%C k=2, f(x) g.f. A000129 and g(x) g.f. 0,A039834(n+1)
%C k=3, f(x) g.f. A006190 and g(x) g.f. 0,A001519(-1)^n
%C k=4, f(x) g.f. A001076 and g(x) g.f. 0,A093129(-1)^n
%C k=5, f(x) g.f. 0,A052918 and g(x) g.f. 0,A192240(-1)^n
%C k=6, f(x) g.f. A005668 and g(x)=(x+5*x^2)/(1+9*x+19*x^2)
%C k=7, f(x) g.f. 0,A054413 and g(x)=(x+6*x^2)/(1+11*x+29*x^2).
%C (End)
%H Colin Barker, <a href="/A270863/b270863.txt">Table of n, a(n) for n = 0..1000</a>
%H Oboifeng Dira, <a href="http://www.seams-bull-math.ynu.edu.cn/downloadfile.jsp?filemenu=_201706&filename=07_41(6).pdf">A Note on Composition and Recursion</a>, Southeast Asian Bulletin of Mathematics (2017), Vol. 41, Issue 6, 849-853.
%H Oboifeng Dira, <a href="/A270863/a270863_1.pdf">Family of composition pairs g(f(x)) generating A270683</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (3,1,-3,-1).
%F a(n) = 3*a(n-1)+a(n-2)-3*a(n-3)-a(n-4) for n > 3, a(0)=0, a(1)=1, a(2)=2, a(3)=6.
%F G.f.: x*(1-x-x^2) / (1-3*x-x^2+3*x^3+x^4). - _Colin Barker_, Mar 24 2016
%F G.f.: B(B(x)) where B(x) is the g.f. of A000045. - _Joerg Arndt_, Mar 25 2016
%F a(n) = (phi*((phi^2 + 5^(1/4)*sqrt(3*phi))^n - (phi^2 - 5^(1/4)*sqrt(3*phi))^n) + (psi^2 + 5^(1/4)*sqrt(3*psi))^n - (psi^2 - 5^(1/4)*sqrt(3*psi))^n)/(2^n * 5^(3/4) * sqrt(3*phi)), where phi = (sqrt(5) + 1)/2 is the golden ratio, and psi = 1/phi = (sqrt(5) - 1)/2. - _Vladimir Reshetnikov_, Aug 01 2019
%F 0 = a(n)*(a(n) +6*a(n+1) -a(n+2)) +a(n+1)*(8*a(n+1) -9*a(n+2) +a(n+3)) +a(n+2)*(-8*a(n+2) +6*a(n+3)) +a(n+3)*(-a(n+3)) if n>=0. - _Michael Somos_, Feb 05 2022
%e a(5) = 3*a(4)+a(3)-3*a(2)-a(1) = 51+6-6-1 = 50.
%p f:= x-> x/(1-x-x^2):
%p a:= n-> coeff(series(f(f(x)), x, n+1), x, n):
%p seq(a(n), n=0..30);
%o (PARI) a(n)=([0,1,0,0; 0,0,1,0; 0,0,0,1; -1,-3,1,3]^(n-1)*[1;2;6;17])[1,1] \\ _Charles R Greathouse IV_, Mar 24 2016
%o (PARI) concat(0, Vec(x*(1-x-x^2)/(1-3*x-x^2+3*x^3+x^4) + O(x^40))) \\ _Colin Barker_, Mar 24 2016
%o (Magma) I:=[0, 1, 2, 6]; [m le 4 select I[m] else 3*Self(m-1)+Self(m-2)-3*Self(m-3)-Self(m-4): m in [1..30]]; // _Marius A. Burtea_, Aug 03 2019
%Y Cf. A000027, A000045, A001622, A030267, A000035, A001519, A000129, A039834.
%K nonn,easy
%O 0,3
%A _Oboifeng Dira_, Mar 24 2016
|
