OFFSET
1,4
COMMENTS
For prime p, a(p) = 1.
If n is not a prime or the square of a prime then n divides a(n).
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
For n = 16 the nontrivial divisors d are 2,4 and 8, so a(16) = (16-2)*(16-4)*(16-8) = 14*12*8 = 1344.
MAPLE
f:= proc(n) local d; mul(n-d, d = numtheory:-divisors(n) minus {1, n}) end proc:
map(f, [$1..50]); # Robert Israel, Dec 30 2024
PROG
(PARI) a(n) = my(d=divisors(n)); prod(j=2, matsize(d)[2]-1, n-d[j]);
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Amarnath Murthy, Jul 28 2002
EXTENSIONS
Edited and extended by Klaus Brockhaus, Jul 31 2002
STATUS
approved