A276083 -id:A276083 - cp's OEIS frontend

A059590 Numbers obtained by reinterpreting base-2 representation of n in the factorial base: a(n) = Sum_{k>=0} A030308(n,k)*A000142(k+1).

0, 1, 2, 3, 6, 7, 8, 9, 24, 25, 26, 27, 30, 31, 32, 33, 120, 121, 122, 123, 126, 127, 128, 129, 144, 145, 146, 147, 150, 151, 152, 153, 720, 721, 722, 723, 726, 727, 728, 729, 744, 745, 746, 747, 750, 751, 752, 753, 840, 841, 842, 843, 846, 847, 848, 849, 864, 865

Offset: 0

Henry Bottomley, Jan 24 2001

Numbers that are sums of distinct factorials (0! and 1! not treated as distinct).
Complement of A115945; A115944(a(n)) > 0; A115647 is a subsequence. - Reinhard Zumkeller, Feb 02 2006
A115944(a(n)) = 1. - Reinhard Zumkeller, Dec 04 2011
From Tilman Piesk, Jun 04 2012: (Start)
The inversion vector (compare A007623) of finite permutation a(n) (compare A055089, A195663) has only zeros and ones. Interpreted as a binary number it is 2*n (or n when the inversion vector is defined without the leading 0).
The inversion set of finite permutation a(n) interpreted as a binary number (compare A211362) is A211364(n).
(End)

			128 is in the sequence since 5! + 3! + 2! = 128.
a(22) = 128. a(22) = a(6) + (1 + floor(log(16) / log(2)))! = 8 + 5! = 128. Also, 22 = 10110_2. Therefore, a(22) = 1 * 5! + 0 * 4! + 1 * 3! + 1 + 2! + 0 * 0! = 128. - _David A. Corneth_, Aug 21 2016

Reinhard Zumkeller (terms 0..500) & Antti Karttunen, Table of n, a(n) for n = 0..8191
Index entries for sequences related to factorial base representation
Index entries for sequences related to factorial numbers

Cf. A007088, A007623, A014597, A051760, A051761, A059589.
Indices of zeros in A257684.
Cf. A275736 (left inverse).
Cf. A025494, A060112 (subsequences).
Cf. A153880, A225901.
Subsequence of A060132, A256450 and A275804.
Other sequences that are built by replacing 2^k in the binary representation with other numbers: A029931 (naturals), A089625 (primes), A022290 (Fibonacci), A197433 (Catalans), A276091 (n*n!), A275959 ((2n)!/2). Cf. also A276082 & A276083.

import Data.List (elemIndices)
a059590 n = a059590_list !! n
a059590_list = elemIndices 1 $ map a115944 [0..]
-- Reinhard Zumkeller, Dec 04 2011

[seq(bin2facbase(j),j=0..64)]; bin2facbase := proc(n) local i; add((floor(n/(2^i)) mod 2)*((i+1)!),i=0..floor_log_2(n)); end;
floor_log_2 := proc(n) local nn,i; nn := n; for i from -1 to n do if(0 = nn) then RETURN(i); fi; nn := floor(nn/2); od; end;
# next Maple program:
a:= n-> (l-> add(l[j]*j!, j=1..nops(l)))(Bits[Split](n)):
seq(a(n), n=0..57);  # Alois P. Heinz, Aug 12 2025

a[n_] :=  Reverse[id = IntegerDigits[n, 2]].Range[Length[id]]!; Table[a[n], {n, 0, 60}] (* Jean-François Alcover, Jun 19 2012, after Philippe Deléham *)

a(n) = if(n>0, a(n-msb(n)) + (1+logint(n,2))!, 0)
msb(n) = 2^#binary(n)>>1
{my(b = binary(n)); sum(i=1,#b,b[i]*(#b+1-i)!)} \\ David A. Corneth, Aug 21 2016

def facbase(k, f):
    return sum(f[i] for i, bi in enumerate(bin(k)[2:][::-1]) if bi == "1")
def auptoN(N): # terms up to N factorial-base digits; 13 generates b-file
    f = [factorial(i) for i in range(1, N+1)]
    return list(facbase(k, f) for k in range(2**N))
print(auptoN(5)) # Michael S. Branicky, Oct 15 2022

G.f. 1/(1-x) * Sum_{k>=0} (k+1)!*x^2^k/(1+x^2^k). - Ralf Stephan, Jun 24 2003
a(n) = Sum_{k>=0} A030308(n,k)*A000142(k+1). - Philippe Deléham, Oct 15 2011
From Antti Karttunen, Aug 19 2016: (Start)
a(0) = 0, a(2n) = A153880(a(n)), a(2n+1) = 1+A153880(a(n)).
a(n) = A225901(A276091(n)).
a(n) = A276075(A019565(n)).
a(A275727(n)) = A276008(n).
A275736(a(n)) = n.
A276076(a(n)) = A019565(n).
A007623(a(n)) = A007088(n).
(End)
a(n) = a(n - mbs(n)) + (1 + floor(log(n) / log(2)))!. - David A. Corneth, Aug 21 2016

Name changed (to emphasize the functional nature of the sequence) with the old definition moved to the comments by Antti Karttunen, Aug 21 2016

A276091 Numbers obtained by reinterpreting base-2 representation of n in A001563-base (A276326): a(n) = Sum_{k>=0} A030308(n,k)*A001563(k+1).

Original entry on oeis.org

0, 1, 4, 5, 18, 19, 22, 23, 96, 97, 100, 101, 114, 115, 118, 119, 600, 601, 604, 605, 618, 619, 622, 623, 696, 697, 700, 701, 714, 715, 718, 719, 4320, 4321, 4324, 4325, 4338, 4339, 4342, 4343, 4416, 4417, 4420, 4421, 4434, 4435, 4438, 4439, 4920, 4921, 4924, 4925, 4938, 4939, 4942, 4943, 5016, 5017, 5020, 5021, 5034, 5035, 5038, 5039, 35280, 35281

Offset: 0

Antti Karttunen, Aug 19 2016

Numbers that are sums of distinct terms of A001563.
A number is included if and only if all the nonzero digits in its factorial base representation (A007623) are maximal allowed in those digit positions, thus this sequence gives all numbers n for which A060130(n) = A260736(n).
Numbers n for which A276328(n) = A276337(n), thus from 1 onward the positions of ones in A276336.
Conjectured also to give all numbers n for which A255411(n) = A276340(n) (thus zeros of A276339).

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

Cf. A000120, A000142, A001563, A030308, A059590, A060130, A260736, A225901, A255411, A275959, A276082, A276083, A276090, A276326, A276328, A276336, A276337, A276339, A276340.

Table[Total[Times @@@ Transpose@ {Map[# #! &, Range@ Length@ #], Reverse@ #}] &@ IntegerDigits[n, 2], {n, 64}] (* Michael De Vlieger, Aug 31 2016 *)

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==0 else sum(int(y[i])*f(i + 1) for i in range(len(y)))
def a(n): return 0 if n==0 else a255411(a(n//2)) if n%2==0 else 1 + a255411(a((n - 1)//2))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 20 2017

;; This is a standalone program:
(define (A276091 n) (let loop ((n n) (s 0) (f 1) (i 2)) (cond ((zero? n) s) ((even? n) (loop (/ n 2) s (* i f) (+ 1 i))) (else (loop (/ (- n 1) 2) (+ s (* (- i 1) f)) (* i f) (+ 1 i))))))
;; This implements one of the given recurrences:
(definec (A276091 n) (cond ((zero? n) n) ((even? n) (A255411 (A276091 (/ n 2)))) (else (+ 1 (A255411 (A276091 (/ (- n 1) 2)))))))
;; Alternatively, we can use A276340 in place of A255411:
(definec (A276091 n) (cond ((zero? n) n) ((even? n) (A276340 (A276091 (/ n 2)))) (else (+ 1 (A276340 (A276091 (/ (- n 1) 2)))))))

a(0) = 0, a(2n) = A255411(a(n)), a(2n+1) = 1+A255411(a(n)).
a(0) = 0, a(2n) = A276340(a(n)), a(2n+1) = 1+A276340(a(n)).
Other identities. For all n >= 0:
a(n) = A225901(A059590(n)).
a(n) = A276090(A275959(n)).
A276328(a(n)) = A276337(a(n)) = A000120(n).

Name changed (to emphasize the functional nature of the sequence) with the original definition moved to the comments by Antti Karttunen, Sep 01 2016

A275959 Sum of distinct terms of A002674: a(0) = 0, a(2n) = A255411(A153880(a(n))), a(2n+1) = 1+A255411(A153880(a(n))).

Original entry on oeis.org

0, 1, 12, 13, 360, 361, 372, 373, 20160, 20161, 20172, 20173, 20520, 20521, 20532, 20533, 1814400, 1814401, 1814412, 1814413, 1814760, 1814761, 1814772, 1814773, 1834560, 1834561, 1834572, 1834573, 1834920, 1834921, 1834932, 1834933, 239500800, 239500801, 239500812, 239500813, 239501160, 239501161, 239501172, 239501173, 239520960, 239520961

Offset: 0

Antti Karttunen, Aug 16 2016

Fixed points of involution A225901.

This can be also viewed as a function that reinterprets base-2 representation of n in base-((2n)!/2) where the digits are multiplied with the successive terms of A002674, thus a(0) = 0.

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

Cf. A002674, A059590, A153880, A255411, A276082, A276083, A276091.
Fixed points of A225901.
Subsequence of A275956 and of A276089.

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==0 else sum([int(y[i])*f(i + 1) for i in range(len(y))])
def a153880(n):
    x=(str(a007623(n)) + '0')[::-1]
    return 0 if n==0 else sum([int(x[i])*f(i + 1) for i in range(len(x))])
def a(n): return 0 if n==0 else a255411(a153880(a(n//2))) if n%2==0 else 1 + a255411(a153880(a((n - 1)//2)))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 20 2017

a(0) = 0, a(2n) = A255411(A153880(a(n))), a(2n+1) = 1+A255411(A153880(a(n))).

a(n) = A276089(A276091(n)).

A276082 a(0) = 0, a(2n) = A153880(a(n)), a(2n+1) = 1+A255411(a(n)).

Original entry on oeis.org

0, 1, 2, 5, 6, 13, 14, 23, 24, 49, 50, 77, 54, 85, 86, 119, 120, 241, 242, 365, 246, 373, 374, 503, 264, 409, 410, 557, 414, 565, 566, 719, 720, 1441, 1442, 2165, 1446, 2173, 2174, 2903, 1464, 2209, 2210, 2957, 2214, 2965, 2966, 3719, 1560, 2401, 2402, 3245, 2406, 3253, 3254, 4103, 2424, 3289, 3290, 4157, 3294, 4165, 4166, 5039, 5040, 10081

Offset: 0

Antti Karttunen, Aug 21 2016

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

Cf. A153880, A255411

Cf. also A059590, A275959, A276091, A225901, A276083.

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==0 else sum([int(y[i])*f(i + 1) for i in range(len(y))])
def a153880(n):
    x=(str(a007623(n)) + '0')[::-1]
    return 0 if n==0 else sum([int(x[i])*f(i + 1) for i in range(len(x))])
def a(n): return 0 if n==0 else a153880(a(n//2)) if n%2==0 else 1 + a255411(a((n - 1)//2))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 20 2017

a(0) = 0, a(2n) = A153880(a(n)), a(2n+1) = 1+A255411(a(n)).
Other identities. For all n >= 0:
a(n) = A225901(A276083(n)).

cp's OEIS Frontend

A059590 Numbers obtained by reinterpreting base-2 representation of n in the factorial base: a(n) = Sum_{k>=0} A030308(n,k)*A000142(k+1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula

Extensions

A276091 Numbers obtained by reinterpreting base-2 representation of n in A001563-base (A276326): a(n) = Sum_{k>=0} A030308(n,k)*A001563(k+1).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Python

Scheme

Formula

Extensions

A275959 Sum of distinct terms of A002674: a(0) = 0, a(2n) = A255411(A153880(a(n))), a(2n+1) = 1+A255411(A153880(a(n))).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Python

Formula

A276082 a(0) = 0, a(2n) = A153880(a(n)), a(2n+1) = 1+A255411(a(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

Python

Formula