login
a(n) = A001065(n) - A003557(n), where A001065(n) = the sum of proper divisors of n, and A003557(n) = n divided by its largest squarefree divisor.
3

%I #11 Jul 12 2019 07:21:39

%S -1,0,0,1,0,5,0,3,1,7,0,14,0,9,8,7,0,18,0,20,10,13,0,32,1,15,4,26,0,

%T 41,0,15,14,19,12,49,0,21,16,46,0,53,0,38,30,25,0,68,1,38,20,44,0,57,

%U 16,60,22,31,0,106,0,33,38,31,18,77,0,56,26,73,0,111,0,39,44,62,18,89,0,98,13,43,0,138,22,45,32,88,0,141,20

%N a(n) = A001065(n) - A003557(n), where A001065(n) = the sum of proper divisors of n, and A003557(n) = n divided by its largest squarefree divisor.

%H Antti Karttunen, <a href="/A326188/b326188.txt">Table of n, a(n) for n = 1..16384</a>

%H Antti Karttunen, <a href="/A326188/a326188.txt">Data supplement: n, a(n) computed for n = 1..65537</a>

%F a(n) = A326187(n) - n = A000203(n) - A003557(n) - n.

%F a(n) = A001065(n) - A003557(n).

%o (PARI)

%o A003557(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 2] = f[i, 2]-1); factorback(f); }; \\ From A003557

%o A326188(n) = ((sigma(n)-A003557(n))-n);

%Y Cf. A000203, A001065, A003557, A010848, A326187, A326189.

%Y Cf. also A326185.

%K sign

%O 1,6

%A _Antti Karttunen_, Jul 11 2019