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 11-16 of 16 results.

A332215 Mersenne-prime fixing variant of A243071: a(n) = A243071(A332213(n)).

Original entry on oeis.org

0, 1, 3, 2, 15, 6, 7, 4, 5, 30, 63, 12, 255, 14, 29, 8, 511, 10, 1023, 60, 13, 126, 2047, 24, 23, 510, 9, 28, 4095, 58, 31, 16, 125, 1022, 27, 20, 16383, 2046, 509, 120, 32767, 26, 65535, 252, 57, 4094, 262143, 48, 11, 46, 1021, 1020, 1048575, 18, 119, 56, 2045, 8190, 2097151, 116, 4194303, 62, 25, 32, 503, 250, 8388607, 2044, 4093, 54, 16777215, 40
Offset: 1

Views

Author

Antti Karttunen, Feb 09 2020

Keywords

Comments

Any Mersenne prime (A000668) times any power of 2 (i.e., 2^k, for k>=0) is fixed by this sequence, including also all even perfect numbers.
From Antti Karttunen, Jul 10 2020: (Start)
This is a "tuned variant" of A243071, and has many of the same properties.
For example, for n > 1, A007814(a(n)) = A007814(n) - A209229(n), that is, this map preserves the 2-adic valuation of n, except when n is a power of two, in which cases that value is decremented by one, and in particular, a(2^k * n) = 2^k * a(n) for all n > 1. Also, like A243071, this bijection maps primes to the terms of A000225 (binary repunits). However, the "tuning" (A332213) has a specific effect that each Mersenne prime (A000668) is mapped to itself. Therefore the terms of A335431 are fixed by this map. Furthermore, I conjecture that there are no other fixed points. For the starters, see the proof in A335879, which shows that at least none of the terms of A335882 are fixed.
(End)

Crossrefs

Cf. A243071, A332210, A332213, A332214 (inverse permutation), A335431 (conjectured to be all the fixed points), A335879.

Programs

Formula

a(n) = A243071(A332213(n)).
For all n >= 1, a(A335431(n)) = A335431(n), a(A335882(n)) = A335879(n). - Antti Karttunen, Jul 10 2020

A332211 Lexicographically earliest permutation of primes such that a(n) = 2^n - 1 when n is one of the Mersenne prime exponents (in A000043).

Original entry on oeis.org

2, 3, 7, 5, 31, 11, 127, 13, 17, 19, 23, 29, 8191, 37, 41, 43, 131071, 47, 524287, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 2147483647, 103, 107, 109, 113, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 2305843009213693951, 269, 271, 277, 281, 283, 293, 307, 311, 313
Offset: 1

Views

Author

Antti Karttunen, Feb 09 2020

Keywords

Comments

Sequence is well-defined also in case there are only a finite number of Mersenne primes.

Examples

			For p in A000043: 2, 3, 5, 7, 13, 17, 19, ..., a(p) = (2^p)-1, thus a(2) = 2^2 - 1 = 3, a(3) = 7, a(5) = 31, a(7) = 127, a(13) = 8191, a(17) = 131071, etc., with the rest of positions filled by the least unused prime:
1, 2, 3, 4,  5,  6,   7,  8,  9, 10, 11, 12,   13, 14, 15, 16, 17, ...
2, 3, 7, 5, 31, 11, 127, 13, 17, 19, 23, 29, 8191, 37, 41, 43, 131071, ...
		

Crossrefs

Cf. A000040, A000043, A000668, A332210 (inverse permutation of primes), A332220.
Used to construct permutations A332212, A332214.

Programs

  • PARI
    up_to = 127;
    A332211list(up_to) = { my(v=vector(up_to), xs=Map(), i=1, q); for(n=1,up_to, if(isprime(q=((2^n)-1)), v[n] = q, while(mapisdefined(xs,prime(i)), i++); v[n] = prime(i)); mapput(xs,v[n],n)); (v); };
    v332211 = A332211list(up_to);
    A332211(n) = v332211[n];
    \\ For faster computing of larger values, use precomputed values of A000043:
    v000043 = [2,3,5,7,13,17,19,31,61,89,107,127,521,607,1279,2203,2281,3217];
    up_to = v000043[#v000043];
    A332211list(up_to) = { my(v=vector(up_to), xs=Map(), i=1, q); for(n=1,up_to, if(vecsearch(v000043,n), q = (2^n)-1, while(mapisdefined(xs,prime(i)), i++); q = prime(i)); v[n] = q; mapput(xs,q,n)); (v); };

Formula

For all applicable n >= 1, a(A000043(n)) = A000668(n).

A332212 Fully multiplicative with a(p) = A332211(A000720(p)).

Original entry on oeis.org

1, 2, 3, 4, 7, 6, 5, 8, 9, 14, 31, 12, 11, 10, 21, 16, 127, 18, 13, 28, 15, 62, 17, 24, 49, 22, 27, 20, 19, 42, 23, 32, 93, 254, 35, 36, 29, 26, 33, 56, 8191, 30, 37, 124, 63, 34, 41, 48, 25, 98, 381, 44, 43, 54, 217, 40, 39, 38, 131071, 84, 47, 46, 45, 64, 77, 186, 524287, 508, 51, 70, 53, 72, 59, 58, 147, 52, 155, 66, 61, 112, 81, 16382, 67, 60
Offset: 1

Views

Author

Antti Karttunen, Feb 09 2020

Keywords

Crossrefs

Cf. A000043, A000668, A000720, A332211, A332213 (inverse permutation), A332214.

Programs

  • PARI
    \\ Needs also code from A332211:
    A332212(n) = { my(f=factor(n)); f[,1] = apply(A332211,apply(primepi,f[,1])); factorback(f); };

Formula

a(1) = 1, a(p^e) = A332211(A000720(p))^e, a(m*n) = a(m)*a(n).

A365802 Numbers k such that A163511(k) is a fifth power.

Original entry on oeis.org

0, 16, 33, 67, 135, 271, 512, 543, 1025, 1056, 1087, 2051, 2113, 2144, 2175, 4103, 4227, 4289, 4320, 4351, 8207, 8455, 8579, 8641, 8672, 8703, 16384, 16415, 16911, 17159, 17283, 17345, 17376, 17407, 32769, 32800, 32831, 33792, 33823, 34319, 34567, 34691, 34753, 34784, 34815, 65539, 65601, 65632, 65663, 67585, 67616
Offset: 1

Views

Author

Antti Karttunen, Oct 01 2023

Keywords

Comments

Equivalently, numbers k for which A332214(k), and also A332817(k) are fifth powers.
The sequence is defined inductively as:
(a) it contains 0 and 16, and
(b) for any nonzero term a(n), (2*a(n)) + 1 and 32*a(n) are also included as terms.
When iterating n -> 2n+1 mod 31, starting from 16 we obtain five distinct remainders 16, 2, 5, 11, 23, before the cycle starts again from 16. (see A153893), while x^5 mod 31 may obtain only these values: 0, 1, 5, 6, 25, 26, 30. The only common element of these sets is 5. We have x^5 == 5 (mod 31) whenever x == 7, 14, 19, 25, 28 mod 31, with all other x leaving a remainder that is not in the set [16, 2, 5, 11, 23].
On the other hand, when iterating n -> 2n+1 mod 33, starting from 16 we obtain ten distinct remainders 16, 0, 1, 3, 7, 15, 31, 30, 28, 24, before the cycle starts again from 16, while x^5 mod 33 obtain only these values: 0, 1, 10, 11, 12, 21, 22, 23, 32. We have x^5 == 0 (mod 33) iff x == 0 (mod 33) and x^5 == 1 (mod 33) whenever x == 1, 4, 16, 25, 31 mod 33. In the n->2n+1 cycles of 5 and 10 elements starting from 16, the 5's (of every second cycle) in the former and the 1's in the latter are aligned with each other.
In any case, this sequence do not contain any fifth powers after the initial zero. See A365805. - Antti Karttunen, Nov 23 2023

Crossrefs

Positions of multiples of 5 in A365805.
Sequence A243071(n^5), n >= 1, sorted into ascending order.
Subsequences: A013825, A198275.

Programs

  • PARI
    A163511(n) = if(!n, 1, my(p=2, t=1); while(n>1, if(!(n%2), (t*=p), p=nextprime(1+p)); n >>= 1); (t*p));
    isA365802(n) = ispower(A163511(n),5);
    
  • PARI
    isA365802(n) = if(n<=16, !(n%16), if(n%2, isA365802((n-1)/2), if(n%32, 0, isA365802(n/32))));

A366287 Numbers k such that A163511(k) is a seventh power.

Original entry on oeis.org

0, 64, 129, 259, 519, 1039, 2079, 4159, 8192, 8319, 16385, 16512, 16639, 32771, 33025, 33152, 33279, 65543, 66051, 66305, 66432, 66559, 131087, 132103, 132611, 132865, 132992, 133119, 262175, 264207, 265223, 265731, 265985, 266112, 266239, 524351, 528415, 530447, 531463, 531971, 532225, 532352, 532479, 1048576, 1048703
Offset: 1

Views

Author

Antti Karttunen, Oct 09 2023

Keywords

Comments

Equivalently, numbers k for which A332214(k), and also A332817(k) are seventh powers.
The sequence is defined inductively as:
(a) it contains 0 and 64,
and
(b) for any nonzero term a(n), (2*a(n)) + 1 and 128*a(n) are also included as terms.
When iterating n -> 2n+1 mod 127, starting from 64 we get 64, 2, 5, 11, 23, 47, 95, and then cycle starts again from 64 (see A153893), while on the other hand, x^7 mod 127 obtains values: 0, 1, 19, 20, 22, 24, 28, 37, 52, 59, 68, 75, 90, 99, 103, 105, 107, 108, 126. These sets have no terms in common, therefore there are no seventh powers in this sequence after the initial 0.

Crossrefs

Positions of multiples of 7 in A365805.
Sequence A243071(n^7), n >= 1, sorted into ascending order.

Programs

  • PARI
    A163511(n) = if(!n, 1, my(p=2, t=1); while(n>1, if(!(n%2), (t*=p), p=nextprime(1+p)); n >>= 1); (t*p));
    isA366287(n) = ispower(A163511(n),7);
    
  • PARI
    isA366287(n) = if(n<=64, !(n%64), if(n%2, isA366287((n-1)/2), if(n%128, 0, isA366287(n>>7))));

A366391 Numbers k such that A163511(k) is an eleventh power.

Original entry on oeis.org

0, 1024, 2049, 4099, 8199, 16399, 32799, 65599, 131199, 262399, 524799, 1049599, 2097152, 2099199, 4194305, 4196352, 4198399, 8388611, 8392705, 8394752, 8396799, 16777223, 16785411, 16789505, 16791552, 16793599, 33554447, 33570823, 33579011, 33583105, 33585152, 33587199, 67108895, 67141647, 67158023, 67166211, 67170305
Offset: 1

Views

Author

Antti Karttunen, Oct 09 2023

Keywords

Comments

Equivalently, numbers k for which A332214(k), and also A332817(k) are eleventh powers.
The sequence is defined inductively as:
(a) it contains 0 and 1024,
and
(b) for any nonzero term a(n), (2*a(n)) + 1 and 2048*a(n) are also included as terms.
When iterating n -> 2n+1 mod 2047, starting from 1024 we get 1024, 2, 5, 11, 23, 47, 95, 191, 383, 767, 1535, and then cycle starts again from 1024 (see A153893), while on the other hand, x^11 mod 2047 obtains values: 0, 1, 230, 322, 344, 368, 390, 482, 622, 712, 713, 942, 967, 1013, 1034, 1080, 1105, 1334, 1335, 1425, 1565, 1657, 1679, 1703, 1725, 1817, 2046. These sets have no terms in common, therefore there are no eleventh powers in this sequence after the initial 0.

Crossrefs

Positions of multiples of 11 in A365805.
Sequence A243071(n^11), n >= 1, sorted into ascending order.

Programs

  • PARI
    A163511(n) = if(!n, 1, my(p=2, t=1); while(n>1, if(!(n%2), (t*=p), p=nextprime(1+p)); n >>= 1); (t*p));
    isA366391v(n) = ispower(A163511(n),11);
    
  • PARI
    isA366391(n) = if(n<=1024, !(n%1024), if(n%2, isA366391((n-1)/2), if(n%2048, 0, isA366391(n>>11))));
Previous Showing 11-16 of 16 results.