%I #43 Mar 11 2018 09:53:48
%S 1,1,1,2,1,1,2,2,2,2,4,3,3,3,3,6,5,5,6,5,5,10,9,9,9,9,9,9,16,16,16,16,
%T 16,16,16,16,30,28,28,29,28,28,29,28,28,52,51,51,51,51,52,51,51,51,51,
%U 94,93,93,93,93,93,93,93,93,93,93,172,170,170,172,170,170,172
%N Triangular array T(n,k) giving coefficients in expansion of Product_{j=1..n} (1+x^j) mod x^(n+1)-1.
%C T(n,k) = number of binary vectors (x_1,...,x_n) satisfying Sum_{i=1..n} i*x_i = k (mod n+1) = size of Varshamov-Tenengolts code VT_k(n).
%D B. D. Ginsburg, On a number theory function applicable in coding theory, Problemy Kibernetiki, No. 19 (1967), pp. 249-252.
%H Seiichi Manyama, <a href="/A053633/b053633.txt">Rows n = 0..139, flattened</a>
%H F. J. van de Bult, D. C. Gijswijt, J. P. Linderman, N. J. A. Sloane and Allan Wilks, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL10/Sloane/sloane55.html">A Slow-Growing Sequence Defined by an Unusual Recurrence</a>, J. Integer Sequences, Vol. 10 (2007), #07.1.2.
%H F. J. van de Bult, D. C. Gijswijt, J. P. Linderman, N. J. A. Sloane and Allan Wilks, A Slow-Growing Sequence Defined by an Unusual Recurrence [<a href="http://neilsloane.com/doc/gijs.pdf">pdf</a>, <a href="http://neilsloane.com/doc/gijs.ps">ps</a>].
%H N. J. A. Sloane, <a href="http://neilsloane.com/doc/dijen.txt">On single-deletion-correcting codes</a>
%H <a href="/index/Su#subsetsums">Index entries for sequences related to subset sums modulo m</a>
%H <a href="/index/Ge#Gijswijt">Index entries for sequences related to Gijswijt's sequence</a>
%F The Maple code gives an explicit formula.
%e Triangle begins:
%e k 0 1 2 3 4 5 6 7 8 9
%e n
%e 0 1;
%e 1 1, 1;
%e 2 2, 1, 1;
%e 3 2, 2, 2, 2;
%e 4 4, 3, 3, 3, 3;
%e 5 6, 5, 5, 6, 5, 5;
%e 6 10, 9, 9, 9, 9, 9, 9;
%e 7 16, 16, 16, 16, 16, 16, 16, 16;
%e 8 30, 28, 28, 29, 28, 28, 29, 28, 28;
%e 9 52, 51, 51, 51, 51, 52, 51, 51, 51, 51;
%e ...
%e [Edited by _Seiichi Manyama_, Mar 11 2018]
%p with(numtheory): A053633 := proc(n,k) local t1,d; t1 := 0; for d from 1 to n do if n mod d = 0 and d mod 2 = 1 then t1 := t1+(1/(2*n))*2^(n/d)*phi(d)*mobius(d/gcd(d,k))/phi(d/gcd(d,k)); fi; od; t1; end;
%t Flatten[ Table[ CoefficientList[ PolynomialMod[ Product[1+x^j, {j,1,n}], x^(n+1)-1], x], {n,0,11}]] (* _Jean-François Alcover_, May 04 2011 *)
%Y Cf. A053632, A063776, A300328, A300628. Leading coefficients give A000016, next column gives A000048.
%K tabl,nonn,easy,nice
%O 0,4
%A _N. J. A. Sloane_, Mar 22 2000