login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A117109 Moebius transform of binomial(n+3, 4). 1
1, 4, 14, 30, 69, 107, 209, 295, 480, 641, 1000, 1209, 1819, 2166, 2976, 3546, 4844, 5379, 7314, 8110, 10402, 11645, 14949, 15890, 20405, 21927, 26910, 29055, 35959, 37108, 46375, 48484, 57890, 61196, 73536, 75027, 91389, 93951, 110096, 114260 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Partial sums of a(n) give A015650(n).

LINKS

Robert Israel, Table of n, a(n) for n = 1..10000

FORMULA

a(n) = |{(x,y,z,w) : 1 <= x <= y <= z <= w <= n, gcd(x,y,z,w,n) = 1}|.

EXAMPLE

a(2)=4 because of the quadruples (1,1,1,1), (1,1,1,2), (1,1,2,2), (1,2,2,2).

MAPLE

b34:= unapply(expand(binomial(n+3, 4)), n):

f:= proc(n) local k; uses numtheory;

add(b34(k)*mobius(n/k), k=divisors(n))

end proc:

map(f, [$1..100]); # Robert Israel, May 24 2019

PROG

(PARI) a(n) = sumdiv(n, k, binomial(k+3, 4)*moebius(n/k)); \\ Michel Marcus, Nov 04 2018

CROSSREFS

Cf. A015650.

Sequence in context: A103779 A049451 A079776 * A317031 A140063 A051409

Adjacent sequences:  A117106 A117107 A117108 * A117110 A117111 A117112

KEYWORD

nonn

AUTHOR

Steve Butler, Apr 18 2006

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 26 04:11 EDT 2019. Contains 324369 sequences. (Running on oeis4.)