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-10 of 22 results. Next

A321508 Expansion of Product_{k>=1} 1/(1 - x^prime(k))^A056768(k).

Original entry on oeis.org

1, 0, 1, 1, 1, 3, 2, 6, 4, 7, 10, 15, 17, 30, 31, 41, 58, 81, 105, 143, 177, 218, 306, 393, 550, 618, 883, 1024, 1395, 1810, 2372, 2985, 3682, 4762, 6077, 7634, 10160, 12517, 15448, 19820, 24754, 32108, 40085, 50851, 62331, 78548, 98505, 125596, 156565
Offset: 0

Views

Author

Ilya Gutkovskiy, Nov 11 2018

Keywords

Comments

a(n) is the number of partitions of n into prime parts prime(k) of A056768(k) kinds.

Examples

			a(7) = 6 because we have [{7}], [{5, 2}], [{5}, {2}], [{3, 2, 2}], [{3, 2}, {2}] and [{3}, {2}, {2}].
		

Crossrefs

Programs

  • Mathematica
    b[n_] := b[n] = SeriesCoefficient[Product[1/(1 - x^Prime[k]), {k, 1, n}], {x, 0, Prime[n]}]; a[n_] := a[n] = SeriesCoefficient[Product[1/(1 - x^Prime[k])^b[k], {k, 1, n}], {x, 0, n}]; Table[a[n], {n, 0, 48}]

Formula

G.f.: Product_{k>=1} 1/(1 - x^A000040(k))^A000607(A000040(k)).

A000607 Number of partitions of n into prime parts.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 9, 10, 12, 14, 17, 19, 23, 26, 30, 35, 40, 46, 52, 60, 67, 77, 87, 98, 111, 124, 140, 157, 175, 197, 219, 244, 272, 302, 336, 372, 413, 456, 504, 557, 614, 677, 744, 819, 899, 987, 1083, 1186, 1298, 1420, 1552, 1695, 1850, 2018, 2198, 2394, 2605, 2833, 3079, 3344
Offset: 0

Views

Author

Keywords

Comments

a(n) gives the number of values of k such that A001414(k) = n. - Howard A. Landman, Sep 25 2001
Let W(n) = {prime p: There is at least one number m whose spf is p, and sopfr(m) = n}. Let V(n,p) = {m: sopfr(m) = n, p belongs to W(n)}. Then a(n) = sigma(|V(n,p)|). E.g.: W(10) = {2,3,5}, V(10,2) = {30,32,36}, V(10,3) = {21}, V(10,5) = {25}, so a(10) = 3+1+1 = 5. - David James Sycamore, Apr 14 2018
From Gus Wiseman, Jan 18 2020: (Start)
Also the number of integer partitions such that the sum of primes indexed by the parts is n. For example, the sum of primes indexed by the parts of the partition (3,2,1,1) is prime(3)+prime(2)+prime(1)+prime(1) = 12, so (3,2,1,1) is counted under a(12). The a(2) = 1 through a(14) = 10 partitions are:
1 2 11 3 22 4 32 41 33 5 43 6 44
21 111 31 221 222 42 322 331 51 52
211 1111 311 321 411 421 332 431
2111 2211 2221 2222 422 3222
11111 3111 3211 3221 3311
21111 22111 4111 4211
111111 22211 22221
31111 32111
211111 221111
1111111
(End)

Examples

			n = 10 has a(10) = 5 partitions into prime parts: 10 = 2 + 2 + 2 + 2 + 2 = 2 + 2 + 3 + 3 = 2 + 3 + 5 = 3 + 7 = 5 + 5.
n = 15 has a(15) = 12 partitions into prime parts: 15 = 2 + 2 + 2 + 2 + 2 + 2 + 3 = 2 + 2 + 2 + 3 + 3 + 3 = 2 + 2 + 2 + 2 + 2 + 5 = 2 + 2 + 2 + 2 + 7 = 2 + 2 + 3 + 3 + 5 = 2 + 3 + 5 + 5 = 2 + 3 + 3 + 7 = 2 + 2 + 11 = 2 + 13 = 3 + 3 + 3 + 3 + 3 = 3 + 5 + 7 = 5 + 5 + 5.
		

References

  • R. Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; see p. 203.
  • Mohammad K. Azarian, A Generalization of the Climbing Stairs Problem, Mathematics and Computer Education, Vol. 31, No. 1, pp. 24-28, Winter 1997. MathEduc Database (Zentralblatt MATH, 1997c.01891).
  • B. C. Berndt and B. M. Wilson, Chapter 5 of Ramanujan's second notebook, pp. 49-78 of Analytic Number Theory (Philadelphia, 1980), Lect. Notes Math. 899, 1981, see Entry 29.
  • D. M. Burton, Elementary Number Theory, 5th ed., McGraw-Hill, 2002.
  • L. M. Chawla and S. A. Shad, On a trio-set of partition functions and their tables, J. Natural Sciences and Mathematics, 9 (1969), 87-96.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

G.f. = 1 / g.f. for A046675. See A046113 for the ordered (compositions) version.
Row sums of array A116865 and of triangle A261013.
Column sums of A331416.
Partitions whose Heinz number is divisible by their sum of primes are A330953.
Partitions of whose sum of primes is divisible by their sum are A331379.

Programs

  • Haskell
    a000607 = p a000040_list where
       p _      0 = 1
       p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
    -- Reinhard Zumkeller, Aug 05 2012
    
  • Magma
    [1] cat [#RestrictedPartitions(n,{p:p in PrimesUpTo(n)}): n in [1..100]]; // Marius A. Burtea, Jan 02 2019
  • Maple
    with(gfun):
    t1:=mul(1/(1-q^ithprime(n)),n=1..51):
    t2:=series(t1,q,50):
    t3:=seriestolist(t2); # fixed by Vaclav Kotesovec, Sep 14 2014
  • Mathematica
    CoefficientList[ Series[1/Product[1 - x^Prime[i], {i, 1, 50}], {x, 0, 50}], x]
    f[n_] := Length@ IntegerPartitions[n, All, Prime@ Range@ PrimePi@ n]; Array[f, 57] (* Robert G. Wilson v, Jul 23 2010 *)
    Table[Length[Select[IntegerPartitions[n],And@@PrimeQ/@#&]],{n,0,60}] (* Harvey P. Dale, Apr 22 2012 *)
    a[n_] := a[n] = If[PrimeQ[n], 1, 0]; c[n_] := c[n] = Plus @@ Map[# a[#] &, Divisors[n]]; b[n_] := b[n] = (c[n] + Sum[c[k] b[n - k], {k, 1, n - 1}])/n; Table[b[n], {n, 1, 20}] (* Thomas Vogler, Dec 10 2015: Uses Euler transform, caches computed values, faster than IntegerPartitions[] function. *)
    nmax = 100; pmax = PrimePi[nmax]; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 0; poly[[3]] = -1; Do[p = Prime[k]; Do[poly[[j + 1]] -= poly[[j + 1 - p]], {j, nmax, p, -1}];, {k, 2, pmax}]; s = Sum[poly[[k + 1]]*x^k, {k, 0, Length[poly] - 1}]; CoefficientList[Series[1/s, {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 11 2021 *)
  • PARI
    N=66;x='x+O('x^N); Vec(1/prod(k=1,N,1-x^prime(k))) \\ Joerg Arndt, Sep 04 2014
    
  • Python
    from sympy import primefactors
    l = [1, 0]
    for n in range(2, 101):
        l.append(sum(sum(primefactors(k)) * l[n - k] for k in range(1, n + 1)) // n)
    l  # Indranil Ghosh, Jul 13 2017
    
  • Sage
    [Partitions(n, parts_in=prime_range(n + 1)).cardinality() for n in range(100)]  # Giuseppe Coppoletta, Jul 11 2016
    

Formula

Asymptotically a(n) ~ exp(2 Pi sqrt(n/log n) / sqrt(3)) (Ayoub).
a(n) = (1/n)*Sum_{k=1..n} A008472(k)*a(n-k). - Vladeta Jovovic, Aug 27 2002
G.f.: 1/Product_{k>=1} (1-x^prime(k)).
See the partition arrays A116864 and A116865.
From Vaclav Kotesovec, Sep 15 2014 [Corrected by Andrey Zabolotskiy, May 26 2017]: (Start)
It is surprising that the ratio of the formula for log(a(n)) to the approximation 2 * Pi * sqrt(n/(3*log(n))) exceeds 1. For n=20000 the ratio is 1.00953, and for n=50000 (using the value from Havermann's tables) the ratio is 1.02458, so the ratio is increasing. See graph above.
A more refined asymptotic formula is found by Vaughan in Ramanujan J. 15 (2008), pp. 109-121, and corrected by Bartel et al. (2017): log(a(n)) = 2*Pi*sqrt(n/(3*log(n))) * (1 - log(log(n))/(2*log(n)) + O(1/log(n))).
See Bartel, Bhaduri, Brack, Murthy (2017) for a more complete asymptotic expansion. (End)
G.f.: 1 + Sum_{i>=1} x^prime(i) / Product_{j=1..i} (1 - x^prime(j)). - Ilya Gutkovskiy, May 07 2017
a(n) = A184198(n) + A184199(n). - Vaclav Kotesovec, Jan 11 2021

A100118 Numbers whose sum of prime factors is prime (counted with multiplicity).

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 12, 13, 17, 19, 22, 23, 28, 29, 31, 34, 37, 40, 41, 43, 45, 47, 48, 52, 53, 54, 56, 58, 59, 61, 63, 67, 71, 73, 75, 76, 79, 80, 82, 83, 88, 89, 90, 96, 97, 99, 101, 103, 104, 107, 108, 109, 113, 117, 118, 127, 131, 136, 137, 139, 142, 147, 148, 149
Offset: 1

Views

Author

Carlos Alves, Dec 26 2004

Keywords

Comments

Numbers n such that integer log of n is a prime number.
As in A001414, denote sopfr(n) the integer log of n. Since sopfr(p)=p, the sequence includes all prime numbers.
See A046363 for the analog excluding prime numbers. - Hieronymus Fischer, Oct 20 2007
These numbers may be arranged in a family of posets of triangles of multiarrows (see link and example). - Gus Wiseman, Sep 14 2016

Examples

			40 = 2^3*5 and 2*3 + 5 = 11 is a prime number.
These numbers correspond to multiarrows in the multiorder of partitions of prime numbers into prime parts. For example: 2:2<=(2), 3:3<=(3), 6:5<=(2,3), 5:5<=(5), 12:7<=(2,2,3), 10:7<=(2,5), 7:7<=(7), 48:11<=(2,2,2,2,3), 52:11<=(2,3,3,3), 40:11<=(2,2,2,5), 45:11<=(3,3,5), 28:11<=(2,2,7), 11:11<=(11). - _Gus Wiseman_, Sep 14 2016
		

Crossrefs

Programs

  • Maple
    for n from 1 to 200 do
        if isprime(A001414(n)) then
            printf("%d,",n);
        end if;
    end do: # R. J. Mathar, Sep 09 2015
  • Mathematica
    L = {}; Do[ww = Transpose[FactorInteger[k]];w = ww[[1]].ww[[2]]; If[PrimeQ[w], AppendTo[L, k]], {k, 2, 500}];L
    Select[Range[150], PrimeQ[Total[Times @@@ FactorInteger[#]]] &] (* Jayanta Basu, Aug 11 2013 *)
  • PARI
    is(n)=my(f=factor(n)); isprime(sum(i=1,#f~,f[i,1]*f[i,2])) \\ Charles R Greathouse IV, Sep 21 2013

A070215 Number of ways to write the n-th prime as a sum of distinct primes.

Original entry on oeis.org

1, 1, 2, 2, 1, 2, 2, 3, 5, 7, 9, 11, 14, 15, 19, 26, 35, 39, 50, 61, 67, 87, 102, 130, 178, 204, 224, 257, 278, 320, 522, 595, 724, 776, 1064, 1136, 1364, 1634, 1836, 2192, 2601, 2761, 3645, 3863, 4294, 4549, 6262, 8558, 9453, 9964, 11001, 12774, 13438
Offset: 1

Views

Author

Lekraj Beedassy, May 07 2002

Keywords

Examples

			With the 10th prime 29, for instance, we have a(10)=7 distinct-prime partitions, viz. 29 = 2 + 3 + 7 + 17 = 2 + 3 + 5 + 19 = 2 + 3 + 11 + 13 = 3 + 7 + 19 = 5 + 7 + 17 = 5 + 11 + 13.
		

Crossrefs

Cf. A000586, A056768 (parts may repeat).

Programs

  • Haskell
    a070215 = a000586 . a000040  -- Reinhard Zumkeller, Aug 05 2012
  • Mathematica
    nn = PrimePi[300]; t = CoefficientList[Series[Product[(1 + x^Prime[k]), {k, nn}], {x, 0, Prime[nn]}], x]; t[[1 + Prime[Range[nn]]]] (* T. D. Noe, Nov 13 2013 *)

Formula

a(n) = A000586(prime(n)). - R. J. Mathar, Apr 30 2007

Extensions

More terms from Naohiro Nomoto and Don Reble, May 11 2002
Offset in b-file corrected by N. J. A. Sloane, Aug 31 2009

A338903 Number of integer partitions of the n-th squarefree semiprime into squarefree semiprimes.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 3, 3, 5, 4, 6, 5, 12, 14, 19, 22, 27, 36, 38, 51, 77, 86, 128, 141, 163, 163, 207, 233, 259, 260, 514, 657, 813, 983, 1010, 1215, 1255, 1720, 2112, 2256, 3171, 3370, 3499, 3864, 4103, 6292, 7313, 7620, 8374, 10650, 17579, 18462, 23034, 25180
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2020

Keywords

Comments

A squarefree semiprime (A006881) is a product of any two distinct prime numbers.

Examples

			The a(n) partitions for n = 1, 5, 7, 9, 10, 11, 13:
  6  21    26       34          35        38           46
     15,6  14,6,6   22,6,6      21,14     26,6,6       34,6,6
           10,10,6  14,14,6     15,14,6   22,10,6      26,14,6
                    14,10,10    15,10,10  14,14,10     21,15,10
                    10,6,6,6,6            14,6,6,6,6   22,14,10
                                          10,10,6,6,6  26,10,10
                                                       15,15,10,6
                                                       22,6,6,6,6
                                                       14,14,6,6,6
                                                       14,10,10,6,6
                                                       10,10,10,10,6
                                                       10,6,6,6,6,6,6
		

Crossrefs

A002100 counts partitions into squarefree semiprimes.
A056768 uses primes instead of squarefree semiprimes.
A101048 counts partitions into semiprimes.
A338902 is the not necessarily squarefree version.
A339113 includes the Heinz numbers of these partitions.
A001358 lists semiprimes, with odd and even terms A046315 and A100484.
A006881 lists squarefree semiprimes, with odd and even terms A046388 and A100484.
A320656 counts factorizations into squarefree semiprimes.
A338898/A338912/A338913 give prime indices of semiprimes, with sum/difference/product A176504/A176506/A087794.
A338899, A270650, and A270652 give the prime indices of squarefree semiprimes.

Programs

  • Mathematica
    nn=100;
    sqs=Select[Range[nn],SquareFreeQ[#]&&PrimeOmega[#]==2&];
    Table[Length[IntegerPartitions[n,All,sqs]],{n,sqs}]

Formula

a(n) = A002100(A006881(n)).

A338902 Number of integer partitions of the n-th semiprime into semiprimes.

Original entry on oeis.org

1, 1, 1, 2, 3, 2, 4, 7, 7, 10, 17, 25, 21, 34, 34, 73, 87, 103, 149, 176, 206, 281, 344, 479, 725, 881, 1311, 1597, 1742, 1841, 2445, 2808, 3052, 3222, 6784, 9298, 11989, 14533, 15384, 17414, 18581, 19680, 28284, 35862, 38125, 57095, 60582, 64010, 71730, 76016
Offset: 1

Views

Author

Gus Wiseman, Nov 24 2020

Keywords

Comments

A semiprime (A001358) is a product of any two prime numbers.

Examples

			The a(1) = 1 through a(33) = 17 partitions of 4, 6, 9, 10, 14, 15, 21, 22, 25, 26, 33, where A-Z = 10-35:
  4  6  9  A   E    F   L     M      P      Q       X
           64  A4   96  F6    994    FA     M4      EA9
               644      966   A66    L4     AA6     F99
                        9444  E44    A96    E66     FE4
                              6664   F64    9944    L66
                              A444   9664   A664    P44
                              64444  94444  E444    9996
                                            66644   AA94
                                            A4444   E964
                                            644444  F666
                                                    FA44
                                                    L444
                                                    96666
                                                    A9644
                                                    F6444
                                                    966444
                                                    9444444
		

Crossrefs

A002100 counts partitions into squarefree semiprimes.
A056768 uses primes instead of semiprimes.
A101048 counts partitions into semiprimes.
A338903 is the squarefree version.
A339112 includes the Heinz numbers of these partitions.
A001358 lists semiprimes, with odd and even terms A046315 and A100484.
A037143 lists primes and semiprimes.
A084126 and A084127 give the prime factors of semiprimes.
A320655 counts factorizations into semiprimes.
A338898/A338912/A338913 give prime indices of semiprimes, with sum/difference/product A176504/A176506/A087794.
A338899/A270650/A270652 give prime indices of squarefree semiprimes.

Programs

  • Mathematica
    nn=100;Table[Length[IntegerPartitions[n,All,Select[Range[nn],PrimeOmega[#]==2&]]],{n,Select[Range[nn],PrimeOmega[#]==2&]}]

Formula

a(n) = A101048(A001358(n)).

A307610 Number of partitions of prime(n) into consecutive primes.

Original entry on oeis.org

1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 3, 1, 1, 2, 2, 1, 2, 2, 1, 1, 3, 1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 3, 3, 2, 3, 1, 1, 2, 1, 1, 3, 1, 2, 2, 2, 1, 3, 1, 1, 1, 5, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 1, 3, 1, 1, 2, 2, 1, 3, 2, 2, 2, 1, 2, 1, 2, 2, 3, 2, 2, 1, 1, 2
Offset: 1

Views

Author

Ilya Gutkovskiy, Apr 18 2019

Keywords

Comments

a(n) - 1 = number of partitions of prime(n) into two or more consecutive primes. - Ray Chandler, Sep 26 2023

Examples

			prime(13) = 41 = 2 + 3 + 5 + 7 + 11 + 13 = 11 + 13 + 17, so a(13) = 3.
		

Crossrefs

Formula

a(n) = [x^prime(n)] Sum_{i>=1} Sum_{j>=i} Product_{k=i..j} x^prime(k).
a(n) = A054845(A000040(n)).

A316092 Heinz numbers of integer partitions of prime numbers into prime parts.

Original entry on oeis.org

3, 5, 11, 15, 17, 31, 33, 41, 45, 59, 67, 83, 93, 109, 127, 153, 157, 177, 179, 191, 211, 241, 275, 277, 283, 297, 327, 331, 353, 367, 369, 375, 401, 405, 425, 431, 459, 461, 509, 537, 547, 563, 587, 599, 603, 605, 617, 709, 739, 773, 775, 797, 825, 831, 837
Offset: 1

Views

Author

Gus Wiseman, Jun 24 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of integer partitions of prime numbers into prime parts together with their Heinz numbers begins:
   3: (2)
   5: (3)
  11: (5)
  15: (2,3)
  17: (7)
  31: (11)
  33: (2,5)
  41: (13)
  45: (2,2,3)
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[900],And[PrimeQ[Total[primeMS[#]]],And@@PrimeQ/@primeMS[#]]&]

A316154 Number of integer partitions of prime(n) into a prime number of prime parts.

Original entry on oeis.org

0, 0, 1, 2, 3, 5, 9, 12, 19, 39, 50, 93, 136, 166, 239, 409, 682, 814, 1314, 1774, 2081, 3231, 4272, 6475, 11077, 14270, 16265, 20810, 23621, 30031, 68251, 85326, 118917, 132815, 226097, 251301, 342448, 463940, 565844, 759873, 1015302, 1117708, 1787452, 1961624
Offset: 1

Views

Author

Gus Wiseman, Jun 25 2018

Keywords

Examples

			The a(7) = 9 partitions of 17 into a prime number of prime parts: (13,2,2), (11,3,3), (7,7,3), (7,5,5), (7,3,3,2,2), (5,5,3,2,2), (5,3,3,3,3), (5,2,2,2,2,2,2), (3,3,3,2,2,2,2).
		

Crossrefs

Programs

  • Maple
    b:= proc(n, p, c) option remember; `if`(n=0 or p=2,
          `if`(n::even and isprime(c+n/2), 1, 0),
          `if`(p>n, 0, b(n-p, p, c+1))+b(n, prevprime(p), c))
        end:
    a:= n-> b(ithprime(n)$2, 0):
    seq(a(n), n=1..50);  # Alois P. Heinz, Jun 26 2018
  • Mathematica
    Table[Length[Select[IntegerPartitions[Prime[n]],And[PrimeQ[Length[#]],And@@PrimeQ/@#]&]],{n,20}]
    (* Second program: *)
    b[n_, p_, c_] := b[n, p, c] = If[n == 0 || p == 2, If[EvenQ[n] && PrimeQ[c + n/2], 1, 0], If[p>n, 0, b[n - p, p, c + 1]] + b[n, NextPrime[p, -1], c]];
    a[n_] := b[Prime[n], Prime[n], 0];
    Array[a, 50] (* Jean-François Alcover, May 20 2021, after Alois P. Heinz *)
  • PARI
    seq(n)={my(p=vector(n,k,prime(k))); my(v=Vec(1/prod(k=1, n, 1 - x^p[k]*y + O(x*x^p[n])))); vector(n, k, sum(i=1, k, polcoeff(v[1+p[k]], p[i])))} \\ Andrew Howroyd, Jun 26 2018

Formula

a(n) = A085755(A000040(n)). - Alois P. Heinz, Jun 26 2018

Extensions

Terms a(21) and beyond from Andrew Howroyd, Jun 26 2018

A316185 Number of strict integer partitions of the n-th prime into a prime number of prime parts.

Original entry on oeis.org

0, 0, 1, 1, 0, 1, 0, 2, 2, 3, 5, 5, 6, 8, 10, 13, 18, 20, 26, 32, 34, 45, 54, 66, 90, 106, 117, 135, 142, 165, 269, 311, 375, 398, 546, 579, 689, 823, 938, 1107, 1301, 1352, 1790, 1850, 2078, 2153, 2878, 3811, 4241, 4338, 4828, 5495, 5637, 7076, 8000, 9032
Offset: 1

Views

Author

Gus Wiseman, Jun 25 2018

Keywords

Examples

			The a(14) = 8 partitions of 43 into a prime number of distinct prime parts: (41,2), (31,7,5), (29,11,3), (23,17,3), (23,13,7), (19,17,7), (19,13,11), (17,11,7,5,3).
		

Crossrefs

Programs

  • Maple
    h:= proc(n) option remember; `if`(n=0, 0,
         `if`(isprime(n), n, h(n-1)))
        end:
    b:= proc(n, i, c) option remember; `if`(n=0,
          `if`(isprime(c), 1, 0), `if`(i<2, 0, b(n, h(i-1), c)+
          `if`(i>n, 0, b(n-i, h(min(n-i, i-1)), c+1))))
        end:
    a:= n-> b(ithprime(n)$2, 0):
    seq(a(n), n=1..56);  # Alois P. Heinz, May 26 2021
  • Mathematica
    Table[Length[Select[IntegerPartitions[Prime[n]],And[UnsameQ@@#,PrimeQ[Length[#]],And@@PrimeQ/@#]&]],{n,10}]
    (* Second program: *)
    h[n_] := h[n] = If[n == 0, 0, If[PrimeQ[n], n, h[n - 1]]];
    b[n_, i_, c_] := b[n, i, c] = If[n == 0,
         If[PrimeQ[c], 1, 0], If[i < 2, 0, b[n, h[i - 1], c] +
         If[i > n, 0, b[n - i, h[Min[n - i, i - 1]], c + 1]]]];
    a[n_] := b[Prime[n], Prime[n], 0];
    Array[a, 56] (* Jean-François Alcover, Jun 11 2021, after Alois P. Heinz *)
  • PARI
    seq(n)={my(p=vector(n, k, prime(k))); my(v=Vec(prod(k=1, n, 1 + x^p[k]*y + O(x*x^p[n])))); vector(n, k, sum(i=1, k, polcoeff(v[1+p[k]], p[i])))} \\ Andrew Howroyd, Jun 26 2018

Formula

a(n) = A045450(A000040(n)).

Extensions

More terms from Alois P. Heinz, Jun 26 2018
Showing 1-10 of 22 results. Next