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.

A109580 Primes in A109579.

This page as a plain text file.
%I A109580 #2 Mar 30 2012 17:34:20
%S A109580 3,3,7,7,11,11,11,11,13,13,17,17,17,17,23,23,29,29,29,29,61,61,71,71,
%T A109580 71,71,83,83,83,83,83,83,107,107,107,107,137,137,137,137,137,137,137,
%U A109580 137,149,149,193,193,227,227,227,227,227,227,241,241,241,241,281,281,307
%N A109580 Primes in A109579.
%p A109580 with(numtheory): a:=proc(n) if isprime(sum((pi(k+1)-pi(k))*(ithprime(k+1)-ithprime(k)),k=1..n))=true then sum((pi(k+1)-pi(k))*(ithprime(k+1)-ithprime(k)),k=1..n) else fi end: seq(a(n),n=1..230);
%t A109580 an = Table[(PrimePi[n + 1] - PrimePi[n])*(Prime[n + 1] - Prime[n]), {n, 1, 400}] a[0] = 0; a[n_] := a[n] = a[n - 1] + an[[n]] aa = Flatten[Table[If[PrimeQ[a[n]] == True, a[n], {}], {n, 0, Length[an]}]]
%Y A109580 Cf. A001223, A010051.
%K A109580 nonn
%O A109580 1,1
%A A109580 _Roger L. Bagula_, Jun 29 2005