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.

A050705 Composite number such that when sum of its prime factors is added or subtracted becomes prime.

This page as a plain text file.
%I A050705 #16 Jul 03 2017 02:05:17
%S A050705 10,12,14,15,20,21,26,33,35,38,44,48,51,65,68,86,93,96,111,112,116,
%T A050705 123,161,188,201,203,206,209,210,215,221,278,297,300,304,306,321,352,
%U A050705 356,371,384,395,398,413,420,441,471,485,524,533,543,545,546,551,570,626
%N A050705 Composite number such that when sum of its prime factors is added or subtracted becomes prime.
%C A050705 Prime factors counted with multiplicity, e.g., 44 = 2*2*11 so the sum of its prime factors is 15 (not 13). - _Harvey P. Dale_, May 30 2012
%H A050705 Harvey P. Dale, <a href="/A050705/b050705.txt">Table of n, a(n) for n = 1..1000</a>
%e A050705 E.g., 545 = 5*109 so 545 +- (5+109) = 545 +- 114 = 659 and 431 and both are primes.
%t A050705 spfQ[n_]:=Module[{s=Total[Times@@@FactorInteger[n]]},!PrimeQ[n] && PrimeQ[ n+s]&&PrimeQ[n-s]]; Select[Range[700],spfQ] (* _Harvey P. Dale_, May 30 2012 *)
%o A050705 (PARI) lista(nn) = {forcomposite(n=2, nn, f = factor(n); sopfr = sum(j=1, #f~, f[j, 1]*f[j, 2]); if (isprime(n+sopfr) && isprime(n-sopfr), print1(n, ", ")););} \\ _Michel Marcus_, Jul 03 2017
%Y A050705 Cf. A050703, A050704, A050706, A050707, A050708, A050709, A050710.
%K A050705 nonn,nice
%O A050705 1,1
%A A050705 _Patrick De Geest_, Aug 15 1999