|
%I #28 Sep 08 2022 08:45:12
%S 2,5,2,17,2,37,2,5,2,101,2,5,2,197,2,257,2,5,2,401,2,5,2,577,2,677,2,
%T 5,2,17,2,5,2,13,2,1297,2,5,2,1601,2,5,2,13,2,29,2,5,2,41,2,5,2,2917,
%U 2,3137,2,5,2,13,2,5,2,17,2,4357,2,5,2,13,2,5,2,5477,2,53,2,5,2,37,2,5,2
%N Smallest prime factor of n^2 + 1.
%C This includes A002496, primes that are of the form n^2+1.
%C Note that a(n) is the smallest prime p such that n^(p+1) == -1 (mod p). - _Thomas Ordowski_, Nov 08 2019
%D H. Rademacher, Lectures on Elementary Number Theory, pp. 33-38.
%H Michael De Vlieger, <a href="/A089120/b089120.txt">Table of n, a(n) for n = 1..10000</a>
%F a(2k+1)=2; a(10k +/- 2)=5, else a(26k +/- 8)=13, else a(34k +/- 4)=17, else a(58k +/- 12)=29, else a(74k +/- 6)=37,... - _M. F. Hasler_, Mar 11 2012
%F A089120(n) = 2 if n is odd, else A089120(n) = min { A002144(k) | n = +/- A209874(k) (mod 2*A002144(k)) }.
%t Array[FactorInteger[#^2 + 1][[1, 1]] &, {83}] (* _Michael De Vlieger_, Sep 08 2015 *)
%o (PARI) smallasqp1(m) = { for(a=1,m, y=a^2 + 1; f = factor(y); v = component(f,1); v1 = v[length(v)]; print1(v[1]",") ) }
%o (PARI) A089120(n)=factor(n^2+1)[1,1] \\ _M. F. Hasler_, Mar 11 2012
%o (Magma) [Min(PrimeDivisors(n^2+1)):n in [1..100]]; // _Marius A. Burtea_, Nov 13 2019
%Y Cf. A002496, A014442, A085722, A144255, A193432.
%K easy,nonn
%O 1,1
%A _Cino Hilliard_, Dec 05 2003
|
