A096113 -id:A096113 - cp's OEIS frontend

A096111 If n = 2^k - 1, then a(n) = k+1, otherwise a(n) = (A000523(n)+1)*a(A053645(n)).

1, 2, 2, 3, 3, 6, 6, 4, 4, 8, 8, 12, 12, 24, 24, 5, 5, 10, 10, 15, 15, 30, 30, 20, 20, 40, 40, 60, 60, 120, 120, 6, 6, 12, 12, 18, 18, 36, 36, 24, 24, 48, 48, 72, 72, 144, 144, 30, 30, 60, 60, 90, 90, 180, 180, 120, 120, 240, 240, 360, 360, 720, 720, 7, 7, 14, 14, 21, 21

Offset: 0

Amarnath Murthy, Jun 29 2004

nonn

Each n > 1 occurs 2*A045778(n) times in the sequence.
f(n+2^k) = (k+1)*f(n) if 2^k > n+1. - Robert Israel, Apr 25 2016
If the binary indices of n (row n of A048793) are the positions 1's in its reversed binary expansion, then a(n) is the product of all binary indices of n + 1. The number of binary indices of n is A000120(n), their sum is A029931(n), and their average is A326699(n)/A326700(n). - Gus Wiseman, Jul 27 2019

Peter Kagey, Table of n, a(n) for n = 0..10000

Permutation of A096115, i.e. a(n) = A096115(A122198(n+1)) [Note the different starting offsets]. Bisection: A121663. Cf. A096113, A052330.
Cf. A029931.
Cf. A000120, A034797, A048793, A070939, A291166, A326031, A326672, A326673.

f:= proc(n) local L;
    L:= convert(2*n+2,base,2);
    convert(subs(0=NULL,zip(`*`,L, [$0..nops(L)-1])),`*`);
end proc:
map(f, [$0..100]); # Robert Israel, Apr 25 2016

CoefficientList[(Product[1 + k x^(2^(k - 1)), {k, 7}] - 1)/x, x] (* Michael De Vlieger, Apr 08 2016 *)
bpe[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];Table[Times@@bpe[n+1],{n,0,100}] (* Gus Wiseman, Jul 26 2019 *)

N=166; q='q+O('q^N);
gf= (prod(n=1,1+ceil(log(N)/log(2)), 1+n*q^(2^(n-1)) ) - 1) / q;
Vec(gf)
/* Joerg Arndt, Oct 06 2012 */

(define (A096111 n) (cond ((pow2? (+ n 1)) (+ 2 (A000523 n))) (else (* (+ 1 (A000523 n)) (A096111 (A053645 n))))))
(define (pow2? n) (and (> n 0) (zero? (A004198bi n (- n 1)))))

G.f.: ( prod(k>=1, 1+k*x^(2^(k-1)) )- 1 ) / x. - Vladeta Jovovic, Nov 08 2005

a(n) is the product of the exponents in the binary expansion of 2*n + 2. - Peter Kagey, Apr 24 2016

Edited, extended and Scheme code added by Antti Karttunen, Aug 25 2006

A052330 Let S_k denote the first 2^k terms of this sequence and let b_k be the smallest positive integer that is not in S_k; then the numbers b_k*S_k are the next 2^k terms.

Original entry on oeis.org

1, 2, 3, 6, 4, 8, 12, 24, 5, 10, 15, 30, 20, 40, 60, 120, 7, 14, 21, 42, 28, 56, 84, 168, 35, 70, 105, 210, 140, 280, 420, 840, 9, 18, 27, 54, 36, 72, 108, 216, 45, 90, 135, 270, 180, 360, 540, 1080, 63, 126, 189, 378, 252, 504, 756, 1512, 315, 630, 945, 1890

Offset: 0

Christian G. Bower, Dec 15 1999

Inverse of sequence A064358 considered as a permutation of the positive integers. - Howard A. Landman, Sep 25 2001
This sequence is not exactly a permutation because it has offset 0 but doesn't contain 0. A052331 is its exact inverse, which has offset 1 and contains 0. See also A064358.
Are there any other values of n besides 4 and 36 with a(n) = n? - Thomas Ordowski, Apr 01 2005
4 = 100 = 4^1 * 3^0 * 2^0, 36 = 100100 = 9^1 * 7^0 * 5^0 * 4^1 * 3^0 * 2^0. - Thomas Ordowski, May 26 2005
Ordering of positive integers by increasing "Fermi-Dirac representation", which is a representation of the "Fermi-Dirac factorization", term implying that each prime power with a power of two as exponent may appear at most once in the "Fermi-Dirac factorization" of n. (Cf. comment in A050376; see also the OEIS Wiki page.) - Daniel Forgues, Feb 11 2011
The subsequence consisting of the squarefree terms is A019565. - Peter Munn, Mar 28 2018
Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. The FDH-number of a strict integer partition (y_1,...,y_k) is f(y_1)*...*f(y_k). A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793. Then a(n) is the number whose binary indices are the parts of the strict integer partition with FDH-number n. - Gus Wiseman, Aug 19 2019
The set of indices of odd-valued terms has asymptotic density 0. In this sense (using the order they appear in this permutation) 100% of numbers are even. - Peter Munn, Aug 26 2019

			Terms following 5 are 10, 15, 30, 20, 40, 60, 120; this is followed by 7 as 6 has already occurred. - _Philippe Deléham_, Jun 03 2015
From _Antti Karttunen_, Apr 13 2018, after also _Philippe Deléham_'s Jun 03 2015 example: (Start)
This sequence can be regarded also as an irregular triangle with rows of lengths 1, 1, 2, 4, 8, 16, ..., that is, it can be represented as a binary tree, where each left hand child contains A300841(k), and each right hand child contains 2*A300841(k), when their parent contains k:
                                     1
                                     |
                  ...................2...................
                 3                                       6
       4......../ \........8                  12......../ \........24
      / \                 / \                 / \                 / \
     /   \               /   \               /   \               /   \
    /     \             /     \             /     \             /     \
   5       10         15       30         20       40         60      120
  7 14   21  42     28  56   84  168    35  70  105  210   140 280  420 840
  etc.
Compare also to trees like A005940 and A283477, and sequences A207901 and A302783.
(End)

Antti Karttunen, Table of n, a(n) for n = 0..8191 (Terms 0..1023 from T. D. Noe)
Michael De Vlieger, Fan style binary tree showing a(n), n = 0..16383, with a color function representing primes in red, perfect powers of primes in gold, squarefree composites in greens, and numbers neither squarefree nor prime powers in blue or purple. Purple represents powerful numbers that are not prime powers. This is a 15 level version of Karttunen's diagram in the example.
OEIS Wiki, Ordering of positive integers by increasing "Fermi-Dirac representation"
Index entries for sequences that are permutations of the natural numbers

Subsequences: A019565 (squarefree terms), A050376 (the left edge from 2 onward), A336882 (odd terms).
Cf. A050030, A052331 (inverse), A096111, A096113, A096114, A096115, A096116, A096118, A096119, A182979, A207901, A300841, A302023, A302783.
Cf. A000120, A029931, A048793, A064547, A070939, A213925, A299755, A299757, A327041.

a = {1}; Do[a = Join[a, a*Min[Complement[Range[Max[a] + 1], a]]], {n, 1, 6}]; a (* Ivan Neretin, May 09 2015 *)

up_to_e = 13; \\ Good for computing up to n = (2^13)-1
v050376 = vector(up_to_e);
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to_e,break));
A050376(n) = v050376[n];
A052330(n) = { my(p=1,i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); }; \\ Antti Karttunen, Apr 12 2018

a(0)=1; a(n+2^k)=a(n)*b(k) for n < 2^k, k = 0, 1, ... where b is A050376. - Thomas Ordowski, Mar 04 2005
The binary representation of n, n = Sum_{i=0..1+floor(log_2(n))} n_i * 2^i, n_i in {0,1}, is taken as the "Fermi-Dirac representation" (A182979) of a(n), a(n) = Product_{i=0..1+floor(log_2(n))} (b_i)^(n_i) where b_i is A050376(i), i.e., the i-th "Fermi-Dirac prime" (prime power with exponent being a power of 2). - Daniel Forgues, Feb 12 2011
From Antti Karttunen, Apr 12 & 17 2018: (Start)
a(0) = 1; a(2n) = A300841(a(n)), a(2n+1) = 2*A300841(a(n)).
a(n) = A207901(A006068(n)) = A302783(A003188(n)) = A302781(A302845(n)).
(End)

Entry revised Mar 17 2005 by N. J. A. Sloane, based on comments from several people, especially David Wasserman and Thomas Ordowski

A096115 If n = (2^k)-1, a(n) = a((n+1)/2) = k, if n = 2^k, a(n) = a(n-1)+1 = k+1, otherwise a(n) = (A000523(n)+1)*a(A035327(n-1)).

Original entry on oeis.org

1, 2, 2, 3, 6, 6, 3, 4, 12, 24, 24, 12, 8, 8, 4, 5, 20, 40, 40, 60, 120, 120, 60, 20, 15, 30, 30, 15, 10, 10, 5, 6, 30, 60, 60, 90, 180, 180, 90, 120, 360, 720, 720, 360, 240, 240, 120, 30, 24, 48, 48, 72, 144, 144, 72, 24, 18, 36, 36, 18, 12, 12, 6, 7, 42, 84, 84, 126

Offset: 1

Amarnath Murthy, Jun 30 2004

nonn

A fractal sequence. For k in range [1,(2^n)-1], a(2^n + k)/a(2^n - k) = n+1. Each n > 1 occurs 2*A045778(n) times in the sequence.

Permutation of A096111, i.e. a(n) = A096111(A122199(n)-1) [Note the different starting offsets]. Cf. A096113, A052330, A096114, A096116.

(define (A096115 n) (cond ((pow2? (+ n 1)) (+ 1 (A000523 n))) ((pow2? n) (+ 1 (A096115 (- n 1)))) (else (* (+ (A000523 n) 1) (A096115 (A035327 (- n 1)))))))
(define (pow2? n) (and (> n 0) (zero? (A004198bi n (- n 1)))))

Edited, extended and Scheme code added by Antti Karttunen, Aug 25 2006

A096114 a(1)=1, a(2)=2, a(32^k) = 32^k, a(32^k + i) = 32^k + a(32^k - i), for i in range [1, 32^k - 1].

Original entry on oeis.org

1, 2, 3, 5, 4, 6, 10, 11, 9, 8, 7, 12, 19, 20, 21, 23, 22, 18, 16, 17, 15, 14, 13, 24, 37, 38, 39, 41, 40, 42, 46, 47, 45, 44, 43, 36, 31, 32, 33, 35, 34, 30, 28, 29, 27, 26, 25, 48, 73, 74, 75, 77, 76, 78, 82, 83, 81, 80, 79, 84, 91, 92, 93, 95, 94, 90, 88, 89, 87, 86, 85, 72

Offset: 1

Amarnath Murthy, Jun 29 2004

nonn

Inverse: A121664. Cf. A096111, A096113, A052330, A096116.

a = {1, 2}; Do[a = Join[a, {3*2^k}, 3*2^k + Reverse[a]], {k, 0, 4}]; a (* Ivan Neretin, Sep 04 2017 *)

Edited, extended and Scheme code added by Antti Karttunen, Aug 25 2006

A096116 a(1)=1, if n=(2^k)+1, a(n) = k+2, otherwise a(n) = 2+A000523(n-1)+a(2+A035327(n-1)).

Original entry on oeis.org

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

Offset: 1

Amarnath Murthy, Jun 30 2004

nonn

Each n > 1 occurs A025147(n) times in the sequence.

Ivan Neretin, Table of n, a(n) for n = 1..10000

Cf. A096111, A050029, A050030, A052330, A096113, A096114, A096115.

a = {1}; Do[AppendTo[a, If[BitAnd[n - 1, n - 2] == 0, Log2[n - 1] + 2, 2 + Floor[Log2[n - 1]] + a[[2 + BitXor[n - 1, 2^Ceiling[Log2[n]] - 1]]]]], {n, 2, 74}]; a (* Ivan Neretin, Jun 24 2016 *)

(define (A096116 n) (cond ((= 1 n) 1) ((pow2? (- n 1)) (+ 2 (A000523 (- n 1)))) (else (+ 2 (A000523 (- n 1)) (A096116 (+ 2 (A035327 (- n 1))))))))
(define (pow2? n) (and (> n 0) (zero? (A004198bi n (- n 1)))))
;; Antti Karttunen, Aug 25 2006

Edited and extended by Antti Karttunen, Aug 25 2006

A123664 a(1) = 1, a(2) = 2; then all new products of subsets of pre-existing terms which include the most recent, then the first integer not present and so on.

Original entry on oeis.org

1, 2, 3, 6, 4, 8, 12, 24, 48, 72, 144, 5, 10, 15, 20, 30, 40, 60, 80, 90, 120, 160, 180, 240, 320, 360, 480, 720, 960, 1080, 1440, 1920, 2160, 2880, 3840, 4320, 5760, 6480, 7680, 8640, 11520, 12960, 15360, 17280, 23040, 25920, 34560, 46080, 51840, 69120, 77760

Offset: 1

Joel B. Lewis, Nov 15 2006

nonn

Similar to A096113. However, each product must include the most recently added singleton. Thus after adding 4, the terms 18 and 36 are not added because they have no representation as a product of earlier terms, including 4. A110797 is similar, but only allows products of pairs (not of subsets).

			After a(1) = 1, a(2) = 2, all products are present, so we add the first integer not included, namely 3. Then we add all products of any subset of {1, 2, 3} which include 3 and are not already present, in this case just 6. Then we add the next integer not already present, 4. Then we add all products of any subset of {1, 2, 3, 6, 4} which include 4 and are not already present, 8 (=2*4), 12 (=3*4), 24 (=2*3*4=6*4), 48 (=2*6*4), 72 (=3*6*4) and 144 (=2*3*6*4). Then we add 5, the next integer not already present. And so on.

Cf. A096113.

M[2]={1,2} M[n_]:= Join[M[n-1], Complement[Union[M[n-1][[ -1]] * Exp[Map[Total, Log[Subsets[Delete[Delete[M[n-1],1],-1]]]]]], M[n-1]],{n}] M[6]

cp's OEIS Frontend

A096111 If n = 2^k - 1, then a(n) = k+1, otherwise a(n) = (A000523(n)+1)*a(A053645(n)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Scheme

Formula

Extensions

A052330 Let S_k denote the first 2^k terms of this sequence and let b_k be the smallest positive integer that is not in S_k; then the numbers b_k*S_k are the next 2^k terms.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A096115 If n = (2^k)-1, a(n) = a((n+1)/2) = k, if n = 2^k, a(n) = a(n-1)+1 = k+1, otherwise a(n) = (A000523(n)+1)*a(A035327(n-1)).

Views

Author

Keywords

Comments

Crossrefs

Programs

Scheme

Extensions

A096114 a(1)=1, a(2)=2, a(3*2^k) = 3*2^k, a(3*2^k + i) = 3*2^k + a(3*2^k - i), for i in range [1, 3*2^k - 1].

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Extensions

A096116 a(1)=1, if n=(2^k)+1, a(n) = k+2, otherwise a(n) = 2+A000523(n-1)+a(2+A035327(n-1)).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Scheme

Extensions

A123664 a(1) = 1, a(2) = 2; then all new products of subsets of pre-existing terms which include the most recent, then the first integer not present and so on.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A096114 a(1)=1, a(2)=2, a(32^k) = 32^k, a(32^k + i) = 32^k + a(32^k - i), for i in range [1, 32^k - 1].