A325317 a(n) = A048250(n) XOR A162296(n), where XOR is the bitwise-XOR, A003987.
1, 3, 4, 7, 6, 12, 8, 15, 13, 18, 12, 28, 14, 24, 24, 31, 18, 23, 20, 10, 32, 36, 24, 60, 31, 42, 32, 56, 30, 72, 32, 63, 48, 54, 48, 67, 38, 60, 56, 90, 42, 96, 44, 20, 46, 72, 48, 124, 57, 89, 72, 18, 54, 96, 72, 120, 80, 90, 60, 40, 62, 96, 104, 127, 84, 144, 68, 126, 96, 144, 72, 187, 74, 114, 124, 108, 96, 168, 80
Offset: 1
Keywords
Links
Programs
-
Mathematica
Array[BitXor @@ Map[Total, {#3, Complement[#2, #3]}] & @@ {#1, #2, Select[#2, SquareFreeQ]} & @@ {#, Divisors[#]} &, 79] (* Michael De Vlieger, Apr 21 2019 *)
-
PARI
A048250(n) = factorback(apply(p -> p+1,factor(n)[,1])); A162296(n) = sumdiv(n, d, d*(1-issquarefree(d))); A325317(n) = bitxor(A048250(n),A162296(n));