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.

A271368 Number of ways to write n as the sum of distinct super-primes (A006450).

This page as a plain text file.
%I A271368 #16 Apr 07 2016 16:20:01
%S A271368 0,0,1,0,1,0,0,1,0,0,1,0,0,1,0,1,1,0,1,1,0,1,0,0,1,0,0,1,0,0,2,0,1,1,
%T A271368 0,2,0,0,1,0,1,1,0,1,1,1,1,1,1,1,1,1,1,0,1,1,1,1,2,1,1,2,1,2,0,1,3,0,
%U A271368 1,2,0,3,1,1,3,1,2,2,1,1,2,0,3,2,0,3,2,2,2,2,2,3,1,2,3,0,2,3,1,4
%N A271368 Number of ways to write n as the sum of distinct super-primes (A006450).
%C A271368 a(n) > 0 for n > 96 (cf. Dressler, Parker, 1975).
%H A271368 R. E. Dressler and S. T. Parker, <a href="http://dx.doi.org/10.1145/321892.321900">Primes with a Prime Subscript</a>, Journal of the ACM, Vol. 22, No. 3 (1975), 380-381.
%H A271368 Wikipedia, <a href="https://en.wikipedia.org/wiki/Super-prime">Super-prime</a>
%F A271368 G.f.: prod(k>=1, 1 + x^A006450(k) ). [_Joerg Arndt_, Apr 06 2016]
%e A271368 There are two ways to write 31 as the sum of distinct super-primes: 31 (a single summand, as 31 is itself a super-prime) and 17 + 11 + 3 (three summands), so a(31) = 2.
%o A271368 (PARI) isokp(pt) = {for (k=1, #pt, if (! isprime(pt[k]) || !isprime(primepi(pt[k])), return (0));); #pt == #Set(pt);}
%o A271368 a(n) = {if (n < 3, return (0)); nb = 0; forpart(pt = n, if (isokp(pt), nb++), [3, n]); nb;} \\ _Michel Marcus_, Apr 06 2016
%Y A271368 Cf. A000586, A006450.
%K A271368 nonn
%O A271368 1,31
%A A271368 _Felix Fröhlich_, Apr 05 2016
%E A271368 More terms from _Michel Marcus_, Apr 06 2016