|
| |
|
|
A077608
|
|
Number of compositions of n into twin primes (i.e. primes that are members of a twin prime pair, like 3,5,7,11,13, but not 2 or 23).
|
|
2
| |
|
|
1, 0, 0, 1, 0, 1, 1, 1, 2, 1, 3, 4, 3, 7, 7, 8, 14, 15, 21, 28, 33, 47, 58, 75, 103, 125, 167, 220, 275, 370, 474, 610, 806, 1028, 1347, 1752, 2253, 2954, 3812, 4944, 6451, 8329, 10841, 14077, 18226, 23720, 30745, 39903, 51857, 67214, 87313, 113340, 147017, 190974
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,9
|
|
|
LINKS
| Alois P. Heinz, Table of n, a(n) for n = 0..1000
P. Flajolet, Publications
|
|
|
FORMULA
| A77608 := proc(n) coeff(series(1/(1-add(z^numtheory[ithprime](j)* subs([true=1, false=0], evalb(isprime(ithprime(j)-2) or isprime(ithprime(j)+2))), j=2..n+2)), z=0, n+1), z, n): end;
|
|
|
EXAMPLE
| a(15) = 8 since 15=11+7=7+11=5+13=13+5=3+5+7=3+7+5=5+3+7=5+7+3=7+3+5=7+5+3 and 3,5,7,11 belong to twin pairs.
|
|
|
MAPLE
| A077608 := proc(n) coeff(series(1/(1-add(z^numtheory[ithprime](j)* subs([true=1, false=0], evalb(isprime(ithprime(j)-2) or isprime(ithprime(j)+2))), j=2..n+2)), z=0, n+1), z, n): end;
|
|
|
CROSSREFS
| Cf. A002124, A023360.
Sequence in context: A022466 A144254 A133310 * A002124 A097564 A128270
Adjacent sequences: A077605 A077606 A077607 * A077609 A077610 A077611
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Philippe Flajolet (Philippe.Flajolet(AT)inria.fr), Nov 11 2002
|
| |
|
|
