cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A373202 Product over all sums of {n, +/-(n-1), +/-(n-2), ..., +/-1}, where "+/-" means here we use all possible combination of signs, this means for n > 0 that we have 2^(n-1) factors.

Original entry on oeis.org

1, 1, 3, 0, 0, -28733079375, -821329806742140930609375, 0, 0
Offset: 0

Views

Author

Thomas Scheuerle, May 28 2024

Keywords

Comments

a(9) is roughly 10^250.72 and is too big for the data section.
Consider the polynomial recurrence: P_{n}(x) = P_{n-1}(x + n)*P_{n-1}(x - n), with P_{1}(x) = 1 + x and P_{0}(x) = 1, then |a(n)| = P_{n}(0).
Each nonzero term divides all nonzero terms with higher index. This can be understood by looking at the tree structure in A365968 as the factors of a(n) are obtained in row n (second half of row). Each factor is only dependent by inheritance from parent and grandparent in this tree. Each row contains either all even numbers from 0 to A000217(n) or all odd numbers from 0 to A000217(n). In the even case 0 will be a factor and thus a(n) = 0. In the odd case some factors are duplicates, but by inheritance we will keep the duplicates from parent rows included.

Examples

			a(4) = (4+3+2+1)*(4+3+2-1)*(4+3-2+1)*(4+3-2-1)*(4-3+2+1)*(4-3+2-1)*(4-3-2+1)*(4-3-2-1) = 10*8*6*4*4*2*0*(-2) = 0.
		

Crossrefs

Cf. A365968 (row n gives factors of a(n)).
Cf. A001787 (result if we would use summation only).

Programs

  • Mathematica
    a[n_]:=Product[Sum[(n+1-m)(2*Part[IntegerDigits[2^(n-1)+k,2],m]-1),{m,n}],{k,0,2^(n-1)-1}]; Array[a,10,0] (* Stefano Spezia, May 29 2024 *)
  • PARI
    a(n) = prod(k=0, 2^(n-1)-1, sum(m=1, n, (n+1-m)*(-1+2*(digits(2^(n-1)+k, 2)[m]))))