Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #9 May 28 2019 19:35:52
%S 27,75,99,203,171,255,219,243,171,335,315,363,539,539,435,575,531,707,
%T 767,579,603,735,651,675,815,507,819,1055,891,675,939,1211,1343,1215,
%U 1035,1419,1083,1155,1179,1691,1547,1455,1715,1859,1371,1395,1443,759,1775,1539,1179,1611,1935,1659,2015,1683,2219,1755,2507,2175,1875,2255
%N a(n) = A325814(A228058(n)).
%C First negative term occurs as a(16307) = -210973, with A228058(16307) = 1289925. The next negative terms occurs as a(20807) = -242901, with A228058(20807) = 1686825.
%H Antti Karttunen, <a href="/A325824/b325824.txt">Table of n, a(n) for n = 1..25000</a>
%F a(n) = A325814(A228058(n)).
%F a(n) = A325379(n) + A325823(n).
%o (PARI)
%o up_to = 10000;
%o isA228058(n) = if(!(n%2)||(omega(n)<2),0,my(f=factor(n),y=0); for(i=1,#f~,if(1==(f[i,2]%4), if((1==y)||(1!=(f[i,1]%4)),return(0),y=1), if(f[i,2]%2, return(0)))); (y));
%o A228058list(up_to) = { my(v=vector(up_to), k=0, n=0); while(k<up_to, n++; if(isA228058(n), k++; v[k] = n)); (v); };
%o v228058 = A228058list(up_to);
%o A228058(n) = v228058[n]; \\ _Antti Karttunen_, May 23 2019
%o A034448(n) = { my(f=factorint(n)); prod(k=1, #f~, 1+(f[k, 1]^f[k, 2])); }; \\ After code in A034448
%o A048146(n) = (sigma(n)-A034448(n));
%o A325814(n) = (n-A048146(n));
%o A325824(n) = A325814(A228058(n));
%Y Cf. A228058, A325379, A325814, A325822, A325823.
%Y Cf. also A325320.
%K sign
%O 1,1
%A _Antti Karttunen_, May 23 2019