|
|
A353976
|
|
Product of the divisors of n whose arithmetic derivative is even.
|
|
1
|
|
|
1, 1, 1, 4, 1, 1, 1, 32, 9, 1, 1, 48, 1, 1, 15, 512, 1, 9, 1, 80, 21, 1, 1, 9216, 25, 1, 9, 112, 1, 15, 1, 16384, 33, 1, 35, 15552, 1, 1, 39, 25600, 1, 21, 1, 176, 135, 1, 1, 7077888, 49, 25, 51, 208, 1, 9, 55, 50176, 57, 1, 1, 864000, 1, 1, 189, 1048576, 65, 33, 1, 272, 69
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Product_{d|n} d^((1+d') mod 2), where d' is the arithmetic derivative of d.
|
|
EXAMPLE
|
a(8) = 32; the product of the divisors of 8 whose arithmetic derivative is even is 4*8 = 32.
|
|
MATHEMATICA
|
d[1] = 0; d[n_] := n*Plus @@ ((Last[#]/First[#]) & /@ FactorInteger[n]); a[n_] := Product[k^Mod[d[k] + 1, 2], {k, Divisors[n]}]; Array[a, 60] (* Amiram Eldar, May 13 2022 *)
|
|
PROG
|
(PARI) ad(n) = vecsum([n/f[1]*f[2]|f<-factor(n+!n)~]); \\ A003415
a(n) = vecprod(select(x->!(ad(x) % 2), divisors(n))); \\ Michel Marcus, May 13 2022
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|