%I #11 Jul 25 2024 08:34:24

%S 1,3,17,120,885,6713,51932,407214,3224845,25733325,206584437,

%T 1666561042,13498994796,109713432390,894291885000,7307812510970,

%U 59847327807597,491062976039618,4036174402666925,33224883837921930,273873806179142545,2260338391869532332

%N a(n) is the (2n)-th term of the n-fold self-convolution of the sum of divisors function sigma.

%H Alois P. Heinz, <a href="/A340993/b340993.txt">Table of n, a(n) for n = 0..1080</a>

%F a(n) = [x^(2n)] (Sum_{j>=1} sigma(j)*x^j)^n.

%F a(n) = A319083(2n,n).

%F a(n) ~ c * d^n / sqrt(n), where d = 8.455610430383829836198938524980234226695900064615457328971640722426861925... and c = 0.352126317954512592610958969393229871240824031029408304023118123356... - _Vaclav Kotesovec_, Jul 25 2024

%p b:= proc(n, k) option remember; `if`(k=0, 1,

%p `if`(k=1, numtheory[sigma](n+1), (q->

%p add(b(j, q)*b(n-j, k-q), j=0..n))(iquo(k, 2))))

%p end:

%p a:= n-> b(n$2):

%p seq(a(n), n=0..23);

%Y Cf. A000203, A319083.

%K nonn

%O 0,2

%A _Alois P. Heinz_, Feb 01 2021