|
|
A074730
|
|
Squares whose sum of anti-divisors is a square.
|
|
1
|
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
See A066272 for definition of anti-divisor.
|
|
LINKS
|
|
|
PROG
|
(Python)
from gmpy2 import is_square
from sympy import divisors
A074730 = [n for n in (x**2 for x in range(1, 10**4)) if is_square(int(sum([2*d for d in divisors(n) if n > 2*d and n%(2*d)] + [d for d in divisors(2*n-1) if n > d >=2 and n%d] + [d for d in divisors(2*n+1) if n > d >=2 and n%d])))]
(PARI) has(n)=my(k=valuation(n, 2)); issquare(sigma(2*n+1)+sigma(2*n-1)+sigma(n>>k)<<(k+1)-6*n-2)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|