|
| |
|
|
A093890
|
|
Number of primes arising as the sum of one or more divisors of n.
|
|
6
| |
|
|
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, 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, 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
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,2
|
|
|
COMMENTS
| a(2^n) = pi(2^(n+1)-1).
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. [From T. D. Noe (noe(AT)sspectra.com), Mar 19 2010]
|
|
|
LINKS
| T. D. Noe, Table of n, a(n) for n=1..10000
|
|
|
EXAMPLE
| a(4) = 4, the divisors of 4 are 1, 2 and 4.
Primes arising are 2, 3= 1+2, 5 = 1+4 and 7 = 1 + 2 + 4.
|
|
|
MATHEMATICA
| Do[l = Subsets[Divisors[n]]; l = Union[Map[Plus @@ #&, l]]; Print[Length[Select[l, PrimeQ]]], {n, 100}] - Ryan Propper (rpropper(AT)stanford.edu), Jun 04 2006
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}] [From T. D. Noe (noe(AT)sspectra.com), Mar 19 2010]
|
|
|
CROSSREFS
| Cf. A093891, A093892.
Cf. A161510 (primes counted with repetition) [From T. D. Noe (noe(AT)sspectra.com), Mar 19 2010]
Sequence in context: A007104 A102627 A088296 * A006306 A083711 A018783
Adjacent sequences: A093887 A093888 A093889 * A093891 A093892 A093893
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Amarnath Murthy (amarnath_murthy(AT)yahoo.com), Apr 23 2004
|
|
|
EXTENSIONS
| Corrected and extended by Ryan Propper (rpropper(AT)stanford.edu), Jun 04 2006
|
| |
|
|
