OFFSET
1,4
PROG
(PARI) a(n)=sum(k=1, n!, if(sum(i=1, n, gcd(i, component(numtoperm(n, k), i)))-sigma(n), 0, 1))
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Benoit Cloitre, Aug 23 2002
EXTENSIONS
a(10) from Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 14 2004
a(11) from Rick L. Shepherd, Mar 15 2004
a(12)-a(13) from Ryan Propper, Oct 25 2006
a(14)-a(20) from Robert Gerbicz, Nov 19 2010
STATUS
approved