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.

A116661 Integers in both sequences A114522 and A063989.

This page as a plain text file.
%I A116661 #18 Sep 08 2022 08:45:24
%S A116661 4,6,8,9,10,12,18,20,22,25,27,32,34,44,48,49,50,58,68,72,80,82,108,
%T A116661 116,118,121,125,128,142,162,164,165,169,176,192,200,202,214,236,242,
%U A116661 243,272,273,274,284,288,289,298,320,343,345,358,361,382,385,394,399,404
%N A116661 Integers in both sequences A114522 and A063989.
%H A116661 Amiram Eldar, <a href="/A116661/b116661.txt">Table of n, a(n) for n = 1..10000</a>
%e A116661 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.
%t A116661 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 *)
%o A116661 (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
%Y A116661 Cf. A114522, A063989, A116660.
%K A116661 nonn
%O A116661 1,1
%A A116661 _Leroy Quet_, Feb 21 2006
%E A116661 More terms from _Robert Gerbicz_, Jun 09 2007