OFFSET
1,3
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
15 is in the sequence because 15_10=1111_2 and 1+1+1+1=4 is composite.
MATHEMATICA
Select[Range[200], !PrimeQ[DigitCount[#, 2, 1]]&] (* Harvey P. Dale, Jan 21 2013 *)
PROG
(PARI) for(n=0, 200, b=binary(n); if(!isprime(sum(m=1, matsize(b)[2], b[m])), print1(n, ", ")))
(Haskell)
a084345 n = a084345_list !! (n-1)
a084345_list = filter ((== 0) . a010051' . a000120) [0..]
-- Reinhard Zumkeller, Aug 28 2013, Nov 16 2012
(Python)
from math import comb
from sympy import isprime, primerange
def A084345(n):
def f(x):
s = bin(x)[-1:1:-1]
m = x.bit_count()
l = x.bit_length()
c = n-1+isprime(m)
for i in range(l):
j = int(s[i])
if j:
m -= 1
for p in primerange(m, l+1):
c += comb(i, p-m)
return c
m, k = n-1, f(n-1)
while m != k: m, k = k, f(k)
return m # Chai Wah Wu, Jun 08 2026
CROSSREFS
KEYWORD
easy,nonn,base,changed
AUTHOR
Zak Seidov, Jun 22 2003
EXTENSIONS
More terms from Rick L. Shepherd, Jun 23 2003
Term 0 added by Michel Marcus, Aug 26 2013
b-file adjusted by Reinhard Zumkeller, Aug 28 2013
STATUS
approved
