cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 51-60 of 62 results. Next

A336882 a(0) = 1; for k >= 0, 0 <= i < 2^k, a(2^k + i) = m_k * a(i), where m_k is the least odd number not in terms 0..2^k - 1.

Original entry on oeis.org

1, 3, 5, 15, 7, 21, 35, 105, 9, 27, 45, 135, 63, 189, 315, 945, 11, 33, 55, 165, 77, 231, 385, 1155, 99, 297, 495, 1485, 693, 2079, 3465, 10395, 13, 39, 65, 195, 91, 273, 455, 1365, 117, 351, 585, 1755, 819, 2457, 4095, 12285, 143, 429, 715, 2145, 1001
Offset: 0

Views

Author

Peter Munn, Aug 16 2020

Keywords

Comments

A permutation of the odd numbers.
Every positive integer, m, is the product of a unique subset of the terms of A050376. The members of the subset are often known as the Fermi-Dirac factors of m. In this sequence, the odd numbers appear lexicographically according to their Fermi-Dirac factors (with those factors listed in decreasing order). The equivalent sequence for all positive integers is A052330.
The sequence has a conditional exponential identity shown in the formula section. This relies on the offset being 0, as in related sequences, notably A019565 and A052330.

Examples

			a(0) = 1, as specified explicitly.
m_0 = 3, the least odd number not in terms 0..0.
So a(1) = a(2^0 + 0) = m_0 * a(0) = 3 * 1 = 3.
m_1 = 5, the least odd number not in terms 0..1.
So a(2) = a(2^1 + 0) = m_1 * a(0) = 5 * 1 = 5;
and a(3) = a(2^1 + 1) = m_1 * a(1) = 5 * 3 = 15.
The initial terms are tabulated below, equated with the product of their Fermi-Dirac factors to exhibit the lexicographic order. We start with 1, since 1 is factored as the empty product and the empty list is first in lexicographic order.
   n     a(n)
   0    1,
   1    3 = 3,
   2    5 = 5,
   3   15 = 5 * 3,
   4    7 = 7,
   5   21 = 7 * 3,
   6   35 = 7 * 5,
   7  105 = 7 * 5 * 3,
   8    9 = 9,
   9   27 = 9 * 3,
  10   45 = 9 * 5,
  11  135 = 9 * 5 * 3,
  12   63 = 9 * 7.
		

Crossrefs

Permutation of A005408.
Subsequence of A052330.
Subsequences: A062090, A332382 (squarefree terms).
A003986, A003987, A004198, A059896, A059897 are used to express relationship between terms of this sequence.

Formula

a(2^k) = min({ 2*m+1 : m >= 0, 2*m+1 <> a(j), 0 <= j < 2^k }) = A062090(k+2).
If x AND y = 0, a(x+y) = a(x) * a(y), where AND denotes the bitwise operation, A004198(.,.).
a(x XOR y) = A059897(a(x), a(y)), where XOR denotes bitwise exclusive-or, A003987(.,.).
a(x OR y) = A059896(a(x), a(y)), where OR denotes the bitwise operation, A003986(.,.).

A344534 For any number n with binary expansion Sum_{k = 1..m} 2^e_k (where 0 <= e_1 < ... < e_m), a(n) = Product_{k = 1..m} prime(1+A002262(e_k))^2^A025581(e_k) (where prime(k) denotes the k-th prime number).

Original entry on oeis.org

1, 2, 4, 8, 3, 6, 12, 24, 16, 32, 64, 128, 48, 96, 192, 384, 9, 18, 36, 72, 27, 54, 108, 216, 144, 288, 576, 1152, 432, 864, 1728, 3456, 5, 10, 20, 40, 15, 30, 60, 120, 80, 160, 320, 640, 240, 480, 960, 1920, 45, 90, 180, 360, 135, 270, 540, 1080, 720, 1440
Offset: 0

Views

Author

Rémy Sigrist, May 22 2021

Keywords

Comments

The ones in the binary expansion of n encode the Fermi-Dirac factors of a(n).
The following table gives the rank of the bit corresponding to the Fermi-Dirac factor p^2^k:
...
7| 9
5| 5 8
3| 2 4 7
2| 0 1 3 6
---+--------
p/k| 0 1 2 3 ...
This sequence is a bijection from the nonnegative integers to the positive integers with inverse A344536.
This sequence establishes a bijection from A261195 to A225547.
This sequence and A344535 each map between two useful choices for encoding sets of elements drawn from a 2-dimensional array. To give a very specific example, each mapping is an isomorphism between two alternative integer representations of the polynomial ring GF2[x,y]. The relevant set is {x^i*y^j : i, j >= 0}. The mappings between the two representations of the ring's addition operation are from XOR (A003987) to A059897(.,.) and for the multiplication operation, they are from A329331(.,.) to A329329(.,.). - Peter Munn, May 31 2021

Examples

			For n = 42:
- 42 = 2^5 + 2^3 + 2^1,
- so we have the following Fermi-Dirac factors p^2^k:
      5| X
      3|
      2|   X X
    ---+------
    p/k| 0 1 2
- a(42) = 2^2^1 * 2^2^2 * 5^2^0 = 320.
		

Crossrefs

Comparable mappings that also use Fermi-Dirac factors: A052330, A059900.
Maps binary operations A003987 to A059897, A329331 to A329329.

Programs

  • PARI
    A002262(n)=n-binomial(round(sqrt(2+2*n)), 2)
    A025581(n)=binomial(1+floor(1/2+sqrt(2+2*n)), 2)-(n+1)
    a(n) = { my (v=1, e); while (n, n-=2^e=valuation(n, 2); v* = prime(1 + A002262(e))^2^A025581(e)); v }

Formula

a(n) = A344535(A344531(n)).
a(n) = A344535(n) iff n belongs to A261195.
A064547(a(n)) = A000120(n).
a(A036442(n)) = prime(n).
a(A006125(n+1)) = 2^2^n for any n >= 0.
a(m + n) = a(m) * a(n) when m AND n = 0 (where AND denotes the bitwise AND operator).
From Peter Munn, Jun 06 2021: (Start)
a(n) = A225546(A344535(n)).
a(n XOR k) = A059897(a(n), a(k)), where XOR denotes bitwise exclusive-or, A003987.
a(A329331(n, k)) = A329329(a(n), a(k)).
(End)

A344535 For any number n with binary expansion Sum_{k = 1..m} 2^e_k (where 0 <= e_1 < ... < e_m), a(n) = Product_{k = 1..m} prime(1+A025581(e_k))^2^A002262(e_k) (where prime(k) denotes the k-th prime number).

Original entry on oeis.org

1, 2, 3, 6, 4, 8, 12, 24, 5, 10, 15, 30, 20, 40, 60, 120, 9, 18, 27, 54, 36, 72, 108, 216, 45, 90, 135, 270, 180, 360, 540, 1080, 16, 32, 48, 96, 64, 128, 192, 384, 80, 160, 240, 480, 320, 640, 960, 1920, 144, 288, 432, 864, 576, 1152, 1728, 3456, 720, 1440
Offset: 0

Views

Author

Rémy Sigrist, May 22 2021

Keywords

Comments

The ones in the binary expansion of n encode the Fermi-Dirac factors of a(n).
The following table gives the rank of the bit corresponding to the Fermi-Dirac factor p^2^k:
...
7| 6
5| 3 7
3| 1 4 8
2| 0 2 5 9
---+--------
p/k| 0 1 2 3 ...
This sequence is a bijection from the nonnegative integers to the positive integers with inverse A344537.
This sequence establishes a bijection from A261195 to A225547.

Examples

			For n = 42:
- 42 = 2^5 + 2^3 + 2^1,
- so we have the following Fermi-Dirac factors p^2^k:
      5| X
      3| X
      2|     X
    ---+------
    p/k| 0 1 2
- a(42) = 3^2^0 * 5^2^0 * 2^2^2 = 240.
		

Crossrefs

Programs

  • PARI
    A002262(n)=n-binomial(round(sqrt(2+2*n)), 2)
    A025581(n)=binomial(1+floor(1/2+sqrt(2+2*n)), 2)-(n+1)
    a(n) = { my (v=1, e); while (n, n-=2^e=valuation(n, 2); v *= prime(1 + A025581(e))^2^A002262(e)); v }

Formula

a(n) = A344534(A344531(n)).
a(n) = A344534(n) iff n belongs to A261195.
A064547(a(n)) = A000120(n).
a(A006125(n)) = prime(n) for any n > 0.
a(A036442(n+1)) = 2^2^n for any n >= 0.
a(m + n) = a(m) * a(n) when m AND n = 0 (where AND denotes the bitwise AND operator).

A303772 Inverse of A303771.

Original entry on oeis.org

0, 1, 3, 5, 9, 2, 17, 6, 33, 10, 65, 4, 129, 18, 12, 257, 513, 34, 1025, 14, 20, 66, 2049, 7, 4097, 130, 36, 22, 8193, 11, 16385, 258, 68, 514, 26, 38, 32769, 1026, 132, 15, 65537, 19, 131073, 70, 42, 2050, 262145, 260
Offset: 1

Views

Author

Antti Karttunen, May 02 2018

Keywords

Crossrefs

Programs

  • PARI
    default(parisizemax,2^31);
    up_to_e = 18;
    up_to = (2 + 2^up_to_e);
    v050376 = vector(2+up_to_e);
    A050376(n) = v050376[n];
    ispow2(n) = (n && !bitand(n,n-1));
    i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == 2+up_to_e,break));
    A052330(n) = { my(p=1,i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); };
    A053669(n) = forprime(p=2, , if (n % p, return(p))); \\ From A053669
    v303760 = vector(up_to);
    m_inverses = Map();
    prev=1; for(n=1,up_to,fordiv(prev,d,if(!mapisdefined(m_inverses,d),v303760[n] = d;mapput(m_inverses,d,n);break)); if(!v303760[n], apu = prev; while(mapisdefined(m_inverses,try = prev*A053669(apu)), apu *= A053669(apu)); v303760[n] = try; mapput(m_inverses,try,n)); prev = v303760[n]);
    A303760(n) = v303760[n+1];
    A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
    A303771(n) = A052330(A048675(A303760(n)));
    m303772 = Map();
    for(n=0,up_to-1,mapput(m303772,A303771(n),n));
    A303772(n) = mapget(m303772,n);

Formula

a(n) = A303768(A052331(n)).

A304534 Inverse of A304533.

Original entry on oeis.org

0, 1, 3, 2, 5, 19, 4, 20, 10, 11, 13, 12, 15, 22, 14, 21, 41, 42, 44, 43, 46, 60, 45, 61, 51, 52, 54, 53, 56, 63, 55, 62, 7, 8, 24, 25, 6, 81, 29, 88, 17, 9, 27, 26, 16, 18, 28, 23, 48, 49, 65, 66, 47, 157, 70, 301, 58, 50, 68, 67, 57, 59, 69, 64, 236, 237, 239, 238, 241, 255, 240, 256, 246, 247, 249, 248, 251, 258, 250, 257, 277, 278, 280, 279, 282, 296, 281
Offset: 0

Views

Author

Antti Karttunen, May 14 2018

Keywords

Crossrefs

Cf. A304533 (inverse).

Formula

For all n >= 0, a(A304533(n)) = n.
a(n) = A304532(A052330(n))-1. [This formula works if A304531 and thus A304533 are indeed permutations, containing all natural numbers.]

A323074 Lexicographically earliest such sequence a that a(i) = a(j) => f(i) = f(j) for all i, j, where f(n) = -(n mod 4) if n is a prime, and f(n) = A300840(n) for any other number.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Jan 04 2019

Keywords

Comments

For all i, j:
A319704(i) = A319704(j) => a(i) = a(j) => A323082(i) = A323082(j).

Crossrefs

Programs

  • PARI
    up_to = 65539;
    rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
    ispow2(n) = (n && !bitand(n,n-1));
    A302777(n) = ispow2(isprimepower(n));
    A050376list(up_to) = { my(v=vector(up_to), i=0); for(n=1,oo,if(A302777(n), i++; v[i] = n); if(i == up_to,return(v))); };
    v050376 = A050376list(up_to);
    A050376(n) = v050376[n];
    A052330(n) = { my(p=1,i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); };
    A052331(n) = { my(s=0,e); while(n > 1, fordiv(n, d, if(((n/d)>1)&&A302777(n/d), e = vecsearch(v050376, n/d); if(!e, print("v050376 too short!"); return(1/0)); s += 2^(e-1); n = d; break))); (s); };
    A300840(n) = A052330(A052331(n)>>1);
    A323074aux(n) = if(isprime(n),-(n%4),A300840(n));
    v323074 = rgs_transform(vector(up_to,n,A323074aux(n)));
    A323074(n) = v323074[n];

A371985 For n a power of 2, a(n) = n. Otherwise a(n) is the smallest novel multiple of a(n - 2^m), where 2^m is the greatest power of 2 not exceeding n.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 9, 8, 7, 10, 12, 20, 15, 18, 27, 16, 11, 14, 21, 24, 25, 30, 36, 40, 28, 50, 48, 60, 45, 54, 81, 32, 13, 22, 33, 44, 35, 42, 63, 56, 49, 70, 72, 80, 75, 90, 108, 96, 55, 84, 105, 120, 100, 150, 144, 160, 112, 200, 192, 180, 135, 162, 243, 64, 17
Offset: 1

Views

Author

David James Sycamore, Apr 15 2024

Keywords

Comments

Reminiscent of the Doudna sequence A005940; also of A052330 and A269848.
All powers of 2 (a(n) = n) are assigned first in order to avoid the second part of the definition giving a(n) = 2^k for some n which is not a power of 2 (see Example for a(12) = 20).
It follows from the definition that all powers of 2, all primes and all multiples of all primes are terms so this sequence is a permutation of the positive integers (A000027), with primes in order.
Each prime power appears before any of its multiples, meaning that this sequence has "property S" as defined in A368900.

Examples

			a(3) = 3, because 2 is the greatest power of 2 not exceeding 3 and 3-2 = 1, so a(3) = 3, the least novel multiple of a(1) = 1.
a(12) is the smallest novel multiple of a(12-8) = a(4) = 4, and at this point in the sequence 4,8,12 are all prior terms and a(16) = 16 is already taken, so a(12) = 20.
		

Crossrefs

Programs

  • Mathematica
    nn = 10; c[] := False; m[] := 1; a[1] = 1; c[1] = True;
    Do[If[i == 0,
       k = 2^j + i,
       (While[Set[k, m[#] #]; Or[c[k], IntegerQ@ Log2[k]], m[#]++]) &@ a[i]];
      Set[{a[2^j + i], c[k]}, {k, True}], {j, nn}, {i, 0, 2^j - 1}];
    Array[a, 2^(nn + 1) - 1] (* Michael De Vlieger, Apr 15 2024 *)
  • PARI
    \\ See PARI link

Formula

a(2^k + 1) = prime(k+1).

A079708 Metaprime binary to standard binary conversion series.

Original entry on oeis.org

0, 1, 2, 3, 6, 12, 20, 28, 140, 260, 64, 11, 30, 420, 7488, 1922800, 11285855256250575, 54979022626732989863421863670075405480
Offset: 0

Views

Author

Will Nicholes, Jan 31 2003

Keywords

Comments

Each term in the series is computed by translating the previous term to binary, then reinterpreting the binary expression as a product of metaprimes. Metaprimes follow the form p^(2^n) where p is a prime number and n is a nonnegative integer. See the link for more detailed explanation.

Examples

			20 (decimal) = 10100 (binary) 10100 (metaprime binary) = 7 * 1 * 4 * 1 * 1 = 28 (decimal).
		

Formula

a(0)=0, a(n) = A052330(a(n-1)). - Thomas Ordowski, Jun 20 2005

Extensions

Link updated by Will Nicholes, Jun 07 2010

A329330 Multiplication operation of a ring over the positive integers that has A059897(.,.) as addition operation and is isomorphic to GF(2)[x] with polynomial x^i mapped to A050376(i+1). Square array read by descending antidiagonals: A(n,k), n >= 1, k >= 1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 4, 4, 1, 1, 5, 5, 5, 5, 1, 1, 6, 7, 7, 7, 6, 1, 1, 7, 12, 9, 9, 12, 7, 1, 1, 8, 9, 20, 11, 20, 9, 8, 1, 1, 9, 15, 11, 35, 35, 11, 15, 9, 1, 1, 10, 11, 28, 13, 8, 13, 28, 11, 10, 1, 1, 11, 21, 13, 45, 63, 63, 45, 13, 21, 11, 1
Offset: 1

Views

Author

Peter Munn, Nov 10 2019

Keywords

Comments

When creating A329329, the author realized it was isomorphic to multiplication in the GF(2)[x,y] polynomial ring. However, A329329 was unusual in having A059897(.,.) as additive operator, whereas the equivalent univariate polynomial ring, GF(2)[x], is more commonly mapped (to integers) with bitwise exclusive-or (A003987) representing polynomial addition (and A048720(.,.) representing polynomial multiplication). This sequence shows how multiplication in GF(2)[x] can look when mapped to integers with A059897(.,.) representing polynomial addition.
The group defined by the binary operation A059897(.,.) over the positive integers is commutative with all elements self-inverse, and isomorphic to the additive group of the polynomial ring GF(2)[x]. There is a unique isomorphism extending each bijective mapping between respective minimal generating sets. The lexicographically earliest minimal generating set for the A059897 group is A050376, often called the Fermi-Dirac primes. The most meaningful generating set for the additive group of GF(2)[x] is {x^i: i >= 0}.
Using f to denote the intended isomorphism from GF(2)[x], we specify f(x^i) = A050376(i+1). This maps minimal generating sets of the additive groups, so the definition of f is completed by specifying f(a+b) = A059897(f(a), f(b)). We then calculate the image under f of polynomial multiplication in GF(2)[x], giving us this sequence as the matching multiplication operation for an isomorphic ring over the positive integers. Using g to denote the inverse of f, A(n,k) = f(g(n) * g(k)).
Note that A050376 is closed with respect to A(.,.).
Recall that GF(2)[x] is more usually mapped to integers with A003987(.,.) as addition and A048720(.,.) as multiplication. With this usual mapping, under which A000079(i) is the image of x^i, A052330(.) is the relevant isomorphism from nonnegative integers under A048720(.,.) and A003987(.,.) to positive integers under A(.,.) and A059897(.,.), with A052331(.) its inverse.

Examples

			Square array A(n,k) begins:
  n\k |  1    2    3    4    5    6    7    8    9   10   11   12
  ----+----------------------------------------------------------
   1  |  1    1    1    1    1    1    1    1    1    1    1    1
   2  |  1    2    3    4    5    6    7    8    9   10   11   12
   3  |  1    3    4    5    7   12    9   15   11   21   13   20
   4  |  1    4    5    7    9   20   11   28   13   36   16   35
   5  |  1    5    7    9   11   35   13   45   16   55   17   63
   6  |  1    6   12   20   35    8   63  120   99  210  143   15
   7  |  1    7    9   11   13   63   16   77   17   91   19   99
   8  |  1    8   15   28   45  120   77   14  117  360  176  420
   9  |  1    9   11   13   16   99   17  117   19  144   23  143
  10  |  1   10   21   36   55  210   91  360  144   22  187  756
  11  |  1   11   13   16   17  143   19  176   23  187   25  208
  12  |  1   12   20   35   63   15   99  420  143  756  208   28
		

Crossrefs

Distributes over A059897, and isomorphic to A048720 over A003987, with A052331 (inverse A052330) as isomorphism.
Row/column 3: A300841.
Row/column k sorted into increasing order: A003159 (k=3), A339690 (k=4), A000379 (k=6).
Subsequences of row/column k: A240521 (k=6), A240522 (k=8), A240536 (k=10), A240524 (k=24), A241025 (k=30), A241024 (k=40).

Formula

A(n,k) = A052330(A048720(A052331(n), A052331(k))), n >= 1, k >= 1.
A059897-based definition: (Start)
A(A050376(i), A050376(j)) = A050376(i+j-1).
A(A059897(n,k), m) = A059897(A(n,m), A(k,m)).
A(m, A059897(n,k)) = A059897(A(m,n), A(m,k)).
(End)
Derived identities: (Start)
A(n,1) = A(1,n) = 1.
A(n,2) = A(2,n) = n.
A(n,k) = A(k,n).
A(n, A(m,k)) = A(A(n,m), k).
(End)
A(A300841(n), k) = A(n, A300841(k)) = A300841(A(n,k)).
A(n,3) = A(3,n) = A300841(n).
A(n,4) = A(4,n) = A300841^2(n).
A(n,5) = A(5,n) = A300841^3(n).
A(A050376(m), 6) = A(6, A050376(m)) = A240521(m).
A(n,7) = A(7,n) = A300841^4(n).
A(A050376(m), 8) = A(8, A050376(m)) = A240522(m).
A(n,9) = A(9,n) = A300841^5(n).
A(A050376(m), 10) = A(10, A050376(m)) = A240536(m).
A(A050376(m), 12) = A(12, A050376(m)) = A300841(A240521(m)).
A(A050376(m), 24) = A(24, A050376(m)) = A240524(m).
A(A050376(m), 30) = A(30, A050376(m)) = A241025(m).
A(A050376(m), 40) = A(40, A050376(m)) = A241024(m).

A064359 Inverse of sequence A052331 considered as a permutation of the natural numbers.

Original entry on oeis.org

1, 2, 5, 3, 7, 11, 23, 4, 9, 14, 29, 19, 39
Offset: 1

Views

Author

Howard A. Landman, Sep 25 2001

Keywords

Comments

Sequence A052331 claims to be the "inverse" of A052330.

Crossrefs

Previous Showing 51-60 of 62 results. Next