OFFSET
1,4
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
David A. Corneth, PARI program
FORMULA
Coefficient of x^(n-1) in expansion of (1/(1-x))*1/Product_{d is prime divisor of n} (1-x^d). - Vladeta Jovovic, Apr 11 2004
PROG
(PARI)
\\ This is for computing just a moderate number of terms:
prime_factors_with1_reversed(n) = vecsort(setunion([1], factor(n)[, 1]~), , 4);
partitions_into_with_trailing_ones(n, parts, from=1) = if(!n, 1, if(#parts<=(from+1), if(#parts == from, 1, (1+(n\parts[from]))), my(s=0); for(i=from, #parts, if(parts[i]<=n, s += partitions_into_with_trailing_ones(n-parts[i], parts, i))); (s)));
A014652(n) = partitions_into_with_trailing_ones(n-1, prime_factors_with1_reversed(n)); \\ Antti Karttunen, Sep 10 2018
(PARI) \\ For an efficient program to compute large numbers of terms, see David A. Corneth's PARI program included in the Links section. - Antti Karttunen, Sep 12 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved