|
|
A324640
|
|
Dirichlet inverse of the Doudna sequence, A005940.
|
|
4
|
|
|
1, -2, -3, 0, -5, 6, -9, 0, 2, 10, -15, 0, -25, 18, 3, 0, -11, -4, -21, 0, 19, 30, -45, 0, -24, 50, -60, 0, -125, -6, -81, 0, 77, 22, 57, 0, -55, 42, 87, 0, -77, -38, -105, 0, -78, 90, -135, 0, -40, 48, -81, 0, -245, 120, -75, 0, -217, 250, -375, 0, -625, 162, -150, 0, 233, -154, -39, 0, 205, -114, -99, 0, -91, 110, 174, 0, -5, -174, -189, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
a(1) = 1; for n > 1, a(n) = -Sum_{d|n, d<n} a(d) * A005940(n/d).
a(p) = -A005940(p) for all primes p.
|
|
PROG
|
(PARI)
up_to = 16384;
DirInverse(v) = { my(u=vector(#v)); u[1] = (1/v[1]); for(n=2, #v, u[n] = -sumdiv(n, d, if(d<n, v[n/d]*u[d], 0))); (u) }; \\ Compute the Dirichlet inverse of the sequence given in input vector v.
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t }; \\ From A005940
v324640 = DirInverse(vector(up_to, n, A005940(n)));
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|