|
%I #9 Feb 06 2018 02:53:05
%S 57,239,1985,10133,9466,11389,27590,51412,153765,344464,107551,296344,
%T 172078,432436,931837,753090,676541,2321221,2027724,3394758,1706203,
%U 4841182,1438398,2947125,398366,5657795,4942017,9400802,11906503
%N Smallest k such that k^2 + 1 is divisible by A002144(n)^3.
%H Klaus Brockhaus, <a href="/A145296/b145296.txt">Table of n, a(n) for n = 1..150</a>
%e a(3) = 1985 since A002144(3) = 17, 1985^2 + 1 = 3940226 = 2*17^3*401 and for no k < 1985 does 17^3 divide k^2+1.
%o (PARI) {m=12000000; pmax=300; z=70; v=vector(z); for(n=1, m, fac=factor(n^2+1); for(j=1, #fac[, 1], if(fac[j, 2]>=3&&fac[j, 1]<=pmax, q=primepi(fac[j, 1]); if(q<=z&&v[q]==0, v[q]=n)))); t=1; j=0; while(t&&j<z, j++; p=prime(j); if(p%4==1, if(v[j]==0, t=0, print1(v[j], ","))))}
%o (PARI) {e=3; forprime(p=2, 300, if(p%4==1, q=p^e; m=q; while(!ispower(m-1,2,&n), m=m+q); print1(n, ",")))} \\ _Klaus Brockhaus_, Oct 09 2008
%Y Cf. A002144 (primes of form 4n+1), A002313 (-1 is a square mod p), A059321, A145297, A145298, A145299.
%K nonn
%O 1,1
%A _Klaus Brockhaus_, Oct 08 2008
|
