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 21-30 of 33 results. Next

A029747 Numbers of the form 2^k times 1, 3 or 5.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 64, 80, 96, 128, 160, 192, 256, 320, 384, 512, 640, 768, 1024, 1280, 1536, 2048, 2560, 3072, 4096, 5120, 6144, 8192, 10240, 12288, 16384, 20480, 24576, 32768, 40960, 49152, 65536, 81920, 98304, 131072, 163840, 196608
Offset: 1

Views

Author

Keywords

Comments

Fixed points of the Doudna sequence: A005940(a(n)) = A005941(a(n)) = a(n). - Reinhard Zumkeller, Aug 23 2006
Subsequence of A103969. - R. J. Mathar, Mar 06 2010
Question: Is there a simple proof that A005940(c) = c would never allow an odd composite c as a solution? See also my comments in A163511 and in A335431 concerning similar problems, also A364551 and A364576. - Antti Karttunen, Jul 28 & Aug 11 2023

Examples

			128 = 2^7 * 1 is in the sequence as well as 160 = 2^5 * 5. - _David A. Corneth_, Sep 18 2020
		

Crossrefs

Subsequence of the following sequences: A103969, A253789, A364541, A364542, A364544, A364546, A364548, A364550, A364560, A364565.
Even terms form a subsequence of A320674.

Programs

  • Mathematica
    m = 200000; Select[Union @ Flatten @ Outer[Times, {1, 3, 5}, 2^Range[0, Floor[Log2[m]]]], # < m &] (* Amiram Eldar, Oct 15 2020 *)
  • PARI
    is(n) = n>>valuation(n, 2) <= 5 \\ David A. Corneth, Sep 18 2020
    
  • Python
    def A029747(n):
        if n<3: return n
        a, b = divmod(n,3)
        return 1<Chai Wah Wu, Apr 02 2025

Formula

a(n) = if n < 6 then n else 2*a(n-3). - Reinhard Zumkeller, Aug 23 2006
G.f.: (1+x+x^2)^2/(1-2*x^3). - R. J. Mathar, Mar 06 2010
Sum_{n>=1} 1/a(n) = 46/15. - Amiram Eldar, Oct 15 2020

Extensions

Edited by David A. Corneth and Peter Munn, Sep 18 2020

A368900 LCM-transform of Doudna sequence.

Original entry on oeis.org

1, 2, 3, 2, 5, 1, 3, 2, 7, 1, 1, 1, 5, 1, 3, 2, 11, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 5, 1, 3, 2, 13, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 5, 1, 3, 2, 17, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 13, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Keywords

Comments

Let's define "property S" for sequences as follows: If s is any sequence of positive natural numbers, normalized to begin with offset 1, then it satisfies the S-property if LCM-transform(s) is equal to the sequence obtained by applying A014963 to sequence s, or in other words, when for all n >= 1, lcm {s(1)..s(n)} / lcm {s(1)..s(n-1)} = A014963(s(n)). This holds if and only if, for all n >= 1, when, either (case A): s(n) is of the form p^k, p prime, then gcd(s(n), lcm {s(1)..s(n-1)}) must be equal to p^(k-1), or (case B): when s(n) is not a prime power, then gcd(s(n), lcm {s(1)..s(n-1)}) must be equal to s(n). Together the cases (A) and (B) reduce to the condition that each prime power should appear in s before any of its multiples do.
Clearly the Doudna-sequence satisfies the property by the way of its construction, as do many of its variants like A356867 (see A369060).
Also, for any base-2 related permutation b that keeps all the numbers of range [2^k, 2^(1+k)[ in the same range, i.e., if for all n >= 1, A000523(b(n)) = A000523(n), then the above property is automatically satisfied.
Furthermore, because in Doudna-sequence no multiple of any term is located on the same row as the term itself (see the tree-illustration in A005940), it follows that any composition of A005940 with any such base-2 related permutation as mentioned above also automatically satisfies the S-property, for example, the permutations A163511, A243353, A253563, A253565, A366260, A366263 and A366275.
Note: Like A005940 itself, also this sequence might be more logical with the starting offset 0 instead of 1, to better align with the underlying mapping from the binary expansion of n to the prime factorization. - Antti Karttunen, Jan 24 2024

Crossrefs

List of LCM-transforms of permutations (permutation given in parentheses):
Cf. A265576 (A064413; note that the EKG sequence permutation does not satisfy the S-property).
In all following cases, the permutation satisfies the S-property:
Cf. A369041 (A003188), A369042 (A006068), A369043 (A193231), A369044 (A057889), A369041 (A054429). [Base-2 related permutations]
Other permutations that have the same property: A303767, (and when used as an offset=1 sequence): A052330.

Programs

  • Mathematica
    nn = 120; Array[Set[{s[#], a[#]}, {#, #}] &, 2]; j = 2;
    Do[If[EvenQ[n],
      Set[s[n], 2 s[n/2]],
      Set[s[n],
        Times @@ Power @@@ Map[{Prime[PrimePi[#1] + 1], #2} & @@ # &,
          FactorInteger[s[(n + 1)/2]]]]];
      k = LCM[j, s[n]]; a[n] = k/j; j = k, {n, 3, nn}];
    Array[a, nn] (* Michael De Vlieger, Mar 24 2024 *)
  • PARI
    up_to = 16384;
    LCMtransform(v) = { my(len = length(v), b = vector(len), g = vector(len)); b[1] = g[1] = 1; for(n=2,len, g[n] = lcm(g[n-1],v[n]); b[n] = g[n]/g[n-1]); (b); };
    A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); (t) };
    v368900 = LCMtransform(vector(up_to,i,A005940(i)));
    A368900(n) = v368900[n];
    
  • PARI
    A000265(n) = (n>>valuation(n,2));
    A209229(n) = (n && !bitand(n,n-1));
    A368900(n)  = if(1==n, 1, my(x=A000265(n-1)); if(A209229(1+x), prime(1+valuation(n-1,2)), 1));

Formula

a(n) = A368901(n) / A368901(n-1) = lcm {1..A005940(n)} / lcm {1..A005940(n-1)}.
a(n) = A005940(n) / gcd(A005940(n), A368901(n-1)).
a(n) = A014963(A005940(n)). [Because A005940 satisfies the property given in the comments]
For n >= 1, Product_{d|n} a(A005941(d)) = n. [Implied by above]
For n >= 1, a(n) = A369030(1+A054429(n-1)).
For n > 1, if n-1 is a number of the form 2^i - 2^j with i >= j, then a(n) = prime(1+j), otherwise a(n) = 1.

A246375 Permutation of natural numbers: a(1) = 1, a(2n) = 2*a(n), a(2n+1) = A003961(1+a(n)). [Where A003961(n) shifts the prime factorization of n one step towards larger primes].

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 9, 8, 7, 10, 15, 12, 11, 18, 21, 16, 25, 14, 27, 20, 13, 30, 81, 24, 17, 22, 45, 36, 23, 42, 39, 32, 19, 50, 51, 28, 35, 54, 99, 40, 55, 26, 33, 60, 37, 162, 129, 48, 49, 34, 75, 44, 29, 90, 87, 72, 41, 46, 135, 84, 47, 78, 189, 64, 65, 38, 63, 100, 95, 102, 153, 56, 31, 70
Offset: 1

Views

Author

Antti Karttunen, Aug 27 2014

Keywords

Comments

This can be viewed as yet another "entanglement permutation" where the two complementary pairs to be interwoven together are even and odd numbers (A005843/A005408) which are entangled with the complementary pair even numbers (taken straight) and odd numbers in the order they appear in A003961: (A005843/A003961). Sequence A163511 has almost the same definition, but its domain starts from 0, which results a different permutation.

Crossrefs

Inverse: A246376.
Similar or related permutations: A005940, A005941, A163511, A245606, A246378, A246379.

Programs

  • PARI
    default(primelimit, (2^31)+(2^30));
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Using code of Michel Marcus
    A246375(n) = if(1==n, 1, if(!(n%2), 2*A246375(n/2), A003961(1+A246375((n-1)/2))));
    for(n=1, 16384, write("b246375.txt", n, " ", A246375(n)));
    (Scheme, with memoizing definec-macro)
    (definec (A246375 n) (cond ((<= n 1) n) ((even? n) (* 2 (A246375 (/ n 2)))) (else (A003961 (+ 1 (A246375 (/ (- n 1) 2)))))))

Formula

a(1) = 1, a(2n) = 2*a(n), a(2n+1) = A003961(1+a(n)). [Where A003961(n) shifts the prime factorization of n one step towards larger primes].
As a composition of related permutations:
a(n) = A246379(A246378(n)).
Other identities. For all n >= 1 the following holds:
A000035(a(n)) = A000035(n). [Like A005940 & A005941, this also preserves the parity].

A246376 Permutation of natural numbers: a(1) = 1, a(2n) = 2 * a(n), a(2n-1) = 1 + (2 * a(A064989(2n-1)-1)).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 9, 8, 7, 10, 13, 12, 21, 18, 11, 16, 25, 14, 33, 20, 15, 26, 29, 24, 17, 42, 19, 36, 53, 22, 73, 32, 43, 50, 37, 28, 45, 66, 31, 40, 57, 30, 81, 52, 27, 58, 61, 48, 49, 34, 35, 84, 117, 38, 41, 72, 87, 106, 169, 44, 213, 146, 67, 64, 65, 86, 89, 100, 91, 74, 173, 56, 149, 90, 51, 132, 101, 62, 113, 80, 23
Offset: 1

Views

Author

Antti Karttunen, Aug 27 2014

Keywords

Crossrefs

Inverse: A246375.
Similar or related permutations: A005940, A005941, A064216, A243071, A245605, A246377, A246380.

Formula

a(1) = 1, a(2n) = 2 * a(n), a(2n-1) = 1 + (2 * a(A064989(2n-1)-1)).
As a composition of related permutations:
a(n) = A246377(A246380(n)).
Other identities. For all n >= 1 the following holds:
A000035(a(n)) = A000035(n). [Like A005940 & A005941, this also preserves the parity].

A364560 Numbers k for which A156552(k) < k.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 9, 10, 12, 15, 16, 18, 20, 21, 24, 25, 27, 30, 32, 35, 36, 40, 42, 45, 48, 49, 50, 54, 55, 60, 63, 64, 70, 72, 75, 77, 80, 81, 84, 90, 91, 96, 98, 99, 100, 105, 108, 110, 120, 121, 125, 126, 128, 135, 140, 143, 144, 147, 150, 154, 160, 162, 165, 168, 169, 175, 180, 182, 187, 189, 192, 195, 196
Offset: 1

Views

Author

Antti Karttunen, Jul 28 2023

Keywords

Comments

Numbers k such that A005941(k) <= k.
Sequence A005940(A364542(.)) sorted into ascending order.
If k is a term, then also 2*k is present in this sequence, and vice versa.

Crossrefs

Positions of nonpositive terms in A364559.
Cf. A005941, A156552, A364542, A364562 (complement).
Subsequences: A029747, A364550, A364561 (odd terms).

Programs

  • PARI
    A156552(n) = { my(f = factor(n), p, p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res };
    isA364560(n) = (A156552(n) < n);

A253551 Square array: A(row,col) = 2^(row-1) * 1+(2*A156552(col)) = A156552(A246278(row,col)), read by antidiagonals A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ...

Original entry on oeis.org

1, 3, 2, 5, 6, 4, 7, 10, 12, 8, 9, 14, 20, 24, 16, 11, 18, 28, 40, 48, 32, 17, 22, 36, 56, 80, 96, 64, 15, 34, 44, 72, 112, 160, 192, 128, 13, 30, 68, 88, 144, 224, 320, 384, 256, 19, 26, 60, 136, 176, 288, 448, 640, 768, 512, 33, 38, 52, 120, 272, 352, 576, 896, 1280, 1536, 1024, 23, 66, 76, 104, 240, 544, 704, 1152, 1792, 2560, 3072, 2048, 65, 46
Offset: 1

Views

Author

Antti Karttunen, Jan 03 2015

Keywords

Comments

Shares with A135764 the property that A001511(n) = k for all terms n on row k and when going downwards in each column, terms grow by doubling.

Examples

			The top left corner of the array:
   1,  3,  5,   7,   9, 11,  17,  15,  13,  19,  33,  23,  65,  35,  21,
   2,  6, 10,  14,  18, 22,  34,  30,  26,  38,  66,  46, 130,  70,  42,
   4, 12, 20,  28,  36, 44,  68,  60,  52,  76, 132,  92, 260, 140,  84,
   8, 24, 40,  56,  72, 88, 136, 120, 104, 152, 264, 184, 520, 280, 168,
  16, 48, 80, 112, 144,176, 272, 240, 208, 304, 528, 368,1040, 560, 336,
...
		

Crossrefs

Inverse: A253552.
Differs from A135764 for the first time at n=22, where a(22) = 17, while A135764(22) = 13.

Formula

A(row,col) = A156552(A246278(row,col)).
A(row,col) = A135764(row, A005941(col)). [Is otherwise the same array as A135764, but the column positions have been permuted by A005941.]
A(row,col) = 2^(row-1) * ((2*A005941(col)) - 1) = 2^(row-1) * A005408(A156552(col)). [The above expands to this.]
As a composition of other permutations:
a(n) = A156552(A246278(n+1)). [When all three sequences are interpreted as one-dimensional sequences.]

A364542 Numbers k for which A005940(k) >= k.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 16, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 71, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95
Offset: 1

Views

Author

Antti Karttunen, Jul 28 2023

Keywords

Comments

Sequence A005941(A364560(.)) sorted into ascending order.
A029747 is included as a subsequence, because it gives the known fixed points of map n -> A005940(n).
Differs from A343107 for the first time at a(22) = 25, which term is not present in A343107. On the other hand, 35 is the first term of A343107 that is not present in this sequence.

Crossrefs

Positions of nonnegative terms in A364499.
Complement of A364540.
Cf. A005940, A005941, A029747 (subsequence), A343107 (not a subsequence), A364560.

Programs

  • Mathematica
    nn = 95; Array[Set[a[#], #] &, 2]; Do[If[EvenQ[n], Set[a[n], 2 a[n/2]], Set[a[n], Times @@ Power @@@ Map[{Prime[PrimePi[#1] + 1], #2} & @@ # &, FactorInteger[a[(n + 1)/2]]]]], {n, 3, nn}]; Select[Range[nn], a[#] >= # &] (* Michael De Vlieger, Jul 28 2023 *)
  • PARI
    A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
    isA364542(n) = (A005940(n)>=n);

A364561 Odd numbers k for which A156552(k) < k.

Original entry on oeis.org

1, 3, 5, 9, 15, 21, 25, 27, 35, 45, 49, 55, 63, 75, 77, 81, 91, 99, 105, 121, 125, 135, 143, 147, 165, 169, 175, 187, 189, 195, 221, 225, 231, 243, 245, 273, 275, 289, 297, 315, 323, 325, 343, 351, 357, 363, 375, 385, 405, 425, 429, 441, 455, 495, 507, 525, 539, 561, 567, 585, 595, 605, 625, 627, 637, 663, 665
Offset: 1

Views

Author

Antti Karttunen, Jul 28 2023

Keywords

Comments

Odd numbers k such that A005941(k) <= k.

Crossrefs

Odd terms in A364560.
Cf. A005940, A005941, A156552, A364545, A364564 (largest prime factor).
Cf. also A364551, A364576 (subsequences).

Programs

  • PARI
    A156552(n) = { my(f = factor(n), p, p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res };
    isA364561(n) = ((n%2)&&(A156552(n) < n));

A364576 Starting from k=1, each subsequent term is the next larger odd k such that A156552(k) < k and the ratio A156552(k)/k is nearer to 1.0 than for any previous k in the sequence.

Original entry on oeis.org

1, 3, 5, 21, 323, 66297, 139965, 263375, 264845, 528581
Offset: 1

Views

Author

Antti Karttunen, Aug 06 2023

Keywords

Comments

All the odd fixed points of map n -> A005940(n) [and its inverse, map n -> A005941(n)] are included in this sequence. This includes both the known odd fixed points, 1, 3 and 5 (see A029747), and any additional hypothetical odd composites that would satisfy the condition n == A005940(n).
This is a subsequence of A364561, so the comments given in A364564 apply also here.

Examples

			       k  A156552(k)    A156552(k)/k  k-(1+A156552(k)) factorization of k
       1:       0         0                0
       3:       2         0.6666667        0
       5:       4         0.8              0
      21:      18         0.8571429        2           (3 * 7)
     323:     320         0.9907121        2           (17 * 19)
   66297:   65714         0.9912062      582           (3 * 7^2 * 11 * 41)
  139965:  139306         0.9952917      658           (3 * 5 * 7 * 31 * 43)
  263375:  262364         0.9961614     1010           (5^3 * 7^2 * 43)
  264845:  264244         0.9977307      600           (5 * 7^2 * 23 * 47)
  528581:  528576         0.9999905        4           (17^2 * 31 * 59).
		

Crossrefs

Subsequence of A364561.
Cf. also A364551, A364564, A364572.

A364546 Numbers k such that k is a multiple of A005940(k).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 10, 12, 16, 20, 24, 32, 40, 48, 64, 80, 96, 128, 160, 192, 256, 320, 384, 512, 640, 768, 1024, 1035, 1280, 1536, 2048, 2070, 2560, 3072, 4096, 4140, 5120, 6144, 8192, 8280, 10240, 12288, 16384, 16560, 20480, 24576, 32768, 33120, 40960, 49152, 65536, 66240, 81920, 98304, 131072, 132480, 163840
Offset: 1

Views

Author

Antti Karttunen, Jul 28 2023

Keywords

Comments

Sequence A005941(A364548(.)) sorted into ascending order.
If k is a term, then also 2*k is present in this sequence, and vice versa.
A029747 is included as a subsequence, because it gives the known fixed points of map n -> A005940(n).

Crossrefs

Positions of 1's in A364502.
Subsequence of A364541.
Subsequences: A029747, A364547 (odd terms).
Cf. also A364496.

Programs

  • Mathematica
    nn = 2^18; Array[Set[a[#], #] &, 2]; Do[If[EvenQ[n], Set[a[n], 2 a[n/2]], Set[a[n], Times @@ Power @@@ Map[{Prime[PrimePi[#1] + 1], #2} & @@ # &, FactorInteger[a[(n + 1)/2]]]]], {n, 3, nn}]; Select[Range[nn], Divisible[#, a[#]] &] (* Michael De Vlieger, Jul 28 2023 *)
  • PARI
    A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
    isA364546(n) = !(n%A005940(n));
Previous Showing 21-30 of 33 results. Next