login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A325378 a(n) = A162296(A228058(n)) - A048250(A228058(n)). 6
30, 70, 90, 246, 150, 266, 190, 210, 678, 342, 270, 310, 654, 574, 370, 570, 450, 738, 930, 490, 510, 722, 550, 570, 798, 1582, 690, 1026, 750, 2034, 790, 1230, 1626, 1178, 870, 1526, 910, 970, 990, 2046, 1558, 1406, 1722, 1962, 1150, 1170, 1210, 4062, 1710, 1290, 3390, 1350, 1862, 1390, 1938, 1410, 2214, 1470, 3030 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
LINKS
FORMULA
a(n) = A162296(A228058(n)) - A048250(A228058(n)).
a(n) = A325319(n) + A325320(n).
PROG
(PARI)
A048250(n) = factorback(apply(p -> p+1, factor(n)[, 1]));
A162296(n) = sumdiv(n, d, d*(1-issquarefree(d)));
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));
k=0; n=0; while(k<100, n++; if(isA228058(n), k++; print1(A162296(n) - A048250(n), ", ")));
CROSSREFS
Sequence in context: A064623 A112343 A182996 * A164596 A295102 A131647
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 22 2019
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 05:39 EDT 2024. Contains 371235 sequences. (Running on oeis4.)