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 47 results. Next

A339216 Numbers k such that k and k+2 are both binary self numbers (A010061).

Original entry on oeis.org

4, 13, 21, 30, 37, 46, 54, 78, 86, 95, 102, 111, 119, 128, 133, 142, 150, 159, 166, 175, 183, 207, 215, 224, 231, 240, 248, 270, 278, 287, 294, 303, 311, 335, 343, 352, 359, 368, 376, 385, 390, 399, 407, 416, 423, 432, 440, 464, 472, 481, 488, 497, 505, 526, 534
Offset: 1

Views

Author

Amiram Eldar, Nov 27 2020

Keywords

Comments

The least difference between consecutive binary self numbers is 2 (see Macris's proof at A010061).

Examples

			4 is a term since 4 and 6 = 4 + 2 are both binary self numbers.
		

Crossrefs

Cf. A010061.

Programs

  • Mathematica
    s[n_] := n + DigitCount[n, 2, 1]; m = 550; c = Complement[Range[m], Array[s, m]]; d = Differences[c]; ind = Position[d, 2] // Flatten; c[[ind]]

A339239 Binary self numbers (A010061) with a record gap to the next binary self number.

Original entry on oeis.org

1, 6, 63, 250, 131070, 1048574, 33554426, 17179869180
Offset: 1

Views

Author

Amiram Eldar, Nov 28 2020

Keywords

Comments

The corresponding gaps are 3, 7, 8, 20, 21, 24, 37, 42, ...
a(9) <= 288230376151711738.
Apparently, the records gaps occur for pairs of consecutive binary self numbers with a power of 2 between them. If this is generally true, then the next terms are 288230376151711738, 147573952589676412923, 37778931862957161709564, 10633823966279326983230456482242756602, 5444517870735015415413993718908291383294, 43556142965880123323311949751266331066367, ..., with the corresponding gaps 70, 77, 83, 135, 136, 137, ...

Examples

			The first 4 binary self numbers are 1, 4, 6 and 13. The gaps between them are 3, 2 and 7. The record gaps are 3 and 7, and the corresponding terms are 1 and 6.
		

Crossrefs

Programs

  • Mathematica
    s[n_] := n + DigitCount[n, 2, 1]; selfQ[n_] := AllTrue[Range[n, n - Floor@Log2[n], -1], s[#] != n &]; dm = 0; seq = {}; n1 = 1; Do[If[selfQ[n], d = n - n1; If[d > dm, dm = d; AppendTo[seq, n1]]; n1 = n], {n, 2, 150000}]; seq

A003052 Self numbers or Colombian numbers (numbers that are not of the form m + sum of digits of m for any m).

Original entry on oeis.org

1, 3, 5, 7, 9, 20, 31, 42, 53, 64, 75, 86, 97, 108, 110, 121, 132, 143, 154, 165, 176, 187, 198, 209, 211, 222, 233, 244, 255, 266, 277, 288, 299, 310, 312, 323, 334, 345, 356, 367, 378, 389, 400, 411, 413, 424, 435, 446, 457, 468, 479, 490, 501, 512, 514, 525
Offset: 1

Views

Author

Keywords

Comments

From Amiram Eldar, Nov 28 2020: (Start)
The term "self numbers" was coined by Kaprekar (1959). The term "Colombian number" was coined by Recamán (1973) of Bogota, Colombia.
The asymptotic density of this sequence is approximately 0.0977778 (Guaraldo, 1978). (End)

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 2.24.
  • Martin Gardner, Time Travel and Other Mathematical Bewilderments. Freeman, NY, 1988, p. 116.
  • V. S. Joshi, A note on self-numbers. Volume dedicated to the memory of V. Ramaswami Aiyar. Math. Student, Vol. 39 (1971), pp. 327-328. MR0330032 (48 #8371).
  • D. R. Kaprekar, Puzzles of the Self-Numbers. 311 Devlali Camp, Devlali, India, 1959.
  • D. R. Kaprekar, The Mathematics of the New Self Numbers, Privately Printed, 311 Devlali Camp, Devlali, India, 1963.
  • D. R. Kaprekar, The Mathematics of the New Self Numbers (Part V). 311 Devlali Camp, Devlali, India, 1967.
  • Bernardo Recamán, The Bogota Puzzles, Dover Publications, Inc., 2020, chapter 36, p. 33.
  • József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, pp. 384-386.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Author?, J. Recreational Math., vol. 23, no. 1, p. 244, 1991.

Crossrefs

For self primes, i.e., self numbers which are primes, see A006378.
Complement of A176995.
See A010061 for the binary version, A283002 for a base-100 version.
Cf. A247104 (subsequence of squarefree terms).
Cf. A377472 for first differences, A377474 for indices where new differences appear.

Programs

  • Haskell
    a003052 n = a003052_list !! (n-1)
    a003052_list = filter ((== 0) . a230093) [1..]
    -- Reinhard Zumkeller, Oct 11 2013, Aug 21 2011
  • Maple
    isA003052 := proc(n) local k ; for k from 0 to n do if k+A007953(k) = n then RETURN(false): fi; od: RETURN(true) ; end:
    A003052 := proc(n) option remember; if n = 1 then 1; else for a from procname(n-1)+1 do if isA003052(a) then RETURN(a) ; fi; od; fi; end: # R. J. Mathar, Jul 27 2009
  • Mathematica
    nn = 525; Complement[Range[nn], Union[Table[n + Total[IntegerDigits[n]], {n, nn}]]] (* T. D. Noe, Mar 31 2013 *)
  • PARI
    is_A003052(n)={for(i=1,min(n\2,9*#digits(n)), sumdigits(n-i)==i && return); n}  \\ M. F. Hasler, Mar 20 2011, updated Nov 08 2018
    
  • PARI
    is(n) = {if(n < 30, return((n < 10 && n%2 == 1) || n == 20)); qd = 1 + logint(n, 10); r = 1 + (n-1)%9; h = (r + 9 * (r%2))/2; ld = 10; while(h + 9*qd >= n % ld, ld*=10); vs = vecsum(digits(n \ ld)); n %= ld; for(i = 0, qd, if(vs + vecsum(digits(n - h - 9*i)) == h + 9*i, return(0))); 1} \\ David A. Corneth, Aug 20 2020
    

Formula

A230093(a(n)) = 0. - Reinhard Zumkeller, Oct 11 2013
In fact this defines the sequence: x is in the sequence iff A230093(x) = 0. - M. F. Hasler, Nov 08 2018

Extensions

More terms from James Sellers, Jul 06 2000

A055938 Integers not generated by b(n) = b(floor(n/2)) + n (complement of A005187).

Original entry on oeis.org

2, 5, 6, 9, 12, 13, 14, 17, 20, 21, 24, 27, 28, 29, 30, 33, 36, 37, 40, 43, 44, 45, 48, 51, 52, 55, 58, 59, 60, 61, 62, 65, 68, 69, 72, 75, 76, 77, 80, 83, 84, 87, 90, 91, 92, 93, 96, 99, 100, 103, 106, 107, 108, 111, 114, 115, 118, 121, 122, 123, 124, 125, 126, 129
Offset: 1

Views

Author

Alford Arnold, Jul 21 2000

Keywords

Comments

Note that the lengths of the consecutive runs in a(n) form sequence A001511.
Integers that are not a sum of distinct integers of the form 2^k-1. - Vladeta Jovovic, Jan 24 2003
Also n! never ends in this many 0's in base 2 - Carl R. White, Jan 21 2008
A079559(a(n)) = 0. - Reinhard Zumkeller, Mar 18 2009
These numbers are dead-end points when trying to apply the iterated process depicted in A071542 in reverse, i.e. these are positive integers i such that there does not exist k with A000120(i+k)=k. See also comments at A179016. - Antti Karttunen, Oct 26 2012
Conjecture: a(n)=b(n) defined as b(1)=2, for n>1, b(n+1)=b(n)+1 if n is already in the sequence, b(n+1)=b(n)+3 otherwise. If so, then see Cloitre comment in A080578. - Ralf Stephan, Dec 27 2013
Numbers n for which A257265(m) = 0. - Reinhard Zumkeller, May 06 2015. Typo corrected by Antti Karttunen, Aug 08 2015
Numbers which have a 2 in their skew-binary representation (cf. A169683). - Allan C. Wechsler, Feb 28 2025

Examples

			Since A005187 begins 0 1 3 4 7 8 10 11 15 16 18 19 22 23 25 26 31... this sequence begins 2 5 6 9 12 13 14 17 20 21
		

Crossrefs

Complement of A005187. Setwise difference of A213713 and A213717.
Row 1 of arrays A257264, A256997 and also of A255557 (when prepended with 1). Equally: column 1 of A256995 and A255555.
Cf. also arrays A254105, A254107 and permutations A233276, A233278.
Left inverses: A234017, A256992.
Gives positions of zeros in A213714, A213723, A213724, A213731, A257265, positions of ones in A213725-A213727 and A256989, positions of nonzeros in A254110.
Cf. also A010061 (integers that are not a sum of distinct integers of the form 2^k+1).
Analogous sequence for factorial base number system: A219658, for Fibonacci number system: A219638, for base-3: A096346. Cf. also A136767-A136774.

Programs

  • Haskell
    a055938 n = a055938_list !! (n-1)
    a055938_list = concat $
       zipWith (\u v -> [u+1..v-1]) a005187_list $ tail a005187_list
    -- Reinhard Zumkeller, Nov 07 2011
    
  • Mathematica
    a[0] = 0; a[1] = 1; a[n_Integer] := a[Floor[n/2]] + n; b = {}; Do[ b = Append[b, a[n]], {n, 0, 105}]; c =Table[n, {n, 0, 200}]; Complement[c, b]
    (* Second program: *)
    t = Table[IntegerExponent[(2n)!, 2], {n, 0, 100}]; Complement[Range[t // Last], t] (* Jean-François Alcover, Nov 15 2016 *)
  • PARI
    L=listcreate();for(n=1,1000,for(k=2*n-hammingweight(n)+1,2*n+1-hammingweight(n+1),listput(L,k)));Vec(L) \\ Ralf Stephan, Dec 27 2013
    
  • Python
    def a053644(n): return 0 if n==0 else 2**(len(bin(n)[2:]) - 1)
    def a043545(n):
        x=bin(n)[2:]
        return int(max(x)) - int(min(x))
    def a079559(n): return 1 if n==0 else a043545(n + 1)*a079559(n + 1 - a053644(n + 1))
    print([n for n in range(1, 201) if a079559(n)==0]) # Indranil Ghosh, Jun 11 2017, after the comment by Reinhard Zumkeller
  • Scheme
    ;; utilizing COMPLEMENT-macro from Antti Karttunen's IntSeq-library)
    (define A055938 (COMPLEMENT 1 A005187))
    ;; Antti Karttunen, Aug 08 2015
    

Formula

a(n) = A080578(n+1) - 2 = A080468(n+1) + 2*n (conjectured). - Ralf Stephan, Dec 27 2013
From Antti Karttunen, Aug 08 2015: (Start)
Other identities. For all n >= 1:
A234017(a(n)) = n.
A256992(a(n)) = n.
A257126(n) = a(n) - A005187(n).
(End)

Extensions

More terms from Robert G. Wilson v, Jul 24 2000

A092391 a(n) = n + wt(n), where wt(n) = A000120(n) = binary weight of n.

Original entry on oeis.org

0, 2, 3, 5, 5, 7, 8, 10, 9, 11, 12, 14, 14, 16, 17, 19, 17, 19, 20, 22, 22, 24, 25, 27, 26, 28, 29, 31, 31, 33, 34, 36, 33, 35, 36, 38, 38, 40, 41, 43, 42, 44, 45, 47, 47, 49, 50, 52, 50, 52, 53, 55, 55, 57, 58, 60, 59, 61, 62, 64, 64, 66, 67, 69, 65, 67, 68, 70, 70, 72, 73, 75
Offset: 0

Views

Author

Reinhard Zumkeller, May 08 2004

Keywords

Crossrefs

A010061 gives the numbers not occurring in this sequence. A228082 gives the terms of this sequence sorted into ascending order, with duplicates removed. A228085(n) gives the number of times n occurs in this sequence.

Programs

Formula

a(n) = n + A000120(n).
A010062(n+1) = a(A010062(n)).
G.f.: (1/(1 - x))*Sum_{k>=0} (2^k + 1)*x^(2^k)/(1 + x^(2^k)). - Ilya Gutkovskiy, Jul 23 2017

A010062 a(0)=1; thereafter a(n+1) = a(n) + number of 1's in binary representation of a(n).

Original entry on oeis.org

1, 2, 3, 5, 7, 10, 12, 14, 17, 19, 22, 25, 28, 31, 36, 38, 41, 44, 47, 52, 55, 60, 64, 65, 67, 70, 73, 76, 79, 84, 87, 92, 96, 98, 101, 105, 109, 114, 118, 123, 129, 131, 134, 137, 140, 143, 148, 151, 156, 160, 162, 165, 169, 173, 178, 182, 187, 193, 196, 199, 204
Offset: 0

Views

Author

Leonid Broukhis, Mar 15 1996

Keywords

Comments

Sequence A230297 (and A157845 without initial term) converted from binary to decimal, cf. formula. - M. F. Hasler, Nov 18 2019

Examples

			a(7) = 14 because a(6) = 12, which is 1100 in binary (having 2 on bits), and 12 + 2 = 14.
a(8) = 17 because a(7) = 14, which is 1110 in binary (having 3 on bits), and 14 + 3 = 17.
		

Crossrefs

First row of A228083.
For the base-10 analog see A004207.
Cf. A000120, A010061, A092391, A229167, A096303, A229743, A229744, A230297 (this sequence written in binary), A230298 (read mod 2).
See A230088 for partial sums.
Equals A028897 o A230297 = A028897 o A157845 (up to offset); see also A007088.

Programs

Formula

a(n) = (n/2)*log n + O(n*sqrt(log n * loglog n)), where log means log_2. In particular, a(n) ~ (n/2)*log n. [Stolarsky]
a(n + 1) = A092391(a(n)) = a(n) + A000120(a(n)). - Reinhard Zumkeller, May 27 2012, May 08 2004; corrected thanks to a notice by Lambert Herrgesell
a(n) = A028897(A230297(n)) = A028897(A157845(n+1)). - M. F. Hasler, Nov 18 2019

Extensions

More terms from Benoit Cloitre, Jun 02 2002
Stolarsky reference from Matthew C. Russell, Oct 08 2013

A010064 Base 4 self or Colombian numbers (not of form k + sum of base 4 digits of k).

Original entry on oeis.org

1, 3, 8, 13, 18, 20, 25, 30, 35, 37, 42, 47, 52, 54, 59, 64, 73, 78, 83, 85, 90, 95, 100, 102, 107, 112, 117, 119, 124, 129, 138, 143, 148, 150, 155, 160, 165, 167, 172, 177, 182, 184, 189, 194, 203, 208, 213, 215, 220, 225, 230, 232, 237, 242, 247, 249, 254
Offset: 1

Views

Author

Keywords

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 2.24, pp. 179-180.
  • József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, pp. 384-386.

Crossrefs

Programs

  • Mathematica
    s[n_] := n + Plus @@ IntegerDigits[n, 4]; m = 250; Complement[Range[m], Array[s, m]] (* Amiram Eldar, Nov 28 2020 *)

A228085 a(n) = number of distinct k which satisfy n = k + wt(k), where wt(k) (A000120) gives the number of 1's in binary representation of a nonnegative integer k.

Original entry on oeis.org

1, 0, 1, 1, 0, 2, 0, 1, 1, 1, 1, 1, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0, 2, 0, 1, 1, 1, 1, 1, 1, 0, 2, 0, 2, 1, 1, 2, 0, 2, 0, 1, 1, 1, 1, 1, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0, 2, 0, 1, 1, 1, 1, 1, 1, 0, 2, 1, 1, 2, 1, 1, 2, 0, 1, 1, 1, 1, 1, 1, 0, 2, 0, 1, 2, 0, 2, 1, 0
Offset: 0

Views

Author

Antti Karttunen, Aug 09 2013

Keywords

Comments

wt(k) is also called bitcount(k).
a(n) = number of times n occurs in A092391.
The first 2 occurs at n = A230303(2) = 5 (as we have two solutions A092391(3) = A092391(4) = 5).
The first 3 occurs at n = A230303(3) = 129 (as we have three solutions A092391(123) = A092391(124) = A092391(128) = 129).
The first 4 occurs at n = A230303(4) = 4102, where we have solutions A092391(4091) = A092391(4092) = A092391(4099) = A092391(4100) = 4102.
For n>=1, a(2^n) = a(n-1) since an integer k = m is a solution to n-1 = m + wt(m) if and only if k = 2^n - 1 - m is a solution to 2^n = k + wt(k). - Max Alekseyev, Feb 23 2021

Crossrefs

A010061 gives the position of zeros, A228082 the positions of nonzeros, A228088 the positions of ones.
Cf. A000120, A010062, A092391, A228086, A228087, A228091 (positions of 2's), A227643, A230058, A230092 (positions of 3's), A230093, A227915 (positions of 4's), A070939, A230303.

Programs

  • Haskell
    a228085 n = length $ filter ((== n) . a092391) [n - a070939 n .. n]
    -- Reinhard Zumkeller, Oct 13 2013
  • Maple
    For Maple code see A230091. - N. J. A. Sloane, Oct 10 2013
    # Find all inverses of m under x -> x + wt(x) - N. J. A. Sloane, Oct 19 2013
    A000120 := proc(n) local w, m, i; w := 0; m := n; while m > 0 do i := m mod 2; w := w+i; m := (m-i)/2; od; w; end: wt := A000120;
    F:=proc(m) local ans,lb,n,i;
    lb:=m-ceil(log(m+1)/log(2)); ans:=[];
    for n from max(1,lb) to m do if (n+wt(n)) = m then ans:=[op(ans),n]; fi; od:
    [seq(ans[i],i=1..nops(ans))];
    end;
  • Mathematica
    nmax = 8191; Clear[a]; a[_] = 0;
    Scan[Set[a[#[[1]]], #[[2]]]&, Tally[Table[n + DigitCount[n, 2, 1], {n, 0, nmax}]]];
    a /@ Range[0, nmax] (* Jean-François Alcover, Oct 29 2019 *)
    a[n_] := Module[{k, cnt = 0}, For[k = n - Floor[Log[2, n]] - 1, k < n, k++, If[n == k + DigitCount[k, 2, 1], cnt++]]; cnt];
    a /@ Range[0, 100] (* Jean-François Alcover, Nov 28 2020 *)

A228082 Numbers that are of the form k + sum of binary digits of k for some nonnegative integer k.

Original entry on oeis.org

0, 2, 3, 5, 7, 8, 9, 10, 11, 12, 14, 16, 17, 19, 20, 22, 24, 25, 26, 27, 28, 29, 31, 33, 34, 35, 36, 38, 40, 41, 42, 43, 44, 45, 47, 49, 50, 52, 53, 55, 57, 58, 59, 60, 61, 62, 64, 65, 66, 67, 68, 69, 70, 72, 73, 74, 75, 76, 77, 79, 81, 82, 84, 85, 87, 89, 90
Offset: 1

Views

Author

Antti Karttunen, Aug 09 2013

Keywords

Comments

Complement of A010061.
Obtained when A092391 is sorted and duplicates are removed.
The asymptotic density of this sequence is 1 - (1/8) * (Sum_{n>=1} 1/2^a(n))^2 = 1 - A242403 = 0.747339... - Amiram Eldar, Nov 28 2020

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, Section 2.24, pp. 179-180.
  • József Sándor and Borislav Crstici, Handbook of Number theory II, Kluwer Academic Publishers, 2004, Chapter 4, p. 384-386.
  • G. Troi and U. Zannier, Note on the density constant in the distribution of self-numbers, Bolletino U. M. I. (7) 9-A (1995), 143-148.

Crossrefs

Numbers that occur to the right of the leftmost column of A228083. Positions of nonzeros in A228085. A superset of A228088.
The even terms are the first row of A350601.

Programs

  • Haskell
    a228082 n = a228082_list !! (n-1)
    a228082_list = 0 : filter ((> 0) . a228085) [1..]
    -- Reinhard Zumkeller, Oct 13 2013
  • Mathematica
    Table[n + Total[IntegerDigits[n, 2]], {n, 0, 100}] // Union (* Jean-François Alcover, Sep 03 2013 *)

A228088 Numbers n for which there is a unique k which satisfies n = k + wt(k), where wt(k) (A000120) gives the number of 1's in binary representation of nonnegative integer k.

Original entry on oeis.org

0, 2, 3, 7, 8, 9, 10, 11, 12, 16, 20, 24, 25, 26, 27, 28, 29, 34, 35, 40, 41, 42, 43, 44, 45, 49, 53, 57, 58, 59, 60, 61, 62, 65, 66, 68, 69, 72, 73, 74, 75, 76, 77, 81, 85, 89, 90, 91, 92, 93, 94, 99, 100, 105, 106, 107, 108, 109, 110, 114, 118, 122, 123, 124
Offset: 1

Views

Author

Antti Karttunen, Aug 09 2013

Keywords

Comments

wt(k) = A000120(k) is also called bitcount(k).
In other words, the positions of ones in A228085.
Numbers that can be expressed as the sum of distinct terms of the form 2^n+1, n=0,1,... in exactly one way. - Matthew C. Russell, Oct 08 2013

Examples

			0 is in this sequence because there is a unique k such that k+A000120(k)=0, in this case k=0.
1 is not in this sequence because there is no such k that k+A000120(k) would be 1. (Instead 1 is in A010061).
2 is in this sequence because there is exactly one k that satisfies k+A000120(k)=2, namely k=1.
3 is in this sequence because there is exactly one k that satisfies k+A000120(k)=3, namely k=2.
4 is not in this sequence because there is no such k that k+A000120(k) would be 4. (Instead 4 is in A010061.)
5 is not in this sequence because there is more than one k that satisfies k+A000120(k)=5, namely k=3 and k=4.
		

Crossrefs

Subset of A228082.
Cf. A228089 (corresponding k's for each a(n)).
Cf. A228090 (the same k's sorted into ascending order).
Cf. A227915.

Programs

Formula

a(n) = A092391(A228089(n)). [Consequence of the definitions of A228088 & A228089. Use the given Scheme-code to actually compute the sequence]
Showing 1-10 of 47 results. Next