login
a(n) = A276085(n) mod bigomega(n), where bigomega (A001222) and A276085 are fully additive with a(p) = 1 and a(p) = p#/p respectively.
4

%I #6 Jun 20 2024 17:50:21

%S 0,0,0,0,1,0,0,0,1,0,1,0,1,0,0,0,2,0,2,0,1,0,1,0,1,0,2,0,0,0,0,0,1,0,

%T 2,0,1,0,1,0,0,0,2,1,1,0,1,0,1,0,2,0,3,0,1,0,1,0,2,0,1,1,0,0,0,0,2,0,

%U 1,0,2,0,1,2,2,0,0,0,0,0,1,0,2,0,1,0,1,0,3,0,2,0,1,0,1,0,1,1,2,0,0,0,1,2,1

%N a(n) = A276085(n) mod bigomega(n), where bigomega (A001222) and A276085 are fully additive with a(p) = 1 and a(p) = p#/p respectively.

%H Antti Karttunen, <a href="/A373850/b373850.txt">Table of n, a(n) for n = 2..100000</a>

%o (PARI)

%o A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*prod(i=1,primepi(f[k, 1]-1),prime(i))); };

%o A373850(n) = (A276085(n)%bigomega(n));

%Y Cf. A001222, A276085, A373832 [= a(A276086(n))], A373851, A373852 (after its initial term, gives the indices of 0's).

%K nonn

%O 2,17

%A _Antti Karttunen_, Jun 20 2024