A334864 a(n) = A064097(A003961(n)) - A064097(n).
0, 1, 1, 2, 1, 2, 1, 3, 2, 2, 0, 3, 0, 2, 2, 4, 1, 3, 1, 3, 2, 1, 0, 4, 2, 1, 3, 3, 0, 3, 0, 5, 1, 2, 2, 4, 0, 2, 1, 4, 1, 3, 1, 2, 3, 1, -1, 5, 2, 3, 2, 2, 1, 4, 1, 4, 2, 1, -1, 4, 1, 1, 3, 6, 1, 2, 0, 3, 1, 3, -1, 5, 1, 1, 3, 3, 1, 2, 0, 5, 4, 2, 0, 4, 2, 2, 1, 3, -1, 4, 1, 2, 1, 0, 2, 6, 1, 3, 2, 4, 0, 3, 1, 3, 3
Offset: 1
Links
- Antti Karttunen, Table of n, a(n) for n = 1..65537
Programs
-
PARI
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); A064097(n) = if(1==n,0,1+A064097(n-(n/vecmin(factor(n)[,1])))); A334864(n) = (A064097(A003961(n))-A064097(n)); \\ Or alternatively as: A334864(n) = { my(f=factor(n)); sum(k=1, #f~, f[k, 2]*(A064097(prime(1+primepi(f[k, 1])))-A064097(f[k,1]))); };
Comments