login
Array read by antidiagonals with T(n,k) = min(n,k).
43

%I #54 Jan 29 2022 09:12:19

%S 1,1,1,1,2,1,1,2,2,1,1,2,3,2,1,1,2,3,3,2,1,1,2,3,4,3,2,1,1,2,3,4,4,3,

%T 2,1,1,2,3,4,5,4,3,2,1,1,2,3,4,5,5,4,3,2,1,1,2,3,4,5,6,5,4,3,2,1,1,2,

%U 3,4,5,6,6,5,4,3,2,1,1,2,3,4,5,6,7,6,5,4,3,2,1,1,2,3,4,5,6,7,7,6,5,4,3,2,1

%N Array read by antidiagonals with T(n,k) = min(n,k).

%C Also, "correlation triangle" for the constant sequence 1. - _Paul Barry_, Jan 16 2006

%C Antidiagonal sums are in A002620.

%C As a triangle, row sums are A002620. T(2n,n)=n+1. Diagonal sums are A001399. Construction: Take antidiagonal triangle of MM^T where M is the sequence array for the constant sequence 1 (lower triangular matrix with all 1's). - _Paul Barry_, Jan 16 2006

%C From _Franklin T. Adams-Watters_, Sep 25 2011: (Start)

%C As a triangle, count up to ceiling(n/2) and back down again (repeating the central term when n is even).

%C When the first two instances of each number are removed from the sequence, the original sequence is recovered.

%C (End)

%H Reinhard Zumkeller, <a href="/A003983/b003983.txt">Rows n=1..100 of triangle, flattened</a>

%F Number triangle T(n, k) = Sum_{j=0..n} [j<=k][j<=n-k]. - _Paul Barry_, Jan 16 2006

%F G.f.: 1/((1-x)*(1-x*y)*(1-x^2*y)). - _Christian G. Bower_, Jan 17 2006

%F a(n) = min(floor( 1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2). - _Leonid Bedratyuk_, Dec 13 2009

%e Triangle version begins

%e 1;

%e 1, 1;

%e 1, 2, 1;

%e 1, 2, 2, 1;

%e 1, 2, 3, 2, 1;

%e 1, 2, 3, 3, 2, 1;

%e 1, 2, 3, 4, 3, 2, 1;

%e 1, 2, 3, 4, 4, 3, 2, 1;

%e 1, 2, 3, 4, 5, 4, 3, 2, 1;

%e ...

%p a(n) = min(floor(1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2) # _Leonid Bedratyuk_, Dec 13 2009

%t Flatten[Table[Min[n-k+1, k], {n, 1, 14}, {k, 1, n}]] (* _Jean-François Alcover_, Feb 23 2012 *)

%o (Haskell)

%o a003983 n k = a003983_tabl !! (n-1) !! (k-1)

%o a003983_tabl = map a003983_row [1..]

%o a003983_row n = hs ++ drop m (reverse hs)

%o where hs = [1..n' + m]

%o (n',m) = divMod n 2

%o -- _Reinhard Zumkeller_, Aug 14 2011

%o (PARI) T(n,k) = min(n,k) \\ _Charles R Greathouse IV_, Feb 06 2017

%Y Cf. A002620, A001399, A087062, A115236, A115237, A124258, A006752, A120268, A173945, A173947, A173948, A173949, A173953, A173954, A173955, A173973, A173982-A173986, A004197.

%K tabl,nonn,easy,nice

%O 1,5

%A _Marc LeBrun_

%E More terms from Larry Reeves (larryr(AT)acm.org), Nov 08 2000

%E Entry revised by _N. J. A. Sloane_, Dec 05 2006