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.

A109315 Numbers n such that prime(n) - n is a prime power.

This page as a plain text file.
%I A109315 #10 May 06 2013 04:14:25
%S A109315 12,15,38,39,118,152,190,258,462,690,746,1396,1632,2119,3370,4522,
%T A109315 4600,7520,15006,24222,33156,34038,51372,52342,64638,77470,90790,
%U A109315 101946,104670,156772,166822,167700,175818,194092,200022,229630,246208,328462,362440,372882
%N A109315 Numbers n such that prime(n) - n is a prime power.
%H A109315 Donovan Johnson, <a href="/A109315/b109315.txt">Table of n, a(n) for n = 1..1000</a>
%F A109315 prime(n) - n = q^k, q is prime and k_Integer >= 2.
%e A109315 690 is OK because prime(690)-690 = 5179-690 = 4489 = 67^2, 67 is prime.
%t A109315 lst = {}; fQ[n_] := Block[{pf=FactorInteger[n]}, (2-Length[pf])(pf[[1, 2]]-1) > 0]; Do[ If[ fQ[Prime[n] - n], Print[n]; AppendTo[lst, n]], {n, 3, 362439}]; lst
%Y A109315 Cf. A025475 = powers of a prime but not prime, also nonprime n such that sigma(n)*phi(n)>(n-1)2; A107712 = values of q, A107713 = values of k; A107714 = values of prime(A109315(n)).
%Y A109315 Cf. A083240.
%K A109315 nonn
%O A109315 1,1
%A A109315 _Zak Seidov_ and _Robert G. Wilson v_, May 22 2005