|
|
A046643
|
|
From square root of Riemann zeta function: form Dirichlet series Sum b_n/n^s whose square is zeta function; sequence gives numerator of b_n.
|
|
23
|
|
|
1, 1, 1, 3, 1, 1, 1, 5, 3, 1, 1, 3, 1, 1, 1, 35, 1, 3, 1, 3, 1, 1, 1, 5, 3, 1, 5, 3, 1, 1, 1, 63, 1, 1, 1, 9, 1, 1, 1, 5, 1, 1, 1, 3, 3, 1, 1, 35, 3, 3, 1, 3, 1, 5, 1, 5, 1, 1, 1, 3, 1, 1, 3, 231, 1, 1, 1, 3, 1, 1, 1, 15, 1, 1, 3, 3, 1, 1, 1, 35, 35, 1, 1, 3, 1, 1, 1, 5, 1, 3
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Multiplicative with a(p^n) = A001790(n).
(End)
|
|
EXAMPLE
|
b_1, b_2, ... = 1, 1/2, 1/2, 3/8, 1/2, 1/4, 1/2, 5/16, 3/8, 1/4, 1/2, 3/16, ...
|
|
MAPLE
|
b := proc(n) option remember; local c, i, t1; if n = 1 then 1 else c := 1; t1 := divisors(n);
for i from 2 to nops(t1)-1 do c := c-b(t1[ i ])*b(n/t1[ i ]); od; c/2; fi; end;
|
|
MATHEMATICA
|
b[1] = 1; b[n_] := b[n] = (dn = Divisors[n]; c = 1;
Do[c = c - b[dn[[i]]]*b[n/dn[[i]]], {i, 2, Length[dn] - 1}]; c/2); a[n_] := Numerator[b[n]]; a /@ Range[90] (* Jean-François Alcover, Apr 04 2011, after Maple version *)
|
|
PROG
|
(PARI)
A046643perA046644(n) = { my(c=1); if(1==n, c, fordiv(n, d, if((d>1)&&(d<n), c -= (A046643perA046644(d)*A046643perA046644(n/d)))); (c/2)); }
(Scheme)
(define (A046643 n) (numerator (A046643perA046644 n)))
(define (A046644 n) (denominator (A046643perA046644 n)))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy,frac,nice,mult
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|