%I #17 Feb 08 2015 14:20:16
%S 1,1,3,11,127,16151,260855055,68045359719085327,
%T 4630170979299719971778494028407039,
%U 21438483297549327871400796194793048411084076762817293736211302918175
%N a(1) = a(2) = 1, a(n+2) = 2*a(n) + a(n+1)^2.
%C In this sequence the primes begin a(3) = 3, a(4) = 11, a(5) = 127, a(9) = 4630170979299719971778494028407039.
%H <a href="/index/Aa#AHSL">Index entries for sequences of form a(n+1)=a(n)^2 + ...</a>
%F a(1) = a(2) = 1, for n>2: a(n) = 2*a(n-2) + a(n-1)^2. a(1) = a(2) = 1, for n>0: a(n+2) = 2*a(n) + a(n+1)^2.
%F a(n) ~ c^(2^n), where c = 1.163464453662702696843453679269882816346479873363677551158525103156732040997... . - _Vaclav Kotesovec_, Dec 18 2014
%e a(1) = 1 by definition.
%e a(2) = 1 by definition.
%e a(3) = 2*1 + 1^2 = 3.
%e a(4) = 2*1 + 3^2 = 11.
%e a(5) = 2*3 + 11^2 = 127.
%e a(6) = 2*11 + 127^2 = 16151.
%t Join[{a=1,b=1},Table[c=1*b^2+2*a;a=b;b=c,{n,10}]] (* _Vladimir Joseph Stephan Orlovsky_, Jan 18 2011 *)
%t RecurrenceTable[{a[1]==1, a[2]==1, a[n] == 2*a[n-2] + a[n-1]^2}, a, {n, 1, 10}] (* _Vaclav Kotesovec_, Dec 18 2014 *)
%Y Cf. A000278, A000283, A014253, A063827, A072878, A112957, A112958, A112959, A112960, A112961, A112969, A113785.
%K easy,nonn
%O 1,3
%A _Jonathan Vos Post_, Jan 24 2006