[go: up one dir, main page]

login
A368919
Number of terms of A368917 less than 10^n, where A368917 lists the numbers k for which there is no prime p such that p^p divides the arithmetic derivative of A276086(k), where A276086 is the primorial base exp-function.
3
8, 88, 863, 8634, 86407, 864150, 8641439, 86414292
OFFSET
1,1
COMMENTS
Value a(n) / 10^n seems to converge to 1 - lim_{n->oo} (A368911(n) / 10^n).
PROG
(PARI)
A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
A359550(n) = { my(f = factor(n)); prod(k=1, #f~, (f[k, 2]<f[k, 1])); };
A368916(n) = ((n>0)&&A359550(A327860(n)));
tp=10; s=0; for(n=1, 10^10, s+=A368916(n); if(1+n==tp, print1(s, ", "), if(n==tp, tp *= 10)));
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Antti Karttunen, Jan 14 2024
STATUS
approved