A003159 -id:A003159 - cp's OEIS frontend

A372588 Numbers k > 1 such that (greatest binary index of k) + (greatest prime index of k) is odd.

2, 6, 7, 8, 10, 11, 15, 18, 19, 21, 24, 26, 27, 28, 29, 32, 33, 34, 40, 41, 44, 45, 46, 47, 50, 51, 55, 59, 60, 62, 65, 70, 71, 72, 74, 76, 78, 79, 81, 84, 86, 87, 89, 91, 95, 96, 98, 101, 104, 105, 106, 107, 108, 111, 112, 113, 114, 116, 117, 122, 126, 128

Offset: 1

Gus Wiseman, May 14 2024

nonn

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The even version is A372589.

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
        {2}   2  (1)
      {2,3}   6  (2,1)
    {1,2,3}   7  (4)
        {4}   8  (1,1,1)
      {2,4}  10  (3,1)
    {1,2,4}  11  (5)
  {1,2,3,4}  15  (3,2)
      {2,5}  18  (2,2,1)
    {1,2,5}  19  (8)
    {1,3,5}  21  (4,2)
      {4,5}  24  (2,1,1,1)
    {2,4,5}  26  (6,1)
  {1,2,4,5}  27  (2,2,2)
    {3,4,5}  28  (4,1,1)
  {1,3,4,5}  29  (10)
        {6}  32  (1,1,1,1,1)
      {1,6}  33  (5,2)
      {2,6}  34  (7,1)
      {4,6}  40  (3,1,1,1)
    {1,4,6}  41  (13)
    {3,4,6}  44  (5,1,1)
  {1,3,4,6}  45  (3,2,2)

For sum (A372428, zeros A372427) we have A372586.
For minimum (A372437) we have A372439, complement A372440.
For length (A372441, zeros A071814) we have A372590, complement A372591.
Positions of odd terms in A372442, zeros A372436.
The complement is A372589.
For just binary indices:
- length: A000069, complement A001969
- sum: A158705, complement A158704
- minimum: A003159, complement A036554
- maximum: A053738, complement A053754
For just prime indices:
- length: A026424 (count A027193), complement A028260 (count A027187)
- sum: A300063 (count A058695), complement A300061 (count A058696)
- minimum: A340932 (count A026804), complement A340933 (count A026805)
- maximum: A244991 (count A027193), complement A244990 (count A027187)
A005408 lists odd numbers.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031368 lists odd-indexed primes, even A031215.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Cf. A000720, A006141, A066208, A160786, A243055, A257991, A300272, A304818, A340604, A341446, A372429-A372433, A372438.

Select[Range[2,100],OddQ[IntegerLength[#,2]+PrimePi[FactorInteger[#][[-1,1]]]]&]

Numbers k such that A070939(k) + A061395(k) is odd.

A329575 Numbers whose smallest Fermi-Dirac factor is 3.

Original entry on oeis.org

3, 12, 15, 21, 27, 33, 39, 48, 51, 57, 60, 69, 75, 84, 87, 93, 105, 108, 111, 123, 129, 132, 135, 141, 147, 156, 159, 165, 177, 183, 189, 192, 195, 201, 204, 213, 219, 228, 231, 237, 240, 243, 249, 255, 267, 273, 276, 285, 291, 297, 300, 303, 309, 321, 327, 336, 339, 345

Offset: 1

Peter Munn, Apr 27 2020

nonn

Every positive integer is the product of a unique subset of the terms of A050376 (sometimes called Fermi-Dirac primes). This sequence lists the numbers where the relevant subset includes 3 but not 2.
Numbers whose squarefree part is divisible by 3 but not 2.
Positive multiples of 3 that survive sieving by the rule: if m appears then 2m, 3m and 6m do not. Asymptotic density is 1/6.

			6 is the product of the following terms of A050376: 2, 3. These terms include 2, so 6 is not in the sequence.
12 is the product of the following terms of A050376: 3, 4. These terms include 3, but not 2, so 12 is in the sequence.
20 is the product of the following terms of A050376: 4, 5. These terms do not include 3, so 20 is not in the sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Jan Snellman, Greedy Regular Convolutions, arXiv:2504.02795 [math.NT], 2025.

Cf. A007913, A036668, A050376, A059897, A223490, A307150, A339690.
Intersection of any 2 of A003159, A145204 and A325424; also subsequence of A028983.
Ordered 3rd quadrisection of A052330.

f[p_, e_] := p^(2^IntegerExponent[e, 2]); fdmin[n_] := Min @@ f @@@ FactorInteger[n]; Select[Range[350], fdmin[#] == 3 &] (* Amiram Eldar, Nov 27 2020 *)

isok(m) = core(m) % 6 == 3; \\ Michel Marcus, May 01 2020

from itertools import count
from sympy import integer_log
def A329575(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        c = n+x
        for i in range(integer_log(x,9)[0]+1):
            i2 = 9**i
            for j in count(0,2):
                k = i2<x:
                    break
                m = x//k
                c -= (m-1)//6+(m-5)//6+2
        return c
    return 3*bisection(f,n,n) # Chai Wah Wu, Apr 10 2025

A223490(a(n)) = 3.
A007913(a(n)) == 3 (mod 6).
A059897(2, a(n)) = 2 * a(n).
A059897(3, a(n)) * 3 = a(n).
{a(n) : n >= 1} = {k : 3 * A307150(k) = 2 * k}.
A003159 = {a(n) / 3 : n >= 1} U {a(n) : n >= 1}.
A036668 = {a(n) / 3 : n >= 1} U {a(n) * 2 : n >= 1}.
A145204 \ {0} = {a(n) : n >= 1} U {a(n) * 2 : n >= 1}.
a(n) = 3*A339690(n). - Chai Wah Wu, Apr 10 2025

A332821 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 1 (mod 3).

Original entry on oeis.org

2, 5, 9, 11, 12, 16, 17, 21, 23, 28, 30, 31, 39, 40, 41, 47, 49, 52, 54, 57, 59, 66, 67, 70, 72, 73, 75, 76, 83, 87, 88, 91, 96, 97, 100, 102, 103, 109, 111, 116, 126, 127, 128, 129, 130, 133, 135, 136, 137, 138, 148, 149, 154, 157, 159, 165, 167, 168, 169, 172, 175, 179, 180, 183, 184, 186, 190, 191, 197, 203, 211, 212

Offset: 1

Antti Karttunen and Peter Munn, Feb 25 2020

nonn

The positive integers are partitioned between A332820, this sequence and A332822.
For each prime p, the terms include exactly one of p and p^2. The primes alternate between this sequence and A332822. This sequence has the primes with odd indexes, those in A031368.
The terms are the even numbers in A332822 halved. The terms are also the numbers m such that 5m is in A332822, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332820, and so on for alternate primes: 7, 13, 19 etc.
The numbers that are half of the even terms of this sequence are in A332820, which consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332822, which consists exactly of those numbers. For larger primes, an alternating pattern applies as described in the previous paragraph.
If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, the resulting number is in A332822, which consists entirely of those numbers.
The product of any 2 terms of this sequence is in A332822, the product of any 3 terms is in A332820, and the product of a term of A332820 and a term of this sequence is in this sequence. So if a number k is present, k^2 is in A332822, k^3 is in A332820, and k^4 is in this sequence.
If k is an even number, exactly one of {k/2, k, 2k} is in the sequence (cf. A191257 / A067368 / A213258); and generally if k is a multiple of a prime p, exactly one of {k/p, k, k*p} is in the sequence.

Positions of ones in A332823; equivalently, numbers in row 3k+1 of A277905 for some k >= 0.
Cf. A048675, A332820, A332822.
Comparable 2 or 3-way classifications: A000379/A000028, A001969/A000069, A003159/A036554, A005843/A005408, A028260/A026424, A191257/A067368/A213258, A325431/A325432, A329609/A329604/A332812.
Subsequences: intersection of A026478 and A066208, A031368 (prime terms), A033431\{0}, A052934\{1}, A069486, A099800, A167747\{1}, A244725\{0}, A244728\{0}, A338911 (semiprime terms).

Select[Range@ 212, Mod[Total@ #, 3] == 1 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]] &] (* Michael De Vlieger, Mar 15 2020 *)

isA332821(n) =  { my(f = factor(n)); (1==((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3)); };

{a(n) : n >= 1} = {2 * A332820(k) : k >= 1} U {A003961(A332822(k)) : k >= 1}.

{a(n) : n >= 1} = {A332822(k)^2 : k >= 1} U {A331590(2, A332820(k)) : k >= 1}.

A332822 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 2 (mod 3).

Original entry on oeis.org

3, 4, 7, 10, 13, 18, 19, 22, 24, 25, 29, 32, 34, 37, 42, 43, 45, 46, 53, 55, 56, 60, 61, 62, 71, 78, 79, 80, 81, 82, 85, 89, 94, 98, 99, 101, 104, 105, 107, 108, 113, 114, 115, 118, 121, 131, 132, 134, 139, 140, 144, 146, 150, 151, 152, 153, 155, 163, 166, 173, 174, 176, 181, 182, 187, 189, 192, 193, 194, 195, 199, 200, 204

Offset: 1

Antti Karttunen and Peter Munn, Feb 25 2020

nonn

The positive integers are partitioned between A332820, A332821 and this sequence.
For each prime p, the terms include exactly one of p and p^2. The primes alternate between this sequence and A332821. This sequence has the primes with even indexes, those in A031215.
The terms are the even numbers in A332820 halved. The terms are also the numbers m such that 5m is in A332820, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332821, and so on for alternate primes: 7, 13, 19, 29 etc.
If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, we get the same set of numbers as we get from halving the even terms of this sequence, and A332821 consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332820, which consists exactly of those numbers. The numbers that are one fifth of the terms that are multiples of 5 constitute A332821, and for larger primes, an alternating pattern applies as described in the previous paragraph.
The product of any 2 terms of this sequence is in A332821, the product of any 3 terms is in A332820, and the product of a term of A332820 and a term of this sequence is in this sequence. So if a number k is present, k^2 is in A332821, k^3 is in A332820, and k^4 is in this sequence.
If k is an even number, exactly one of {k/2, k, 2k} is in the sequence (cf. A191257 / A067368 / A213258); and generally if k is a multiple of a prime p, exactly one of {k/p, k, k*p} is in the sequence.

Positions of terms valued -1 in A332823; equivalently, numbers in row 3k-1 of A277905 for some k >= 1.
Cf. A048675, A332820, A332821.
Comparable 2 or 3-way classifications: A000379/A000028, A001969/A000069, A003159/A036554, A005843/A005408, A028260/A026424, A191257/A067368/A213258, A325431/A325432, A329609/A329604/A332812.
Subsequences: intersection of A026478 and A066207, A031215 (prime terms), A033430\{0}, A117642\{0}, A169604, A244727\{0}, A244729\{0}, A338910 (semiprime terms).

Select[Range@ 204, Mod[Total@ #, 3] == 2 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]] &] (* Michael De Vlieger, Mar 15 2020 *)

isA332822(n) =  { my(f = factor(n)); (2==((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3)); };

{a(n) : n >= 1} = {2 * A332821(k) : k >= 1} U {A003961(A332821(k)) : k >= 1}.

{a(n) : n >= 1} = {A332821(k)^2 : k >= 1} U {A331590(2, A332821(k)) : k >= 1}.

A372586 Numbers k such that (sum of binary indices of k) + (sum of prime indices of k) is odd.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 9, 12, 15, 16, 17, 20, 21, 29, 32, 36, 42, 43, 45, 46, 47, 48, 51, 53, 54, 55, 59, 60, 61, 63, 64, 65, 66, 67, 68, 71, 73, 78, 79, 80, 81, 84, 89, 91, 93, 94, 95, 97, 99, 101, 105, 110, 111, 113, 114, 115, 116, 118, 119, 121, 122, 125, 127

Offset: 1

Gus Wiseman, May 14 2024

nonn

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The even version is A372587.

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
            {1}   1  ()
            {2}   2  (1)
          {1,2}   3  (2)
            {3}   4  (1,1)
          {1,3}   5  (3)
            {4}   8  (1,1,1)
          {1,4}   9  (2,2)
          {3,4}  12  (2,1,1)
      {1,2,3,4}  15  (3,2)
            {5}  16  (1,1,1,1)
          {1,5}  17  (7)
          {3,5}  20  (3,1,1)
        {1,3,5}  21  (4,2)
      {1,3,4,5}  29  (10)
            {6}  32  (1,1,1,1,1)
          {3,6}  36  (2,2,1,1)
        {2,4,6}  42  (4,2,1)
      {1,2,4,6}  43  (14)
      {1,3,4,6}  45  (3,2,2)
      {2,3,4,6}  46  (9,1)
    {1,2,3,4,6}  47  (15)
          {5,6}  48  (2,1,1,1,1)

Positions of odd terms in A372428, zeros A372427.
For minimum (A372437) we have A372439, complement A372440.
For length (A372441, zeros A071814) we have A372590, complement A372591.
For maximum (A372442, zeros A372436) we have A372588, complement A372589.
The complement is A372587.
For just binary indices:
- length: A000069, complement A001969
- sum: A158705, complement A158704
- minimum: A003159, complement A036554
- maximum: A053738, complement A053754
For just prime indices:
- length: A026424 (count A027193), complement A028260 (count A027187)
- sum: A300063 (count A058695), complement A300061 (count A058696)
- minimum: A340932 (count A026804), complement A340933 (count A026805)
- maximum: A244991 (count A027193), complement A244990 (count A027187)
A005408 lists odd numbers.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031368 lists odd-indexed primes, even A031215.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Cf. A000720, A066208, A160786, A257991, A300272, A304818, A340604, A341446, A372429-A372433, A372438.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
Select[Range[100],OddQ[Total[bix[#]]+Total[prix[#]]]&]

Numbers k such that A029931(k) + A056239(k) is odd.

A372589 Numbers k > 1 such that (greatest binary index of k) + (greatest prime index of k) is even.

Original entry on oeis.org

3, 4, 5, 9, 12, 13, 14, 16, 17, 20, 22, 23, 25, 30, 31, 35, 36, 37, 38, 39, 42, 43, 48, 49, 52, 53, 54, 56, 57, 58, 61, 63, 64, 66, 67, 68, 69, 73, 75, 77, 80, 82, 83, 85, 88, 90, 92, 93, 94, 97, 99, 100, 102, 103, 109, 110, 115, 118, 119, 120, 121, 123, 124

Offset: 1

Gus Wiseman, May 14 2024

nonn

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The odd version is A372588.

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
        {1,2}   3  (2)
          {3}   4  (1,1)
        {1,3}   5  (3)
        {1,4}   9  (2,2)
        {3,4}  12  (2,1,1)
      {1,3,4}  13  (6)
      {2,3,4}  14  (4,1)
          {5}  16  (1,1,1,1)
        {1,5}  17  (7)
        {3,5}  20  (3,1,1)
      {2,3,5}  22  (5,1)
    {1,2,3,5}  23  (9)
      {1,4,5}  25  (3,3)
    {2,3,4,5}  30  (3,2,1)
  {1,2,3,4,5}  31  (11)
      {1,2,6}  35  (4,3)
        {3,6}  36  (2,2,1,1)
      {1,3,6}  37  (12)
      {2,3,6}  38  (8,1)
    {1,2,3,6}  39  (6,2)
      {2,4,6}  42  (4,2,1)
    {1,2,4,6}  43  (14)

For sum (A372428, zeros A372427) we have A372587, complement A372586.
For minimum (A372437) we have A372440, complement A372439.
For length (A372441, zeros A071814) we have A372591, complement A372590.
Positions of even terms in A372442, zeros A372436.
The complement is A372588.
For just binary indices:
- length: A001969, complement A000069
- sum: A158704, complement A158705
- minimum:  A036554, complement A003159
- maximum: A053754, complement A053738
For just prime indices:
- length: A026424 A028260 (count A027187), complement (count A027193)
- sum: A300061 (count A058696), complement A300063 (count A058695)
- minimum: A340933 (count A026805), complement A340932 (count A026804)
- maximum: A244990 (count A027187), complement A244991 (count A027193)
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031215 lists even-indexed primes, odd A031368.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Cf. A000720, A006141, A066207, A243055, A257991, A300272, A304818, A340604, A341446, A372429-A372433, A372438.

Select[Range[2,100],EvenQ[IntegerLength[#,2]+PrimePi[FactorInteger[#][[-1,1]]]]&]

Numbers k such that A070939(k) + A061395(k) is even.

A372590 Numbers whose binary weight (A000120) plus bigomega (A001222) is odd.

Original entry on oeis.org

1, 3, 4, 5, 12, 14, 16, 17, 18, 20, 21, 22, 23, 25, 26, 27, 29, 30, 35, 38, 43, 45, 48, 49, 53, 55, 56, 62, 63, 64, 66, 68, 69, 71, 72, 74, 75, 78, 80, 81, 82, 83, 84, 87, 88, 89, 91, 92, 93, 94, 99, 100, 101, 102, 104, 105, 108, 113, 114, 115, 116, 118, 120

Offset: 1

Gus Wiseman, May 14 2024

nonn

The even version is A372591.

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
        {1}   1  ()
      {1,2}   3  (2)
        {3}   4  (1,1)
      {1,3}   5  (3)
      {3,4}  12  (2,1,1)
    {2,3,4}  14  (4,1)
        {5}  16  (1,1,1,1)
      {1,5}  17  (7)
      {2,5}  18  (2,2,1)
      {3,5}  20  (3,1,1)
    {1,3,5}  21  (4,2)
    {2,3,5}  22  (5,1)
  {1,2,3,5}  23  (9)
    {1,4,5}  25  (3,3)
    {2,4,5}  26  (6,1)
  {1,2,4,5}  27  (2,2,2)
  {1,3,4,5}  29  (10)
  {2,3,4,5}  30  (3,2,1)
    {1,2,6}  35  (4,3)
    {2,3,6}  38  (8,1)
  {1,2,4,6}  43  (14)
  {1,3,4,6}  45  (3,2,2)

For sum (A372428, zeros A372427) we have A372586, complement A372587.
For minimum (A372437) we have A372439, complement A372440.
Positions of odd terms in A372441, zeros A071814.
For maximum (A372442, zeros A372436) we have A372588, complement A372589.
The complement is A372591.
For just binary indices:
- length: A000069, complement A001969
- sum: A158705, complement A158704
- minimum: A003159, complement A036554
- maximum: A053738, complement A053754
For just prime indices:
- length: A026424 (count A027193), complement A028260 (count A027187)
- sum: A300063 (count A058695), complement A300061 (count A058696)
- minimum: A340932 (count A026804), complement A340933 (count A026805)
- maximum: A244991 (count A027193), complement A244990 (count A027187)
A005408 lists odd numbers.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031368 lists odd-indexed primes, even A031215.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Cf. A000720, A066208, A160786, A257991, A300272, A304818, A340604, A341446, A372429-A372433, A372438.

Select[Range[100],OddQ[DigitCount[#,2,1]+PrimeOmega[#]]&]

A036552 List of pairs (m,2m) where m is the least unused positive number.

Original entry on oeis.org

1, 2, 3, 6, 4, 8, 5, 10, 7, 14, 9, 18, 11, 22, 12, 24, 13, 26, 15, 30, 16, 32, 17, 34, 19, 38, 20, 40, 21, 42, 23, 46, 25, 50, 27, 54, 28, 56, 29, 58, 31, 62, 33, 66, 35, 70, 36, 72, 37, 74, 39, 78, 41, 82, 43, 86, 44, 88, 45, 90, 47, 94, 48, 96, 49, 98, 51, 102

Offset: 1

Tom Verhoeff

A permutation of the natural numbers. Inverse permutation is A065037.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
P. Erdős, R. Freud, and N. Hegyvári, Arithmetical properties of permutations of integers, Acta Mathematica Hungarica 41:1-2 (1983), pp. 169-176.
Dana G. Korssjoen, Biyao Li, Stefan Steinerberger, Raghavendra Tripathi, and Ruimin Zhang, Finding structure in sequences of real numbers via graph theory: a problem list, arXiv:2012.04625[math.CO], 2020-2021.
Index entries for sequences that are permutations of the natural numbers

Alternating merge of A003159 and A036554. Cf. A064736, A065037.

import Data.List (delete)
a036552 n = a036552_list !! (n-1)
a036552_list = g [1..] where
   g (x:xs) = x : (2*x) : (g $ delete (2*x) xs)
-- Reinhard Zumkeller, Feb 07 2011

w = {}; Do[ w = If[ FreeQ[w, k], w = Join[w, {k, 2k}], w], {k, 100}]; w
(* Jean-François Alcover, Nov 04 2011, after Wouter Meeussen *)

apairs(N) = my(m=0, r=List(), i=0); while(#rRuud H.G. van Tol, May 09 2024

A225838 Numbers of form 2^i3^j(6k+5), i, j, k >= 0.

Original entry on oeis.org

5, 10, 11, 15, 17, 20, 22, 23, 29, 30, 33, 34, 35, 40, 41, 44, 45, 46, 47, 51, 53, 58, 59, 60, 65, 66, 68, 69, 70, 71, 77, 80, 82, 83, 87, 88, 89, 90, 92, 94, 95, 99, 101, 102, 105, 106, 107, 113, 116, 118, 119, 120, 123, 125, 130, 131, 132, 135, 136, 137, 138

Offset: 1

Ralf Stephan, May 16 2013

Are a(n) > A225837(n) for all n? - Zak Seidov, May 17 2013
Yes. Imagine every 3-smooth number, m, visits you regularly, depositing a gold coin for safe keeping at each epoch (6k+1)*m and collecting it at epoch (6k+5)*m. If you run out of coins, you are doing something other than keeping them in a vault! - Peter Munn, Nov 13 2023
The asymptotic density of this sequence is 1/2. - Amiram Eldar, Apr 03 2022

Zak Seidov, Table of n, a(n) for n = 1..10000
Zak Seidov, Graph of A225838(n) - A225837(n) for n=1..126454.

Complement of A225837.

Symmetric difference of A003159 and A026225; of A189716 and A325424.

[n: n in [1..200] | d mod 6 eq 5 where d is n div (2^Valuation(n,2)*3^Valuation(n,3))]; // Bruno Berselli, May 16 2013

mx = 153; t = {}; Do[n = 2^i*3^j (6 k + 5); If[n <= mx, AppendTo[t, n]], {i, 0, Log[2, mx]}, {j, 0, Log[3, mx]}, {k, 0, mx/6}]; Union[t] (* T. D. Noe, May 16 2013 *)

for(n=1,200,t=n/(2^valuation(n,2)*3^valuation(n,3));if((t%6==5),print1(n,",")))

from sympy import integer_log
def A225838(n):
    def f(x): return n+sum(((x//3**i>>j)+5)//6 for i in range(integer_log(x,3)[0]+1) for j in range((x//3**i).bit_length()))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Feb 02 2025

A280049 Canonical representation of n as a sum of distinct Jacobsthal numbers J(n) (A001045) (see Comments for details); also binary numbers that end in an even number of zeros.

Original entry on oeis.org

1, 11, 100, 101, 111, 1001, 1011, 1100, 1101, 1111, 10000, 10001, 10011, 10100, 10101, 10111, 11001, 11011, 11100, 11101, 11111, 100001, 100011, 100100, 100101, 100111, 101001, 101011, 101100, 101101, 101111, 110000, 110001, 110011, 110100, 110101, 110111

Offset: 1

N. J. A. Sloane, Dec 31 2016

Every positive integer has a unique expression as a sum of distinct Jacobsthal numbers in which the index of the smallest summand is odd, with J(1) = 1 and J(2) = 1 both allowed. [Carlitz-Scoville-Hoggatt, 1972]. - Based on a comment in A001045 from Ira M. Gessel, Dec 31 2016.

The highest-order bits are on the left. Interpreting these as binary numbers we get A003159.

			9 = 5+3+1 = J(4)+J(3)+J(1) = 1101.

Lars Blomberg, Table of n, a(n) for n = 1..10000
L. Carlitz, R. Scoville, and V. E. Hoggatt, Jr., Representations for a special sequence, Fibonacci Quarterly 10.5 (1972), 499-518, 550.

Cf. A001045, A003159, A007088.

FromDigits[IntegerDigits[#, 2]] & /@ Select[Range[100], EvenQ[IntegerExponent[#, 2]] &] (* Amiram Eldar, Jul 14 2023 *)

lista(kmax) = for(k = 1, kmax, if(!(valuation(k, 2)%2), print1(fromdigits(binary(k), 10), ", "))); \\ Amiram Eldar, Jul 14 2023

from itertools import count, islice
def A280049_gen(): # generator of terms
    return map(lambda n:int(bin(n)[2:]),filter(lambda n:(n&-n).bit_length()&1,count(1)))
A280049_list = list(islice(A280049_gen(),20)) # Chai Wah Wu, Mar 19 2024

def A280049(n):
    def f(x):
        c, s = n+x, bin(x)[2:]
        l = len(s)
        for i in range(l&1^1,l,2):
            c -= int(s[i])+int('0'+s[:i],2)
        return c
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return int(bin(m)[2:]) # Chai Wah Wu, Jan 30 2025

a(n) = A007088(A003159(n)). - Amiram Eldar, Jul 14 2023

Corrected a(5), a(16) and more terms from Lars Blomberg, Jan 02 2017

cp's OEIS Frontend

A372588 Numbers k > 1 such that (greatest binary index of k) + (greatest prime index of k) is odd.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A329575 Numbers whose smallest Fermi-Dirac factor is 3.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A332821 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 1 (mod 3).

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

PARI

Formula

A332822 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 2 (mod 3).

Views

Author

Keywords

Comments

Crossrefs

Programs

Mathematica

PARI

Formula

A372586 Numbers k such that (sum of binary indices of k) + (sum of prime indices of k) is odd.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A372589 Numbers k > 1 such that (greatest binary index of k) + (greatest prime index of k) is even.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A372590 Numbers whose binary weight (A000120) plus bigomega (A001222) is odd.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A036552 List of pairs (m,2m) where m is the least unused positive number.

Views

Author

Keywords

Comments

Links

Crossrefs

A225838 Numbers of form 2^i3^j(6k+5), i, j, k >= 0.