cp's OEIS Frontend

A158786 Irregular triangle T(n, k) = A000032(n-2*k+1) if (n-2*k) mod 2 = 0, otherwise 25*A000032(n-2*k), read by rows.

%I A158786 #11 Apr 25 2024 09:31:16
%S A158786 1,25,4,1,100,25,11,4,1,275,100,25,29,11,4,1,725,275,100,25,76,29,11,
%T A158786 4,1,1900,725,275,100,25,199,76,29,11,4,1,4975,1900,725,275,100,25,
%U A158786 521,199,76,29,11,4,1,13025,4975,1900,725,275,100,25,1364,521,199,76,29,11,4,1,34100,13025,4975,1900,725,275,100,25
%N A158786 Irregular triangle T(n, k) = A000032(n-2*k+1) if (n-2*k) mod 2 = 0, otherwise 25*A000032(n-2*k), read by rows.
%D A158786 H. S. M. Coxeter, Regular Polytopes, 3rd ed., Dover, NY, 1973, pp. 159-162.
%H A158786 G. C. Greubel, <a href="/A158786/b158786.txt">Rows n = 0..50 of the irregular triangle, flattened</a>
%F A158786 T(n, k) = 5*e(n, k), where e(n,k) = (e(n-1, k)*e(n, k-1) + 1)/e(n-1, k-1), and e(n, 0) = sqrt(5)*(GoldenRatio^(n) + GoldenRatio^(-n)).
%F A158786 From _G. C. Greubel_, Dec 06 2021: (Start)
%F A158786 T(n, k) = A000032(n-2*k+1) if (n-2*k) mod 2 = 0, otherwise 25*A000032(n-2*k).
%F A158786 Sum_{k=0..floor(n/2)} T(n, k) = A000032(n) - 2 if (n mod 2 = 0), otherwise 25*(A000032(n-1) - 2). (End)
%e A158786 Irregular triangle begins as:
%e A158786      1;
%e A158786     25;
%e A158786      4,    1;
%e A158786    100,   25;
%e A158786     11,    4,   1;
%e A158786    275,  100,  25;
%e A158786     29,   11,   4,   1;
%e A158786    725,  275, 100,  25;
%e A158786     76,   29,  11,   4,  1;
%e A158786   1900,  725, 275, 100, 25;
%e A158786    199,   76,  29,  11,  4,   1;
%e A158786   4975, 1900, 725, 275, 100, 25;
%t A158786 (* First program *)
%t A158786 e[n_, 0]:= Sqrt[5]*(GoldenRatio^(n) + GoldenRatio^(-n));
%t A158786 e[n_, k_]:= If[k>n-1, 0, (e[n-1, k]*e[n, k-1] +1)/e[n-1, k-1]];
%t A158786 T[n_,k_]:= 5*Rationalize[N[e[n, k]]];
%t A158786 Table[T[n, k], {n, 2, 16}, {k, Mod[n, 2] +1, n-1,2}]//Flatten
%t A158786 (* Second program *)
%t A158786 f[n_]:= f[n]= If[EvenQ[n], LucasL[n-1], 25*LucasL[n-2]];
%t A158786 T[n_, k_]:= f[n-2*k];
%t A158786 Table[T[n, k], {n, 2, 16}, {k, 0, (n-2)/2}]//Flatten (* _G. C. Greubel_, Dec 06 2021 *)
%o A158786 (Sage)
%o A158786 def A158786(n,k): return lucas_number2(n-2*k-1,1,-1) if ((n-2*k)%2==0) else 25*lucas_number2(n-2*k-2,1,-1)
%o A158786 flatten([[A158786(n,k) for k in (0..((n-2)//2))] for n in (2..16)]) # _G. C. Greubel_, Dec 06 2021
%Y A158786 Cf. A000032, A002878, A004146, A158753.
%K A158786 nonn,tabf
%O A158786 2,2
%A A158786 _Roger L. Bagula_ and _Gary W. Adamson_, Mar 26 2009
%E A158786 Edited by _G. C. Greubel_, Dec 06 2021