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.

Showing 1-2 of 2 results.

A136443 Numbers m such that A102863(m) = 1.

Original entry on oeis.org

1, 3, 5, 7, 8, 9, 10, 11, 13, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78
Offset: 1

Views

Author

Giovanni Teofilatto, Apr 03 2008

Keywords

Comments

Except first term, numbers n such that the sum of the first n primes is composite: Essentially the complement of A013916.

Crossrefs

Programs

A105783 Number of terms among the first n primes that are divisors of the sum of the first n primes.

Original entry on oeis.org

1, 0, 2, 0, 2, 0, 1, 2, 2, 1, 2, 0, 3, 0, 2, 1, 3, 1, 1, 2, 1, 1, 3, 1, 3, 2, 2, 1, 3, 2, 3, 1, 3, 1, 1, 1, 3, 2, 3, 2, 3, 1, 3, 1, 3, 1, 2, 2, 3, 3, 3, 2, 4, 1, 1, 3, 4, 2, 1, 0, 2, 1, 2, 0, 1, 2, 2, 3, 2, 3, 3, 1, 3, 1, 1, 2, 4, 1, 3, 3, 1, 1, 1, 4, 3, 2, 4, 3, 3, 3, 4, 1, 1, 2, 1, 0, 2, 3, 2, 0, 2, 0, 4, 1, 4
Offset: 1

Views

Author

Emeric Deutsch, Apr 19 2005

Keywords

Comments

Sequence inspired by A102863 (Giovanni Teofilatto).

Examples

			a(2)=0 because neither 2 nor 3 is a divisor of 5;
a(5)=2 because exactly two terms from {2,3,5,7,11} are divisors of 2+3+5+7+11=28.
		

Crossrefs

Cf. A102863.

Programs

  • Maple
    with(numtheory): a:=n->nops(factorset(sum(ithprime(k),k=1..n)) intersect {seq(ithprime(j),j=1..n)}): seq(a(n),n=1..130);
    # second Maple program:
    s:= proc(n) option remember; `if`(n<1, 0, ithprime(n)+s(n-1)) end:
    a:= n-> nops(select(x-> x <= ithprime(n), numtheory[factorset](s(n)))):
    seq(a(n), n=1..100);  # Alois P. Heinz, Apr 11 2018
  • Mathematica
    a[n_] := Module[{pp = Prime[Range[n]], s}, s = Total[pp]; Count[pp, p_ /; Divisible[s, p]]];
    Array[a, 105] (* Jean-François Alcover, Jun 19 2018 *)
  • PARI
    a(n) = #select(x->(x <= prime(n)), factor(sum(k=1, n, prime(k)))[,1]); \\ Michel Marcus, Apr 11 2018
Showing 1-2 of 2 results.