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.

A272713 Prime powers (p^k, k>=2) that are the sum of consecutive prime numbers.

This page as a plain text file.
%I A272713 #23 May 06 2016 15:13:52
%S A272713 8,49,121,128,169,243,625,841,961,1331,1369,1681,1849,2209,3125,5329,
%T A272713 6241,6859,6889,8192,10201,11449,11881,12167,12769,16384,18769,22801,
%U A272713 24649,26569,32768,36481,39601,44521
%N A272713 Prime powers (p^k, k>=2) that are the sum of consecutive prime numbers.
%C A272713 In other words, prime powers (p^k, k>=2) that are the sum of two or more consecutive prime numbers.
%C A272713 Intersection of A025475 and A034707.
%C A272713 Terms of this sequence are 2^3, 7^2, 11^2, 2^7, 13^2, 3^5, 5^4, 29^2, ...
%H A272713 Charles R Greathouse IV, <a href="/A272713/b272713.txt">Table of n, a(n) for n = 1..4469</a>
%e A272713 8 is a term because 8 = 2^3 = 3 + 5.
%e A272713 49 is a term because 49 = 7^2 = 13 + 17 + 19.
%e A272713 121 is a term because 121 = 11^2 = 37 + 41 + 43.
%o A272713 (PARI) list(lim)=my(v=List(),n=1,p,q,t,s); while(1, t=primes(n++); p=2; q=t[n]; s=vecsum(t); if(s>lim, return(Set(v))); while(s<=lim, if(isprimepower(s)>1, listput(v,s)); q=nextprime(q+1); s+=q-p; p=nextprime(p+1))) \\ _Charles R Greathouse IV_, May 05 2016
%Y A272713 Cf. A025475, A034707, A067377, A050936.
%K A272713 nonn
%O A272713 1,1
%A A272713 _Altug Alkan_, May 05 2016
%E A272713 a(9)-a(34) from _Charles R Greathouse IV_, May 05 2016