|
%I #15 Mar 21 2024 08:31:39
%S 1,3,3,4,8,4,6,11,11,6,8,16,15,16,8,9,21,22,22,21,9,11,24,29,32,29,24,
%T 11,12,29,33,42,42,33,29,12,14,32,40,48,55,48,40,32,14,16,37,44,58,63,
%U 63,58,44,37,16,17,42,51,64,76,72,76,64,51,42,17
%N Array T(n, k) is the number x such that frac(x*phi) + frac(n*phi)*frac(k*phi) = 1 where phi is the golden ratio A001622 and frac(y) is the fractional part of y, read by antidiagonals.
%H Paolo Xausa, <a href="/A340429/b340429.txt">Table of n, a(n) for n = 1..11325</a> (first 150 antidiagonals, flattened).
%F T(n, k) = 2*n*k + A339765(n)*k/2 + A339765(k)*n/2.
%F T(n, k) = T(k, n), array is symmetric.
%F T(n, k) = 3*n*k - n*h(k) - k*h(n) where h(n) = ceiling(2*n / (sqrt(5) + 3)) = A189663(n + 1). - _Peter Luschny_, Mar 21 2024
%e Array begins:
%e 1 3 4 6 8 ...
%e 3 8 11 16 21 ...
%e 4 11 15 22 29 ...
%e 6 16 22 32 42 ...
%e 8 21 29 42 55 ...
%e ...
%p h := n -> ceil(2*n / (sqrt(5) + 3)):
%p T := (n, k) -> 3*n*k - n*h(k) - k*h(n):
%p seq(lprint(seq(T(n, k), k = 1..9)), n = 1..7); # _Peter Luschny_, Mar 21 2024
%t A340429[n_, k_] := Floor[n * GoldenRatio] * k + Floor[k * GoldenRatio] * n - n * k;
%t Table[A340429[n - k + 1, k], {n, 15}, {k, n}] (* _Paolo Xausa_, Mar 21 2024 *)
%o (PARI) f(n) = 2*floor(n*(1+sqrt(5))/2) - 3*n; \\ A339765
%o T(n, k) = 2*n*k + f(n)*k/2 + f(k)*n/2;
%Y Cf. A001622, A101330, A189663, A339765.
%Y Cf. A000201 (row 1), A003623 (row 2), A190509 (row 3), A371388 (main diagonal).
%K nonn,tabl
%O 1,2
%A _Michel Marcus_, Jan 07 2021
|
