login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A116661
Integers in both sequences A114522 and A063989.
2
4, 6, 8, 9, 10, 12, 18, 20, 22, 25, 27, 32, 34, 44, 48, 49, 50, 58, 68, 72, 80, 82, 108, 116, 118, 121, 125, 128, 142, 162, 164, 165, 169, 176, 192, 200, 202, 214, 236, 242, 243, 272, 273, 274, 284, 288, 289, 298, 320, 343, 345, 358, 361, 382, 385, 394, 399, 404
OFFSET
1,1
LINKS
EXAMPLE
20 = 2^2 *5^1. Both the number of prime divisors (counted with multiplicity), 2+1 = 3 and the sum of the distinct prime divisors, 2+5 = 7, are primes. So 20 is in the sequence.
MATHEMATICA
Select[Range[500], AllTrue[{PrimeOmega[#], Total[FactorInteger[#][[All, 1]]]}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Feb 03 2019 *)
PROG
(Magma) f:=func<n|&+[p[2]: p in Factorization(n)]>; [k:k in [2..450]| IsPrime(f(k)) and IsPrime(&+PrimeDivisors(k))]; // Marius A. Burtea, Nov 14 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Feb 21 2006
EXTENSIONS
More terms from Robert Gerbicz, Jun 09 2007
STATUS
approved