login
A267114
Numbers n for which A001222(n) = A267115(n) + A267116(n).
5
1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 33, 34, 35, 36, 38, 39, 40, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 57, 58, 62, 63, 65, 68, 69, 72, 74, 75, 76, 77, 80, 82, 85, 86, 87, 88, 91, 92, 93, 94, 95, 96, 98, 99, 100, 104, 106, 108, 111, 112, 115, 116, 117, 118, 119, 122, 123, 124, 129, 133, 134, 135, 136, 141, 142, 143, 144
OFFSET
1,2
LINKS
EXAMPLE
6 = 2^1 * 3^1 is included as bitwise-or of its exponents is 1 and likewise bitwise-and(1,1) = 1 and 1+1 = A001222(6) = 2, the number of the prime factors of 6 when counted with multiplicity.
12 = 2^2 * 3^1 is included as bitwise-or of its exponents ("10" and "01" in binary) is 3 ("11"), bitwise-and(1,2) = 0 and 3+0 = A001222(12).
60 = 2^2 * 3^1 * 5^1 is NOT included as bitwise-or(2,1,1) = 3, bitwise-and(2,1,1) = 0 and 3+0 < 4 = A001222(60).
MATHEMATICA
{1}~Join~Select[Range@ 144, Function[n, PrimeOmega@ n == BitAnd @@ # + BitOr @@ # &@ Map[Last, FactorInteger@ n]]] (* Michael De Vlieger, Feb 07 2016 *)
PROG
(Scheme, with Antti Karttunen's IntSeq-library)
(define A267114 (ZERO-POS 1 1 (lambda (n) (- (A001222 n) (+ (A267115 n) (A267116 n))))))
(PARI) is(n)=if(n>1, my(f=factor(n)[, 2]); fold(bitand, f) + fold(bitor, f) == vecsum(f), 1) \\ Charles R Greathouse IV, Aug 04 2016
CROSSREFS
Differs from A030231 for the first time at n=118, where A030231(118) = 210 (= 2*3*5*7), which term is missing from this sequence.
a(n+1) differs from A007774 for the first time at n=375, as a(376) = 720 = 2^4 * 3^2 * 5^1.
Cf. A007774 (subsequence).
Sequence in context: A327476 A007774 A030231 * A275665 A056760 A328956
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 03 2016
EXTENSIONS
Erroneous claim corrected by Antti Karttunen, Feb 07 2016
STATUS
approved