login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A335088 a(n) = n! * [x^n] exp(Sum_{k=1..n, gcd(n,k) = 1} x^k / k). 3
1, 1, 1, 4, 9, 96, 145, 4320, 11025, 179200, 693441, 36288000, 62610625, 5748019200, 15259154625, 378002200576, 4108830350625, 334764638208000, 643373556450625, 115242726703104000, 309281197713168681, 14870534508544000000, 168541849294187630625 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,4
COMMENTS
Number of permutations of [n] whose cycle lengths are relatively prime to n.
LINKS
MAPLE
b:= proc(n, m) option remember; `if`(n=0, 1, add(`if`(
igcd(j, m)=1, b(n-j, m), 0)*(n-1)!/(n-j)!, j=1..n))
end:
a:= n-> b(n$2):
seq(a(n), n=0..23); # Alois P. Heinz, Oct 12 2020
MATHEMATICA
Table[n! SeriesCoefficient[Exp[Sum[Boole[GCD[n, k] == 1] x^k/k, {k, 1, n}]], {x, 0, n}], {n, 0, 22}]
CROSSREFS
Sequence in context: A227419 A242097 A220972 * A247009 A245241 A115690
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 12 2020
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 16:58 EDT 2024. Contains 371254 sequences. (Running on oeis4.)