%I #9 Apr 09 2014 10:13:40
%S 0,1,5,2,34,324,2988,28944,300816,3371040,793440,49438080,1121696640,
%T 21930773760,416139897600,7965912729600,156343306291200,
%U 3169866941337600,66629845591142400,1454242012231680000
%N n! * (fractional part of n-th harmonic number).
%F a(n) = n! * fractional part ( harmonic number (n)).
%F a(n) also equals n!*(sum{k=1 to n} 1/k) (mod n!). - _Leroy Quet_, Mar 20 2008
%t f[n_] := n!*FractionalPart[HarmonicNumber[n]]; Table[ f[n], {n, 22}]
%K nonn
%O 1,3
%A _Robert G. Wilson v_, Jul 19 2004, based on a suggestion from _Leroy Quet_, Jun 18 2004