A110376 a(n) = Sum_{r < n, gcd(r,n)=1} n!/r.

1, 2, 9, 32, 250, 864, 12348, 67584, 804816, 5760000, 116915040, 686776320, 19323757440, 157991178240, 2951575200000, 42301705420800, 1202482800691200, 10048607738265600, 425162773111910400, 4541227794432000000

Offset: 1

Author

Amarnath Murthy, Jul 24 2005

Keywords

• easy
• nonn

Examples

			a(6) = 6!(1/1 + 1/5) = 864.
		

Crossrefs

Cf. A038048, A110373, A110374.

Programs

Maple

a:=proc(n) local s,r: s:=0: for r from 1 to n do if gcd(r,n)=1 then s:=s+1/r else s:=s: fi: od: n!*s end: seq(a(n),n=1..23); # Emeric Deutsch, Jul 25 2005

Mathematica

Do[Print[n! * Plus @@ Map[(1/#)&, Select[Range[n], GCD[ #, n] == 1 &]]], {n, 1, 30}] (* Ryan Propper, Jul 25 2005 *)

Extensions

More terms from Emeric Deutsch, Ryan Propper and Reinhard Zumkeller, Jul 25 2005