OFFSET
0,2
COMMENTS
Doudna-gram for A099308 (numbers whose k-th arithmetic derivative is zero for some k).
The creases in the scatter plot graph (these are more easily seen in Michael De Vlieger's attached png image) occur because in general the expected value of A328308(n) decreases if the maximal exponent (A051903) in the prime factorization of n grows, and because the Doudna sequence (A005940) reorders the natural numbers in such a way that numbers with a high value of A051903 occur towards the end of each range [2^n, 2^(n-1)[. The folding effect is even more pronounced in A328307. - Antti Karttunen, Jan 12 2024
LINKS
Antti Karttunen, Table of n, a(n) for n = 0..65537
Michael De Vlieger, Scatterplot of a(n), n = 0..65537 (x axis labeled instead n-1).
Michael De Vlieger, Log log scatterplot of a(n) n = 0..65537 (x axis labeled instead n-1).
FORMULA
a(0) = 1, and for n > 0, a(n) = a(n-1) + A368907(n).
PROG
(PARI)
up_to = 65537;
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); (t); };
A003415checked(n) = if(n<=1, 0, my(f=factor(n), s=0); for(i=1, #f~, if(f[i, 2]>=f[i, 1], return(0), s += f[i, 2]/f[i, 1])); (n*s));
A328308(n) = if(!n, 1, while(n>1, n = A003415checked(n)); (n));
A368908list(up_to) = { my(v=vector(up_to), s=A368907(0)); for(i=1, up_to, s += A368907(i); v[i] = s); (v); };
v368908 = A368908list(up_to);
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 11 2024
STATUS
approved