|
%I #10 Sep 07 2023 07:14:35
%S 7,247,21756,3613701,974243088,388409565699,214946329538304,
%T 157727064375306153,148245464311769260800,173696139110375108022159,
%U 248243987235370949531025408,425095516929076538387157860013
%N v(n+1)/v(n), where v=A203012.
%C See A093883 for a discussion and guide to related sequences.
%F a(n) ~ 3^(3*n/2 + 1) * exp((n+1)*Pi/(2*sqrt(3)) - 2*n) * n^(2*n). - _Vaclav Kotesovec_, Sep 07 2023
%t f[j_] := j; z = 12;
%t v[n_] := Product[Product[f[j]^2 + f[j] f[k] + f[k]^2,
%t {j, 1, k - 1}], {k, 2, n}]
%t Table[v[n], {n, 1, z}] (* A203012 *)
%t Table[v[n + 1]/v[n], {n, 1, z}] (* A203158 *)
%t Table[Product[k^2 + k*(n+1) + (n+1)^2, {k, 1, n}], {n, 1, 15}] (* _Vaclav Kotesovec_, Sep 07 2023 *)
%K nonn
%O 1,1
%A _Clark Kimberling_, Jan 04 2012
|
