cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A227890 Primes of the form prime(k)^2 - k.

This page as a plain text file.
%I A227890 #15 Jan 22 2017 17:46:47
%S A227890 3,7,163,353,5021,12739,32719,49681,52391,78901,113501,252913,361091,
%T A227890 452807,551917,993841,1559797,1956979,2193127,3463037,4067983,5003837,
%U A227890 5138953,6115363,6723271,7251857,7447043,7578607,8426989,9479801,11295847,12186593,12439237
%N A227890 Primes of the form prime(k)^2 - k.
%H A227890 K. D. Bajpai, <a href="/A227890/b227890.txt">Table of n, a(n) for n = 1..7200</a>
%e A227890 a(3)= 163: prime(6)^2 - 6= 13^2 - 6= 169 - 6= 163 which is prime.
%e A227890 a(4)= 353: prime(8)^2 - 8= 19^2 - 8= 361 - 8= 353 which is prime.
%p A227890 with(numtheory):KD := proc() local a; a:= (ithprime(k)^2-k); if isprime(a) then RETURN (a); fi; end: seq(KD(), k=1..1000);
%t A227890 Select[Table[Prime[k]^2-k,{k,1000}],PrimeQ]
%o A227890 (PARI) for(k=1, 10^5, if(ispseudoprime(KD=((prime(k)^2-k))), print1(KD", ")));
%Y A227890 Cf. A000040 (prime numbers).
%Y A227890 Cf. A064713 (for the integers k).
%Y A227890 Cf. A184935 (primes: k^2 + prime(k)).
%Y A227890 Cf. A188831 (primes: k^2 - prime(k)).
%Y A227890 Cf. A229203 (primes: k^3 - prime(k)).
%K A227890 nonn
%O A227890 1,1
%A A227890 _K. D. Bajpai_, Oct 26 2013