OFFSET
0,6
COMMENTS
k is strongly prime to n iff k is relatively prime to n and k does not divide n-1.
For 0 we have the empty product, giving 1. - Daniel Forgues, Aug 03 2012
From Robert G. Wilson v, Aug 04 2012: (Start)
Records appear at positions 0, 5, 7, 9, 11, 13, 17, 19, 23, 29, 31, ....
Except for 0 and 9, all records appear at prime positions and beginning with the sixth term, are == 0 (mod 100).
There are some primes which are not records: 2, 3, 61, 73, 109, 151, 181, 193, 229, 241, 271, 313, 349, 421, 433, 463, ....
Anti-records appear at positions 6, 10, 12, 14, 15, 18, 20, 24, 30, 36, 42, 48, 60, 66, 70, 78, 84, 90, 96, ..., and their values are odd. (End)
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 0..1000
Peter Luschny, Strong coprimality.
EXAMPLE
a(11) = 3 * 4 * 6 * 7 * 8 * 9 = 36288.
MAPLE
with(numtheory):
StrongCoprimes := n -> select(k->igcd(k, n)=1, {$1..n}) minus divisors(n-1):
A181832 := proc(n) local i; mul(i, i=StrongCoprimes(n)) end:
coprimorial := proc(n) local i; mul(i, i=select(k->igcd(k, n)=1, [$1..n])) end:
divisorial := proc(n) local i; mul(i, i=divisors(n)) end:
A181832a := n -> `if`(n=0, 1, coprimorial(n)/divisorial(n-1)):
MATHEMATICA
f[n_] := Times @@ Select[ Range@ n, GCD[#, n] == 1 && Mod[n - 1, #] != 0 &]; Array[f, 27, 0] (* Robert G. Wilson v, Aug 03 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Nov 17 2010
STATUS
approved