OFFSET
1,4
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(p) = 1, for prime p.
MAPLE
f:= proc(n) convert(select(t -> t^2 <= n, NumberTheory:-PrimeFactors(n)), `*`) end proc:
map(f, [$1..100]); # Robert Israel, Feb 17 2026
MATHEMATICA
Table[Times @@ Select[Divisors[n], PrimeQ[#] && # <= Sqrt[n] &], {n, 1, 85}]
PROG
(PARI) a(n) = vecprod(select(x->x<=sqrt(n), factor(n)[, 1])); \\ Michel Marcus, Apr 17 2025
(Python)
from math import isqrt
from sympy import primefactors
def A382486(n):
m = isqrt(n)
return prod(p for p in primefactors(n) if p<=m) # Chai Wah Wu, May 20 2026
CROSSREFS
KEYWORD
AUTHOR
Ilya Gutkovskiy, Apr 10 2025
STATUS
approved
