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-3 of 3 results.

A255411 Shift factorial base representation of n one digit left (with 0 added to right), increment all nonzero digits by one, then convert back to decimal; Numbers with no digit 1 in their factorial base representation.

Original entry on oeis.org

0, 4, 12, 16, 18, 22, 48, 52, 60, 64, 66, 70, 72, 76, 84, 88, 90, 94, 96, 100, 108, 112, 114, 118, 240, 244, 252, 256, 258, 262, 288, 292, 300, 304, 306, 310, 312, 316, 324, 328, 330, 334, 336, 340, 348, 352, 354, 358, 360, 364, 372, 376, 378, 382, 408, 412, 420, 424, 426, 430, 432, 436, 444
Offset: 0

Views

Author

Antti Karttunen, Apr 16 2015

Keywords

Comments

Nonnegative integers such that the number of ones (A257511) in their factorial base representation (A007623) is zero.
Nonnegative integers such that the least missing nonzero digit (A257079) in their factorial base representation is one.
a(n) can be also directly computed from n by "shifting left" its factorial base representation (that is, by appending one zero to the right, see A153880) and then incrementing all nonzero digits by one, and then converting the resulting (still valid) factorial base number back to decimal. See the examples.
The sequences A227130 and A227132 are closed under a(n), in other words, permutation listed as the a(n)-th entry in tables A060117 & A060118 has the same parity as the n-th entry in those same tables.

Examples

			Factorial base representation (A007623) of 1 is "1", shifting it left yields "10", and when we increment all nonzero digits by one, we get "20", which is the factorial base representation of 4 (as 4 = 2*2! + 0*1!), thus a(1) = 4.
F.b.r. of 2 is "10", shifting it left yields "100", and "200" is f.b.r. of 12, thus a(2) = 12.
F.b.r. of 43 is "1301", shifting it left and incrementing all nonzeros by one yields "24020", which is f.b.r of 340, thus a(43) = 340.

Links

Crossrefs

Complement: A256450.
Positions of ones in A257079, fixed points of A257080, positions of zeros in A257511, A257081 and A257261.
Cf. A007623, A227157, A153880, A231720.
Cf. A255341, A255342, A255343.
Cf. A257262, A257263.
Cf. also A227130/A227132, A060117/A060118 and also arrays A257503 & A257505.

Programs

  • Mathematica
    factBaseIntDs[n_] := Module[{m, i, len, dList, currDigit}, i = 1; While[n > i!, i++]; m = n; len = i; dList = Table[0, {len}]; Do[currDigit = 0; While[m >= j!, m = m - j!; currDigit++]; dList[[len - j + 1]] = currDigit, {j, i, 1, -1}]; If[dList[[1]] == 0, dList = Drop[dList, 1]]; dList]; s = Table[FromDigits[factBaseIntDs[n]], {n, 500}]; {0}~Join~Flatten@ Position[s, x_ /; DigitCount[x][[1]] == 0](* Michael De Vlieger, Apr 27 2015, after Alonso del Arte at A007623 *)
Select[Range[0, 444], ! MemberQ[IntegerDigits[#, MixedRadix[Reverse@ Range@ 12]], 1] &] (* Michael De Vlieger, May 30 2016, Version 10.2 *)
r = MixedRadix[Reverse@Range[2, 12]]; Table[FromDigits[Map[If[# == 0, 0, # + 1] &, IntegerDigits[n, r]]~Join~{0}, r], {n, 0, 60}] (* Michael De Vlieger, Aug 14 2016, Version 10.2 *)
  • Python
    from sympy import factorial as f
def a007623(n, p=2): return n if n
    0 else '0' for i in x)[::-1]
    return 0 if n==0 else sum(int(y[i])*f(i + 1) for i in range(len(y)))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 20 2017

A227157 Numbers k whose factorial base representation A007623(k) does not contain any nonleading zeros.

Original entry on oeis.org

1, 3, 5, 9, 11, 15, 17, 21, 23, 33, 35, 39, 41, 45, 47, 57, 59, 63, 65, 69, 71, 81, 83, 87, 89, 93, 95, 105, 107, 111, 113, 117, 119, 153, 155, 159, 161, 165, 167, 177, 179, 183, 185, 189, 191, 201, 203, 207, 209, 213, 215, 225, 227, 231, 233, 237, 239, 273
Offset: 1

Views

Author

Antti Karttunen, Jul 04 2013

Keywords

Comments

a(A003422(n)) = A007489(n).
a(A007489(n)) = (n+1)!-1 thus A007489(n) gives the number of terms less than (n+1)! in this sequence.
Equivalently, there are n! terms in the sequence with their magnitude in range n!..(n+1)!.
Also numbers k such that A304036(k) = 1 for k > 0. - Seiichi Manyama, May 06 2018

Links

Crossrefs

The sequence gives all n for which A208575(n) is not zero. Complement of A227187. Subsets: A071156 (apart from zero), A231716, A231720.
Cf. also A003422, A007489, A007623, A153880, A227130, A227132, A304036.

Programs

  • Mathematica
    q[n_] := Module[{k = n, m = 2, c = 0, r}, While[{k, r} = QuotientRemainder[k, m]; k != 0 || r != 0, If[r == 0, c++]; m++]; c == 0]; Select[Range[300], q] (* Amiram Eldar, Jan 23 2024 *)

A220655 For n with a unique factorial base representation n = du*u! + ... + d2*2! + d1*1! (each di in range 0..i, cf. A007623), a(n) = (du+1)*u! + ... + (d2+1)*2! + (d1+1)*1!; a(n) = n + A007489(A084558(n)).

Original entry on oeis.org

2, 5, 6, 7, 8, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99
Offset: 1

Views

Author

Antti Karttunen, Dec 17 2012

Keywords

Comments

Used for computing A107346.
Term a(n) can be obtained by adding one to each digit of factorial base representation of n (A007623(n)) and then reinterpreting it as a kind of pseudo-factorial base representation, ignoring the fact that now some of the digits might be over the maximum allowed in that position. Please see the example section. - Antti Karttunen, Nov 29 2013

Examples

			1 has a factorial base representation A007623(1) = '1', as 1 = 1*1!. Incrementing the digit 1 with 1, we get 2*1! = 2, thus a(1) = 2. (Note that although '2' is not a valid factorial base representation, it doesn't matter here.)
2 has a factorial base representation '10', as 2 = 1*2! + 0*1!. Incrementing the digits by one, we get 2*2! + 1*1! = 5, thus a(2) = 5.
3 has a factorial base representation '11', as 3 = 1*2! + 1*1!. Incrementing the digits by one, we get 2*2! + 2*1! = 6, thus a(3) = 6.

Links

Crossrefs

Complement: A220695.
One less than A220656.
Cf. A000142, A003422, A007489, A007623, A084558, A212598, A153880, A231720.

Programs

  • Mathematica
    Block[{nn = 66, m = 1}, While[Factorial@ m < nn, m++]; m = MixedRadix[Reverse@ Range[2, m]]; Array[FromDigits[1 + IntegerDigits[#, m], m] &, nn]] (* Michael De Vlieger, Jan 20 2020 *)
  • Scheme
    ;; Standalone iterative implementation (Nov 29 2013):
(define (A220655 n) (let loop ((n n) (z 0) (i 2) (f 1)) (cond ((zero? n) z) (else (loop (quotient n i) (+ (* f (+ 1 (remainder n i))) z) (+ 1 i) (* f i))))))
;; Alternative implementation:
(define (A220655 n) (+ n (A007489 (A084558 n))))

Formula

a(n) = A220656(n)-1 = A003422(A084558(n)+1) + A000142(A084558(n)) + A212598(n) - 1. [The original definition]
a(n) = n + A007489(A084558(n)). [The above formula reduces to this, which proves that the original Dec 17 2012 description and the new main description produce the same sequence. Essentially, we are adding to n a factorial base repunit '1...111' with as many fact.base digits as n has.] - Antti Karttunen, Nov 29 2013
For n >= 1, A231720(n) = a(A153880(n)).

Extensions

Name changed by Antti Karttunen, Nov 29 2013
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.