|
|
A327858
|
|
Greatest common divisor of the arithmetic derivative and the primorial base exp-function: a(n) = gcd(A003415(n), A276086(n)).
|
|
40
|
|
|
1, 2, 1, 1, 1, 1, 5, 1, 3, 6, 1, 1, 1, 1, 3, 2, 1, 1, 1, 1, 3, 10, 1, 1, 1, 10, 15, 3, 1, 1, 1, 1, 1, 14, 1, 6, 5, 1, 21, 2, 1, 1, 1, 1, 3, 3, 25, 1, 7, 14, 15, 10, 7, 1, 1, 2, 1, 2, 1, 1, 1, 1, 3, 3, 3, 18, 1, 1, 3, 2, 1, 1, 1, 1, 3, 5, 5, 18, 1, 1, 1, 6, 1, 1, 1, 2, 15, 2, 35, 1, 1, 2, 3, 2, 49, 6, 1, 1, 7, 15, 35, 1, 7, 1, 1, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
COMMENTS
|
Sequence contains only terms of A048103.
|
|
LINKS
|
|
|
FORMULA
|
a(p) = 1 for all primes p.
(End)
|
|
MATHEMATICA
|
Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12], f, g}, f[n_] := If[Abs@ n < 2, 0, n Total[#2/#1 & @@@ FactorInteger[Abs@ n]]]; g[n_] := Times @@ Power @@@ # &@ Transpose@ {Prime@ Range@ Length@ #, Reverse@ #} &@ IntegerDigits[n, b]; Array[GCD[f@ #, g@ #] &, 105]] (* Michael De Vlieger, Sep 30 2019 *)
|
|
PROG
|
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
|
|
CROSSREFS
|
Cf. A003415, A048103, A235992, A276086, A327859, A328382, A351234, A354348, A356299, A358669, A359423, A359589 (Dirichlet inverse of a(n)-1).
Cf. A046337 (positions of even terms), A356311 (positions of 1's), A356310 (their characteristic function).
|
|
KEYWORD
|
nonn,base,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|