login
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.
0

%I #33 Jun 17 2024 15:54:37

%S 1,1,3,0,0,-28733079375,-821329806742140930609375,0,0

%N 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.

%C a(9) is roughly 10^250.72 and is too big for the data section.

%C 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).

%C 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.

%e 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.

%t 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 *)

%o (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]))))

%Y Cf. A365968 (row n gives factors of a(n)).

%Y Cf. A001787 (result if we would use summation only).

%K sign

%O 0,3

%A _Thomas Scheuerle_, May 28 2024