OFFSET
1,1
COMMENTS
This sequence is a subsequence of A008846. - Ray Chandler, Jan 27 2017
FORMULA
A002144(n)^A002365(n) and A002144(n)^A002366(n) are terms of the sequence for all n. - Ray Chandler, Jan 27 2017
EXAMPLE
n = 23125, n' = 19125 and sqrt(n^2-n'^2) = 13000.
MAPLE
with(numtheory); ListA211176:= proc(q)local a, n, p;
for n from 2 to q do a:=n*add(op(2, p)/op(1, p), p=ifactors(n)[2]);
if n<>a and type(sqrt(n^2-a^2), integer) then print(n); fi;
od; end: ListA211176(10^9);
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Paolo P. Lava, Feb 01 2013
EXTENSIONS
Name and Maple program corrected by Paolo P. Lava, Sep 30 2013
a(12)-a(16) from Donovan Johnson, Sep 30 2013
a(17)-a(18) from Ray Chandler, Jan 25 2017
a(19)-a(20) from Ray Chandler, Jan 27 2017
STATUS
approved