login
Squarefree part of n^2+1.
5

%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