A323082 -id:A323082 - 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

A322809 Lexicographically earliest such sequence a that a(i) = a(j) => f(i) = f(j) for all i, j, where f(n) = -1 if n is an odd prime, and f(n) = floor(n/2) for all other numbers.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Dec 26 2018

nonn

This sequence is a restricted growth sequence transform of a function f which is defined as f(n) = A004526(n), unless n is an odd prime, in which case f(n) = -1, which is a constant not in range of A004526. See the Crossrefs section for a list of similar sequences.
For all i, j:
  A305801(i) = A305801(j) => a(i) = a(j),
  a(i) = a(j) => A039636(i) = A039636(j).
For all i, j: a(i) = a(j) <=> A323161(i+1) = A323161(j+1).
The shifted version of this filter, A323161, has a remarkable ability to find many sequences related to primes and prime chains. - Antti Karttunen, Jan 06 2019

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences related to binary expansion of n

Cf. A004526, A039636, A305801, A323161.
A list of few similarly constructed sequences follows, where each sequence is an rgs-transform of such function f, for which the value of f(n) is the n-th term of the sequence whose A-number follows after a parenthesis, unless n is of the form ..., in which case f(n) is given a constant value outside of the range of that sequence:
  A322809 (A004526, unless an odd prime) [This sequence],
  A322589 (A007913, unless an odd prime),
  A322591 (A007947, unless an odd prime),
  A322805 (A252463, unless an odd prime),
  A323082 (A300840, unless an odd prime),
  A322822 (A300840, unless n > 2 and a Fermi-Dirac prime, A050376),
  A322988 (A322990, unless a prime power > 2),
  A323078 (A097246, unless an odd prime),
  A322808 (A097246, unless a squarefree number > 2),
  A322816 (A048675, unless an odd prime),
  A322807 (A285330, unless an odd prime),
  A322814 (A286621, unless an odd prime),
  A322824 (A242424, unless an odd prime),
  A322973 (A006370, unless an odd prime),
  A322974 (A049820, unless n > 1 and n is in A046642),
  A323079 (A060681, unless an odd prime),
  A322587 (A295887, unless an odd prime),
  A322588 (A291751, unless an odd prime),
  A322592 (A289625, unless an odd prime),
  A323369 (A323368, unless an odd prime),
  A323371 (A295886, unless an odd prime),
  A323374 (A323373, unless an odd prime),
  A323401 (A323372, unless an odd prime),
  A323405 (A323404, unless an odd prime).

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; };
A322809aux(n) = if((n>2)&&isprime(n),-1,(n>>1));
v322809 = rgs_transform(vector(up_to,n,A322809aux(n)));
A322809(n) = v322809[n];

a(n) = A323161(n+1) - 1.

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

Antti Karttunen, Jan 04 2019

nonn

For all i, j:

A319704(i) = A319704(j) => a(i) = a(j) => A323082(i) = A323082(j).

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

Cf. A052330, A300840, A319704, A323081, A323082.

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];

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

A322809 Lexicographically earliest such sequence a that a(i) = a(j) => f(i) = f(j) for all i, j, where f(n) = -1 if n is an odd prime, and f(n) = floor(n/2) for all other numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

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.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI