

A305901


Filter sequence for all such sequences b, for which b(A006254(k)) = constant for all k >= 3.


6



1, 2, 3, 4, 5, 4, 4, 6, 4, 4, 7, 4, 8, 9, 4, 4, 10, 11, 4, 12, 4, 4, 13, 4, 14, 15, 4, 16, 17, 4, 4, 18, 19, 4, 20, 4, 4, 21, 22, 4, 23, 4, 24, 25, 4, 26, 27, 28, 4, 29, 4, 4, 30, 4, 4, 31, 4, 32, 33, 34, 35, 36, 37, 4, 38, 4, 39, 40, 4, 4, 41, 42, 43, 44, 4, 4, 45, 46, 4, 47, 48, 4, 49, 4, 50, 51, 4, 52, 53, 4, 4, 54, 55, 56, 57, 4, 4, 58, 4, 4, 59, 60, 61
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

Restricted growth sequence transform of A305900(A064216(n)).
For all i, j:
a(i) = a(j) => A278223(i) = A278223(j).
a(i) = a(j) => A253786(i) = A253786(j).


LINKS

Antti Karttunen, Table of n, a(n) for n = 1..100000


FORMULA

For n <= 3, a(n) = n, and for n >= 4, a(n) = 4 if 2n1 is a prime (for all n in A006254[3..] = 4, 6, 7, 9, 10, 12, 15, ...), and for all other n (numbers n such that 2n1 is composite), a(n) = running count from 5 onward.


PROG

(PARI)
up_to = 1000;
partialsums(f, up_to) = { my(v = vector(up_to), s=0); for(i=1, up_to, s += f(i); v[i] = s); (v); }
v_partsums = partialsums(x > isprime(x+x1), up_to);
A305901(n) = if(n<=3, n, if(isprime(n+n1), 4, 3+nv_partsums[n]));


CROSSREFS

Cf. A006254, A064216, A305900.
Cf. also A305902.
Sequence in context: A270434 A204982 A017850 * A291520 A317646 A305748
Adjacent sequences: A305898 A305899 A305900 * A305902 A305903 A305904


KEYWORD

nonn


AUTHOR

Antti Karttunen, Jun 14 2018


STATUS

approved



