OFFSET
1,1
COMMENTS
For n to qualify, the sum of these reciprocals must be less than one. For example, 350 = 2*5*5*7, with sum of reciprocals 1/2 + 1/5 + 1/5 + 1/7 = 73/70, which is greater than 1, so 350 is not in the sequence.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
Take 238 = 2*7*17 with sum of reciprocals 1/2 + 1/7 + 1/17 = 167/238, giving 238*(167/238) = 167 and 238-167=71. Both 167 and 71 are prime.
MATHEMATICA
srpfQ[n_]:=Module[{s=Total[Flatten[1/Table[First[#], {Last[#]}]&/@ FactorInteger[ n]]], ns}, ns=n*s; And@@PrimeQ[{ns, n-ns}]]; Select[Range[ 1600], srpfQ] (* Harvey P. Dale, Jul 01 2012 *)
PROG
(Magma) T:=[ Integers()!(n*&+[ d[2]/d[1]: d in Factorization(n) ]): n in [2..1600] ]; [ k+1: k in [1..#T] | IsPrime(s) and IsPrime((k+1-s)) where s is T[k] ]; // Klaus Brockhaus, Apr 27 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
J. M. Bergot, Apr 23 2009
EXTENSIONS
Edited, 282 and 390 inserted, and extended beyond 418 by Klaus Brockhaus, Apr 27 2009
STATUS
approved