A319347 -id:A319347 - cp's OEIS frontend

A305801 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = 0 if n is an odd prime, with f(n) = n for all other n.

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

Offset: 1

Antti Karttunen, Jun 14 2018

nonn

The original name was: "Filter sequence for a(odd prime) = constant sequences", which stemmed from the fact that for all i, j, a(i) = a(j) => b(i) = b(j) for any sequence b that obtains a constant value for all odd primes A065091.
For example, we have for all i, j:
  a(i) = a(j) => A305800(i) = A305800(j),
  a(i) = a(j) => A007814(i) = A007814(j),
  a(i) = a(j) => A305891(i) = A305891(j) => A291761(i) = A291761(j).
There are several filter sequences "above" this one (meaning that they have finer equivalence class partitioning), for example, we have, for all i, j:
[where odd primes are further distinguished by]
  A305900(i) = A305900(j) => a(i) = a(j),   [whether p = 3 or > 3]
  A319350(i) = A319350(j) => a(i) = a(j),   [A007733(p)]
  A319704(i) = A319704(j) => a(i) = a(j),   [p mod 4]
  A319705(i) = A319705(j) => a(i) = a(j),   [A286622(p)]
  A331304(i) = A331304(j) => a(i) = a(j),   [parity of A000720(p)]
  A336855(i) = A336855(j) => a(i) = a(j).   [distance to the next larger prime]

Antti Karttunen, Table of n, a(n) for n = 1..100000

Cf. A000720, A065091, A305800, A305900.
Cf. A305900, A319350, A319704, A319705, A331304, A336855 (sequences with finer equivalence class partitioning).
Cf. A305800, A305890, A305891, A305896, A318500, A318888, A319346, A319347, A319349, A319701, A322591, A322809, A322810, A323078, A323367, A323082, A323369, A323370, A323371, A323374, A323400, A324401, A326199, A326201, A326203, A326203, A328470, A329608, A331174, A331730, A331301 (sequences with coarser equivalence class partitioning).
Cf. also A003602, A103391, A295300, A305795, A324400, A331300, A336460 (for similar constructions or similarly useful sequences).

Array[If[# <= 2, #, If[PrimeQ[#], 3, 2 + # - PrimePi[#]]] &, 105] (* Michael De Vlieger, Oct 18 2021 *)

A305801(n) = if(n<=2,n,if(isprime(n),3,2+n-primepi(n)));

a(1) = 1, a(2) = 2; for n > 2, a(n) = 3 for odd primes, and a(n) = 2+n-A000720(n) for composite n.

For n > 2, a(n) = 1 + A305800(n).

Name changed and Comment section rewritten by Antti Karttunen, Oct 17 2021

A294877 Lexicographically earliest such sequence a that a(i) = a(j) => A003557(i) = A003557(j) and A046523(i) = A046523(j), for all i, j.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 5, 6, 4, 2, 7, 2, 4, 4, 8, 2, 9, 2, 7, 4, 4, 2, 10, 11, 4, 12, 7, 2, 13, 2, 14, 4, 4, 4, 15, 2, 4, 4, 10, 2, 13, 2, 7, 9, 4, 2, 16, 17, 18, 4, 7, 2, 19, 4, 10, 4, 4, 2, 20, 2, 4, 9, 21, 4, 13, 2, 7, 4, 13, 2, 22, 2, 4, 18, 7, 4, 13, 2, 16, 23, 4, 2, 20, 4, 4, 4, 10, 2, 24, 4, 7, 4, 4, 4, 25, 2, 26, 9, 27, 2, 13, 2, 10, 13, 4, 2, 28, 2, 13

Offset: 1

Antti Karttunen, Nov 11 2017

nonn

Restricted growth sequence transform of A291757, which means that this is the lexicographically least sequence a, such that for all i, j: a(i) = a(j) <=> A291757(i) = A291757(j) <=> A003557(i) = A003557(j) and A046523(i) = A046523(j). That this is equal to the definition given in the title follows because any such lexicographically least sequence satisfying relation <=> is also the least sequence satisfying relation => with the same parameters.
Also the restricted growth sequence transform of A294876, Product_{d|n, d>1} prime(gcd(d,n/d)). (This was the original definition).
For all i, j:
  A295300(i) = A295300(j) => a(i) = a(j),
  A319347(i) = A319347(j) => a(i) = a(j),
  a(i) = a(j) => A055155(i) = A055155(j).

Antti Karttunen, Table of n, a(n) for n = 1..65537

Cf. also A291751, A291757, A294876, A295300, A295888, A319347, A322021.

Cf. A000188, A055155, A294897, A295666, A322020 (a few of the matched sequences).

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; };
A294876(n) = { my(m=1); fordiv(n,d,if(d>1, m *= prime(gcd(d,n/d)))); m; };
v294877 = rgs_transform(vector(up_to,n,A294876(n)));
A294877(n) = v294877[n];

A003557(n) = n/factorback(factor(n)[, 1]); \\ From A003557
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
v294877 = rgs_transform(vector(up_to,n,[A003557(n),A046523(n)]));
A294877(n) = v294877[n]; \\ Antti Karttunen, Nov 28 2018

Name changed and comments added by Antti Karttunen, Nov 28 2018

A291757 a(n) = (1/2)(2 + ((A003557(n)+A046523(n))^2) - A003557(n) - 3A046523(n)).

Original entry on oeis.org

1, 2, 2, 12, 2, 16, 2, 59, 18, 16, 2, 80, 2, 16, 16, 261, 2, 94, 2, 80, 16, 16, 2, 355, 33, 16, 129, 80, 2, 436, 2, 1097, 16, 16, 16, 826, 2, 16, 16, 355, 2, 436, 2, 80, 94, 16, 2, 1493, 52, 125, 16, 80, 2, 505, 16, 355, 16, 16, 2, 1832, 2, 16, 94, 4497, 16, 436, 2, 80, 16, 436, 2, 3415, 2, 16, 125, 80, 16, 436, 2, 1493, 888, 16, 2, 1832, 16, 16, 16, 355, 2

Offset: 1

Antti Karttunen, Sep 10 2017

nonn

Antti Karttunen, Table of n, a(n) for n = 1..16385

Cf. A000027, A003557, A046523, A291750, A291752, A291756, A291758, A294877 (rgs-transform), A319347.

A003557(n) = n/factorback(factor(n)[, 1]); \\ From A003557
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
A291757(n) = (1/2)*(2 + ((A003557(n)+A046523(n))^2) - A003557(n) - 3*A046523(n));

a(n) = (1/2)*(2 + ((A003557(n)+A046523(n))^2) - A003557(n) - 3*A046523(n)).

Name changed by Antti Karttunen, Nov 28 2018

cp's OEIS Frontend

A305801 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = 0 if n is an odd prime, with f(n) = n for all other n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A294877 Lexicographically earliest such sequence a that a(i) = a(j) => A003557(i) = A003557(j) and A046523(i) = A046523(j), for all i, j.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Extensions

A291757 a(n) = (1/2)(2 + ((A003557(n)+A046523(n))^2) - A003557(n) - 3A046523(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

Extensions

cp's OEIS Frontend

A305801 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = 0 if n is an odd prime, with f(n) = n for all other n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A294877 Lexicographically earliest such sequence a that a(i) = a(j) => A003557(i) = A003557(j) and A046523(i) = A046523(j), for all i, j.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Extensions

A291757 a(n) = (1/2)*(2 + ((A003557(n)+A046523(n))^2) - A003557(n) - 3*A046523(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

Extensions

A291757 a(n) = (1/2)(2 + ((A003557(n)+A046523(n))^2) - A003557(n) - 3A046523(n)).