login
A073343
Number of permutations p of (1,2,3,...,n) such that Sum_{k=1..n} gcd(k,p(k)) = sigma(n).
0
1, 1, 1, 6, 33, 75, 877, 2075, 38692, 239794, 1250340, 18495201, 43023264, 4554211553, 70490483945, 911526001064, 182477118528, 219638126087356, 13070186026272, 79249480638985626
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
Sequence in context: A140521 A171141 A069065 * A157872 A153127 A274218
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