login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A186814 a(n) = smallest number m such that A002144(n) divides gcd(A002314(n)^2+1,(A002314(n)+m)^2+1). 2
1, 3, 9, 5, 25, 23, 7, 39, 19, 21, 53, 81, 43, 83, 63, 61, 101, 13, 143, 31, 169, 15, 55, 113, 225, 105, 157, 175, 17, 263, 89, 41, 77, 269, 165, 159, 271, 361, 123, 363, 75, 315, 239, 365, 93, 51, 437, 321, 397, 529, 439, 351, 543, 229, 333, 355, 449, 557, 625, 431, 517, 27, 583 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Sequence A002314 gives the minimal integer square root of -1 modulo p(n),where p(n) = n-th prime of form 4k+1.
LINKS
EXAMPLE
for n=1, k = A002314(1) = 2 => a(1) = 1, because 2^2+1 = 5 and (2+1)^2+1 = 2*5 ;
for n=2, k = A002314(2) = 5 => a(2) = 3, because 5^2+1 = 2*13 and (5+3)^2+1 = 5*13 ;
for n=3, k = A002314(3) = 4 => a(3) = 9, because 4^2+1 = 17 and (4+9)^2+1 = 2*5*17;
for n=4, k = A002314(4)= 12 => a(4)= 5, because 12^2+1
= 5*29 and (12+5)^2+1 = 2*5*29, and 29 divides
GCD(5*29, 2*5*29)=145.
MAPLE
with(numtheory):T:=array(1..90):j:=1:for i from 1 to 250 do:x:=4*i+1:if type(x, prime)=true
then T[j]:=x:j:=j+1:else fi:od:for p from 1 to j do:u:=T[p]:id:=0: for m from
1 to 1000 while(id=0) do: z:=m^2+1:for d from 1 to u while(id=0) do: z1:=(m+d)^2+1:zz:=
gcd(z, z1):if irem(zz, u)=0 then id:=1:printf(`%d, `, d):else fi:od:od:od:
CROSSREFS
Sequence in context: A252117 A103934 A365203 * A077384 A122943 A276148
KEYWORD
nonn
AUTHOR
Michel Lagneau, Feb 27 2011
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 09:44 EDT 2024. Contains 371268 sequences. (Running on oeis4.)