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.

Showing 1-10 of 12 results. Next

A096111 If n = 2^k - 1, then a(n) = k+1, otherwise a(n) = (A000523(n)+1)*a(A053645(n)).

Original entry on oeis.org

1, 2, 2, 3, 3, 6, 6, 4, 4, 8, 8, 12, 12, 24, 24, 5, 5, 10, 10, 15, 15, 30, 30, 20, 20, 40, 40, 60, 60, 120, 120, 6, 6, 12, 12, 18, 18, 36, 36, 24, 24, 48, 48, 72, 72, 144, 144, 30, 30, 60, 60, 90, 90, 180, 180, 120, 120, 240, 240, 360, 360, 720, 720, 7, 7, 14, 14, 21, 21
Offset: 0

Views

Author

Amarnath Murthy, Jun 29 2004

Keywords

Comments

Each n > 1 occurs 2*A045778(n) times in the sequence.
f(n+2^k) = (k+1)*f(n) if 2^k > n+1. - Robert Israel, Apr 25 2016
If the binary indices of n (row n of A048793) are the positions 1's in its reversed binary expansion, then a(n) is the product of all binary indices of n + 1. The number of binary indices of n is A000120(n), their sum is A029931(n), and their average is A326699(n)/A326700(n). - Gus Wiseman, Jul 27 2019

Crossrefs

Permutation of A096115, i.e. a(n) = A096115(A122198(n+1)) [Note the different starting offsets]. Bisection: A121663. Cf. A096113, A052330.
Cf. A029931.

Programs

  • Maple
    f:= proc(n) local L;
        L:= convert(2*n+2,base,2);
        convert(subs(0=NULL,zip(`*`,L, [$0..nops(L)-1])),`*`);
    end proc:
    map(f, [$0..100]); # Robert Israel, Apr 25 2016
  • Mathematica
    CoefficientList[(Product[1 + k x^(2^(k - 1)), {k, 7}] - 1)/x, x] (* Michael De Vlieger, Apr 08 2016 *)
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];Table[Times@@bpe[n+1],{n,0,100}] (* Gus Wiseman, Jul 26 2019 *)
  • PARI
    N=166; q='q+O('q^N);
    gf= (prod(n=1,1+ceil(log(N)/log(2)), 1+n*q^(2^(n-1)) ) - 1) / q;
    Vec(gf)
    /* Joerg Arndt, Oct 06 2012 */
  • Scheme
    (define (A096111 n) (cond ((pow2? (+ n 1)) (+ 2 (A000523 n))) (else (* (+ 1 (A000523 n)) (A096111 (A053645 n))))))
    (define (pow2? n) (and (> n 0) (zero? (A004198bi n (- n 1)))))
    

Formula

G.f.: ( prod(k>=1, 1+k*x^(2^(k-1)) )- 1 ) / x. - Vladeta Jovovic, Nov 08 2005
a(n) is the product of the exponents in the binary expansion of 2*n + 2. - Peter Kagey, Apr 24 2016

Extensions

Edited, extended and Scheme code added by Antti Karttunen, Aug 25 2006

A291166 Connected Haar graph numbers.

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 33, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81
Offset: 1

Views

Author

Eric W. Weisstein, Aug 19 2017

Keywords

Comments

Complement of A291165.
These appear to be numbers whose positions of 1's in their reversed binary expansion are relatively prime. If so, this sequence lists all positions of 1's in A326674. Numbers whose positions of 1's in their reversed binary expansion are pairwise coprime (as opposed to relatively prime) are A326675. - Gus Wiseman, Jul 19 2019

Crossrefs

A326675 The positions of 1's in the reversed binary expansion of n are pairwise coprime, where a singleton is not coprime unless it is {1}.

Original entry on oeis.org

1, 3, 5, 6, 7, 9, 12, 13, 17, 18, 19, 20, 21, 22, 23, 24, 25, 28, 29, 33, 48, 49, 65, 66, 67, 68, 69, 70, 71, 72, 73, 76, 77, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 92, 93, 96, 97, 112, 113, 129, 132, 133, 144, 145, 148, 149, 192, 193, 196, 197, 208, 209, 212
Offset: 1

Views

Author

Gus Wiseman, Jul 17 2019

Keywords

Examples

			41 has reversed binary expansion (1,0,0,1,0,1) with positions of 1's being {1,4,6}, which are not pairwise coprime, so 41 is not in the sequence.
		

Crossrefs

Equals the complement of A131577 in A087087.
Numbers whose prime indices are pairwise coprime are A302696.
Taking relatively prime instead of pairwise coprime gives A291166.

Programs

  • Maple
    extend:= proc(L) local C,c;
      C:= select(t -> andmap(s -> igcd(s,t)=1, L), [$1..L[-1]-1]);
      L, seq(procname([op(L),c]),c=C)
    end proc:
    g:= proc(L) local i;
      add(2^(i-1),i=L)
    end proc:
    map(g, [[1],seq(extend([k])[2..-1], k=2..10)]); # Robert Israel, Jul 19 2019
  • Mathematica
    Select[Range[100],CoprimeQ@@Join@@Position[Reverse[IntegerDigits[#,2]],1]&]
  • PARI
    is(n) = my (p=1); while (n, my (o=1+valuation(n,2)); if (gcd(p,o)>1, return (0), n-=2^(o-1); p*=o)); return (1) \\ Rémy Sigrist, Jul 19 2019

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

A326669 Numbers k such that the average position of the ones in the binary expansion of k is an integer.

Original entry on oeis.org

1, 2, 4, 5, 7, 8, 10, 14, 16, 17, 20, 21, 27, 28, 31, 32, 34, 35, 39, 40, 42, 49, 54, 56, 57, 62, 64, 65, 68, 70, 73, 78, 80, 84, 85, 93, 98, 99, 107, 108, 112, 114, 119, 124, 127, 128, 130, 133, 136, 140, 141, 146, 147, 155, 156, 160, 161, 167, 168, 170, 175
Offset: 1

Views

Author

Gus Wiseman, Jul 17 2019

Keywords

Comments

These are numbers whose exponents in their representation as a sum of distinct powers of 2 have integer average.

Examples

			42 is in the sequence because 42 = 2^1 + 2^3 + 2^5 and the average of {1,3,5} is 3, an integer.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[100],IntegerQ[Mean[Join@@Position[IntegerDigits[#,2],1]]]&]
  • PARI
    isok(m) = my(b=binary(m)); denominator(vecsum(Vec(select(x->(x==1), b, 1)))/hammingweight(m)) == 1; \\ Michel Marcus, Jul 02 2021

A359402 Numbers whose binary expansion and reversed binary expansion have the same sum of positions of 1's, where positions in a sequence are read starting with 1 from the left.

Original entry on oeis.org

0, 1, 3, 5, 7, 9, 15, 17, 21, 27, 31, 33, 45, 51, 63, 65, 70, 73, 78, 85, 93, 99, 107, 119, 127, 129, 150, 153, 165, 189, 195, 219, 231, 255, 257, 266, 273, 282, 294, 297, 310, 313, 325, 334, 341, 350, 355, 365, 371, 381, 387, 397, 403, 413, 427, 443, 455, 471
Offset: 1

Views

Author

Gus Wiseman, Jan 05 2023

Keywords

Comments

Also numbers whose binary expansion and reversed binary expansion have the same sum of partial sums.
Also numbers whose average position of a 1 in their binary expansion is (c+1)/2, where c is the number of digits.
Conjecture: Also numbers whose binary expansion has as least squares fit a line of zero slope, counted by A222955.

Examples

			The binary expansion of 70 is (1,0,0,0,1,1,0), with positions of 1's {1,5,6}, while the reverse positions are {2,3,7}. Both sum to 12, so 70 is in the sequence.
		

Crossrefs

Binary words of this type appear to be counted by A222955.
For greater instead of equal sums we have A359401.
These are the indices of 0's in A359495.
A030190 gives binary expansion, reverse A030308.
A048793 lists partial sums of reversed standard compositions, sums A029931.
A070939 counts binary digits, 1's A000120.
A326669 lists numbers with integer mean position of a 1 in binary expansion.

Programs

  • Mathematica
    Select[Range[0,100],#==0||Mean[Join@@Position[IntegerDigits[#,2],1]]==(IntegerLength[#,2]+1)/2&]
  • Python
    from functools import reduce
    from itertools import count, islice
    def A359402_gen(startvalue=0): # generator of terms
        return filter(lambda n:(r:=reduce(lambda c, d:(c[0]+d[0]*(e:=int(d[1])),c[1]+e),enumerate(bin(n)[2:],start=1),(0,0)))[0]<<1==(n.bit_length()+1)*r[1],count(max(startvalue,0)))
    A359402_list = list(islice(A359402_gen(),30)) # Chai Wah Wu, Jan 08 2023

Formula

A230877(a(n)) = A029931(a(n)).

A359495 Sum of positions of 1's in binary expansion minus sum of positions of 1's in reversed binary expansion, where positions in a sequence are read starting with 1 from the left.

Original entry on oeis.org

0, 0, -1, 0, -2, 0, -2, 0, -3, 0, -2, 1, -4, -1, -3, 0, -4, 0, -2, 2, -4, 0, -2, 2, -6, -2, -4, 0, -6, -2, -4, 0, -5, 0, -2, 3, -4, 1, -1, 4, -6, -1, -3, 2, -5, 0, -2, 3, -8, -3, -5, 0, -7, -2, -4, 1, -9, -4, -6, -1, -8, -3, -5, 0, -6, 0, -2, 4, -4, 2, 0, 6
Offset: 0

Views

Author

Gus Wiseman, Jan 05 2023

Keywords

Comments

Also the sum of partial sums of reversed binary expansion minus sum of partial sums of binary expansion.

Examples

			The binary expansion of 158 is (1,0,0,1,1,1,1,0), with positions of 1's {1,4,5,6,7} with sum 23, reversed {2,3,4,5,8} with sum 22, so a(158) = 1.
		

Crossrefs

Indices of positive terms are A359401.
Indices of 0's are A359402.
A030190 gives binary expansion, reverse A030308.
A070939 counts binary digits.
A230877 adds up positions of 1's in binary expansion, reverse A029931.

Programs

  • Maple
    a:= n-> (l-> add(i*(l[-i]-l[i]), i=1..nops(l)))(Bits[Split](n)):
    seq(a(n), n=0..127);  # Alois P. Heinz, Jan 09 2023
  • Mathematica
    sap[q_]:=Sum[q[[i]]*(2i-Length[q]-1),{i,Length[q]}];
    Table[sap[IntegerDigits[n,2]],{n,0,100}]
  • Python
    def A359495(n):
        k = n.bit_length()-1
        return sum((i<<1)-k for i, j in enumerate(bin(n)[2:]) if j=='1') # Chai Wah Wu, Jan 09 2023

Formula

a(n) = A029931(n) - A230877(n).
If n = Sum_{i=1..k} q_i * 2^(i-1), then a(n) = Sum_{i=1..k} q_i * (2i-k-1).

A326672 The positions of ones in the binary expansion of n have integer geometric mean.

Original entry on oeis.org

1, 2, 4, 8, 9, 13, 16, 18, 26, 32, 36, 52, 64, 72, 104, 128, 144, 208, 256, 257, 288, 321, 416, 512, 514, 576, 642, 832, 1024, 1028, 1152, 1284, 1664, 2048, 2056, 2304, 2568, 3328, 4096, 4112, 4608, 5136, 6656, 8192, 8224, 9216, 10272, 13312, 16384, 16448
Offset: 1

Views

Author

Gus Wiseman, Jul 17 2019

Keywords

Crossrefs

Partitions with integer geometric mean are A067539.
Subsets with integer geometric mean are A326027.
Factorizations with integer geometric mean are A326028.
Numbers whose binary expansion positions have integer mean are A326669.
Numbers whose binary expansion positions are relatively prime are A326674.
Numbers whose reversed binary expansion positions have integer geometric mean are A326673.

Programs

  • Mathematica
    Select[Range[100],IntegerQ[GeometricMean[Join@@Position[IntegerDigits[#,2],1]]]&]

A359401 Nonnegative integers whose sum of positions of 1's in their binary expansion is greater than the sum of positions of 1's in their reversed binary expansion, where positions in a sequence are read starting with 1 from the left.

Original entry on oeis.org

11, 19, 23, 35, 37, 39, 43, 47, 55, 67, 69, 71, 75, 77, 79, 83, 87, 91, 95, 103, 111, 131, 133, 134, 135, 137, 139, 141, 142, 143, 147, 149, 151, 155, 157, 158, 159, 163, 167, 171, 173, 175, 179, 183, 187, 191, 199, 203, 207, 215, 223, 239, 259, 261, 262, 263
Offset: 1

Views

Author

Gus Wiseman, Jan 05 2023

Keywords

Comments

First differs from A161601 in having 134, with binary expansion (1,0,0,0,0,1,1,0), positions of 1's 1 + 6 + 7 = 14, reversed 2 + 3 + 8 = 13.

Crossrefs

Indices of positive terms in A359495; indices of 0's are A359402.
A030190 gives binary expansion, reverse A030308.
A070939 counts binary digits.
A230877 adds up positions of 1's in binary expansion, reverse A029931.
A326669 lists numbers with integer mean position of a 1 in binary expansion.

Programs

  • Mathematica
    sap[q_]:=Sum[q[[i]]*(2i-Length[q]-1),{i,Length[q]}];
    Select[Range[0,100],sap[IntegerDigits[#,2]]>0&]

Formula

A230877(a(n)) > A029931(a(n)).

A327368 The positions of ones in the reversed binary expansion of n have integer mean and integer geometric mean.

Original entry on oeis.org

1, 2, 4, 8, 16, 32, 64, 128, 130, 256, 257, 512, 1024, 2048, 2084, 2316, 4096, 8192, 16384, 32768, 32776, 32777, 65536, 131072, 131074, 131200, 131457, 131462, 133390, 165920, 262144, 524288, 1048576, 2097152, 4194304, 8388608, 8388640, 8388897, 8390688, 8519840
Offset: 1

Views

Author

Gus Wiseman, Sep 27 2019

Keywords

Examples

			The sequence of terms together with their binary indices begins:
  2      {2}
  4      {3}
  8      {4}
  16     {5}
  32     {6}
  64     {7}
  128    {8}
  130    {2,8}
  256    {9}
  257    {1,9}
  512    {10}
  1024   {11}
  2048   {12}
  2084   {3,6,12}
  2316   {3,4,9,12}
  4096   {13}
  8192   {14}
  16384  {15}
  32768  {16}
  32776  {4,16}
		

Crossrefs

A superset of A327777.
Numbers whose binary indices have integer mean: A326669
Numbers whose binary indices have integer geometric mean: A326673

Programs

  • Mathematica
    bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[1000],IntegerQ[Mean[bpe[#]]]&&IntegerQ[GeometricMean[bpe[#]]]&]
  • PARI
    ok(n)={my(s=0,p=1,k=0); for(i=0, logint(n,2), if(bittest(n,i), s+=i+1; p*=i+1; k++)); s%k==0 && ispower(p,k)}
    { for(n=1, 10^7, if(ok(n), print1(n, ", "))) } \\ Andrew Howroyd, Sep 29 2019

Extensions

a(33)-a(40) from Andrew Howroyd, Sep 29 2019
Showing 1-10 of 12 results. Next