Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A305817
Number of terms of A091206 <= n; Partial sums of A305816.
3
0, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 11, 11, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14
OFFSET
1,3
FORMULA
a(1) = 0; for n > 1, a(n) = A305816(n-1) + a(n).
For all n >= 1, a(A091206(n)) = n.
PROG
(PARI)
up_to = 65537;
A305816(n) = (isprime(n)&&polisirreducible(Pol(binary(n))*Mod(1, 2)));
partialsums(f, up_to) = { my(v = vector(up_to), s=0); for(i=1, up_to, s += f(i); v[i] = s); (v); }
v305817 = partialsums(A305816, up_to);
A305817(n) = v305817[n];
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 15 2018
STATUS
approved