A248527 - OEIS

%I #17 Sep 08 2022 08:46:10

%S 34,44,60,70,86,96,164,174,190,200,216,226,294,304,320,330,346,356,

%T 424,434,450,460,476,486,554,564,580,590,606,616,684,694,710,720,736,

%U 746,814,824,840,850,866,876,944,954,970,980,996,1006,1074,1084,1100,1110

%N Numbers n such that the smallest prime divisor of n^2+1 is 13.

%C Or numbers n such that the smallest prime divisor of A002522(n) is A002313(3).

%C a(n) == 8 (mod 26) if n is odd and a(n) == 18 (mod 26) if n is even.

%C It is interesting to observe that a(n) is given by a linear formula (see the formula below).

%H Amiram Eldar, <a href="/A248527/b248527.txt">Table of n, a(n) for n = 1..10000</a>

%F {a(n)} = {8+(k + m)*26} union {18+(k + m)*26} for m = 0, 5, 10,...,5p,... and k = 1, 2, 3 (values in increasing order).

%e 34 is in the sequence because 34^2+1= 13*89.

%p * first program *

%p with(numtheory):p:=13:

%p    for n from 1 to 1000 do:

%p     if factorset(n^2+1)[1] = p then printf(`%d, `, n):

%p     else

%p     fi:

%p    od:

%p * second program using the formula*

%p for n from 0 to 100 by 5 do:

%p    for k from 1 to 3 do:

%p      x:=8+(k+n)*26:y:=18+(k+n)*26:

%p      printf(`%d, `,x):printf(`%d, `,y):

%p    od:

%p   od:

%t lst={};Do[If[FactorInteger[n^2+1][[1,1]]==13,AppendTo[lst,n]],{n,2,2000}];lst

%t p = 13; ps = Select[Range[p - 1], Mod[#, 4] != 3 && PrimeQ[#] &]; Select[Range[1200], Divisible[(nn = #^2 + 1), p] && ! Or @@ Divisible[nn, ps] &] (* _Amiram Eldar_, Aug 16 2019 *)

%o (PARI) isok(n) = factor(n^2+1)[1,1] == 13; \\ _Michel Marcus_, Oct 08 2014

%o (Magma) [n: n in [2..3000] | PrimeDivisors(n^2+1)[1] eq 13]; // _Bruno Berselli_, Oct 08 2014

%Y Cf. A002522, A089120, A002313.

%K nonn

%O 1,1

%A _Michel Lagneau_, Oct 08 2014