A293226 -id:A293226 - cp's OEIS frontend

A293215 Restricted growth sequence transform of A293214, where A293214(n) = Product_{d|n, dA019565(d).

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

Offset: 1

Antti Karttunen, Oct 03 2017

nonn

For all i, j: a(i) = a(j) => A001065(i) = A001065(j).

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

Cf. A001065, A019565, A048675, A293214, A293217 (a variant), A293223, A293224, A293226.

Differs from related A293232 for the first time at n=55, where a(55) = 28, while A293232(55) = 39.

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; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A293214(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(d))); m; };
write_to_bfile(1,rgs_transform(vector(16384,n,A293214(n))),"b293215.txt");

A293221 a(n) = Product_{d|n, dA019565(A289813(d)); a product obtained from the 1-digits present in ternary expansions of proper divisors of n.

Original entry on oeis.org

1, 2, 2, 2, 2, 6, 2, 12, 6, 6, 2, 36, 2, 4, 18, 12, 2, 30, 2, 360, 12, 10, 2, 540, 6, 60, 30, 360, 2, 900, 2, 120, 30, 10, 12, 2700, 2, 4, 180, 360, 2, 540, 2, 360, 450, 6, 2, 5400, 4, 120, 30, 360, 2, 210, 30, 5040, 12, 14, 2, 1701000, 2, 84, 180, 2520, 180, 1260, 2, 840, 18, 12600, 2, 94500, 2, 140, 180, 840, 20, 18900, 2, 756000, 210, 210, 2, 23814000, 30

Offset: 1

Antti Karttunen, Oct 03 2017

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

Cf. A019565, A289813, A293214, A293222, A293223 (restricted growth sequence transform), A293226.

Cf. also A290091.

A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A289813(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); } \\ From Remy Sigrist
A293221(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289813(d)))); m; };

a(n) = Product_{d|n, dA019565(A289813(d)).

For all n >= 0, a(3^n) = A002110(n).

A293222 a(n) = Product_{d|n, dA019565(A289814(d)); a product obtained from the 2-digits present in ternary expansions of proper divisors of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 1, 4, 1, 6, 1, 6, 2, 12, 1, 6, 1, 4, 3, 4, 1, 36, 2, 2, 1, 12, 1, 36, 1, 36, 2, 12, 6, 30, 1, 10, 1, 240, 1, 180, 1, 20, 6, 20, 1, 1620, 3, 60, 6, 60, 1, 30, 4, 72, 5, 4, 1, 360, 1, 2, 15, 72, 2, 180, 1, 36, 10, 144, 1, 2700, 1, 2, 90, 20, 6, 180, 1, 720, 1, 4, 1, 540, 12, 6, 2, 720, 1, 900, 3, 100, 1, 20, 10, 16200, 1, 60, 6

Offset: 1

Antti Karttunen, Oct 03 2017

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

Cf. A019565, A289814, A293221, A293224 (restricted growth sequence transform), A293226.

A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A289814(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==2, 1, 0)), 2); } \\ From Remy Sigrist
A293222(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289814(d)))); m; };

a(n) = Product_{d|n, dA019565(A289814(d)).

A300833 Filter sequence combining A300830(n), A300831(n) and A300832(n), three products formed from such proper divisors d of n for which mu(n/d) = 0, +1 or -1 respectively, where mu is Möbius mu function (A008683).

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Mar 16 2018

nonn

Restricted growth sequence transform of triple [A300830(n), A300831(n), A300832(n)].
For all i, j:
  a(i) = a(j) => A293215(i) = A293215(j) => A001065(i) = A001065(j).
  a(i) = a(j) => A051953(i) = A051953(j).
  a(i) = a(j) => A295885(i) = A295885(j).
Apparently this is also the restricted growth sequence transform of ordered pair [A300831(n), A300832(n)], which is true if it holds that whenever we have A300831(i) = A300831(j) and A300832(i) = A300832(j) for any i, j, then also A300830(i) = A300830(j). This has been checked for the first 65537 terms.

			a(39) = a(55) (= 28) as A300830(39) = A300830(55) = 1, A300831(39) = A300831(55) = 2 and A300832(39) = A300832(55) = 420.

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

Cf. A019565, A051953, A300831, A300832.

Cf. also A293214, A293215, A293226, A295885, A300825.

allocatemem(2^30);
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; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A300830(n) = { my(m=1); fordiv(n,d,if(!moebius(n/d),m *= A019565(d))); m; };
A300831(n) = { my(m=1); fordiv(n,d,if((d < n)&&(1==moebius(n/d)), m *= A019565(d))); m; };
A300832(n) = { my(m=1); fordiv(n,d,if(-1==moebius(n/d), m *= A019565(d))); m; };
Aux300833(n) = [A300830(n), A300831(n), A300832(n)];
write_to_bfile(1,rgs_transform(vector(up_to,n,Aux300833(n))),"b300833.txt");

A293223 Restricted growth sequence transform of A293221, a product formed from the 1-digits present in the ternary expansion of proper divisors of n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Oct 03 2017

nonn

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

Cf. A293221.

Cf. also A290091, A293224, A293225, A293226, A293215, A293217, A293232.

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; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A289813(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); } \\ From Remy Sigrist
A293221(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289813(d)))); m; };
write_to_bfile(1,rgs_transform(vector(19683,n,A293221(n))),"b293223.txt");

A293224 Restricted growth sequence transform of A293222, a product formed from the 2-digits present in the ternary expansions of proper divisors of n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Oct 03 2017

nonn

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

Cf. A293222.

Cf. also A290092, A293223, A293215, A293217, A293232, A293225, A293226.

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; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A289814(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==2, 1, 0)), 2); } \\ From Remy Sigrist
A293222(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289814(d)))); m; };
write_to_bfile(1,rgs_transform(vector(19683,n,A293222(n))),"b293224.txt");

A293225 Compound filter: a(n) = P(A293224(n), A293223(n)), where P(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

1, 2, 2, 5, 2, 8, 2, 12, 4, 13, 2, 32, 2, 40, 30, 33, 2, 59, 2, 58, 42, 69, 2, 143, 8, 80, 29, 83, 2, 178, 2, 197, 38, 96, 25, 239, 2, 100, 121, 163, 2, 221, 2, 202, 194, 103, 2, 448, 61, 365, 59, 245, 2, 333, 48, 576, 187, 256, 2, 720, 2, 278, 546, 718, 138, 606, 2, 503, 114, 1009, 2, 1101, 2, 437, 651, 678, 532, 831, 2, 1400, 172, 213, 2, 1508, 71, 500, 597

Offset: 1

Antti Karttunen, Oct 03 2017

nonn

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

Cf. A000027, A019565, A293221, A293222, A293223, A293224, A293226 (rgs-version of this filter).

Cf. also A289813, A289814, A290093, A290094.

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; };
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A289813(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); };
A289814(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==2, 1, 0)), 2); };
A293221(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289813(d)))); m; };
A293222(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289814(d)))); m; };
v293223 = rgs_transform(vector(19683,n,A293221(n)));
A293223(n) = v293223[n];
v293224 = rgs_transform(vector(19683,n,A293222(n)));
A293224(n) = v293224[n];
A293225(n) = (1/2)*(2 + ((A293224(n) + A293223(n))^2) - A293224(n) - 3*A293223(n));

(define (A293225 n) (* 1/2 (+ (expt (+ (A293224 n) (A293223 n)) 2) (- (A293224 n)) (- (* 3 (A293223 n))) 2)))

a(n) = (1/2)*(2 + ((A293224(n) + A293223(n))^2) - A293224(n) - 3*A293223(n)).

A296073 Filter combining A296071(n) and A296072(n), related to the deficiencies of proper divisors of n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Dec 04 2017

nonn

Construction: Pack the values of A296071(n) and A296072(n) to a single value with any injective N x N -> N packing function, like for example as f(n) = (1/2)*(2 + ((A296071(n)+A296072(n))^2) - A296071(n) - 3*A296072(n)) (the packing function here is the two-argument form of A000027). Then apply the restricted growth sequence transform to the sequence f(1), f(2), f(3), ... The transform assigns a unique increasing number for each newly encountered term of the sequence, and for any subsequent occurrences of the same term it gives the same number that term obtained for the first time.
For all i, j: a(i) = a(j) => A296074(i) = A296074(j).
Note that this is NOT restricted growth transform of A239968, which is A305800. Apart from 2's that occur at every prime, there are other duplicates also, first at a(125) = a(46) = 33.

			To see that a(46) and a(125) have the same value (33), consider the proper divisors of 46 = 1, 2, 23 and of 125 = 1, 5, 25. Their deficiencies are 1, 1, 22 and 1, 4, 19 respectively. When we look at their balanced ternary representations [as here all elements are positive, it can be obtained as A007089(A117967(n)) with 2's standing for -1's]:
   1 =    1
   1 =    1
  22 = 1211 (as 22 = 1*(3^3) + -1*(3^2) + 1*(3^1) + 1*(3^0))
and
   1 =    1
   4 =   11
  19 = 1201 (as 19 = 1*(3^3) + -1*(3^2) + 0*(3^1) + 1*(3^0)).
we see that in each column there is an equal number of 1's and an equal number of 2's. Moreover, this then implies also that the sums of those two sequences of deficiencies {1, 1, 22} and {1, 4, 19} are equal, as A296074(n) is a function of (can be computed from) a(n).

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

Cf. A019565, A033879, A117967, A117968, A295882, A296071, A296072, A296074.
Cf. also A293226.
Differs from A305800 for the first time at n=125.

up_to = 65536;
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; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ This function from M. F. Hasler
A117967(n) = if(n<=1,n,if(!(n%3),3*A117967(n/3),if(1==(n%3),1+3*A117967((n-1)/3),2+3*A117967((n+1)/3))));
A117968(n) = if(1==n,2,if(!(n%3),3*A117968(n/3),if(1==(n%3),2+3*A117968((n-1)/3),1+3*A117968((n+1)/3))));
A289813(n) = { my (d=digits(n, 3)); from digits(vector(#d, i, if (d[i]==1, 1, 0)), 2); } \\ From Rémy Sigrist
A289814(n) = { my (d=digits(n, 3)); from digits(vector(#d, i, if (d[i]==2, 1, 0)), 2); } \\ From Rémy Sigrist
A295882(n) = { my(x = (2*n)-sigma(n)); if(x >= 0,A117967(x),A117968(-x)); };
A296071(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289813(A295882(d))))); m; };
A296072(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A289814(A295882(d))))); m; };
Anotsubmitted3(n) = (1/2)*(2 + ((A296071(n)+A296072(n))^2) - A296071(n) - 3*A296072(n));
write_to_bfile(1,rgs_transform(vector(up_to,n,Anotsubmitted3(n))),"b296073.txt");

Data section extended up to a(125) by Antti Karttunen, Jun 14 2018

A351030 Lexicographically earliest infinite sequence such that a(i) = a(j) => A351031(i) = A351031(j) and A351032(i) = A351032(j), for all i, j >= 1.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Jan 29 2022

nonn

Restricted growth sequence transform of the ordered pair [A351031(n), A351032(n)], or equally, of the ordered pair [A351033(n), A351034(n)].

For all i, j: a(i) = a(j) => A349910(i) = A349910(j).

Antti Karttunen, Table of n, a(n) for n = 1..65537
Index entries for sequences computed from indices in prime factorization

Cf. A019565, A048673, A289813, A289814, A291759, A304759, A351031, A351032 A351033, A351034.

Cf. also A293226, A349910.

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; };
A019565(n) = { my(m=1, p=1); while(n>0, p = nextprime(1+p); if(n%2, m *= p); n >>= 1); (m); };
A048673(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); (1+factorback(f))/2; };
A289813(n) = { my(d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); };
A289814(n) = { my(d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==2, 1, 0)), 2); };
A291759(n) = A289814(A048673(n));
A304759(n) = A289813(A048673(n));
A351031(n) = { my(m=1); fordiv(n,d,if(dA019565(A304759(d)))); (m); };
A351032(n) = { my(m=1); fordiv(n,d,if(dA019565(A291759(d)))); (m); };
Aux351030(n) = [A351031(n),A351032(n)];
v351030 = rgs_transform(vector(up_to, n, Aux351030(n)));
A351030(n) = v351030[n];

A351090 Lexicographically earliest infinite sequence such that a(i) = a(j) => A351091(i) = A351091(j) and A351092(i) = A351092(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, 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, 31, 16, 32, 1, 33, 17, 34, 9, 35, 18, 36, 5, 37, 19, 38, 10, 39, 20, 40, 3, 41, 21, 42, 11, 43, 22, 44, 6, 45, 23

Offset: 1

Antti Karttunen, Jan 31 2022

Restricted growth sequence transform of the ordered pair [A351091(n), A351092(n)], or equally, of the ordered pair [A351093(n), A351094(n)].

For all i, j: A003602(i) = A003602(j) => a(i) = a(j) => A000593(i) = A000593(j).

			Consider two odd semiprimes, 689 and 697. The divisors of 689 are 1, 13, 53, 689, and the divisors of 697 are 1, 17, 41, 697. Applying A019565(A289813(x)) to the former gives [2, 30, 7, 105], while with the latter it gives [2, 5, 105, 42], and the product of both sequences is 44100. Applying A019565(A289814(x)) to the former gives [1, 1, 30, 286], while with the latter it gives [1, 6, 2, 715]. Product of both sequences is 8580. Therefore, because A351091(689) = A351091(697) and A351092(689) = A351092(697), also a(689) = a(697).

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

Cf. A000593, A351037, A351091, A351092, A351093, A351094.
Differs from A003602 for the first time at n=697, where a(697) = 345 while A003602(697) = 349.
Cf. also A293226, A351030.

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; };
A019565(n) = { my(m=1, p=1); while(n>0, p = nextprime(1+p); if(n%2, m *= p); n >>= 1); (m); };
A289813(n) = { my(d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==1, 1, 0)), 2); }; \\ From A289813
A289814(n) = { my (d=digits(n, 3)); fromdigits(vector(#d, i, if (d[i]==2, 1, 0)), 2); }; \\ From A289814
A351091(n) = { my(m=1); fordiv(n>>valuation(n,2),d,m *= A019565(A289813(d))); (m); };
A351092(n) = { my(m=1); fordiv(n>>valuation(n,2),d,m *= A019565(A289814(d))); (m); };
Aux351090(n) = [A351091(n),A351092(n)];
v351090 = rgs_transform(vector(up_to, n, Aux351090(n)));
A351090(n) = v351090[n];

cp's OEIS Frontend

A293215 Restricted growth sequence transform of A293214, where A293214(n) = Product_{d|n, dA019565(d).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A293221 a(n) = Product_{d|n, dA019565(A289813(d)); a product obtained from the 1-digits present in ternary expansions of proper divisors of n.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A293222 a(n) = Product_{d|n, dA019565(A289814(d)); a product obtained from the 2-digits present in ternary expansions of proper divisors of n.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A300833 Filter sequence combining A300830(n), A300831(n) and A300832(n), three products formed from such proper divisors d of n for which mu(n/d) = 0, +1 or -1 respectively, where mu is Möbius mu function (A008683).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A293223 Restricted growth sequence transform of A293221, a product formed from the 1-digits present in the ternary expansion of proper divisors of n.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

A293224 Restricted growth sequence transform of A293222, a product formed from the 2-digits present in the ternary expansions of proper divisors of n.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

A293225 Compound filter: a(n) = P(A293224(n), A293223(n)), where P(n,k) is sequence A000027 used as a pairing function.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Scheme

Formula

A296073 Filter combining A296071(n) and A296072(n), related to the deficiencies of proper divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Extensions

A351030 Lexicographically earliest infinite sequence such that a(i) = a(j) => A351031(i) = A351031(j) and A351032(i) = A351032(j), for all i, j >= 1.

Views

Author

Keywords

Comments