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.

Showing 1-3 of 3 results.

A364297 a(n) = A348717(A163511(n)).

Original entry on oeis.org

1, 2, 4, 2, 8, 4, 6, 2, 16, 8, 18, 4, 12, 6, 10, 2, 32, 16, 54, 8, 36, 18, 50, 4, 24, 12, 30, 6, 20, 10, 14, 2, 64, 32, 162, 16, 108, 54, 250, 8, 72, 36, 150, 18, 100, 50, 98, 4, 48, 24, 90, 12, 60, 30, 70, 6, 40, 20, 42, 10, 28, 14, 22, 2, 128, 64, 486, 32, 324, 162, 1250, 16, 216, 108, 750, 54, 500, 250, 686, 8, 144
Offset: 0

Views

Author

Antti Karttunen, Aug 15 2023

Keywords

Comments

For all i, j: a(i) = a(j) => A278531(i) = A278531(j).
As the underlying sequence A163511 can be represented as a binary tree, so can this be also:
1
|
...................2...................
4 2
8......../ \........4 6......../ \........2
/ \ / \ / \ / \
/ \ / \ / \ / \
/ \ / \ / \ / \
16 8 18 4 12 6 10 2
32 16 54 8 36 18 50 4 24 12 30 6 20 10 14 2
etc.
Each rightward leaning branch stays constant, because a(2n+1) = a(n).
Conjecture: Mersenne primes (A000668) gives all such odd numbers k for which a(k) = A348717(k). If true, then it immediately implies that map n -> A163511(n) [or equally: map n -> A243071(n)] has no other fixed points than those given by A007283. But see also A364959. - Edited Sep 03 2023

Crossrefs

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));
    A348717(n) = if(1==n, 1, my(f = factor(n), k = primepi(f[1, 1])-1); for (i=1, #f~, f[i, 1] = prime(primepi(f[i, 1])-k)); factorback(f));
    A364297(n) = A348717(A163511(n));

Formula

a(0) = 1, a(1) = 2, a(2n) = A163511(2n) = 2*A163511(n), and for n > 0, a(2n+1) = a(n).

A365422 Numbers k for which k and A163511(k) have the same prime signature.

Original entry on oeis.org

3, 6, 7, 12, 14, 24, 28, 31, 45, 48, 55, 56, 62, 90, 96, 110, 111, 112, 119, 123, 124, 127, 175, 180, 192, 207, 220, 222, 224, 238, 246, 247, 248, 253, 254, 350, 360, 384, 414, 440, 444, 447, 448, 476, 492, 494, 496, 506, 508, 567, 700, 720, 768, 828, 880, 888, 894, 895, 896, 927, 945, 952, 957, 959, 984, 987, 988, 992
Offset: 1

Views

Author

Antti Karttunen, Sep 04 2023

Keywords

Comments

If k is present, then 2*k is also a term, and vice versa.

Crossrefs

Cf. A046523, A163511, A278531, A365421 (characteristic function).
Subsequences: A007283, A335431, A365423 (odd terms).

Programs

A374478 Lexicographically earliest infinite sequence such that a(i) = a(j) => A348717(i) = A348717(j) and A364255(i) = A364255(j), for all i, j >= 1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 5, 7, 8, 9, 5, 10, 5, 11, 12, 13, 5, 14, 5, 15, 16, 17, 5, 18, 19, 20, 21, 22, 5, 23, 5, 24, 25, 26, 27, 28, 5, 29, 30, 31, 5, 32, 5, 33, 34, 35, 5, 36, 37, 38, 39, 40, 5, 41, 42, 43, 44, 45, 5, 46, 5, 47, 48, 49, 50, 51, 5, 52, 53, 54, 5, 55, 5, 56, 57, 58, 59, 60, 5, 61, 62, 63, 5, 64, 65, 66, 67, 68, 5, 69, 70, 71, 72, 73, 39, 74, 5, 75, 76, 77, 5
Offset: 1

Views

Author

Antti Karttunen, Aug 07 2024

Keywords

Comments

Restricted growth sequence transform of the ordered pair [A348717(n), A364255(n)].
For all i, j >= 1:
A305900(i) = A305900(j) => a(i) = a(j),
a(i) = a(j) => A305891(i) = A305891(j),
a(i) = a(j) => A374477(i) = A374477(j).

Crossrefs

Differs from A374040 first at n=77, where a(77) = 59, while A374040(77) = 50.
Differs from A305900 first at n=95, where a(95) = 39, while A305900(95) = 74.

Programs

  • PARI
    up_to = 65537;
    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; };
    A348717(n) = if(1==n, 1, my(f = factor(n), k = primepi(f[1, 1])-1); for (i=1, #f~, f[i, 1] = prime(primepi(f[i, 1])-k)); factorback(f));
    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));
    A364255(n) = gcd(n, A163511(n));
    Aux374478(n) = [A348717(n), A364255(n)];
    v374478 = rgs_transform(vector(up_to, n, Aux374478(n)));
    A374478(n) = v374478[n];
Showing 1-3 of 3 results.