login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A130767
a(n) = product{k=1 to n} (k-th integer from among those positive integers which are coprime to (n+1-k)).
1
1, 2, 9, 40, 420, 2700, 56595, 419328, 8820900, 88488400, 2327925600, 38767286880, 1912404574080, 21612951360000, 644047087612500, 10985391056640000, 634391869996684800, 14046187624838328960, 764077915447610400000, 15840110879873280000000, 755098009918296312668400
OFFSET
1,2
COMMENTS
a(n) is the product of the terms in the n-th antidiagonal of the A126572 array. - Michel Marcus, Mar 14 2018
EXAMPLE
The integers coprime to 1 are: 1,2,3,4,5,6,... The 5th of these is 5. The integers coprime to 2 are: 1,3,5,7,9,... The 4th of these is 7. The integers coprime to 3 are: 1,2,4,5,7,... The 3rd of these is 4. The integers coprime to 4 are: 1,3,5,... The 2nd of these is 3. And the integers coprime to 5 are: 1,2,3,4,6,... The first of these is 1. So a(5) = 5 * 7 * 4 * 3 * 1 = 420.
PROG
(PARI) cop(k, j) = {my(nbc = 0, i = 0); while (nbc != j, i++; if (gcd(i, k)==1, nbc++)); i; }
a(n) = {my(vc = vector(n, k, cop(k, n-k+1))); prod(k=1, n, vc[k]); } \\ Michel Marcus, Mar 14 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Aug 18 2007
EXTENSIONS
More terms from Michel Marcus, Mar 14 2018
STATUS
approved