login
A378880
a(n) = A378879(n) - A083025(n).
1
0, 1, 1, 2, -1, 2, 1, 3, 2, 0, 1, 3, -1, 2, 0, 4, -1, 3, 1, 1, 2, 2, 1, 4, -2, 0, 3, 3, -1, 1, 1, 5, 2, 0, 0, 4, -1, 2, 0, 2, -1, 3, 1, 3, 1, 2, 1, 5, 2, -1, 0, 1, -1, 4, 0, 4, 2, 0, 1, 2, -1, 2, 3, 6, -2, 3, 1, 1, 2, 1, 1, 5, -1, 0, -1, 3, 2, 1, 1, 3, 4, 0
OFFSET
1,4
EXAMPLE
a(10) = 0 because the factorization 2*5 has 1 each of a Pythagorean prime, 5, and a non-Pythagorean prime, 2.
MATHEMATICA
f[{x_, y_}] := If[Mod[x, 4] == 1, y, -y];
s[n_] := Map[f, FactorInteger[n]];
p[n_] := {Total[Select[s[n], # > 0 &]], -Total[Select[s[n], # < 0 &]]};
p[1] = {0, 0};
t = Table[p[n], {n, 1, 135}]
u = Map[First, t] (* A083025 *)
v = Map[Last, t] (* A378879 *)
v - u (* A377625 *)
CROSSREFS
KEYWORD
sign
AUTHOR
Clark Kimberling, Jan 14 2025
STATUS
approved