A093688 -id:A093688 - cp's OEIS frontend

A129771 Evil odd numbers.

3, 5, 9, 15, 17, 23, 27, 29, 33, 39, 43, 45, 51, 53, 57, 63, 65, 71, 75, 77, 83, 85, 89, 95, 99, 101, 105, 111, 113, 119, 123, 125, 129, 135, 139, 141, 147, 149, 153, 159, 163, 165, 169, 175, 177, 183, 187, 189, 195, 197, 201, 207, 209, 215, 219, 221, 225, 231, 235

Offset: 1

Tanya Khovanova, May 16 2007

A heuristic argument suggests that, as n tends to infinity, a(n)/n converges to 4. - Stefan Steinerberger, May 17 2007
These numbers may be called primitive evil numbers because every evil number is a power of 2 multiplied by one of these numbers. Note that the difference between consecutive terms is either 2, 4, or 6. - T. D. Noe, Jun 06 2007
If m is in the sequence, then so is 2m-1 because in binary, m is x1 and 2m-1 is x01. Presumably the numbers that generate the whole sequence by application of n -> 2n-1 are the evil numbers times 4 plus 3. - Ralf Stephan, May 25 2013

Francisco J. Muñoz, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Francisco J. Muñoz and Juan Carlos Nuño, Rule-based Generation of de Bruijn Sequences: Memory and Learning, arXiv:2507.09764 [cs.FL], 2025. See p. 9.

Intersection of A001969 and A005408.
Supersequence of A093688.
Cf. A092246 (odd odious numbers).
Column 2 of A277880, positions of 1's in A277808 (2's in A277822).
Cf. A000069, A128309, A277823, A277902.

Select[Range[300], OddQ[ # ] && EvenQ[DigitCount[ #, 2, 1]] &] (* Stefan Steinerberger, May 17 2007 *)
Select[Range[300], EvenQ[Plus @@ IntegerDigits[ #, 2]] && OddQ[ # ] &]

is(n)=n%2 && hammingweight(n)%2==0 \\ Charles R Greathouse IV, Mar 21 2013

a(n)=4*n-if(hammingweight(n-1)%2,3,1) \\ Charles R Greathouse IV, Mar 21 2013

def A129771(n): return (((m:=n-1)<<1)+(m.bit_count()&1^1)<<1)+1 # Chai Wah Wu, Mar 09 2023

a(n) = 2*A000069(n) + 1. a(n) is 1 plus twice odious numbers.
a(n) = A128309(n) + 1. a(n) is 1 plus odious even numbers.
A132680(a(n)) = A132680((a(n)-1)/2) + 2. - Reinhard Zumkeller, Aug 26 2007
a(n) = 4n + O(1). - Charles R Greathouse IV, Mar 21 2013
a(n) = A001969(1+A000069(n)) = A277902(A277823(n)). - Antti Karttunen, Nov 05 2016

More terms from Stefan Steinerberger, May 17 2007

A227872 Number of odious divisors (A000069) of n.

Original entry on oeis.org

1, 2, 1, 3, 1, 2, 2, 4, 1, 2, 2, 3, 2, 4, 1, 5, 1, 2, 2, 3, 3, 4, 1, 4, 2, 4, 1, 6, 1, 2, 2, 6, 2, 2, 3, 3, 2, 4, 2, 4, 2, 6, 1, 6, 1, 2, 2, 5, 3, 4, 1, 6, 1, 2, 3, 8, 2, 2, 2, 3, 2, 4, 3, 7, 2, 4, 2, 3, 2, 6, 1, 4, 2, 4, 2, 6, 3, 4, 2, 5, 2, 4, 1, 9, 1, 2, 2

Offset: 1

Vladimir Shevelev, Oct 25 2013

Peter J. C. Moses, Table of n, a(n) for n = 1..10000

Cf. A000005, A000069, A093688, A093696, A129771, A330289, A355968, A355969, A227873 (sum odious divs.).

Cf. A010060.

A227872 := proc(n)
    option remember ;
    local a,d ;
    a := 0 ;
    for d in numtheory[divisors](n) do
        if not isA001969(d) then
            a := a+1 ;
        end if;
    end do:
    a ;
end proc:
seq(A227872(n),n=1..200) ; # R. J. Mathar, Aug 07 2022

a[n_] := DivisorSum[n, 1 &, OddQ[DigitCount[#, 2, 1]] &]; Array[a, 100] (* Amiram Eldar, Jul 23 2022 *)

a(n) = sumdiv(n, d, hammingweight(d) % 2); \\ Michel Marcus, Feb 06 2016

isod(n) = hammingweight(n) % 2; \\ A000069
a(n) = my(v=valuation(n, 2)); n >>= v; sumdiv(n,d,isod(d)) * (v+1) \\ David A. Corneth, Jul 23 2022

from sympy import divisors
def c(n): return bin(n).count("1")&1
def a(n): return sum(1 for d in divisors(n, generator=True) if c(d))
print([a(n) for n in range(1, 101)]) # Michael S. Branicky, Jul 23 2022

a(n) + A356018(n) = A000005(n).
a(2^n) = n+1. - Bernard Schott, Jul 22 2022
a(n) = 1 iff n is in A093688. - Bernard Schott, Jul 23 2022
a(n) = Sum_{d|n} A010060(d). - Ridouane Oudra, Apr 12 2025

More terms from Peter J. C. Moses, Oct 25 2013

A356018 a(n) is the number of evil divisors (A001969) of n.

Original entry on oeis.org

0, 0, 1, 0, 1, 2, 0, 0, 2, 2, 0, 3, 0, 0, 3, 0, 1, 4, 0, 3, 1, 0, 1, 4, 1, 0, 3, 0, 1, 6, 0, 0, 2, 2, 1, 6, 0, 0, 2, 4, 0, 2, 1, 0, 5, 2, 0, 5, 0, 2, 3, 0, 1, 6, 1, 0, 2, 2, 0, 9, 0, 0, 3, 0, 2, 4, 0, 3, 2, 2, 1, 8, 0, 0, 4, 0, 1, 4, 0, 5, 3, 0, 1, 3, 3, 2

Offset: 1

Bernard Schott, Jul 23 2022

a(n) = 0 iff n is in A093696.

			12 has 6 divisors: {1, 2, 3, 4, 6, 12} of which three {3, 6, 12} have an even number of 1's in their binary expansion with 11, 110 and 11100 respectively; hence a(12) = 3.

Cf. A000005, A001969, A093688, A093696 (location of 0s), A227872, A356019, A356020.

Similar sequences: A083230, A087990, A087991, A332268, A355302.

A356018 := proc(n)
    local a,d ;
    a := 0 ;
    for d in numtheory[divisors](n) do
        if isA001969(d) then
            a := a+1 ;
        end if;
    end do:
    a ;
end proc:
seq(A356018(n),n=1..200) ;  # R. J. Mathar, Aug 07 2022

a[n_] := DivisorSum[n, 1 &, EvenQ[DigitCount[#, 2, 1]] &]; Array[a, 100] (* Amiram Eldar, Jul 23 2022 *)

a(n) = my(v = valuation(n, 2)); n>>=v; d=divisors(n); sum(i=1, #d, bitand(hammingweight(d[i]), 1) == 0) * (v+1) \\ David A. Corneth, Jul 23 2022

from sympy import divisors
def c(n): return bin(n).count("1")&1 == 0
def a(n): return sum(1 for d in divisors(n, generator=True) if c(d))
print([a(n) for n in range(1, 101)]) # Michael S. Branicky, Jul 23 2022

a(n) = A000005(n) - A227872(n).

More terms from David A. Corneth, Jul 23 2022

A356020 Positions of records in A356018, i.e., integers whose number of evil divisors sets a new record.

Original entry on oeis.org

1, 3, 6, 12, 18, 30, 60, 90, 120, 180, 360, 540, 720, 1080, 1440, 2160, 3780, 4320, 6120, 7560, 8640, 12240, 15120, 24480, 27720, 30240, 36720, 48960, 50400, 55440, 73440, 83160, 110880, 128520, 138600, 166320, 221760, 257040, 277200, 332640, 471240, 514080, 554400

Offset: 1

Bernard Schott, Jul 24 2022

Corresponding records of number of evil divisors are 0, 1, 2, 3, 4, 6, 9, 10, 12, 15, ...

			60 is in the sequence because A356018(60) = 9 is larger than any earlier value in A356018.

David A. Corneth, Table of n, a(n) for n = 1..189

Cf. A001969, A093688, A093696, A356018, A356019.

Similar sequences: A093036, A093037, A330815, A350756, A355969.

f[n_] := DivisorSum[n, 1 &, EvenQ[DigitCount[#, 2, 1]] &]; fm = -1; s = {}; Do[If[(fn = f[n]) > fm, fm = fn; AppendTo[s, n]], {n, 1, 10^5}]; s (* Amiram Eldar, Jul 24 2022 *)

upto(n) = my(res = List(), r=-1); forfactored(i=1, n, if(numdiv(i[2]) > r, d = divisors(i[2]); c=sum(j=1, #d, isevil(d[j])); if(c>r, r=c; listput(res,i[1])))); res
isevil(n) = bitand(hammingweight(n), 1)==0 \\ David A. Corneth, Jul 24 2022

from sympy import divisors
from itertools import count, islice
def c(n): return bin(n).count("1")&1 == 0
def f(n): return sum(1 for d in divisors(n, generator=True) if c(d))
def agen(record=-1):
    for k in count(1):
        if f(k) > record: record = f(k); yield k
print(list(islice(agen(), 40))) # Michael S. Branicky, Jul 24 2022

More terms from Amiram Eldar, Jul 24 2022

A356019 a(n) is the smallest number that has exactly n evil divisors (A001969).

Original entry on oeis.org

1, 3, 6, 12, 18, 45, 30, 135, 72, 60, 90, 765, 120, 1575, 270, 180, 600, 3465, 480, 13545, 360, 540, 1530, 10395, 1260, 720, 3150, 1980, 1080, 49725, 1440, 45045, 2520, 3060, 6930, 2160, 3780, 58905, 27090, 6300, 5040, 184275, 4320, 135135, 6120, 7920, 20790, 329175, 7560, 8640

Offset: 0

Bernard Schott, Jul 23 2022

Differs from A327328 since a(7) = 135 while A327328(7) = 105.

			a(4) = 18 since 18 has six divisors: {1, 2, 3, 6, 9, 18} of which four {3, 6, 9, 18} have an even number of 1's in their binary expansion: 11, 110, 1001 and 10010 respectively; also, no positive integer smaller than 18 has exactly four divisors that are evil.

David A. Corneth, Table of n, a(n) for n = 0..396
David A. Corneth, Upperbounds on a(n), terms <= 8*10^9 are certain

Cf. A001969, A093688, A227872, A327328, A356018, A356020, A356040.

Similar sequences: A087997, A333456, A355303, A355699, A355968.

# output in unsorted b-file style
A356019_list := [seq(0,i=1..1000)] ;
for n from 1 do
    evd := A356018(n) ;
    if evd < nops(A356019_list) then
        if op(evd+1,A356019_list) <= 0 then
            printf("%d %d\n",evd,n) ;
            A356019_list := subsop(evd+1=n,A356019_list) ;
        end if;
    end if;
end do:  # R. J. Mathar, Aug 07 2022

f[n_] := DivisorSum[n, 1 &, EvenQ[DigitCount[#, 2, 1]] &]; seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s]; seq[50, 10^6] (* Amiram Eldar, Jul 23 2022 *)

from sympy import divisors
from itertools import count, islice
def c(n): return bin(n).count("1")&1 == 0
def f(n): return sum(1 for d in divisors(n, generator=True) if c(d))
def agen():
    n, adict = 0, dict()
    for k in count(1):
        fk = f(k)
        if fk not in adict: adict[fk] = k
        while n in adict: yield adict[n]; n += 1
print(list(islice(agen(), 50))) # Michael S. Branicky, Jul 23 2022

a(n) <= A356040(n). - David A. Corneth, Jul 26 2022

More terms from Amiram Eldar, Jul 23 2022

A237417 Numbers that are the product of an odiousfree number and an evilfree number.

Original entry on oeis.org

3, 5, 6, 9, 10, 12, 15, 17, 18, 20, 21, 23, 24, 27, 29, 30, 33, 34, 35, 36, 39, 40, 42, 43, 45, 46, 48, 51, 53, 54, 55, 57, 58, 60, 63, 65, 66, 68, 70, 71, 72, 78, 80, 83, 84, 85, 86, 89, 90, 92, 93, 95, 96, 99, 101, 102, 105, 106, 108, 110, 111, 113, 114, 116, 117, 119, 120, 123, 126, 129

Offset: 1

Irina Gerasimova, Feb 23 2014, following a suggestion from Juri-Stepan Gerasimov

Odiousfree*evilfree numbers: numbers of the form odiousfree*evilfree.
Subsequence of this sequence (A237417): numbers that are not the products of two odious numbers or the products of two evil numbers: 3, 5, 6, 10, 12, 17, 20, 23, 24, 29, 33, 34, 39, 40, 43, 46, 48, 57, 58, 63, 65, 66, 68, 71, 78, 80, 83, 86, 89, 92, 95, 101, 105, 106, 111, 113, 114, 116, 119,...
Putting the 1 aside in A093688, these could be called odiousfree numbers, and are a subsequence of A001969. A093696 would be the evilfree numbers then, and are a subsequence of A000069.

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A093688, A093696.

N:= 200: # to get all terms <= N
Ofree:= {$2..N}: Efree:= {$1..N/3}:
for n from 2 to N do
  t:= convert(convert(n,base,2),`+`) mod 2;
  if t = 0 then Efree:= Efree minus {seq(i,i=n..N/3,n)}
  else Ofree:= Ofree minus {seq(i,i=n..N,n)}
  fi
od:
sort(convert(select(`<=`,{seq(seq(s*t,s=Ofree),t=Efree)},N),list)); # Robert Israel, May 09 2019

odFreeQ[n_] := AllTrue[Rest @ Divisors[n], EvenQ[DigitCount[#, 2, 1]] &]; evFreeQ[n_] := AllTrue[Divisors[n], OddQ[DigitCount[#, 2, 1]] &]; m = 100; o = Select[Range[2, m], odFreeQ]; e = Select[Range[m], evFreeQ]; Union @ Select[Times @@@ Tuples[{o, e}], # <= m &] (* Amiram Eldar, Oct 16 2020 *)

isA093696(n)= fordiv(n, d, if(hammingweight(d)%2==0, return(0))); 1;
isA093688(n)= if (n==1, 0, sumdiv(n, d, hammingweight(d)%2)==1);
lista(nn) = {vn = vector(2*nn, i, i); vof = select(n->isA093696(n), vn); vef = select(n->isA093688(n), vn); vp = []; for (i = 1, #vof, for (j = 1, #vef, vp = Set(concat(vp, vof[i]*vef[j])););); for (i = 1, #vp, if (vp[i] <= nn, print1(vp[i], ", ")););} \\ Michel Marcus, Mar 05 2014

a(n) = A093688(k+1)*A093696(m).

Definition corrected by Jon E. Schoenfield, Feb 26 2014

A238989 Numbers that are not odiousfree*evilfree. The complement of A237417.

Original entry on oeis.org

1, 2, 4, 7, 8, 11, 13, 14, 16, 19, 22, 25, 26, 28, 31, 32, 37, 38, 41, 44, 47, 49, 50, 52, 56, 59, 61, 62, 64, 67, 69, 73, 74, 75, 76, 77, 79, 81, 82, 87, 88, 91, 94, 97, 98, 100, 103, 104, 107, 109, 112, 115, 118, 121, 122, 124, 125, 127, 128

Offset: 1

Juri-Stepan Gerasimov, Mar 08 2014

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A000069 (odious numbers), A001969 (evil numbers), A093688 (1 together with odiousfree numbers), A093696 (evilfree numbers), A237417 (odiousfree*evilfree numbers).

odFreeQ[n_] := AllTrue[Rest @ Divisors[n], EvenQ[DigitCount[#, 2, 1]] &]; evFreeQ[n_] := AllTrue[Divisors[n], OddQ[DigitCount[#, 2, 1]] &]; m = 100; o = Select[Range[2, m], odFreeQ]; e = Select[Range[m], evFreeQ]; Complement[Range[m], Union @ Select[Times @@@ Tuples[{o, e}], # <= m &]] (* Amiram Eldar, Oct 16 2020 *)

A331833 Numbers k such that all the divisors of k, excluding 1, have an even number of 1's in their negabinary representations.

Original entry on oeis.org

1, 2, 5, 7, 10, 13, 14, 17, 19, 25, 31, 34, 37, 49, 61, 62, 65, 67, 73, 79, 85, 97, 107, 127, 133, 155, 167, 170, 173, 179, 193, 214, 217, 223, 229, 241, 247, 254, 257, 259, 271, 277, 289, 310, 313, 325, 334, 337, 347, 359, 365, 395, 419, 425, 427, 431, 434, 443

Offset: 1

Amiram Eldar, Jan 28 2020

			10 is a term since all of its divisors exclusing 1, i.e., 2, 5 and 10, or 110, 101, and 11110 in negabinary representation, have an even number of 1's.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Subsequence of A268272.

Cf. A027615, A039724, A093688.

negaBinWt[n_] := negaBinWt[n] = If[n==0, 0, negaBinWt[Quotient[n-1, -2]] + Mod[n, 2]]; eveNegaBinQ[n_] := EvenQ[negaBinWt[n]]; seqQ[n_] := AllTrue[Rest @ Divisors[n], eveNegaBinQ]; Select[Range[401],seqQ]

A338611 Numbers all of whose divisors, excluding the divisor 1, are evil numbers (A001969) with a record number of divisors.

Original entry on oeis.org

1, 3, 9, 15, 45, 135, 765, 2295, 196605, 589815, 12884901885, 38654705655

Offset: 1

Amiram Eldar, Nov 03 2020

A number m is in this sequence if it is in A093688, and d(m) > d(k) for all terms k < m in A093688, where d(m) is the number of divisors of m (A000005).
The corresponding record numbers of divisors are 1, 2, 3, 4, 6, 8, 12, 16, 24, 32, 48, 64, ...
Apparently, all the terms except for 1 are products of powers of Fermat primes (A019434). 3 seems to be the only prime with multiplicity larger than 1 in some of the terms. There are no other terms in this sequence that are products of powers of the 5 known Fermat primes.

			The first 4 terms of A093688 are 1, 3, 5, 9, and their numbers of divisors are 1, 2, 2, 3. The record values 1, 2 and 3 occur at 1, 3 and 9 that are the first 3 terms of this sequence.

Subsequence of A001969 and A093688.
Similar sequence with odious numbers: A330289.
Cf. A000005, A019434.

evilQ[n_] := EvenQ @ DigitCount[n, 2, 1]; allDivEvilQ[n_] := AllTrue[Rest @ Divisors[n], evilQ]; divNumMax = 0; seq={}; Do[If[allDivEvilQ[n] && (divNum = DivisorSigma[0, n]) > divNumMax, divNumMax = divNum; AppendTo[seq, n]], {n, 1, 6*10^5}]; seq

cp's OEIS Frontend

A129771 Evil odd numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Python

Formula

Extensions

A227872 Number of odious divisors (A000069) of n.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A356018 a(n) is the number of evil divisors (A001969) of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A356020 Positions of records in A356018, i.e., integers whose number of evil divisors sets a new record.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Extensions

A356019 a(n) is the smallest number that has exactly n evil divisors (A001969).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Python

Formula

Extensions

A237417 Numbers that are the product of an odiousfree number and an evilfree number.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica