OFFSET
1,5
COMMENTS
A refinement of Stirling numbers of the first kind.
FORMULA
If n = Product prime(x_i)^y_i is the prime factorization of n, then a(n) = (-1)^(Sum x_i * y_i - Sum y_i) (Sum x_i * y_i)! / (Product x_i^y_i * Product y_i!).
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
numPermsOfType[ptn_]:=Total[ptn]!/Times@@ptn/Times@@Factorial/@Length/@Split[ptn];
Table[(-1)^(Total[primeMS[n]]-PrimeOmega[n])*numPermsOfType[primeMS[n]], {n, 100}]
CROSSREFS
KEYWORD
sign
AUTHOR
Gus Wiseman, Sep 13 2018
STATUS
approved