|
0, 0, 0, 2, 0, 4, 2, 0, 0, 0, 0, 0, 0, 0, 8, 6, 0, 0, 2, 0, 4, 0, 16, 0, 0, 0, 0, 2, 8, 0, 6, 0, 0, 0, 0, 0, 0, 0, 0, 34, 0, 0, 10, 0, 0, 4, 0, 0, 0, 0, 0, 2, 32, 0, 2, 20, 8, 0, 0, 0, 4, 0, 0, 10, 0, 0, 2, 0, 64, 0, 4, 0, 0, 0, 0, 2, 0, 8, 2, 0, 0, 0, 0, 0, 0, 68, 0, 2, 16, 0, 2, 8, 0, 0, 0, 4, 0, 0, 0, 2, 4, 0, 66
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,4
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
Table[If[# == 1, 0, Total[#2*2^PrimePi[#1] & @@@ FactorInteger[#]]/2] &@ GCD[n, Times @@ Prime@ Flatten@ Position[Reverse@ IntegerDigits[n, 2], 1]], {n, 102}] (* Michael De Vlieger, Feb 20 2022 *)
|
|
PROG
|
(PARI)
A019565(n) = { my(m=1, p=1); while(n>0, p = nextprime(1+p); if(n%2, m *= p); n >>= 1); (m); };
A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
|
|
CROSSREFS
|
|
|
KEYWORD
|
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|