cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-20 of 36 results. Next

A337461 Number of pairwise coprime ordered triples of positive integers summing to n.

Original entry on oeis.org

0, 0, 0, 1, 3, 3, 9, 3, 15, 9, 21, 9, 39, 9, 45, 21, 45, 21, 87, 21, 93, 39, 87, 39, 153, 39, 135, 63, 153, 57, 255, 51, 207, 93, 225, 93, 321, 81, 291, 135, 321, 105, 471, 105, 393, 183, 381, 147, 597, 147, 531, 213, 507, 183, 759, 207, 621, 273, 621, 231
Offset: 0

Views

Author

Gus Wiseman, Sep 02 2020

Keywords

Examples

			The a(3) = 1 through a(9) = 9 triples:
  (1,1,1)  (1,1,2)  (1,1,3)  (1,1,4)  (1,1,5)  (1,1,6)  (1,1,7)
           (1,2,1)  (1,3,1)  (1,2,3)  (1,5,1)  (1,2,5)  (1,3,5)
           (2,1,1)  (3,1,1)  (1,3,2)  (5,1,1)  (1,3,4)  (1,5,3)
                             (1,4,1)           (1,4,3)  (1,7,1)
                             (2,1,3)           (1,5,2)  (3,1,5)
                             (2,3,1)           (1,6,1)  (3,5,1)
                             (3,1,2)           (2,1,5)  (5,1,3)
                             (3,2,1)           (2,5,1)  (5,3,1)
                             (4,1,1)           (3,1,4)  (7,1,1)
                                               (3,4,1)
                                               (4,1,3)
                                               (4,3,1)
                                               (5,1,2)
                                               (5,2,1)
                                               (6,1,1)
		

Crossrefs

A000212 counts the unimodal instead of coprime version.
A220377*6 is the strict case.
A307719 is the unordered version.
A337462 counts these compositions of any length.
A337563 counts the case of partitions with no 1's.
A337603 only requires the *distinct* parts to be pairwise coprime.
A337604 is the intersecting instead of coprime version.
A014612 ranks 3-part partitions.
A302696 ranks pairwise coprime partitions.
A327516 counts pairwise coprime partitions.
A333228 ranks compositions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n,{3}],CoprimeQ@@#&]],{n,0,30}]

A326701 BII-numbers of set partitions.

Original entry on oeis.org

0, 1, 2, 3, 4, 8, 9, 10, 11, 12, 16, 18, 32, 33, 64, 128, 129, 130, 131, 132, 136, 137, 138, 139, 140, 144, 146, 160, 161, 192, 256, 258, 264, 266, 288, 512, 513, 520, 521, 528, 1024, 1032, 2048, 2049, 2050, 2051, 2052, 4096, 4098, 8192, 8193, 16384, 32768, 32769
Offset: 1

Views

Author

Gus Wiseman, Jul 21 2019

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. We define the set-system with BII-number n to be obtained by taking the binary indices of each binary index of n. For example, 18 has reversed binary expansion (0,1,0,0,1), and since the binary indices of 2 and 5 are {2} and {1,3} respectively, and {{2},{1,3}} is a set partition, it follows that 18 belongs to the sequence.

Examples

			The sequence of all set partitions together with their BII numbers begins:
    0: {}
    1: {{1}}
    2: {{2}}
    3: {{1},{2}}
    4: {{1,2}}
    8: {{3}}
    9: {{1},{3}}
   10: {{2},{3}}
   11: {{1},{2},{3}}
   12: {{1,2},{3}}
   16: {{1,3}}
   18: {{2},{1,3}}
   32: {{2,3}}
   33: {{1},{2,3}}
   64: {{1,2,3}}
  128: {{4}}
  129: {{1},{4}}
  130: {{2},{4}}
  131: {{1},{2},{4}}
  132: {{1,2},{4}}
  136: {{3},{4}}
		

Crossrefs

MM-numbers of set partitions are A302521.
BII-numbers of chains of nonempty sets are A326703.
BII-numbers of antichains of nonempty sets are A326704.

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[0,1000],UnsameQ@@Join@@bpe/@bpe[#]&]
  • Python
    from itertools import chain, count, combinations, islice
    from sympy.utilities.iterables import multiset_partitions
    def a_gen():
        yield 0
        for n in count(1):
            t = []
            for i in chain.from_iterable(combinations(range(1,n+1),r) for r in range(n+1)):
                if n in i:
                    for j in multiset_partitions(i):
                        t.append(sum(2**(sum(2**(m-1) for m in k)-1) for k in j))
            yield from sorted(t)
    A326701_list = list(islice(a_gen(), 100)) # John Tyler Rascoe, May 24 2024

A337561 Number of pairwise coprime strict compositions of n, where a singleton is not considered coprime unless it is (1).

Original entry on oeis.org

1, 1, 0, 2, 2, 4, 8, 6, 16, 12, 22, 40, 40, 66, 48, 74, 74, 154, 210, 228, 242, 240, 286, 394, 806, 536, 840, 654, 1146, 1618, 2036, 2550, 2212, 2006, 2662, 4578, 4170, 7122, 4842, 6012, 6214, 11638, 13560, 16488, 14738, 15444, 16528, 25006, 41002, 32802
Offset: 0

Views

Author

Gus Wiseman, Sep 18 2020

Keywords

Examples

			The a(1) = 1 through a(9) = 12 compositions (empty column shown as dot):
   (1)  .  (1,2)  (1,3)  (1,4)  (1,5)    (1,6)  (1,7)    (1,8)
           (2,1)  (3,1)  (2,3)  (5,1)    (2,5)  (3,5)    (2,7)
                         (3,2)  (1,2,3)  (3,4)  (5,3)    (4,5)
                         (4,1)  (1,3,2)  (4,3)  (7,1)    (5,4)
                                (2,1,3)  (5,2)  (1,2,5)  (7,2)
                                (2,3,1)  (6,1)  (1,3,4)  (8,1)
                                (3,1,2)         (1,4,3)  (1,3,5)
                                (3,2,1)         (1,5,2)  (1,5,3)
                                                (2,1,5)  (3,1,5)
                                                (2,5,1)  (3,5,1)
                                                (3,1,4)  (5,1,3)
                                                (3,4,1)  (5,3,1)
                                                (4,1,3)
                                                (4,3,1)
                                                (5,1,2)
                                                (5,2,1)
		

Crossrefs

A072706 counts unimodal strict compositions.
A220377*6 counts these compositions of length 3.
A305713 is the unordered version.
A337462 is the not necessarily strict version.
A000740 counts relatively prime compositions, with strict case A332004.
A051424 counts pairwise coprime or singleton partitions.
A101268 considers all singletons to be coprime, with strict case A337562.
A178472 counts compositions with a common factor > 1.
A327516 counts pairwise coprime partitions, with strict case A305713.
A328673 counts pairwise non-coprime partitions.
A333228 ranks compositions whose distinct parts are pairwise coprime.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],#=={}||UnsameQ@@#&&CoprimeQ@@#&]],{n,0,10}]

Formula

a(n) = A337562(n) - 1 for n > 1.

A367908 Numbers n such that there is only one way to choose a different binary index of each binary index of n.

Original entry on oeis.org

1, 2, 3, 5, 6, 8, 9, 10, 11, 13, 14, 17, 19, 21, 22, 24, 26, 28, 34, 35, 37, 38, 40, 41, 44, 49, 50, 56, 67, 69, 70, 73, 74, 81, 88, 98, 104, 128, 129, 130, 131, 133, 134, 136, 137, 138, 139, 141, 142, 145, 147, 149, 150, 152, 154, 156, 162, 163, 165, 166, 168
Offset: 1

Views

Author

Gus Wiseman, Dec 11 2023

Keywords

Comments

Also BII-numbers of set-systems (sets of nonempty sets) satisfying a strict version of the axiom of choice in exactly one way.
A binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion. A set-system is a finite set of finite nonempty sets. We define the set-system with BII-number n to be obtained by taking the binary indices of each binary index of n. Every finite set of finite nonempty sets has a different BII-number. For example, 18 has reversed binary digits (0,1,0,0,1), and since the binary indices of 2 and 5 are {2} and {1,3} respectively, the BII-number of {{2},{1,3}} is 18.
The axiom of choice says that, given any set of nonempty sets Y, it is possible to choose a set containing an element from each. The strict version requires this set to have the same cardinality as Y, meaning no element is chosen more than once.

Examples

			The set-system {{1},{1,2},{1,3}} with BII-number 21 satisfies the axiom in exactly one way, namely (1,2,3), so 21 is in the sequence.
The terms together with the corresponding set-systems begin:
   1: {{1}}
   2: {{2}}
   3: {{1},{2}}
   5: {{1},{1,2}}
   6: {{2},{1,2}}
   8: {{3}}
   9: {{1},{3}}
  10: {{2},{3}}
  11: {{1},{2},{3}}
  13: {{1},{1,2},{3}}
  14: {{2},{1,2},{3}}
  17: {{1},{1,3}}
  19: {{1},{2},{1,3}}
  21: {{1},{1,2},{1,3}}
  22: {{2},{1,2},{1,3}}
		

Crossrefs

These set-systems are counted by A367904.
Positions of 1's in A367905, firsts A367910, sorted firsts A367911.
If there is at least one choice we get A367906, counted by A367902.
If there are no choices we get A367907, counted by A367903.
If there are multiple choices we get A367909, counted by A367772.
The version for MM-numbers of multiset partitions is A368101.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A058891 counts set-systems, covering A003465, connected A323818.
A059201 counts covering T_0 set-systems.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.
A326031 gives weight of the set-system with BII-number n.
A368098 counts unlabeled multiset partitions for axiom, complement A368097.
BII-numbers: A309314 (hyperforests), A326701 (set partitions), A326703 (chains), A326704 (antichains), A326749 (connected), A326750 (clutters), A326751 (blobs), A326752 (hypertrees), A326754 (covers), A326783 (uniform), A326784 (regular), A326788 (simple), A330217 (achiral).

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[100], Length[Select[Tuples[bpe/@bpe[#]], UnsameQ@@#&]]==1&]
  • Python
    from itertools import count, islice, product
    def bin_i(n): #binary indices
        return([(i+1) for i, x in enumerate(bin(n)[2:][::-1]) if x =='1'])
    def a_gen(): #generator of terms
        for n in count(1):
            p = list(product(*[bin_i(k) for k in bin_i(n)]))
            x,c = len(p),0
            for j in range(x):
                if len(set(p[j])) == len(p[j]): c += 1
                if j+1 == x and c == 1: yield(n)
    A367908_list = list(islice(a_gen(), 100)) # John Tyler Rascoe, Feb 10 2024

Formula

A326674 GCD of the set of positions of 1's in the reversed binary expansion of n.

Original entry on oeis.org

1, 2, 1, 3, 1, 1, 1, 4, 1, 2, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 6, 1, 2, 1, 3, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Jul 17 2019

Keywords

Comments

a(n) is even if and only if n is in A062880. - Robert Israel, Oct 13 2020

Examples

			The reversed binary expansion of 40 is (0,0,0,1,0,1), with positions of 1's being {4,6}, so a(40) = GCD(4,6) = 2.
		

Crossrefs

Positions of 1's are A291166, and non-1's are A291165.
GCDs of prime indices are A289508.
GCDs of strict partitions encoded by FDH numbers are A319826.
Numbers whose binary positions are pairwise coprime are A326675.

Programs

  • Maple
    f:= proc(n) local B;
      B:= convert(n,base,2);
      igcd(op(select(t -> B[t]=1, [$1..ilog2(n)+1])))
    end proc:
    map(f, [$1..100]); # Robert Israel, Oct 13 2020
  • Mathematica
    Table[GCD@@Join@@Position[Reverse[IntegerDigits[n,2]],1],{n,100}]

Formula

Trivially, a(n) <= log_2(n). - Charles R Greathouse IV, Nov 15 2022

A337462 Number of pairwise coprime compositions of n, where a singleton is not considered coprime unless it is (1).

Original entry on oeis.org

1, 1, 1, 3, 6, 12, 21, 37, 62, 100, 159, 253, 402, 634, 983, 1491, 2224, 3280, 4813, 7043, 10270, 14888, 21415, 30585, 43400, 61204, 85747, 119295, 164834, 226422, 309663, 422301, 574826, 781236, 1060181, 1436367, 1942588, 2622078, 3531151, 4742315, 6348410
Offset: 0

Views

Author

Gus Wiseman, Sep 18 2020

Keywords

Comments

A composition of n is a finite sequence of positive integers summing to n.

Examples

			The a(1) = 1 through a(5) = 12 compositions:
  (1)  (1,1)  (1,2)    (1,3)      (1,4)
              (2,1)    (3,1)      (2,3)
              (1,1,1)  (1,1,2)    (3,2)
                       (1,2,1)    (4,1)
                       (2,1,1)    (1,1,3)
                       (1,1,1,1)  (1,3,1)
                                  (3,1,1)
                                  (1,1,1,2)
                                  (1,1,2,1)
                                  (1,2,1,1)
                                  (2,1,1,1)
                                  (1,1,1,1,1)
		

Crossrefs

A000740 counts the relatively prime instead of pairwise coprime version.
A101268 considers all singletons to be coprime, with strict case A337562.
A327516 is the unordered version.
A333227 ranks these compositions, with complement A335239.
A337461 counts these compositions of length 3.
A337561 is the strict case.
A051424 counts pairwise coprime or singleton partitions.
A101268 counts pairwise coprime or singleton compositions.
A178472 counts compositions with a common factor.
A305713 counts strict pairwise coprime partitions.
A328673 counts pairwise non-coprime partitions.
A333228 ranks compositions whose distinct parts are pairwise coprime.
A337667 counts pairwise non-coprime compositions.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],#=={}||CoprimeQ@@#&]],{n,0,10}]

Formula

For n > 1, a(n) = A101268(n) - 1.

A335235 Numbers k such that the k-th composition in standard order (A066099) is pairwise coprime, where a singleton is always considered coprime.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 23, 24, 25, 27, 28, 29, 30, 31, 32, 33, 35, 37, 38, 39, 41, 44, 47, 48, 49, 50, 51, 52, 55, 56, 57, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 71, 72, 75, 77, 78, 79, 80, 83, 89, 92, 95, 96, 97
Offset: 1

Views

Author

Gus Wiseman, May 28 2020

Keywords

Comments

The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence together with the corresponding compositions begins:
   1: (1)          20: (2,3)          48: (1,5)
   2: (2)          23: (2,1,1,1)      49: (1,4,1)
   3: (1,1)        24: (1,4)          50: (1,3,2)
   4: (3)          25: (1,3,1)        51: (1,3,1,1)
   5: (2,1)        27: (1,2,1,1)      52: (1,2,3)
   6: (1,2)        28: (1,1,3)        55: (1,2,1,1,1)
   7: (1,1,1)      29: (1,1,2,1)      56: (1,1,4)
   8: (4)          30: (1,1,1,2)      57: (1,1,3,1)
   9: (3,1)        31: (1,1,1,1,1)    59: (1,1,2,1,1)
  11: (2,1,1)      32: (6)            60: (1,1,1,3)
  12: (1,3)        33: (5,1)          61: (1,1,1,2,1)
  13: (1,2,1)      35: (4,1,1)        62: (1,1,1,1,2)
  14: (1,1,2)      37: (3,2,1)        63: (1,1,1,1,1,1)
  15: (1,1,1,1)    38: (3,1,2)        64: (7)
  16: (5)          39: (3,1,1,1)      65: (6,1)
  17: (4,1)        41: (2,3,1)        66: (5,2)
  18: (3,2)        44: (2,1,3)        67: (5,1,1)
  19: (3,1,1)      47: (2,1,1,1,1)    68: (4,3)
		

Crossrefs

The version counting partitions is A051424, with strict case A007360.
The version for binary indices is A087087.
The version counting compositions is A101268.
The version for prime indices is A302569.
The case without singletons is A333227.
The complement is A335236.
Numbers whose binary indices are pairwise coprime are A326675.
Coprime partitions are counted by A327516.
All of the following pertain to compositions in standard order:
- Length is A000120.
- The parts are row k of A066099.
- Sum is A070939.
- Product is A124758.
- Reverse is A228351
- GCD is A326674.
- Heinz number is A333219.
- LCM is A333226.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,100],Length[stc[#]]==1||CoprimeQ@@stc[#]&]

A335238 Numbers k such that the distinct parts of the k-th composition in standard order (A066099) are not pairwise coprime, where a singleton is not coprime unless it is (1).

Original entry on oeis.org

0, 2, 4, 8, 10, 16, 32, 34, 36, 40, 42, 64, 69, 70, 81, 88, 98, 104, 128, 130, 136, 138, 139, 141, 142, 160, 162, 163, 168, 170, 177, 184, 197, 198, 209, 216, 226, 232, 256, 260, 261, 262, 274, 276, 277, 278, 279, 282, 283, 285, 286, 288, 290, 292, 296, 321
Offset: 1

Views

Author

Gus Wiseman, May 28 2020

Keywords

Comments

We use the Mathematica definition for CoprimeQ, so a singleton is not considered coprime unless it is (1).
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence together with the corresponding compositions begins:
    0: ()          88: (2,1,4)      177: (2,1,4,1)
    2: (2)         98: (1,4,2)      184: (2,1,1,4)
    4: (3)        104: (1,2,4)      197: (1,4,2,1)
    8: (4)        128: (8)          198: (1,4,1,2)
   10: (2,2)      130: (6,2)        209: (1,2,4,1)
   16: (5)        136: (4,4)        216: (1,2,1,4)
   32: (6)        138: (4,2,2)      226: (1,1,4,2)
   34: (4,2)      139: (4,2,1,1)    232: (1,1,2,4)
   36: (3,3)      141: (4,1,2,1)    256: (9)
   40: (2,4)      142: (4,1,1,2)    260: (6,3)
   42: (2,2,2)    160: (2,6)        261: (6,2,1)
   64: (7)        162: (2,4,2)      262: (6,1,2)
   69: (4,2,1)    163: (2,4,1,1)    274: (4,3,2)
   70: (4,1,2)    168: (2,2,4)      276: (4,2,3)
   81: (2,4,1)    170: (2,2,2,2)    277: (4,2,2,1)
		

Crossrefs

The complement is A333228.
Not ignoring repeated parts gives A335239.
Singleton or pairwise coprime partitions are counted by A051424.
Singleton or pairwise coprime sets are ranked by A087087.
Coprime partitions are counted by A327516.
Non-coprime partitions are counted by A335240.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Sum is A070939.
- Product is A124758.
- Reverse is A228351
- GCD is A326674.
- Heinz number is A333219.
- LCM is A333226.
- Coprime compositions are A333227.
- Compositions whose distinct parts are coprime are A333228.
- Number of distinct parts is A334028.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,100],!CoprimeQ@@Union[stc[#]]&]

A253317 Indices in A261283 where records occur.

Original entry on oeis.org

0, 1, 2, 3, 8, 9, 10, 11, 128, 129, 130, 131, 136, 137, 138, 139, 32768, 32769, 32770, 32771, 32776, 32777, 32778, 32779, 32896, 32897, 32898, 32899, 32904, 32905, 32906, 32907, 2147483648, 2147483649, 2147483650, 2147483651, 2147483656, 2147483657
Offset: 1

Views

Author

Philippe Beaudoin, Dec 30 2014

Keywords

Comments

From Gus Wiseman, Dec 29 2023: (Start)
These are numbers whose binary indices are all powers of 2, where a binary index of n (row n of A048793) is any position of a 1 in its reversed binary expansion. For example, the terms together with their binary expansions and binary indices begin:
0: 0 ~ {}
1: 1 ~ {1}
2: 10 ~ {2}
3: 11 ~ {1,2}
8: 1000 ~ {4}
9: 1001 ~ {1,4}
10: 1010 ~ {2,4}
11: 1011 ~ {1,2,4}
128: 10000000 ~ {8}
129: 10000001 ~ {1,8}
130: 10000010 ~ {2,8}
131: 10000011 ~ {1,2,8}
136: 10001000 ~ {4,8}
137: 10001001 ~ {1,4,8}
138: 10001010 ~ {2,4,8}
139: 10001011 ~ {1,2,4,8}
For powers of 3 we have A368531.
(End)

Crossrefs

Cf. A053644 (most significant bit).
A048793 lists binary indices, length A000120, sum A029931.
A070939 gives length of binary expansion.
A096111 gives product of binary indices.

Programs

  • Maple
    a := proc(n) local k, A:
    A := [seq(0,i=1..n)]: A[1]:=0:
    for k from 1 to n-1 do
       A[k+1] := A[k-2^ilog2(k)+1]+2^(2^ilog2(k)-1): od:
    return A[n]: end proc: # Lorenzo Sauras Altuzarra, Dec 18 2019
    # second Maple program:
    a:= n-> (l-> add(l[i+1]*2^(2^i-1), i=0..nops(l)-1))(Bits[Split](n-1)):
    seq(a(n), n=1..38);  # Alois P. Heinz, Dec 13 2023
  • Mathematica
    Nest[Append[#1, #1[[-#2]] + 2^(#2 - 1)] & @@ {#, 2^(IntegerLength[Length[#], 2] - 1)} &, {0, 1}, 36] (* Michael De Vlieger, May 08 2020 *)
  • PARI
    a(n)={if(n<=1, 0, my(t=1<Andrew Howroyd, Dec 20 2019

Formula

a(1) = 0 and a(n) = a(n-A053644(n-1)) + 2^(A053644(n-1)-1). - Lorenzo Sauras Altuzarra, Dec 18 2019
a(n) = A358126(n-1) / 2. - Tilman Piesk, Dec 18 2022
a(2^n+1) = 2^(2^n-1) = A058891(n+1). - Gus Wiseman, Dec 29 2023
a(2^n) = A072639(n). - Gus Wiseman, Dec 29 2023
G.f.: 1/(1-x) * Sum_{k>=0} (2^(-1+2^k))*x^2^k/(1+x^2^k). - John Tyler Rascoe, May 22 2024

Extensions

Corrected reference in name from A253315 to A261283. - Tilman Piesk, Dec 18 2022

A295235 Numbers k such that the positions of the ones in the binary representation of k are in arithmetic progression.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 15, 16, 17, 18, 20, 21, 24, 28, 30, 31, 32, 33, 34, 36, 40, 42, 48, 56, 60, 62, 63, 64, 65, 66, 68, 72, 73, 80, 84, 85, 96, 112, 120, 124, 126, 127, 128, 129, 130, 132, 136, 144, 146, 160, 168, 170, 192, 224, 240, 248
Offset: 1

Views

Author

Rémy Sigrist, Nov 18 2017

Keywords

Comments

Also numbers k of the form Sum_{b=0..h-1} 2^(i+j*b) for some h >= 0, i >= 0, j > 0 (in fact, h = A000120(k), and if k > 0, i = A007814(k)).
There is a simple bijection between the finite sets of nonnegative integers in arithmetic progression and the terms of this sequence: s -> Sum_{i in s} 2^i; the term 0 corresponds to the empty set.
For any n > 0, A054519(n) gives the numbers of terms with n+1 digits in binary representation.
For any n >= 0, n is in the sequence iff 2*n is in the sequence.
For any n > 0, A000695(a(n)) is in the sequence.
The first prime numbers in the sequence are: 2, 3, 5, 7, 17, 31, 73, 127, 257, 8191, 65537, 131071, 262657, 524287, ...
This sequence contains the following sequences: A000051, A000079, A000225, A000668, A002450, A019434, A023001, A048645.
For any k > 0, 2^k - 2, 2^k - 1, 2^k, 2^k + 1 and 2^k + 2 are in the sequence (e.g., 14, 15, 16, 17, and 18).
Every odd term is a binary palindrome (and thus belongs to A006995).
Odd terms are A064896. - Robert Israel, Nov 20 2017

Examples

			The binary representation of the number 42 is "101010" and has ones evenly spaced, hence 42 appears in the sequence.
The first terms, alongside their binary representations, are:
   n  a(n)  a(n) in binary
  --  ----  --------------
   1    0           0
   2    1           1
   3    2          10
   4    3          11
   5    4         100
   6    5         101
   7    6         110
   8    7         111
   9    8        1000
  10    9        1001
  11   10        1010
  12   12        1100
  13   14        1110
  14   15        1111
  15   16       10000
  16   17       10001
  17   18       10010
  18   20       10100
  19   21       10101
  20   24       11000
		

Crossrefs

Cf. A029931, A048793 (binary indices triangle), A070939, A291166, A325328 (prime indices rather than binary indices), A326669, A326675.

Programs

  • Maple
    f:= proc(d) local i,j,k;
      op(sort([seq(seq(add(2^(d-j*k),k=0..m),m=1..d/j),j=1..d),2^(d+1)]))
    end proc:
    0,1,seq(f(d),d=0..10); # Robert Israel, Nov 20 2017
  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[100],SameQ@@Differences[bpe[#]]&] (* Gus Wiseman, Jul 22 2019 *)
  • PARI
    is(n) = my(h=hammingweight(n)); if(h<3, return(1), my(i=valuation(n,2),w=#binary(n)); if((w-i-1)%(h-1)==0, my(j=(w-i-1)/(h-1)); return(sum(k=0,h-1,2^(i+j*k))==n), return(0)))
Previous Showing 11-20 of 36 results. Next