|
%I #15 Jun 20 2015 08:54:22
%S 1,2,5,10,17,26,37,2,65,82,101,122,145,170,197,226,257,290,13,362,401,
%T 442,485,530,577,626,677,730,785,842,901,962,41,1090,1157,1226,1297,
%U 1370,5,1522,1601,2,1765,74,1937,2026,2117,2210,2305,2402,2501,2602
%N Squarefree part of n^2+1.
%C Related to period-1 continued fractions [z,z,z,...].
%H Enrique Pérez Herrero, <a href="/A059591/b059591.txt">Table of n, a(n) for n = 0..5000</a>
%F a(n) = A007913(n^2+1).
%F a(n)*A059592(n)^2 = n^2+1.
%e a(7)=2 since 7^2+1 = 50 = 25*2 = (5^2)*2.
%p A:= proc(n)
%p local F;
%p F:= select(t -> t[2]::odd, ifactors(n^2+1)[2]);
%p mul(t[1],t=F)
%p end proc:
%p map(A, [$0..100]); # _Robert Israel_, Jun 18 2015
%t Table[Times @@ Flatten[Table[#1, {#2}] & @@@ Select[FactorInteger[n^2 + 1], OddQ@ Last@ # &]], {n, 120}] (* _Michael De Vlieger_, Jun 19 2015 *)
%o (PARI) a(n) = core(n^2+1); \\ _Michel Marcus_, Jun 18 2015
%Y Cf. A007913, A059592, A000188.
%K easy,nonn
%O 0,2
%A _Marc LeBrun_, Jan 25 2001
|
