|
|
|
|
3, 1, 9, 1, 3, 1, 1, 3, 3, 1, 9, 1, 1, 7, 1, 5, 9, 1, 1, 1, 3, 1, 1, 3, 1, 1, 3, 1, 3, 9, 1, 1, 1, 1, 3, 1, 1, 1, 9, 1, 1, 1, 7, 1, 1, 9, 1, 3, 5, 3, 3, 9, 1, 1, 1, 3, 1, 3, 1, 5, 1, 1, 3, 1, 1, 1, 1, 1, 1, 3, 3, 1, 1, 3, 1, 1, 1, 3, 3, 1, 25, 1, 1, 9, 1, 1, 9, 1, 3, 1, 27, 1, 1, 1, 1, 3, 9, 1, 49, 1, 1, 1, 1, 1, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
|
|
PROG
|
(PARI)
up_to = 25000;
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));
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); };
v228058 = A228058list(up_to);
A048250(n) = factorback(apply(p -> p+1, factor(n)[, 1]));
A162296(n) = sumdiv(n, d, d*(1-issquarefree(d)));
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|