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

A225901 Write n in factorial base, then replace each nonzero digit d of radix k with k-d.

Original entry on oeis.org

0, 1, 4, 5, 2, 3, 18, 19, 22, 23, 20, 21, 12, 13, 16, 17, 14, 15, 6, 7, 10, 11, 8, 9, 96, 97, 100, 101, 98, 99, 114, 115, 118, 119, 116, 117, 108, 109, 112, 113, 110, 111, 102, 103, 106, 107, 104, 105, 72, 73, 76, 77, 74, 75, 90, 91, 94, 95, 92, 93, 84, 85, 88, 89, 86, 87, 78, 79, 82, 83, 80, 81, 48, 49, 52, 53, 50, 51, 66, 67, 70, 71, 68
Offset: 0

Views

Author

Paul Tek, May 20 2013

Keywords

Comments

Analogous to A004488 or A048647 for the factorial base.
A self-inverse permutation of the natural numbers.
From Antti Karttunen, Aug 16-29 2016: (Start)
Consider the following way to view a factorial base representation of nonnegative integer n. For each nonzero digit d_i present in the factorial base representation of n (where i is the radix = 2.. = one more than 1-based position from the right), we place a pebble to the level (height) d_i at the corresponding column i of the triangular diagram like below, while for any zeros the corresponding columns are left empty:
.
Level
6 o
─ ─
5 . .
─ ─ ─
4 . . .
─ ─ ─ ─
3 . . . .
─ ─ ─ ─ ─
2 . . o . .
─ ─ ─ ─ ─ ─
1 . o . . o o
─ ─ ─ ─ ─ ─ ─
Radix: 7 6 5 4 3 2
Digits: 6 1 2 0 1 1 = A007623(4491)
Instead of levels, we can observe on which "slope" each pebble (nonzero digit) is located at. Formally, the slope of nonzero digit d_i with radix i is (i - d_i). Thus in above example, both the most significant digit (6) and the least significant 1 are on slope 1 (called "maximal slope", because it contains digits that are maximal allowed in those positions), while the second 1 from the right is on slope 2 ("submaximal slope").
This involution (A225901) sends each nonzero digit at level k to the slope k (and vice versa) by flipping such a diagram by the shallow diagonal axis that originates from the bottom right corner. Thus, from above diagram we obtain:
Slope (= digit's radix - digit's value)
1
2 .
3 . .╲
4 . .╲o╲
5 . .╲.╲.╲
6 . .╲.╲o╲.╲
. .╲.╲.╲.╲o╲
o╲.╲.╲.╲.╲o╲
-----------------
1 5 3 0 2 1 = A007623(1397)
and indeed, a(4491) = 1397 and a(1397) = 4491.
Thus this permutation maps between polynomial encodings A275734 & A275735 and all the respective sequences obtained from them, where the former set of sequences are concerned with the "slopes" and the latter set with the "levels" of the factorial base representation. See the Crossrefs section.
Sequences A231716 and A275956 are closed with respect to this sequence, in other words, for all n, a(A231716(n)) is a term of A231716 and a(A275956(n)) is a term of A275956.
(End)

Examples

			a(1000) = a(1*6! + 2*5! + 1*4! + 2*3! + 2*2!) = (7-1)*6! + (6-2)*5! + (5-1)*4! + (4-2)*3! + (3-2)*2! = 4910.
a(1397) = a(1*6! + 5*5! + 3*4! + 0*3! + 2*2! + 1*1!) = (7-1)*6! + (6-5)*5! + (5-3)*4! + (3-2)*2! + (2-1)*1! = 4491.

Links

Crossrefs

Cf. A000142, A007623, A004488, A048647, A001563, A007489, A257684, A257687, A276091, A275736, A276149.
Cf. A275959 (fixed points), A231716, A275956.
Cf. A153880 & A255411.
Cf. also A275734 & A275735, A275952 & A275954.
This involution maps between the following sequences related to "levels" and "slopes" (see comments): A275806 <--> A060502, A257511 <--> A260736, A264990 <--> A275811, A275729 <--> A275728, A275948 <--> A275946, A275949 <--> A275947, A275964 <--> A275962, A059590 <--> A276091.
Related permutations: A275957, A275958, A275835, A275836, A275837, A275838, A276011, A276012.

Programs

  • Mathematica
    b = MixedRadix[Reverse@ Range[2, 12]]; Table[FromDigits[Map[Boole[# > 0] &, #] (Reverse@ Range[2, Length@ # + 1] - #), b] &@ IntegerDigits[n, b], {n, 0, 82}] (* Version 10.2, or *)
f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Range[# + 1] <= n &]; Most@ Rest[a][[All, 1]] /. {} -> {0}]; g[w_List] := Total[Times @@@ Transpose@ {Map[Times @@ # &, Range@ Range[0, Length@ w]], Reverse@ Append[w, 0]}]; Table[g[Map[Boole[# > 0] &, #] (Reverse@ Range[2, Length@ # + 1] - #)] &@ f@ n, {n, 0, 82}] (* Michael De Vlieger, Aug 29 2016 *)
  • PARI
    a(n)=my(s=0,d,k=2);while(n,d=n%k;n=n\k;if(d,s=s+(k-d)*(k-1)!);k=k+1);return(s)
  • Python
    from sympy import factorial as f
def a(n):
    s=0
    k=2
    while(n):
        d=n%k
        n=(n//k)
        if d: s=s+(k - d)*f(k - 1)
        k+=1
    return s
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 19 2017
  • Scheme
    (define (A225901 n) (let loop ((n n) (z 0) (m 2) (f 1)) (cond ((zero? n) z) (else (loop (quotient n m) (if (zero? (modulo n m)) z (+ z (* f (- m (modulo n m))))) (+ 1 m) (* f m))))))
;; One implementing the first recurrence, with memoization-macro definec:
(definec (A225901 n) (if (zero? n) n (+ (A276091 (A275736 n)) (A153880 (A225901 (A257684 n))))))
;; Antti Karttunen, Aug 29 2016

Formula

From Antti Karttunen, Aug 29 2016: (Start)
a(0) = 0; for n >= 1, a(n) = A276091(A275736(n)) + A153880(a(A257684(n))).
or, for n >= 1, a(n) = A276149(n) + a(A257687(n)).
(End)
Other identities. For n >= 0:
a(n!) = A001563(n).
a(n!-1) = A007489(n-1).
From Antti Karttunen, Aug 16 2016: (Start)
A275734(a(n)) = A275735(n) and vice versa, A275735(a(n)) = A275734(n).
A060130(a(n)) = A060130(n). [The flip preserves the number of nonzero digits.]
A153880(n) = a(A255411(a(n))) and A255411(n) = a(A153880(a(n))). [This involution conjugates between the two fundamental factorial base shifts.]
a(n) = A257684(a(A153880(n))) = A266193(a(A255411(n))). [Follows from above.]
A276011(n) = A273662(a(A273670(n))).
A276012(n) = A273663(a(A256450(n))).
(End)

A275734 Prime-factorization representations of "factorial base slope polynomials": a(0) = 1; for n >= 1, a(n) = A275732(n) * a(A257684(n)).

Original entry on oeis.org

1, 2, 3, 6, 2, 4, 5, 10, 15, 30, 10, 20, 3, 6, 9, 18, 6, 12, 2, 4, 6, 12, 4, 8, 7, 14, 21, 42, 14, 28, 35, 70, 105, 210, 70, 140, 21, 42, 63, 126, 42, 84, 14, 28, 42, 84, 28, 56, 5, 10, 15, 30, 10, 20, 25, 50, 75, 150, 50, 100, 15, 30, 45, 90, 30, 60, 10, 20, 30, 60, 20, 40, 3, 6, 9, 18, 6, 12, 15, 30, 45, 90, 30, 60, 9, 18, 27
Offset: 0

Views

Author

Antti Karttunen, Aug 08 2016

Keywords

Comments

These are prime-factorization representations of single-variable polynomials where the coefficient of term x^(k-1) (encoded as the exponent of prime(k) in the factorization of n) is equal to the number of nonzero digits that occur on the slope (k-1) levels below the "maximal slope" in the factorial base representation of n. See A275811 for the definition of the "digit slopes" in this context.

Examples

			For n=23 ("321" in factorial base representation, A007623), all three nonzero digits are maximal for their positions (they all occur on "maximal slope"), thus a(23) = prime(1)^3 = 2^3 = 8.
For n=29 ("1021"), there are three nonzero digits, where both 2 and the rightmost 1 are on the "maximal slope", while the most significant 1 is on the "sub-sub-sub-maximal", thus a(29) = prime(1)^2 * prime(4)^1 = 2*7 = 28.
For n=37 ("1201"), there are three nonzero digits, where the rightmost 1 is on the maximal slope, 2 is on the sub-maximal, and the most significant 1 is on the "sub-sub-sub-maximal", thus a(37) = prime(1) * prime(2) * prime(4) = 2*3*7 = 42.
For n=55 ("2101"), the least significant 1 is on the maximal slope, and the digits "21" at the beginning are together on the sub-sub-maximal slope (as they are both two less than the maximal digit values 4 and 3 allowed in those positions), thus a(55) = prime(1)^1 * prime(3)^2 = 2*25 = 50.

Links

Crossrefs

Cf. A001221, A001222, A002110, A007489, A007814, A048675, A051903, A056169, A056170, A060130, A060502, A225901.
Cf. A257684, A275732.
Cf. A275811.
Cf. A260736, A275728, A275808, A275812, A275946, A275947, A275962.
Cf. A275804 (indices of squarefree terms), A275805 (of terms not squarefree).
Cf. also A275725, A275733, A275735, A276076 for other such prime factorization encodings of A060117/A060118-related polynomials.

Programs

  • Python
    from operator import mul
from sympy import prime, factorial as f
def a007623(n, p=2): return n if n
    0 else '0' for i in x)[::-1]
    return 0 if n==1 else sum(int(y[i])*f(i + 1) for i in range(len(y)))
def a(n): return 1 if n==0 else a275732(n)*a(a257684(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 19 2017

Formula

a(0) = 1; for n >= 1, a(n) = A275732(n) * a(A257684(n)).
Other identities and observations. For all n >= 0:
a(n) = A275735(A225901(n)).
a(A007489(n)) = A002110(n).
A001221(a(n)) = A060502(n).
A001222(a(n)) = A060130(n).
A007814(a(n)) = A260736(n).
A051903(a(n)) = A275811(n).
A048675(a(n)) = A275728(n).
A248663(a(n)) = A275808(n).
A056169(a(n)) = A275946(n).
A056170(a(n)) = A275947(n).
A275812(a(n)) = A275962(n).

A275727 a(0) = 0, for n >= 1, a(n) = A275736(n) + 2*a(A257684(n)).

Original entry on oeis.org

0, 1, 2, 3, 2, 3, 4, 5, 6, 7, 6, 7, 4, 5, 6, 7, 6, 7, 4, 5, 6, 7, 6, 7, 8, 9, 10, 11, 10, 11, 12, 13, 14, 15, 14, 15, 12, 13, 14, 15, 14, 15, 12, 13, 14, 15, 14, 15, 8, 9, 10, 11, 10, 11, 12, 13, 14, 15, 14, 15, 12, 13, 14, 15, 14, 15, 12, 13, 14, 15, 14, 15, 8, 9, 10, 11, 10, 11, 12, 13, 14, 15, 14, 15, 12, 13, 14, 15
Offset: 0

Views

Author

Antti Karttunen, Aug 09 2016

Keywords

Comments

a(n) has ones in those positions of its base-2 representation where n has nonzero digits in its factorial base representation.

Examples

			For n=19, A007623(19) = 301, thus a(19) = 5 because A007088(5) = 101.

Links

Crossrefs

Cf. A000120, A007088, A007623, A048675, A060130, A257684, A275733, A275736.
Cf. also A275728.

Formula

a(0) = 0, for n >= 1, a(n) = A275736(n) + 2*a(A257684(n)).
a(n) = A048675(A275733(n)).
Other identities and observations. For all n >= 0:
A000120(a(n)) = A060130(n).

A275808 a(0) = 0, for n >= 1, a(n) = A275736(n) XOR a(A257684(n)), where XOR is given by A003987.

Original entry on oeis.org

0, 1, 2, 3, 1, 0, 4, 5, 6, 7, 5, 4, 2, 3, 0, 1, 3, 2, 1, 0, 3, 2, 0, 1, 8, 9, 10, 11, 9, 8, 12, 13, 14, 15, 13, 12, 10, 11, 8, 9, 11, 10, 9, 8, 11, 10, 8, 9, 4, 5, 6, 7, 5, 4, 0, 1, 2, 3, 1, 0, 6, 7, 4, 5, 7, 6, 5, 4, 7, 6, 4, 5, 2, 3, 0, 1, 3, 2, 6, 7, 4, 5, 7, 6, 0, 1, 2, 3, 1, 0, 3, 2, 1, 0, 2, 3, 1, 0, 3, 2, 0
Offset: 0

Views

Author

Antti Karttunen, Aug 09 2016

Keywords

Links

Crossrefs

Cf. A003987, A248663, A257684, A275734, A275736.
Cf. A275809 (positions of zeros), A275810 (and their first differences).
Cf. also A275728.

Formula

a(0) = 0, for n >= 1, a(n) = A275736(n) XOR a(A257684(n)), where 2-argument function XOR is given by A003987.
a(n) = A248663(A275734(n)).

A276010 a(0) = 0, for n >= 1, a(n) = A275736(n) OR a(A257684(n)), where OR is given by A003986.

Original entry on oeis.org

0, 1, 2, 3, 1, 1, 4, 5, 6, 7, 5, 5, 2, 3, 2, 3, 3, 3, 1, 1, 3, 3, 1, 1, 8, 9, 10, 11, 9, 9, 12, 13, 14, 15, 13, 13, 10, 11, 10, 11, 11, 11, 9, 9, 11, 11, 9, 9, 4, 5, 6, 7, 5, 5, 4, 5, 6, 7, 5, 5, 6, 7, 6, 7, 7, 7, 5, 5, 7, 7, 5, 5, 2, 3, 2, 3, 3, 3, 6, 7, 6, 7, 7, 7, 2, 3, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 1, 1, 3, 3, 1, 1, 5, 5, 7, 7, 5, 5, 3
Offset: 0

Views

Author

Antti Karttunen, Aug 17 2016

Keywords

Links

Crossrefs

Cf. A003986, A087207, A257684, A275734.
Cf. also A000120, A060502, A275736, A275728, A275808.

Formula

a(0) = 0, for n >= 1, a(n) = A275736(n) OR a(A257684(n)), where OR is given by A003986.
a(n) = A087207(A275734(n)).
Other identities. For all n >= 1:
A000120(a(n)) = A060502(n).
Showing 1-5 of 5 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.