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.

A093890 Number of primes arising as the sum of one or more divisors of n.

This page as a plain text file.
%I A093890 #10 Nov 17 2018 21:50:39
%S A093890 0,2,1,4,1,5,1,6,2,7,1,9,1,5,4,11,1,12,1,13,5,5,1,17,2,5,4,16,1,20,1,
%T A093890 18,4,6,6,24,1,5,5,24,1,24,1,18,11,5,1,30,1,15,3,18,1,30,6,30,5,7,1,
%U A093890 39,1,3,18,31,6,34,1,16,3,34,1,44,1,4,13,16,4,39,1,42,5,5,1,48,5,5,2,41,1,51,2
%N A093890 Number of primes arising as the sum of one or more divisors of n.
%C A093890 a(2^n) = pi(2^(n+1)-1).
%C A093890 Except for n=3 and n=42, it appears that the records occur at the highly abundant numbers A002093. The record values appear to be pi(sigma(n)) for n in A002093, which means that these n are members of A093891. [_T. D. Noe_, Mar 19 2010]
%H A093890 T. D. Noe, <a href="/A093890/b093890.txt">Table of n, a(n) for n = 1..10000</a>
%e A093890 a(4) = 4, the divisors of 4 are 1, 2 and 4.
%e A093890 Primes arising are 2, 3 = 1 + 2, 5 = 1 + 4 and 7 = 1 + 2 + 4.
%t A093890 Do[l = Subsets[Divisors[n]]; l = Union[Map[Plus @@ #&, l]]; Print[Length[Select[l, PrimeQ]]], {n, 100}] (* _Ryan Propper_, Jun 04 2006 *)
%t A093890 CountPrimes[n_] := Module[{d=Divisors[n],t,lim,x}, t=CoefficientList[Product[1+x^i, {i,d}], x]; lim=PrimePi[Length[t]-1]; Count[t[[1+Prime[Range[lim]]]], _?(#>0 &)]]; Table[CountPrimes[n], {n,100}] (* _T. D. Noe_, Mar 19 2010 *)
%Y A093890 Cf. A093891, A093892.
%Y A093890 Cf. A161510 (primes counted with repetition). [_T. D. Noe_, Mar 19 2010]
%K A093890 nonn
%O A093890 1,2
%A A093890 _Amarnath Murthy_, Apr 23 2004
%E A093890 Corrected and extended by _Ryan Propper_, Jun 04 2006