A351461 -id:A351461 - cp's OEIS frontend

A336467 Fully multiplicative with a(2) = 1 and a(p) = A000265(p+1) for odd primes p, with A000265(k) giving the odd part of k.

1, 1, 1, 1, 3, 1, 1, 1, 1, 3, 3, 1, 7, 1, 3, 1, 9, 1, 5, 3, 1, 3, 3, 1, 9, 7, 1, 1, 15, 3, 1, 1, 3, 9, 3, 1, 19, 5, 7, 3, 21, 1, 11, 3, 3, 3, 3, 1, 1, 9, 9, 7, 27, 1, 9, 1, 5, 15, 15, 3, 31, 1, 1, 1, 21, 3, 17, 9, 3, 3, 9, 1, 37, 19, 9, 5, 3, 7, 5, 3, 1, 21, 21, 1, 27, 11, 15, 3, 45, 3, 7, 3, 1, 3, 15, 1, 49, 1, 3, 9, 51, 9, 13, 7, 3

Offset: 1

Antti Karttunen, Jul 22 2020

For the comment here, we extend the definition of the first kind of Cunningham chain (see Wikipedia-article) so that also isolated primes for which neither (p-1)/2 nor 2p+1 is a prime are considered to be in singular chains, that is, in chains of the length one. If we replace one or more instances of any particular odd prime factor p in n with any odd prime q of the same Cunningham chain, so that m = (q^k)*n / p^(e-k), where e is the exponent of p of n, and k <= e is the number of instances of p replaced with q, then it holds that a(m) = a(n), and by induction, the value stays invariant for any number of such replacements. Note also that A001222, but not necessarily A001221 will stay invariant in such changes.

For example, if some of the odd prime divisors p of n are Sophie Germain primes (in A005384), then replacing any of them with 2p+1 ("safe primes", i.e., the corresponding terms of A005385), gives a new number m, for which a(m) = a(n). And vice versa, the same is true for any safe prime factors > 5 of n (that are in A005385), then replacing any one of them with (p-1)/2 will not affect the result. For example, a(5*11*23*47) = a(11*11*23*23) = a(5^4) = a(11^4) = a(23^4) = 81, as 5, 11, 23 and 47 are in the same Cunningham chain of the first kind.

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Wikipedia, Cunningham chain

Cf. A000265, A003959, A206787, A336390, A336392, A336651, A336921, A336922, A351461.

Cf. also A335915, A336466 (similar sequences).

A000265(n) = (n>>valuation(n,2));
A336467(n) = { my(f=factor(n)); prod(k=1,#f~,if(2==f[k,1],1,(A000265(f[k,1]+1))^f[k,2])); };

For all n >= 1, A331410(a(n)) = A336921(n).
From Antti Karttunen, Nov 21 2023: (Start)
a(n) = A335915(n) / A336466(n).
a(1) = 1, and for n > 1, a(n) = A000265(A206787(n)) * a(A336651(n)).
(End)

A351460 Lexicographically earliest infinite sequence such that a(i) = a(j) => A006530(i) = A006530(j), A206787(i) = A206787(j) and A336651(i) = A336651(j) for all i, j >= 1.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Feb 11 2022

nonn

Restricted growth sequence transform of the ordered triplet [A006530(n), A206787(n), A336651(n)].
For all i, j >= 1:
  A324400(i) = A324400(j) => a(i) = a(j),
  a(i) = a(j) => A347241(i) = A347241(j),
  a(i) = a(j) => A351461(i) = A351461(j) => A347240(i) = A347240(j).

			a(429) = a(455) because 429 = 3*11*13 and 455 = 5*7*13, so they have equal largest prime factor (A006530), and they also agree on A206787(429) = A206787(455) = 672 and on A336651(429) = A336651(455) = 1 (because both are squarefree), therefore they get equal value (which is 216) allotted to them by the restricted growth sequence transform. - _Antti Karttunen_, Feb 14 2022

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

Cf. A006530, A206787, A336651, A347241, A351461.
Cf. also A324400, A351452.
Differs from A351454 for the first time at n=121, where a(121) = 62, while A351454(121) = 51.
Differs from A103391(1+n) for the first time after n=1 at n=455, where a(455) = 216, while A103391(456) = 229.

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; };
A006530(n) = if(1==n, n, my(f=factor(n)); f[#f~, 1]);
A206787(n) = sumdiv(n, d, d*(d % 2)*issquarefree(d)); \\ From A206787
A336651(n) = { my(f=factor(n)); prod(i=1, #f~, if(2==f[i,1],1,f[i,1]^(f[i,2]-1))); };
Aux351460(n) = [A006530(n), A206787(n), A336651(n)];
v351460 = rgs_transform(vector(up_to, n, Aux351460(n)));
A351460(n) = v351460[n];

A351040 Lexicographically earliest infinite sequence such that a(i) = a(j) => A336158(i) = A336158(j), A206787(i) = A206787(j) and A336651(i) = A336651(j) for all i, j >= 1.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 4, 1, 5, 3, 6, 2, 7, 4, 8, 1, 9, 5, 10, 3, 11, 6, 12, 2, 13, 7, 14, 4, 15, 8, 16, 1, 17, 9, 17, 5, 18, 10, 19, 3, 20, 11, 21, 6, 22, 12, 23, 2, 24, 13, 25, 7, 26, 14, 25, 4, 27, 15, 28, 8, 29, 16, 30, 1, 31, 17, 32, 9, 33, 17, 34, 5, 35, 18, 36, 10, 33, 19, 37, 3, 38, 20, 39, 11, 40, 21, 41, 6, 42

Offset: 1

Antti Karttunen, Jan 31 2022

Restricted growth sequence transform of the ordered triplet [A336158(n), A206787(n), A336651(n)].
For all i, j >= 1:
  A003602(i) = A003602(j) => a(i) = a(j),
  a(i) = a(j) => A336390(i) = A336390(j) => A336391(i) = A336391(j),
  a(i) = a(j) => A347374(i) = A347374(j),
  a(i) = a(j) => A351036(i) = A351036(j) => A113415(i) = A113415(j),
  a(i) = a(j) => A351461(i) = A351461(j).
From Antti Karttunen, Nov 23 2023: (Start)
Conjectured to be equal to the lexicographically earliest infinite sequence such that b(i) = b(j) => A000593(i) = A000593(j), A336158(i) = A336158(j) and A336467(i) = A336467(j), for all i, j >= 1 (this was the original definition). In any case it holds that a(i) = a(j) => b(i) = b(j) for all i, j >= 1. See comment in A351461.
(End)

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

Cf. A000265, A000593, A336158, A336467.
Cf. also A003602, A336390, A336391, A351037, A366891.
Differs from A347374 for the first time at n=103, where a(103) = 48, while A347374(103) = 30.
Differs from A351035 for the first time at n=175, where a(175) = 80, while A351035(175) = 78.
Differs from A351036 for the first time at n=637, where a(637) = 272, while A351036(637) = 261.

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; };
A000265(n) = (n>>valuation(n,2));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };
A336158(n) = A046523(A000265(n));
A206787(n) = sumdiv(n, d, d*(d % 2)*issquarefree(d));
A336651(n) = { my(f=factor(n)); prod(i=1, #f~, if(2==f[i,1],1,f[i,1]^(f[i,2]-1))); };
Aux351040(n) = [A336158(n), A206787(n), A336651(n)];
v351040 = rgs_transform(vector(up_to, n, Aux351040(n)));
A351040(n) = v351040[n];

Original definition moved to the comment section and replaced with a definition that is at least as encompassing, and conjectured to be equal to the original one. - Antti Karttunen, Nov 23 2023

A366881 Lexicographically earliest infinite sequence such that a(i) = a(j) => A206787(A163511(i)) = A206787(A163511(j)) and A336651(A163511(n)) = A336651(A163511(j)) for all i, j >= 0.

Original entry on oeis.org

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

Offset: 0

Antti Karttunen, Nov 04 2023

nonn

Restricted growth sequence transform of the ordered pair [A206787(A163511(n)), A336651(A163511(n))].
Restricted growth sequence transform of sequence b(n) = A351461(A163511(n)).
For all i, j >= 0:
  a(i) = a(j) => A324186(i) = A324186(j), (similarly for A366806)
  a(i) = a(j) => A366885(i) = A366885(j). (similarly for A366886).

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

Cf. A000265, A000593, A163511, A347385, A351461, A366887.
Cf. also A324186, A366806, A366885, A366886, A366888.
Differs from A366806 for the first time at n=105, where a(105) = 52, while A366806(105) = 19.

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; };
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));
A206787(n) = sumdiv(n, d, d*issquarefree(2*d));
A336651(n) = { my(f=factor(n>>valuation(n,2))); prod(i=1, #f~, f[i,1]^(f[i,2]-1)); };
A366881aux(n) = [A206787(A163511(n)), A336651(A163511(n))];
v366881 = rgs_transform(vector(1+up_to,n,A366881aux(n-1)));
A366881(n) = v366881[1+n];

For all n >= 1, a(n) = a(2*n) = a(A000265(n)).

A351441 Lexicographically earliest infinite sequence such that a(i) = a(j) => A003557(i) = A003557(j) and A351450(i) = A351450(j) for all i, j >= 1.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Feb 12 2022

nonn

Restricted growth sequence transform of the ordered pair [A003557(n), A351450(n)].
For all i, j >= 1:
   a(i) = a(j) => A326042(i) = A326042(j),
   a(i) = a(j) => A351455(i) = A351455(j).

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

Cf. A003557, A003961, A048250, A064989, A326042, A351450, A351455.

Cf. also A291751, A351461.

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; };
A003557(n) = (n/factorback(factorint(n)[, 1]));
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A048250(n) = factorback(apply(p -> p+1,factor(n)[,1]));
A064989(n) = { my(f = factor(n>>valuation(n,2))); for(i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f); };
A351450(n) = A064989(A048250(A003961(n)));
Aux351441(n) = [A003557(n), A351450(n)];
v351441 = rgs_transform(vector(up_to, n, Aux351441(n)));
A351441(n) = v351441[n];

A366891 Lexicographically earliest infinite sequence such that a(i) = a(j) => A365425(i) = A365425(j), A206787(A163511(i)) = A206787(A163511(j)) and A336651(A163511(n)) = A336651(A163511(j)) for all i, j >= 0.

Original entry on oeis.org

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

Offset: 0

Antti Karttunen, Nov 04 2023

nonn

Restricted growth sequence transform of the triplet [A365425(n), A206787(A163511(n)), A336651(A163511(n))], and also by conjecture, of sequence b(n) = A351040(A163511(n)).
For all i, j >= 0:
  a(i) = a(j) => A365395(i) = A365395(j),
  a(i) = a(j) => A366874(i) = A366874(j),
  a(i) = a(j) => A366881(i) = A366881(j).

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

Cf. A000265, A000593, A163511, A347385, A351040, A351461, A365425.
Differs from A366806 for the first time at n=105, where a(105) = 52, while A366806(105) = 19.
Differs from A366881 for the first time at n=511, where a(511) = 249, while A366881(511) = 7.

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; };
A000265(n) = (n>>valuation(n,2));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ From A046523
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));
A365425(n) = A046523(A000265(A163511(n)));
A206787(n) = sumdiv(n, d, d*issquarefree(2*d));
A336651(n) = { my(f=factor(n>>valuation(n,2))); prod(i=1, #f~, f[i,1]^(f[i,2]-1)); };
A366891aux(n) = [A365425(n), A206787(A163511(n)), A336651(A163511(n))];
v366891 = rgs_transform(vector(1+up_to,n,A366891aux(n-1)));
A366891(n) = v366891[1+n];

For all n >= 1, a(n) = a(2*n) = a(A000265(n)).

cp's OEIS Frontend

A336467 Fully multiplicative with a(2) = 1 and a(p) = A000265(p+1) for odd primes p, with A000265(k) giving the odd part of k.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A351460 Lexicographically earliest infinite sequence such that a(i) = a(j) => A006530(i) = A006530(j), A206787(i) = A206787(j) and A336651(i) = A336651(j) for all i, j >= 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A351040 Lexicographically earliest infinite sequence such that a(i) = a(j) => A336158(i) = A336158(j), A206787(i) = A206787(j) and A336651(i) = A336651(j) for all i, j >= 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Extensions

A366881 Lexicographically earliest infinite sequence such that a(i) = a(j) => A206787(A163511(i)) = A206787(A163511(j)) and A336651(A163511(n)) = A336651(A163511(j)) for all i, j >= 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A351441 Lexicographically earliest infinite sequence such that a(i) = a(j) => A003557(i) = A003557(j) and A351450(i) = A351450(j) for all i, j >= 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A366891 Lexicographically earliest infinite sequence such that a(i) = a(j) => A365425(i) = A365425(j), A206787(A163511(i)) = A206787(A163511(j)) and A336651(A163511(n)) = A336651(A163511(j)) for all i, j >= 0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula