OFFSET
1,1
COMMENTS
The asymptotic density of this sequence is Sum_{n>=1} f(A356871(n)) = 9.1348...*10^(-6), where f(n) = (4/(Pi^2*n))*Product_{prime p|n}(p/(p+1)) if n is odd and 0 otherwise. - Amiram Eldar, Sep 02 2022
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
99225 is a term since it is odd and the sum of its coreful divisors is A057723(99225) = 201600 > 2 * 99225.
MATHEMATICA
f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; s[1] = 1; s[n_] := Times @@ (f @@@ FactorInteger[n]); Select[Range[1, 10^6, 2], s[#] > 2*# &]
CROSSREFS
Subsequence of A321147.
KEYWORD
nonn
AUTHOR
Amiram Eldar, Dec 23 2020
STATUS
approved