OFFSET
1,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 687 terms from Matthew Campbell)
EXAMPLE
a(2) = hyperfactorial(2) mod 3 = (2^2*1^1) mod 3 = 4 mod 3 = 1.
MAPLE
a:= proc(n) option remember; local i, p, r, v;
p, r, v:= ithprime(n), 1$2;
for i from p-1 to 1 by -1 do
v:= v*i mod p; r:= r*v mod p
od; r
end:
seq(a(n), n=1..100); # Alois P. Heinz, Jul 21 2015
MATHEMATICA
Table[Mod[Hyperfactorial[Prime[n] - 1], Prime[n]], {n, 1, 200}]
PROG
(PARI) a(n, p=prime(n))=lift(prod(k=2, p-1, Mod(k, p)^k)) \\ Charles R Greathouse IV, Jul 23 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Matthew Campbell, Jul 17 2015
STATUS
approved