%I #11 Oct 21 2023 16:01:54
%S 3,3,4,8,31,133,832,5113,41044,388800,3958704,39916802,518682390,
%T 6302045232,90968651712,1332614649600,22844265373440,356226551466344,
%U 7504470340300800,123358411682195904,2432902126073962432,52279222588118377280,1175121515279802150144
%N a(n) is the sum of the divisors of n!+1.
%F a(n) = sigma(n!+1) = A000203(A038507(n)).
%e a(5) = 133 because the divisors of 5!+1 are {1, 11, 121}.
%p a:=n->numtheory[sigma](n!+1):
%p seq(a(n), n=0..30);
%t DivisorSigma[1,Range[0,25]!+1] (* _Paolo Xausa_, Oct 21 2023 *)
%o (Python)
%o from math import factorial
%o from sympy import divisor_sigma
%o def A366758(n): return divisor_sigma(factorial(n)+1) # _Chai Wah Wu_, Oct 20 2023
%Y Cf. A038507, A000203, A064144, A062569, A366757.
%K nonn
%O 0,1
%A _Sean A. Irvine_, Oct 20 2023
|