[go: up one dir, main page]

login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
Search: a096115 -id:a096115
Displaying 1-5 of 5 results found. page 1
     Sort: relevance | references | number | modified | created      Format: long | short | data
A096111 If n = 2^k - 1, then a(n) = k+1, otherwise a(n) = (A000523(n)+1)*a(A053645(n)). +10
94
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 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
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
LINKS
FORMULA
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
MAPLE
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
MATHEMATICA
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 *)
PROG
(Scheme:) (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)))))
(PARI)
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 */
CROSSREFS
Permutation of A096115, i.e. a(n) = A096115(A122198(n+1)) [Note the different starting offsets]. Bisection: A121663. Cf. A096113, A052330.
Cf. A029931.
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Jun 29 2004
EXTENSIONS
Edited, extended and Scheme code added by Antti Karttunen, Aug 25 2006
STATUS
approved
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. +10
63
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 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
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
LINKS
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.
FORMULA
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)
EXAMPLE
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)
MATHEMATICA
a = {1}; Do[a = Join[a, a*Min[Complement[Range[Max[a] + 1], a]]], {n, 1, 6}]; a (* Ivan Neretin, May 09 2015 *)
PROG
(PARI)
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
CROSSREFS
Subsequences: A019565 (squarefree terms), A050376 (the left edge from 2 onward), A336882 (odd terms).
KEYWORD
nonn,look,tabf
AUTHOR
Christian G. Bower, Dec 15 1999
EXTENSIONS
Entry revised Mar 17 2005 by N. J. A. Sloane, based on comments from several people, especially David Wasserman and Thomas Ordowski
STATUS
approved
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)). +10
4
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 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Each n > 1 occurs A025147(n) times in the sequence.
LINKS
MATHEMATICA
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 *)
PROG
(Scheme) (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
CROSSREFS
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Jun 30 2004
EXTENSIONS
Edited and extended by Antti Karttunen, Aug 25 2006
STATUS
approved
A122198 Permutation of natural numbers: a recursed variant of A122155. +10
4
0, 1, 2, 3, 4, 7, 6, 5, 8, 15, 14, 13, 12, 9, 10, 11, 16, 31, 30, 29, 28, 25, 26, 27, 24, 17, 18, 19, 20, 23, 22, 21, 32, 63, 62, 61, 60, 57, 58, 59, 56, 49, 50, 51, 52, 55, 54, 53, 48, 33, 34, 35, 36, 39, 38, 37, 40, 47, 46, 45, 44, 41, 42, 43, 64, 127, 126, 125, 124, 121 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Maps between A096115 and A096111.
LINKS
FORMULA
a(0)=0, otherwise a(n) = A122155(A053644(n)+a(A053645(n))).
PROG
(Scheme:) (define (A122198 n) (if (< n 1) n (A122155 (+ (A053644 n) (A122198 (A053645 n))))))
CROSSREFS
Inverse: A122199.
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 25 2006
STATUS
approved
A122199 Permutation of natural numbers: a recursed variant of A122155. +10
4
0, 1, 2, 3, 4, 7, 6, 5, 8, 13, 14, 15, 12, 11, 10, 9, 16, 25, 26, 27, 28, 31, 30, 29, 24, 21, 22, 23, 20, 19, 18, 17, 32, 49, 50, 51, 52, 55, 54, 53, 56, 61, 62, 63, 60, 59, 58, 57, 48, 41, 42, 43, 44, 47, 46, 45, 40, 37, 38, 39, 36, 35, 34, 33, 64, 97, 98, 99, 100, 103, 102 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
COMMENTS
Maps between A096111 and A096115.
LINKS
FORMULA
a(0)=0, otherwise a(n) = A053644(A122155(n)) + a(A053645(A122155(n))).
PROG
(Scheme:) (define (A122199 n) (if (< n 1) n (let ((m (A122155 n))) (+ (A053644 m) (A122199 (A053645 m))))))
CROSSREFS
Inverse: A122198.
KEYWORD
nonn
AUTHOR
Antti Karttunen, Aug 25 2006
STATUS
approved
page 1

Search completed in 0.007 seconds

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 29 17:19 EDT 2024. Contains 375518 sequences. (Running on oeis4.)