|
%I #33 Apr 11 2020 22:06:26
%S 1,37,3367,333667,33336667,3333366667,333333666667,33333336666667,
%T 3333333366666667,333333333666666667,33333333336666666667,
%U 3333333333366666666667,333333333333666666666667,33333333333336666666666667,3333333333333366666666666667,333333333333333666666666666667
%N a(n) = (10^(2*n)+10^n+1)/3.
%C Apart from the initial 1, common difference of the arithmetic progression pertaining to the sequence A074991.
%C This is also a root sequence pertaining to the patterned perfect square sequence 1369, 11336689,111333666889,... i.e., k ones, k threes and k sixes followed by (k-1) 8's and a 9. (37^2 = 1369).
%C This is a self-complementing sequence: each term has even number of digits (the first one has to be read 01, the leading zero is important). If you add the first half to the second half of any term, you get the sequence A011557, powers of 10. Furthermore, the reciprocals of the sequence terms, except the first one, give a sequence of periodic terms with period sequence as in A008585, a(n)=3n, and value given by A086574, a(n)=k where R(k+3)=3. - _Rodolfo A. Fiorini_, Jul 14 2016
%H Rodolfo A. Fiorini, <a href="http://dx.doi.org/10.3233/FI-2015-1267">Computerized tomography noise reduction by CICT optimized exponential cyclic sequences (OECS) co-domain</a>, Fundamenta Informaticae, 141(2015), 115-134.
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (111, -1110, 1000).
%F a(n) = 111*a(n-1) - 1110*a(n-2) + 1000*a(n-3), for n>2, a(0)=1, a(1)=37, a(2)=3367.
%F G.f.: (1 - 74*x + 370*x^2)/((1-x)(1-10*x)(1-100*x)). - _Colin Barker_, Sep 23 2013 and _Robert Israel_, Jul 22 2016
%p A074992 := proc(n)
%p (10^(2*n)+10^n+1)/3 ;
%p end proc:
%p seq(A074992(n),n=0..15) ; # _R. J. Mathar_, May 06 2017
%t {01}~Join~Table[FromDigits@ Flatten@ Map[IntegerDigits, {#, 10^n - #}] &@ Floor[10^n/3], {n, 12}] (* _Michael De Vlieger_, Jul 22 2016 *)
%o (PARI) a(n) = (10^(2*n) + 10^n + 1)/3; \\ _Michel Marcus_, Sep 14 2013
%o (PARI) Vec(-x*(1000*x^2-740*x+37)/((x-1)*(10*x-1)*(100*x-1)) + O(x^100)) \\ _Colin Barker_, Sep 23 2013
%o (PARI) a(n)=my(x=10^n); (x^2+x+1)/3 \\ _Charles R Greathouse IV_, Jul 22 2016
%Y Cf. A074991, A011557, A008585, A086574, A274766.
%K nonn,easy
%O 0,2
%A _Amarnath Murthy_, Aug 31 2002
%E Entry revised (new definition, new offset, new initial term, etc.) by _N. J. A. Sloane_, Jul 27 2016 (Some of the old programs may need slight modifications.)
|
