|
%I #37 Feb 05 2019 16:26:53
%S 1,3,7,37,121,1201,5041,62161,423361,5473441,39916801,818959681,
%T 6227020801,130784734081,1536517382401,32256486662401,355687428096001,
%U 10679532671808001,121645100408832001,3770998783116364801,59616236292028416001,1686001119824999577601
%N a(n) = n! * sum 1/k! where the sum is over all positive integers k that divide n.
%C Sets of lists of equal size, cf. A000262. - _Vladeta Jovovic_, Nov 02 2003
%C From _Gus Wiseman_, Jan 10 2019: (Start)
%C Number of matrices whose entries are 1,...,n, up to column permutations. For example, inequivalent representatives of the a(4) = 37 matrices are:
%C One 1 X 4 matrix:
%C [1234]
%C 12 2 X 2 matrices:
%C [12] [12] [13] [13] [14] [14] [23] [23] [24] [24] [34] [34]
%C [34] [43] [24] [42] [23] [32] [14] [41] [13] [31] [12] [21]
%C and 24 4 X 1 matrices:
%C [1][1][1][1][1][1][2][2][2][2][2][2][3][3][3][3][3][3][4][4][4][4][4][4]
%C [2][2][3][3][4][4][1][1][3][3][4][4][1][1][2][2][4][4][1][1][2][2][3][3]
%C [3][4][2][4][2][3][3][4][1][4][1][3][2][4][1][4][1][2][2][3][1][3][1][2]
%C [4][3][4][2][3][2][4][3][4][1][3][1][4][2][4][1][2][1][3][2][3][1][2][1]
%C in total 1+12+24 = 37.
%C (End)
%H Seiichi Manyama, <a href="/A057625/b057625.txt">Table of n, a(n) for n = 1..449</a>
%F E.g.f.: Sum_{n>0} (exp(x^n)-1). - _Vladeta Jovovic_, Dec 30 2001
%F E.g.f.: Sum_{k>0} x^k/k!/(1-x^k). - _Vladeta Jovovic_, Oct 14 2003
%F Equals the logarithmic derivative of A209903. - _Paul D. Hanna_, Jul 26 2012
%e a(4) = 4! (1 + 1/2! + 1/4!) = 24 (1 + 1/2 + 1/24) = 37.
%t a[n_] := n! DivisorSum[n, 1/#! &]; Array[a, 22] (* _Jean-François Alcover_, Dec 23 2015 *)
%o (PARI) a(n)=n! * sumdiv(n, d, 1/d! ); /* _Joerg Arndt_, Oct 07 2012 */
%Y Cf. A000005, A005225, A038041, A038048, A061095, A121860, A209903 (exp), A236696, A320444.
%Y Row sums of A066387.
%K nonn
%O 1,2
%A _Leroy Quet_, Oct 09 2000
| 