%I #24 Aug 13 2022 10:08:43
%S 4,32,3072,2752512,21139292160,1342091380654080,692647993190048071680,
%T 2882479558988139892026900480,96342151992701835341576224427212800,
%U 25811138467998276182105365247324712232550400
%N Order of unitary 2^n X 2^n group H_{n,4} acting on Siegel modular forms.
%D B. Runge, On Siegel modular forms I, J. Reine Angew. Math., 436 (1993), 57-85.
%H G. C. Greubel, <a href="/A027639/b027639.txt">Table of n, a(n) for n = 0..45</a>
%H B. Runge, <a href="https://doi.org/10.1016/0012-365X(94)00271-J">Codes and Siegel modular forms</a>, Discrete Math. 148 (1996), 175-204.
%H <a href="/index/Gre#groups_modular">Index entries for sequences related to modular groups</a>
%F a(n) = 2^(n^2 + 2*n + 2) * Product_{j=1..n} (2^j - 1).
%F a(n) = (-1)^n * 2^(n^2 + 2*n + 2) * (2, 2)_{n}, where (q, q)_{n} is the q-Pochhammer symbol. - _G. C. Greubel_, Aug 04 2022
%p seq(2^(n^2+2*n+2)*mul(2^i-1, i=1..n), n=0..10);
%t a[n_]:= (-1)^n*2^(n^2 +2*n+2)*QPochhammer[2,2,n];
%t Table[a[n], {n, 0, 15}] (* _G. C. Greubel_, Aug 04 2022 *)
%o (Magma)
%o A027639:= func< n | n eq 0 select 4 else 2^(n^2+2*n+2)*(&*[2^j-1: j in [1..n]]) >;
%o [A027639(n): n in [0..15]]; // _G. C. Greubel_, Aug 04 2022
%o (SageMath)
%o from sage.combinat.q_analogues import q_pochhammer
%o def A027639(n): return (-1)^n*2^(n^2 + 2*n + 2)*q_pochhammer(n, 2, 2)
%o [A027639(n) for n in (0..15)] # _G. C. Greubel_, Aug 04 2022
%o (PARI) a(n) = my(ret=1); for(i=1,n, ret = ret<<i-ret); ret << (n^2 + 2*n + 2); \\ _Kevin Ryde_, Aug 13 2022
%Y Cf. A005329.
%K nonn,nice,easy
%O 0,1
%A _N. J. A. Sloane_