%I #16 Mar 17 2019 11:23:10
%S 2,4,16,112,448,3136,12544,87808,351232,2458624,9834496,68841472,
%T 275365888,1927561216,7710244864,53971714048,215886856192,
%U 1511207993344,6044831973376,42313823813632,169255295254528,1184787066781696,4739148267126784,33174037869887488
%N a(1)=2, a(n+1) = a(n)*A010888(a(n)).
%C From a(2) onwards, the digital root follows the pattern alternately 4,7,4,7,4,7,...
%H Colin Barker, <a href="/A110365/b110365.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (0,28).
%F a(1) = 2, a(2) = 4, a(3) = 16. a(2*n) = 4*a(2*n-1), a(2*n+1) = 7*a(2*n) for n > 1.
%F From _Colin Barker_, May 05 2016: (Start)
%F a(n) = 2^(-1+n)*(7^(1/2*(-3+n))*(2-2*(-1)^n + sqrt(7) + (-1)^n*sqrt(7))) for n > 1.
%F a(n) = 2^n*7^(n/2-1) for n > 1 and even.
%F a(n) = 2^(n+1)*7^((n-3)/2) for n > 1 and odd.
%F a(n) = 28*a(n-2) for n > 3.
%F G.f.: 2*x*(1+2*x-20*x^2) / (1-28*x^2).
%F (End)
%F E.g.f.: (-7 + 70*x + 7*cosh(2*Sqrt(7)*x) + 2*sqrt(7)*sinh(2*sqrt(7)*x))/49. - _Ilya Gutkovskiy_, May 05 2016
%t k = 2; Do[Print[k]; k *= Mod[Plus @@ IntegerDigits[k], 9], {n, 1, 30}] (* _Ryan Propper_, Oct 13 2005 *)
%t LinearRecurrence[{0,28},{2,4,16},30] (* _Harvey P. Dale_, Mar 17 2019 *)
%o (PARI) Vec(2*x*(1+2*x-20*x^2)/(1-28*x^2) + O(x^50)) \\ _Colin Barker_, May 05 2016
%Y Cf. A010888, A047892.
%K base,easy,nonn
%O 1,1
%A _Amarnath Murthy_, Jul 24 2005
%E More terms from _Ryan Propper_, Oct 13 2005
%E Name clarified by _Robert Israel_, May 05 2016