cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A181864 a(1) = 1, a(2) = 2. For n >= 3, a(n) is found by concatenating the squares of the first n-1 terms of the sequence and then dividing the resulting number by a(n-1).

This page as a plain text file.
%I A181864 #7 Mar 30 2012 18:40:13
%S A181864 1,2,7,207,700207,207000000700207,70020700000000000000207000000700207,
%T A181864 2070000007002070000000000000000000000000000000000070020700000000000000207000000700207
%N A181864 a(1) = 1, a(2) = 2. For n >= 3, a(n) is found by concatenating the squares of the first n-1 terms of the sequence and then dividing the resulting number by a(n-1).
%C A181864 The calculations for the first few values of the sequence are
%C A181864 ... 2^2 = 4 so a(3) = 14/2 = 7
%C A181864 ... 7^2 = 49 so a(4) = 1449/7 = 207
%C A181864 ... 207^2 = 42849 so a(5) = 144942849/207 = 700207.
%C A181864 For similarly defined sequences see A181754 through A181756 and A181865 through
%C A181864  A181870.
%F A181864 DEFINITION
%F A181864 a(1) = 1, a(2) = 2, and for n >= 3
%F A181864 (1)... a(n) = concatenate(a(1)^2,a(2)^2,...,a(n-1)^2)/a(n-1).
%F A181864 RECURRENCE RELATION
%F A181864 For n >= 2
%F A181864 (2)...a(n+2) = a(n+1) + 10^F(n,2)*a(n) = a(n+1) + 10^Pell(n)*a(n),
%F A181864 where F(n,2) is the Fibonacci polynomial F(n,x) evaluated at x = 2
%F A181864 and where Pell(n) = A000129(n).
%F A181864 RELATION WITH OTHER SEQUENCES
%F A181864 a(n) has A113225(n-2) digits.
%F A181864 a(n)^2 has Pell(n-1) digits.
%p A181864 #A181864
%p A181864 M:=8: a:=array(1..M):s:=array(1..M):
%p A181864 a[1]:=1:a[2]:=2:
%p A181864 s[1]:=convert(a[1]^2,string):
%p A181864 s[2]:=cat(s[1],convert(a[2]^2,string)):
%p A181864 for n from 3 to M do
%p A181864 a[n] := parse(s[n-1])/a[n-1];
%p A181864 s[n]:= cat(s[n-1],convert(a[n]^2,string));
%p A181864 end do:
%p A181864 seq(a[n],n = 1..M);
%Y A181864 A000129, A113225, A181754, A181755, A181756, A181865, A181866,A181867, A181868, A181869, A181870
%K A181864 nonn,easy
%O A181864 1,2
%A A181864 _Peter Bala_, Nov 28 2010