Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.
%I #29 Feb 18 2021 00:19:09
%S 4,27,1808,3125,12204,12707,82377,269827,823543,1412500,7089739,
%T 9534375,46873785,78192979,372241436
%N Numbers that form a Pythagorean 5-tuple with their first three arithmetic derivatives.
%e First three arithmetic derivatives of 1808 are 3632, 7280, 17616 and sqrt(1808^2 + 3632^2 + 7280^2 + 17616^2) = 19488.
%p with(numtheory);
%p Dr:=proc(w) local x,p; x:=w*add(op(2,p)/op(1,p),p=ifactors(w)[2]); end:
%p P:=proc(q,h) local a,b,k,n; for n from 2 to q do a:=n; b:=n^2;
%p for k from 1 to h do a:=Dr(a); b:=b+a^2; od; if type(sqrt(b),integer) then print(n);
%p fi; od; end: P(10^9,3);
%Y A051674 is a subsequence.
%Y Cf. A003415, A210503, A230543, A249106, A249107, A249110.
%K nonn,more
%O 1,1
%A _Paolo P. Lava_, Oct 21 2014
%E a(13) from _Ray Chandler_, Dec 23 2016
%E a(14) from _Ray Chandler_, Dec 24 2016
%E a(15) from _Ray Chandler_, Jan 08 2017