OFFSET
1,1
COMMENTS
m = k^2, where k satisfies k' = (k+1)/2. - Charlie Neder, Mar 08 2019
It follows that 2k is a Giuga number satisfying (2k)' = 2k + 1. - Max Alekseyev, Sep 09 2025
FORMULA
For n <= 12, a(n) = (A007850(n)/2)^2. - Max Alekseyev, Sep 09 2025
EXAMPLE
For m = 225 we have that m' = 240, sqrt(225) = 15 and 240 = 225 + 15.
MAPLE
with(numtheory); P:= proc(q) local n, p, x;
for n from 1 to q do x:=n^2;
if x*add(op(2, p)/op(1, p), p=ifactors(x)[2])=n^2+n then print(n^2);
fi; od; end: P(10^9);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Mar 20 2014
EXTENSIONS
Terms a(5) onward from Max Alekseyev, Sep 09 2025
STATUS
approved
