%I #19 Dec 12 2015 13:15:09
%S 4,8,14,23,37,57,87,133,201,304,457,688,1034,1553,2331,3499,5251,7878,
%T 11819,17731,26598,39899,59851,89778,134669,202005,303010,454517,
%U 681778,1022668,1534004,2301009,3451515,5177275,7765914
%N The nearest integer of perimeter of T-square (fractal) after n-iterations, starting with a unit square.
%C The round (nearest) function uses the convention that values halfway between two integers are rounded up, so a(4) = 37 = round(36 + 1/2).
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/T-square_(fractal)">T-square (fractal)</a>
%H Kival Ngaokrajang, <a href="/A227621/a227621.jpg">Illustration for n = 0..4</a>
%F a(0) = 4; for n >= 1, a(n) = round(A083313(n)/2^(n-3)).
%e The central unit square has perimeter 4. At n=1, the additional 4 squares hide 4*1/2 = 2 units of length and add 4*3*1/2 = 6 units, to give a(1) = 4 - 2 + 6 = 8.
%e At n=2, the additional 12 squares hide 12*1/(2*2) units of length and add 12*3*1/(2*2) to give a(2) = 8 - 12/4 + 36/4 = 14.
%p A227621 := proc(n)
%p if n = 0 then 4
%p else
%p round(A083313(n)/2^(n-3)) ;
%p end if;
%p end proc: # _R. J. Mathar_, Aug 01 2013
%o (Small Basic)
%o a[0]=4
%o For n = 0 To 51
%o a[n+1] = Math.Round((Math.Power(3,n+1) - Math.Power(2,n))/Math.Power(2,n-2))
%o TextWindow.Write(a[n]+", ")
%o EndFor
%Y Cf. A083313.
%K nonn
%O 0,1
%A _Kival Ngaokrajang_, Jul 17 2013
|