OFFSET
1,1
COMMENTS
Contains primes of A095078(n) as a subset. Intersection of a(n) with A049445(n) is A117891(n). - R. J. Mathar, Apr 03 2006
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n)<=A117891(n) - R. J. Mathar, Apr 03 2006
a(n) mod A023416(a(n)) = 0. - Reinhard Zumkeller, Nov 22 2007
EXAMPLE
24 is 11000 in binary. This binary representation has three 0's and 3 divides 24. So 24 is in the sequence.
PROG
(C) #include <stdio.h> int main(int argc, char *argv[]) { for(int n=1; n< 500; n++) { int digs=0; int nshifted=n; while(nshifted) { digs += 1- nshifted & 1; nshifted >>= 1; } if ( digs) if( n % digs == 0 ) printf("%d, ", n); } } - R. J. Mathar, Apr 03 2006
(Haskell)
a117890 n = a117890_list !! (n-1)
a117890_list = [x | x <- [1..], let z = a023416 x, z > 0, mod x z == 0]
-- Reinhard Zumkeller, Mar 31 2015
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Leroy Quet, Mar 30 2006
EXTENSIONS
More terms from R. J. Mathar, Apr 03 2006
STATUS
approved