A257684 -id:A257684 - cp's OEIS frontend

A246359 Maximum digit in the factorial base expansion of n (A007623).

0, 1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4

Offset: 0

Antti Karttunen, Oct 20 2014

Maximum entry in n-th row of A108731.

			Factorial base representation of 46 is "1320" as 46 = 1*4! + 3*3! + 2*2! + 0*1!, and the largest of these digits is 3, thus a(46) = 3.

Antti Karttunen, Table of n, a(n) for n = 0..10080
Index entries for sequences related to factorial base representation

Cf. A007623, A034968, A051903, A060130, A084558, A099563, A257684, A257687, A276076.

Cf. also A249070.

nn = 96; m = 1; While[Factorial@ m < nn, m++]; m; Table[Max@ IntegerDigits[n, MixedRadix[Reverse@ Range[2, m]]], {n, 0, nn}] (* 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}]; Table[Max@ f@ n, {n, 0, 96}] (* Michael De Vlieger, Aug 29 2016 *)

def a007623(n, p=2): return n if nIndranil Ghosh, Jun 21 2017

From Antti Karttunen, Aug 29 2016: (Start)
a(0) = 0; for n >= 1, a(n) = 1 + a(A257684(n)).
a(0) = 0; for n >= 1, a(n) = max(A099563(n), a(A257687(n))).
a(n) = A051903(A276076(n)).
(End)

A275804 Numbers with at most one nonzero digit on each digit slope of the factorial base representation of n.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 9, 10, 12, 13, 16, 18, 20, 24, 25, 26, 27, 28, 30, 31, 32, 33, 34, 36, 37, 40, 42, 44, 48, 49, 50, 51, 52, 60, 61, 64, 66, 68, 72, 73, 76, 78, 79, 82, 90, 96, 98, 102, 104, 108, 120, 121, 122, 123, 124, 126, 127, 128, 129, 130, 132, 133, 136, 138, 140, 144, 145, 146, 147, 148, 150, 151, 152, 153, 154, 156, 157, 160

Offset: 0

Antti Karttunen, Aug 10 2016

Indexing starts from zero, because a(0) = 0 is a special case in this sequence.
Numbers n for which A275947(n) = 0 or equally, for which A275811(n) <= 1.
Numbers n for which A008683(A275734(n)) <> 0, that is, indices of squarefree terms in A275734.
Numbers n for which A060130(n) = A060502(n).
Numbers with at most one nonzero digit on each digit slope of the factorial base representation of n (see A275811 and A060502 for the definition of slopes in this context). More exactly: numbers n in whose factorial base representation (A007623) there does not exist a pair of digit positions i_1 and i_2 with nonzero digits d_1 and d_2, such that (i_1 - d_1) = (i_2 - d_2).

Antti Karttunen, Table of n, a(n) for n = 0..10000
Index entries for sequences related to factorial base representation

Complement: A275805.
Indices of zeros in A275947 and A275962.
Intersection with A276005 gives A261220.
Cf. A059590 (a subsequence).
Cf. A007623, A008683, A060130, A060502, A275734, A275811.

from operator import mul
from sympy import prime, factorial as f
from sympy.ntheory.factor_ import core
def a007623(n, p=2): return n if n0 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))
def ok(n): return 1 if n==0 else core(a(n))==a(n)
print([n for n in range(201) if ok(n)]) # Indranil Ghosh, Jun 19 2017

A264990 a(n) = number of occurrences of a most frequent nonzero digit in factorial base representation (A007623) of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 2, 2, 3, 1, 2, 1, 1, 1, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 3, 1, 2, 2, 3, 3, 4, 2, 3, 1, 2, 2, 3, 2, 2, 1, 2, 2, 3, 1, 2, 1, 1, 1, 2, 2, 2, 1, 2, 2, 3, 2, 2, 2, 2, 2, 2, 3, 3, 1, 1, 1, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 3, 1, 2, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 3, 1, 2, 1, 1, 1, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2

Offset: 0

Antti Karttunen, Dec 22 2015

			   n  A007623(n)   a(n) [highest number of times any nonzero digit occurs].
   0 =   0           0 (because no nonzero digits present)
   1 =   1           1
   2 =  10           1
   3 =  11           2
   4 =  20           1
   5 =  21           1
   6 = 100           1
   7 = 101           2
   8 = 110           2
   9 = 111           3
  10 = 120           1
  11 = 121           2
  12 = 200           1
  13 = 201           1
  14 = 210           1
  15 = 211           2
  16 = 220           2
  17 = 221           2
  18 = 300           1
and for n=63 we have:
  63 = 2211          2.

Antti Karttunen, Table of n, a(n) for n = 0..10080
Index entries for sequences related to factorial base representation.

Cf. A007623, A051903, A225901, A257511, A257684, A275735, A275811.
Cf. A265349 (positions of terms <= 1), A265350 (positions of term > 1).
Cf. also A266117, A266118.

a[n_] := Module[{k = n, m = 2, r, s = {}}, While[{k, r} = QuotientRemainder[k, m]; k != 0|| r != 0, AppendTo[s, r]; m++]; Max[Tally[Select[s, # > 0 &]][[;;,2]]]]; a[0] = 0; Array[a, 100, 0] (* Amiram Eldar, Jan 24 2024 *)

from sympy import prime, factorint
from operator import mul
import collections
def a007623(n, p=2): return n if nIndranil Ghosh, Jun 20 2017

a(0) = 0; for n >= 1, a(n) = max(A257511(n), a(A257684(n))).
Other identities. For all n >= 0:
From Antti Karttunen, Aug 15 2016: (Start)
a(n) = A275811(A225901(n)).
a(n) = A051903(A275735(n)).
(End)

Name changed by Antti Karttunen, Aug 15 2016

A273668 Permutation of nonnegative integers: a(0) = 0, a(A255411(n)) = A153880(a(n)), a(A256450(n)) = A273670(a(n)).

Original entry on oeis.org

0, 1, 3, 5, 2, 9, 4, 15, 7, 21, 11, 29, 8, 17, 41, 13, 14, 23, 6, 57, 19, 20, 32, 33, 10, 77, 27, 28, 44, 45, 16, 101, 39, 40, 61, 63, 22, 129, 53, 55, 83, 87, 31, 165, 71, 75, 107, 111, 12, 43, 213, 95, 56, 99, 137, 141, 18, 59, 269, 119, 26, 76, 125, 177, 80, 183, 38, 25, 81, 341, 134, 153, 30, 37, 100, 161, 62, 225, 104, 231, 52, 35

Offset: 0

Antti Karttunen, May 30 2016

Inverse: A273667.
Cf. A153880, A255411, A256450, A257680, A257684, A273662, A273670.
Similar or related permutations: A255565, A273666.

a(0) = 0; for n >= 1: if A257680(n) = 0 [when n is one of the terms of A255411] then a(n) = A153880(a(A257684(n))), otherwise [when n is one of the terms of A256450], a(n) = A273670(a(A273662(n))).
As a composition of other permutations:
a(n) = A273666(A255565(n)).

A275811 Number of nonzero digits on a maximally occupied slope of factorial base representation of n: a(n) = A051903(A275734(n)). See comments for the definition.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 2, 2, 3, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 2, 2, 3, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 2, 2, 3, 1, 1, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 2, 2, 3, 3, 2, 2, 1, 2, 2, 2, 2, 3, 1, 2, 1, 2, 2, 3, 1, 2, 1, 2, 2, 3, 1, 2, 2, 2, 2, 3, 2, 3, 2, 3, 3, 4, 1

Offset: 0

Antti Karttunen, Aug 10 2016

Digit slopes are called "maximal", "sub-maximal", "sub-sub-maximal", etc. For digit-positions we employ one-based indexing, thus we say that the least significant digit of factorial base expansion of n is in position 1, etc. The maximal digit slope is occupied when there is at least one digit-position k that contains digit k (maximal digit allowed in that position), so that A260736(n) > 0, and n is thus a term of A273670. The sub-maximal digit slope is occupied when there is at least one nonzero digit k in a digit-position k+1. The sub-sub-maximal slope is occupied when there is at least one nonzero digit k in a digit-position k+2, etc. This sequence gives the number of nonzero digits on a slope (of possibly several) for which there exists no other slopes with more nonzero digits. See the examples.
In other words: a(n) gives the number of occurrences of a most common element in the multiset [(i_x - d_x) | where d_x ranges over each nonzero digit present in factorial base representation of n and i_x is that digit's position from the right].
Involution A225901 maps this metric to another metric A264990 which gives the maximal number of equal nonzero digits occurring in factorial base representation (A007623) of n. See also A060502.

			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 the "maximal slope" is also the "maximally occupied slope" (as there are no other occupied slopes present), and a(23) = 3.
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 here the "maximal slope" is also the "maximally occupied slope" (with 2 nonzero digits present), and a(29) = 2.
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 there are three occupied slopes in total, all with just one nonzero digit present, and a(37) = 1.
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 here the sub-sub-maximal slope is the "maximally occupied slope" with its two occupiers, and a(55) = 2.

Antti Karttunen, Table of n, a(n) for n = 0..40320
Index entries for sequences related to factorial base representation

Cf. A051903, A060502, A225901, A264990, A275734.

Cf. A275804 (gives the indices of 0 and 1's), A275805 (gives the indices of terms > 1).

from sympy import prime, factorint
from operator import mul
from functools import reduce
from sympy import factorial as f
def a051903(n): return 0 if n==1 else max(factorint(n).values())
def a007623(n, p=2): return n if n0 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 a275734(n): return 1 if n==0 else a275732(n)*a275734(a257684(n))
def a(n): return 0 if n==0 else a051903(a275734(n))
print([a(n) for n in range(201)]) # Indranil Ghosh, Jun 20 2017

a(n) = A051903(A275734(n)).

a(n) = A264990(A225901(n)).

Signs in comment corrected and clarification added by Antti Karttunen, Aug 16 2016

A275805 Indices of nonsquarefree terms in A275734; numbers with at least one digit slope (in their factorial base representation) with multiple nonzero digits. (See comments for the exact definition).

Original entry on oeis.org

5, 11, 14, 15, 17, 19, 21, 22, 23, 29, 35, 38, 39, 41, 43, 45, 46, 47, 53, 54, 55, 56, 57, 58, 59, 62, 63, 65, 67, 69, 70, 71, 74, 75, 77, 80, 81, 83, 84, 85, 86, 87, 88, 89, 91, 92, 93, 94, 95, 97, 99, 100, 101, 103, 105, 106, 107, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 125, 131, 134, 135, 137, 139, 141, 142, 143, 149, 155

Offset: 1

Antti Karttunen, Aug 10 2016

Numbers n for which A008683(A275734(n)) = 0.
Numbers n for which A275811(n) > 1.
Numbers n in whose factorial base representation (A007623) there exists at least one pair of digit positions i_1 and i_2 with nonzero digits d_1 and d_2 such that (i_1 - d_1) = (i_2 - d_2).

			For n=5, "21" in factorial base (A007623), the pair 2 (in position 2) and 1 (in position 1) satisfies the condition, as (2-2) = (1-1), thus 5 is included.
For n=55, "2101" in factorial base, the pair 2 (in position 4) and 1 (in position 3) satisfies the condition, as (4-2) = (3-1), thus 55 is included.
For n=67, "2301" in factorial base, the pair 3 (in position 3) and 1 (in position 1) satisfies the condition, as (3-3) = (1-1), thus 67 is included in the sequence.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences related to factorial base representation

Complement: A275804.
Cf. A007623, A008683, A275734, A275811.
Cf. A275809 (a subsequence apart from its initial 0-term).
Subsequence of A115945.

from operator import mul
from sympy import prime, factorial as f, mobius
from functools import reduce
def a007623(n, p=2): return n if n0 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([n for n in range(201) if mobius(a(n))==0]) # Indranil Ghosh, Jun 19 2017

A275848 Permutation of natural numbers: a(0) = 0, a(A255411(n)) = A153880(a(n)), a(A256450(n)) = A273670(n).

Original entry on oeis.org

0, 1, 3, 4, 2, 5, 7, 9, 10, 11, 13, 15, 8, 16, 17, 18, 12, 19, 6, 20, 21, 22, 14, 23, 25, 27, 28, 29, 31, 33, 34, 35, 37, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 51, 52, 53, 55, 57, 26, 58, 59, 61, 32, 63, 64, 65, 66, 67, 68, 69, 36, 70, 71, 73, 38, 75, 50, 76, 77, 79, 56, 81, 30, 82, 83, 85, 60, 87, 88, 89, 90, 91, 92, 93, 62, 94, 95, 96, 72, 97, 48

Offset: 0

Antti Karttunen, Aug 13 2016

Inverse: A275847.
Cf. A153880, A255411, A256450, A257680, A257684, A273662, A273670.
Similar permutations: A273668 (a more recursed variant), A275845, A275846.

a(0) = 0; for n >= 1: if A257680(n) = 0 [when n is one of the terms of A255411] then a(n) = A153880(a(A257684(n))), otherwise [when n is one of the terms of A256450], a(n) = A273670(A273662(n)).

A257685 Left inverse for injection A255411: a(0) = 0, after which, if n = A255411(k) for some k, then a(n) = k, otherwise a(n) = 0.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 3, 0, 4, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0, 8, 0, 0, 0, 9, 0, 10, 0, 0, 0, 11, 0, 12, 0, 0, 0, 13, 0, 0, 0, 0, 0, 0, 0, 14, 0, 0, 0, 15, 0, 16, 0, 0, 0, 17, 0, 18, 0, 0, 0, 19, 0, 0, 0, 0, 0, 0, 0, 20, 0, 0, 0, 21, 0, 22, 0, 0, 0, 23, 0, 0

Offset: 0

Antti Karttunen, May 04 2015

nonn

Antti Karttunen, Table of n, a(n) for n = 0..10080

Cf. A255411, A257680, A257681, A257684.

Position[Select[Range[0, 120], ! MemberQ[IntegerDigits[#, MixedRadix[Reverse@ Range@ 12]], 1] &], #] - 1 & /@ Range[0, 120] /. {} -> 0 // Flatten (* Michael De Vlieger, May 30 2016, Version 10.2 *)

from sympy import factorial as f
def a007623(n, p=2):
    return n if n0 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 a257680(n): return 1 if '1' in str(a007623(n)) else 0
def a(n): return (1 - a257680(n))*a257684(n)
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 21 2017

(define (A257685 n) (* (- 1 (A257680 n)) (A257684 n)))

a(0) = 0, after which, if n = A255411(k) for some k, then a(n) = k, otherwise a(n) = 0.
a(n) = (1-A257680(n)) * A257684(n).
Other identities:
For all n >= 0, a(A255411(n)) = n. [This sequence works as a left inverse of A255411.]

A275835 Permutation of nonnegative integers: a(n) = A273667(A225901(n)).

Original entry on oeis.org

0, 1, 6, 3, 4, 2, 56, 20, 36, 17, 21, 9, 48, 15, 30, 13, 16, 7, 18, 8, 24, 10, 12, 5, 495, 135, 74, 31, 132, 53, 582, 401, 147, 59, 157, 158, 361, 116, 216, 173, 117, 47, 136, 155, 380, 46, 78, 82, 420, 111, 61, 25, 108, 45, 490, 347, 123, 51, 133, 134, 312, 93, 192, 149, 94, 41, 112, 131, 327, 40, 64, 68, 360, 270, 80, 38, 88, 89, 416, 303, 99, 44, 109, 110, 288, 34

Offset: 0

Antti Karttunen, Aug 13 2016

Inverse: A275836.
Cf. A225901, A255411, A257684, A273667.
Cf. also A275837, A275838.

(define (A275835 n) (A273667 (A225901 n)))

a(n) = A273667(A225901(n)).
Other identities. For all n >= 0:
a(n) = A257684(a(A255411(n))). [Restriction to A255411 induces the same permutation.]

A275836 Permutation of nonnegative integers: a(n) = A225901(A273668(n)).

Original entry on oeis.org

0, 1, 5, 3, 4, 23, 2, 17, 19, 11, 21, 99, 22, 15, 111, 13, 16, 9, 18, 95, 7, 10, 118, 119, 20, 51, 101, 98, 106, 107, 14, 27, 113, 110, 85, 89, 8, 623, 75, 91, 69, 65, 115, 707, 81, 53, 45, 41, 12, 103, 659, 57, 94, 29, 615, 611, 6, 93, 579, 33, 100, 50, 603, 695, 70, 689, 112, 97, 71, 507, 616, 719, 114, 109, 26, 711, 88, 647, 46, 641, 74, 117, 47, 441, 688

Offset: 0

Antti Karttunen, Aug 13 2016

Inverse: A275835.
Cf. A225901,  A255411, A257684, A273668.
Cf. also A275837, A275838.

(define (A275836 n) (A225901 (A273668 n)))

a(n) = A225901(A273668(n)).
Other identities. For all n >= 0:
a(n) = A257684(a(A255411(n))). [Restriction to A255411 induces the same permutation.]

cp's OEIS Frontend

A246359 Maximum digit in the factorial base expansion of n (A007623).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

A275804 Numbers with at most one nonzero digit on each digit slope of the factorial base representation of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Python

A264990 a(n) = number of occurrences of a most frequent nonzero digit in factorial base representation (A007623) of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Python

Formula

Extensions

A273668 Permutation of nonnegative integers: a(0) = 0, a(A255411(n)) = A153880(a(n)), a(A256450(n)) = A273670(a(n)).

Views

Author

Keywords

Links

Crossrefs

Formula

A275811 Number of nonzero digits on a maximally occupied slope of factorial base representation of n: a(n) = A051903(A275734(n)). See comments for the definition.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Python

Formula

Extensions

A275805 Indices of nonsquarefree terms in A275734; numbers with at least one digit slope (in their factorial base representation) with multiple nonzero digits. (See comments for the exact definition).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Python

A275848 Permutation of natural numbers: a(0) = 0, a(A255411(n)) = A153880(a(n)), a(A256450(n)) = A273670(n).

Views

Author

Keywords

Links

Crossrefs

Formula

A257685 Left inverse for injection A255411: a(0) = 0, after which, if n = A255411(k) for some k, then a(n) = k, otherwise a(n) = 0.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Python

Scheme