A302845 -id:A302845 - cp's OEIS frontend

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.

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

A163355 Permutation of integers for constructing Hilbert curve in N x N grid.

Original entry on oeis.org

0, 1, 3, 2, 14, 15, 13, 12, 4, 7, 5, 6, 8, 11, 9, 10, 16, 19, 17, 18, 20, 21, 23, 22, 30, 29, 31, 28, 24, 25, 27, 26, 58, 57, 59, 56, 54, 53, 55, 52, 60, 61, 63, 62, 50, 51, 49, 48, 32, 35, 33, 34, 36, 37, 39, 38, 46, 45, 47, 44, 40, 41, 43, 42, 234, 235, 233, 232, 236, 239

Offset: 0

Antti Karttunen, Jul 29 2009

nonn

Antti Karttunen, Table of n, a(n) for n = 0..262143
Index entries for sequences that are permutations of the natural numbers

Inverse: A163356. A163357 & A163359 give two variants of Hilbert curve in N x N grid. Cf. also A163332.
Second and third "powers": A163905, A163915.
In range [A000302(n-1)..A024036(n)] of this permutation, the number of cycles is given by A163910, number of fixed points seems to be given by A147600(n-1) (fixed points themselves: A163901). Max. cycle sizes is given by A163911 and LCM's of all cycle sizes by A163912.
See also: A000695, A163890, A163894, A163902-A163903, A163914, A163485, A302843, A302845.

A057300 := proc(n)
    option remember;
    `if`(n=0, 0, procname(iquo(n, 4, 'r'))*4+[0, 2, 1, 3][r+1])
end proc:
A163355 := proc(n)
    option remember ;
    local d,base4,i,r ;
    if n <= 1 then
        return n ;
    end if;
    base4 := convert(n,base,4) ;
    d := op(-1,base4) ;
    i := nops(base4)-1 ;
    r := n-d*4^i ;
    if ( d=1 and type(i,even) ) or ( d=2 and type(i,odd)) then
        4^i+procname(A057300(r)) ;
    elif d= 3 then
        2*4^i+procname(A057300(r)) ;
    else
        3*4^i+procname(4^i-1-r) ;
    end if;
end proc:
seq(A163355(n),n=0..100) ; # R. J. Mathar, Nov 22 2023

A057300(n) = { my(t=1, s=0); while(n>0,  if(1==(n%4),n++,if(2==(n%4),n--)); s += (n%4)*t; n >>= 2; t <<= 2); (s); };
A163355(n) = if(!n,n,my(i = (#binary(n)-1)\2, f = 4^i, d = (n\f)%4, r = (n%f)); if(((1==d)&&!(i%2))||((2==d)&&(i%2)), f+A163355(A057300(r)), if(3==d,f+f+A163355(A057300(r)), (3*f)+A163355(f-1-r)))); \\ Antti Karttunen, Apr 14 2018

a(0) = 0, and given d=1, 2 or 3, then a((d*(4^i))+r)
  = (4^i) + a(A057300(r)), if d=1 and i is even, or if d=2 and i is odd
  = 2*(4^i) + a(A057300(r)), if d=3,
  = 3*(4^i) + a((4^i)-1-r) in other cases.
From Alan Michael Gómez Calderón, May 06 2025: (Start)
a(3*A000695(n)) = 2*A000695(n);
a(3*(A000695(n) + 2^A000695(2*m))) = 2*(A000695(n) + 2^A000695(2*m)) for m >= 2;
a((2 + 16^n)*2^(-1 + 4*m)) = 4^(2*(n + m) - 1) + (11*16^m - 2)/3. (End)

Links to further derived sequences added by Antti Karttunen, Sep 21 2009

A302846 Interleave the Gray-coded X and Y-coordinates of 2-dimensional Hilbert's curve in alternate bit-positions: a(n) = A000695(A003188(A059253(n))) + 2*A000695(A003188(A059252(n))).

Original entry on oeis.org

0, 1, 3, 2, 10, 8, 9, 11, 15, 13, 12, 14, 6, 7, 5, 4, 20, 22, 23, 21, 17, 16, 18, 19, 27, 26, 24, 25, 29, 31, 30, 28, 60, 62, 63, 61, 57, 56, 58, 59, 51, 50, 48, 49, 53, 55, 54, 52, 36, 37, 39, 38, 46, 44, 45, 47, 43, 41, 40, 42, 34, 35, 33, 32, 160, 162, 163, 161, 165, 164, 166, 167, 175, 174, 172, 173, 169, 171, 170, 168, 136

Offset: 0

Antti Karttunen, Apr 14 2018

Like in binary Gray code A003188, also in this permutation the binary expansions of a(n) and a(n+1) differ always by just a single bit-position, that is, A000120(A003987(a(n),a(n+1))) = 1 for all n >= 0. Here A003987 computes bitwise-XOR of its two arguments.

When composed with A052330 this gives A302781.

Cf. A302845 (inverse permutation).
Cf. A000695, A302844, A057300, A059252, A059253, A064706, A163356, A302781.
Cf. also A003188, A163252, A300838 for other permutations satisfying the same condition.

A064706(n) = bitxor(n, n>>2);
A057300(n) = { my(t=1,s=0); while(n>0, if(1==(n%4),n++,if(2==(n%4),n--)); s += (n%4)*t; n >>= 2; t <<= 2); (s); };
A163356(n) = if(!n,n,my(i = (#binary(n)-1)\2, f = 4^i, d = (n\f)%4, r = (n%f)); (((((2+(i%2))^d)%5)-1)*f) + if(3==d,f-1-A163356(r),A057300(A163356(r))));
A302846(n) = A064706(A163356(n));

a(n) = A000695(A003188(A059253(n))) + 2*A000695(A003188(A059252(n))).

a(n) = A064706(A163356(n)) = A003188(A302844(n)).

A302843 Permutation of nonnegative integers: a(n) = A163355(A006068(n)).

Original entry on oeis.org

0, 1, 2, 3, 12, 13, 14, 15, 10, 9, 8, 11, 4, 7, 6, 5, 26, 27, 24, 25, 30, 29, 28, 31, 16, 19, 18, 17, 22, 23, 20, 21, 42, 43, 40, 41, 46, 45, 44, 47, 32, 35, 34, 33, 38, 39, 36, 37, 58, 57, 56, 59, 52, 55, 54, 53, 48, 49, 50, 51, 60, 61, 62, 63, 192, 193, 194, 195, 204, 205, 206, 207, 202, 201, 200, 203, 196, 199, 198, 197, 218, 219, 216, 217, 222

Offset: 0

Antti Karttunen, Apr 14 2018

nonn

Antti Karttunen, Table of n, a(n) for n = 0..16383
Index entries for sequences that are permutations of the natural numbers

Cf. A302844 (inverse).

Cf. A003188, A006068, A057300, A163355, A302845.

A006068(n)= { my(s=1, ns); while(1, ns = n >> s; if(0==ns, break()); n = bitxor(n, ns); s <<= 1; ); return (n); } \\ From A006068
A057300(n) = { my(t=1, s=0); while(n>0,  if(1==(n%4),n++,if(2==(n%4),n--)); s += (n%4)*t; n >>= 2; t <<= 2); (s); };
A163355(n) = if(!n,n,my(i = (#binary(n)-1)\2, f = 4^i, d = (n\f)%4, r = (n%f)); if(((1==d)&&!(i%2))||((2==d)&&(i%2)), f+A163355(A057300(r)), if(3==d,f+f+A163355(A057300(r)), (3*f)+A163355(f-1-r))));
A302843(n) = A163355(A006068(n));

a(n) = A163355(A006068(n)).

a(n) = A302845(A003188(n)).

A302782 Inverse permutation to A302781.

Original entry on oeis.org

0, 1, 3, 15, 5, 2, 21, 14, 63, 6, 255, 12, 85, 20, 4, 341, 1023, 62, 4095, 10, 22, 254, 1365, 13, 5461, 86, 60, 16, 16383, 7, 65535, 340, 252, 1022, 26, 48, 21845, 4094, 84, 9, 87381, 23, 262143, 240, 58, 1366, 1048575, 342, 349525, 5460, 1020, 90, 1398101, 61, 250, 19, 4092, 16382, 4194303, 11, 16777215, 65534, 42

Offset: 1

Antti Karttunen, Apr 14 2018

nonn

Antti Karttunen, Table of n, a(n) for n = 1..8192
Index entries for sequences that are permutations of the natural numbers

Cf. A302781 (inverse).

Cf. A006068, A052331, A163355, A302845.

up_to = 8192;
v050376 = vector(up_to);
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,break));
A052331(n) = { my(s=0,e); while(n > 1, fordiv(n, d, if(((n/d)>1)&&ispow2(isprimepower(n/d)), e = vecsearch(v050376, n/d); if(!e, print("v050376 too short!"); return(1/0)); s += 2^(e-1); n = d; break))); (s); };
A006068(n)= { my(s=1, ns); while(1, ns = n >> s; if(0==ns, break()); n = bitxor(n, ns); s <<= 1; ); return (n); } \\ From A006068
A057300(n) = { my(t=1,s=0); while(n>0, if(1==(n%4),n++,if(2==(n%4),n--)); s += (n%4)*t; n >>= 2; t <<= 2); (s); };
A163355(n) = if(!n,n,my(i = (#binary(n)-1)\2, f = 4^i, d = (n\f)%4, r = (n%f)); if(((1==d)&&!(i%2))||((2==d)&&(i%2)), f+A163355(A057300(r)), if(3==d,f+f+A163355(A057300(r)), (3*f)+A163355(f-1-r))));
A302845(n) = A163355(A006068(A006068(n)));
A302782(n) = A302845(A052331(n));

a(n) = A302845(A052331(n)).

cp's OEIS Frontend

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

A163355 Permutation of integers for constructing Hilbert curve in N x N grid.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

PARI

Formula

Extensions

A302846 Interleave the Gray-coded X and Y-coordinates of 2-dimensional Hilbert's curve in alternate bit-positions: a(n) = A000695(A003188(A059253(n))) + 2*A000695(A003188(A059252(n))).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A302843 Permutation of nonnegative integers: a(n) = A163355(A006068(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A302782 Inverse permutation to A302781.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula