OFFSET
1,5
LINKS
Paolo Xausa, Table of n, a(n) for n = 1..1000
FORMULA
a(p) = A102647(p) if p prime.
MATHEMATICA
A383752[n_] := Times @@ DeleteCases[Mod[n, Prime[Range[PrimePi[n - 2]]]], 0];
Array[A383752, 50] (* Paolo Xausa, Jun 05 2025 *)
PROG
(Python)
from sympy import primerange
def a(n):
s = 1
for p in primerange(0, n):
if p > (n >> 1): s *= (n-p)
elif (x:= n % p) > 0: s*= x
return s
print([a(n) for n in range(1, 41)])
(PARI) a(n) = vecprod(select(x->(x!=0), apply(lift, apply(x->Mod(n, x), primes([2, n-1]))))); \\ Michel Marcus, May 28 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Darío Clavijo, May 28 2025
STATUS
approved
