OFFSET
1,3
COMMENTS
Sequence continues: a(29) = ?, 40320, ?, 1893397524480, 32006016000000, 478511953920, ?, 209018880, ?, 122634632171520, ?, 665481277440, ?, 399168000. - Franklin T. Adams-Watters, Sep 25 2006
EXAMPLE
The integers coprime to 10 and <= 10 are (1,3,7,9).
Of the permutations of these integers, there are 12 permutations where all adjacent terms are coprime, i.e., where 3 and 9 are not adjacent. So the 10th term of the sequence is 12.
PROG
(PARI) {A109813(x)=local(v, A, d, r, M, i, j); n=eulerphi(x); v=vector(n, i, 0); j=0; for(i=1, x, if(gcd(i, x)==1, j++; v[j]=i)); A=matrix(n, n, i, j, gcd(v[i], v[j])==1); r=0; for(s=1, 2^n-1, M=vecextract(A, s, s)^(n-1); d=matsize(M)[1]; r+=(-1)^(n-d)*sum(i=1, d, sum(j=1, d, M[i, j]))); r} \\ Franklin T. Adams-Watters, Sep 25 2006
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Leroy Quet, Aug 16 2005
EXTENSIONS
More terms from Franklin T. Adams-Watters, Sep 25 2006
STATUS
approved