VOOZH
about
URL: https://oeis.org/A068333
⇱ A068333 - OEIS
login
A068333
Product(n/k - k) where the product is over the divisors k of n and where 1 <= k <= sqrt(n).
2
0, 1, 2, 0, 4, 5, 6, 14, 0, 27, 10, 44, 12, 65, 28, 0, 16, 357, 18, 152, 80, 189, 22, 2300, 0, 275, 156, 972, 28, 2639, 30, 1736, 256, 495, 68, 0, 36, 629, 380, 12636, 40, 8569, 42, 6020, 2112, 945, 46, 215072, 0, 5635, 700, 11016, 52, 59625
(
list
;
graph
;
refs
;
listen
;
history
;
text
;
internal format
)
OFFSET
1,3
COMMENTS
From
Robert Israel
, Jun 02 2019: (Start)
a(n) is divisible by n-1.
a(n) = 0 if and only if n is a square.
a(n) = n-1 if n is prime. (End)
LINKS
Robert Israel,
Table of n, a(n) for n = 1..10000
EXAMPLE
a(8) = (8 - 1) (4 - 2) = 14 because 1 and 2 are the divisors of 8 which are <= sqrt(8).
MAPLE
f:= proc(n) local D, k;
D:= select(t -> t^2 <= n, numtheory:-divisors(n));
mul(n/k-k, k=D)
end proc:
map(f, [$1..100]); #
Robert Israel
, Jun 02 2019
MATHEMATICA
a[n_] := Product[If[1 <= k <= Sqrt[n], (n/k - k), 1], {k, Divisors[n]}];
Array[a, 100] (*
Jean-François Alcover
, Aug 16 2020 *)
PROG
(PARI) a(n) = my(p=1); fordiv(n, d, if (d^2 <= n, p *= n/d - d)); p; \\
Michel Marcus
, Jun 02 2019
CROSSREFS
Sequence in context:
A049271
A004178
A385712
*
A378748
A394040
A121451
Adjacent sequences:
A068330
A068331
A068332
*
A068334
A068335
A068336
KEYWORD
nonn
,
look
AUTHOR
Leroy Quet
, Feb 27 2002
STATUS
approved