OFFSET

0,4

COMMENTS

a(n) = (n^2 + 1)^(-1 + Sum_{k=1..n} phi(k)) - f(n) where phi(n) is Euler's totient function, and f(n) is the number of trivial solutions which do not satisfy the equation q_1*x_1 + q_2*x_2 + ... + q_m*x_m = n. Each coefficient is a rational number satisfying the criteria given in the definition, and m = -1 + Sum_{k=1..n} phi(k).

EXAMPLE

For n=3, the number of partitions is equal to the number of nonnegative integer solutions for the equation: (1/2)*x_1 + (1/3)*x_2 + (2/3)*x_3 = 3. The set S of solutions is {[0,1,4], [0,3,3], [0,5,2], [0,7,1], [0,9,0], [2,0,3], [2,2,2], [2,4,1], [2,6,0], [4,1,1], [4,3,0], [6,0,0]}. Therefore, |S| = a(3) = 12.

PROG

(PARI) s(v, n, t) = {if(t==#v, f = n\v[t]; v[t]*f == n, sum(i=0, n\v[t], s(v, n-v[t]*i, t+1)))}

a(n) = {if(n<=2, return(n-1)); my(fractions = List(), q = 0); for(i=2, n, for(j=1, i-1, if(gcd(i, j)==1, listput(fractions, j/i)))); s(fractions, n, 1)} \\ David A. Corneth, Aug 03 2017

CROSSREFS

KEYWORD

nonn

AUTHOR

Joseph Wheat, Aug 03 2017

EXTENSIONS

a(7)-a(17) from Alois P. Heinz, Aug 03 2017

STATUS

approved