|
%I #25 Aug 05 2018 08:19:57
%S 1,1,1,0,0,1,1,1,0,1,0,0,0,0,1,0,0,1,0,0,1,0,0,0,0,0,0,1,1,1,0,1,0,0,
%T 0,1,0,0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,1,0,0,
%U 1,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,1
%N Inverse of matrix (1,x)-(x,x^2) (expressed in Riordan array notation).
%C Row sums are the 'ruler function' A001511. Columns are stretched Fredholm-Rueppel sequences (A036987). Inverse is A115359.
%C Eigensequence of triangle A115361 = A018819 starting with offset 1: (1, 2, 2, 4, 4, 6, 6, 10, 10, 14, 14, 20, 20, ...). - _Gary W. Adamson_, Nov 21 2009
%C From _Gary W. Adamson_, Nov 27 2009: (Start)
%C A115361 * [1, 2, 3, ...] = A129527 = (1, 3, 3, 7, 5, 9, 7, 15, ...).
%C (A115361)^(-1) * [1, 2, 3, ...] = A115359 * [1, 2, 3, ...] = A026741 starting /Q (1, 1, 3, 2, 5, 3, 7, 4, 9, ...). (End)
%C This is the lower-left triangular part of the inverse of the infinite matrix A_{ij} = [i=j] - [i=2j], its upper-right part (above / right to the diagonal) being zero. The n-th row has 1 in column n/2^i, i = 0, 1, ... as long as this is an integer. - _M. F. Hasler_, May 13 2018
%F Number triangle whose k-th column has g.f. x^k*sum{j>=0} x^((2^j-1)*(k+1)).
%F T(n,k) = A209229((n+1)/(k+1)) for k+1 divides n+1, T(n,k) = 0 otherwise. - _Andrew Howroyd_, Aug 05 2018
%e Triangle begins:
%e 1;
%e 1,1;
%e 0,0,1;
%e 1,1,0,1;
%e 0,0,0,0,1;
%e 0,0,1,0,0,1;
%e 0,0,0,0,0,0,1;
%e 1,1,0,1,0,0,0,1;
%e 0,0,0,0,0,0,0,0,1;
%e 0,0,0,0,1,0,0,0,0,1;
%e 0,0,0,0,0,0,0,0,0,0,1;
%p A115361 := proc(n,k)
%p for j from 0 do
%p if k+(2*j-1)*(k+1) > n then
%p return 0 ;
%p elif k+(2^j-1)*(k+1) = n then
%p return 1 ;
%p end if;
%p end do;
%p end proc: # _R. J. Mathar_, Jul 14 2012
%t (*recurrence*)
%t Clear[t]
%t t[1, 1] = 1;
%t t[n_, k_] :=
%t t[n, k] =
%t If[k == 1, Sum[t[n, k + i], {i, 1, 2 - 1}],
%t If[Mod[n, k] == 0, t[n/k, 1], 0], 0]
%t Flatten[Table[Table[t[n, k], {k, 1, n}], {n, 14}]] (* _Mats Granvik_, Jun 26 2014 *)
%o (PARI) tabl(nn) = {T = matrix(nn, nn, n, k, n--; k--; if ((n==k), 1, if (n==2*k+1, -1, 0))); Ti = T^(-1); for (n=1, nn, for (k=1, n, print1(Ti[n, k], ", ");); print(););} \\ _Michel Marcus_, Mar 28 2015
%o (PARI) A115361_row(n,v=vector(n))={until(bittest(n,0)||!n\=2,v[n]=1);v} \\ Yields the n-th row (of length n). - _M. F. Hasler_, May 13 2018
%o (PARI) T(n,k)={if(n%k, 0, my(e=valuation(n/k,2)); n/k==1<<e)}
%o for(n=1, 10, for(k=1, n, print1(T(n,k), ", ")); print) \\ _Andrew Howroyd_, Aug 03 2018
%Y Cf. A016741, A018819, A129527. - _Gary W. Adamson_, Nov 21 2009
%Y Cf. A036987, A209229.
%K easy,nonn,tabl
%O 0,1
%A _Paul Barry_, Jan 21 2006
|
