OFFSET
1,2
COMMENTS
From Robert Israel, Sep 05 2024: (Start)
According to A000203, sigma_1(m) <= (6/Pi^2)*m^(3/2) for m >= 12. Thus since sigma_2(m) > m^2, sigma_2(m)/sigma_1(m) > (Pi^2/6) * m^(1/2).
This would suggest that to find all terms <= K of this sequence we should look at sigma_2(m)/sigma_1(m) for m <= 36 * K^2/Pi^4. But using the b-file for A004394 we may get a good upper bound for sigma_1(m)/m for m in this interval, resulting in a much smaller search region. (End)
LINKS
Robert Israel, Table of n, a(n) for n = 1..1009
MAPLE
# This uses the b-file for A004394
K:= 10000: # to get terms <= K
M:= 36 * K^2/Pi^4:
for i from 1 while A004394[i] < M do od:
V:= Vector(K):
for m from 1 to r*K do
F:= numtheory:-divisors(m);
v:= add(d^2, d=F)/add(d, d=F);
if v::integer and v <= K and V[v] = 0 then V[v]:= m fi;
od:
select(v -> V[v] > 0, [$1..K]); # Robert Israel, Sep 05 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Apr 26 2010
EXTENSIONS
More terms from Robert Israel, Sep 05 2024
STATUS
approved