A293221 -id:A293221 - cp's OEIS frontend

A293226 Restricted growth sequence transform of A293225, a filter combining two products, the other formed from the 1-digits (A293221) and the other from the 2-digits (A293222) present in the ternary expansions of proper divisors of n.

1, 2, 2, 3, 2, 4, 2, 5, 6, 7, 2, 8, 2, 9, 10, 11, 2, 12, 2, 13, 14, 15, 2, 16, 4, 17, 18, 19, 2, 20, 2, 21, 22, 23, 24, 25, 2, 26, 27, 28, 2, 29, 2, 30, 31, 32, 2, 33, 34, 35, 12, 36, 2, 37, 38, 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

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..19683

Cf. A001065, A293221, A293222, A293223, A293224, A293225.

Cf. also A290093, A290094, 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); };
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; };
Anot_submitted(n) = (1/2)*(2 + ((A293222(n) + A293221(n))^2) - A293222(n) - 3*A293221(n)); \\ Eq.class-wise equal to A293225.
write_to_bfile(1,rgs_transform(vector(19683,n,Anot_submitted(n))),"b293226.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");

A293450 Restricted growth sequence transform of (3*A293225(n) + A010872(n)), a filter combining (n mod 3) with two products, the other formed from the 1-digits (A293221) and the other from the 2-digits (A293222) present in the ternary expansions of proper divisors of n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Nov 06 2017

nonn

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

Cf. A002324, A010872, A293221, A293222, 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
A289813(n) = { my (d=digits(n, 3)); from digits(vector(#d, i, if (d[i]==1, 1, 0)), 2); };
A289814(n) = { my (d=digits(n, 3)); from digits(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; };
Anot_submitted(n) = (1/2)*(2 + ((A293222(n) + A293221(n))^2) - A293222(n) - 3*A293221(n)); \\ Eq.class-wise equal to A293225.
Anot2submitted(n) = ((3*Anot_submitted(n))+(n%3));
write_to_bfile(1,rgs_transform(vector(59049,n,Anot2submitted(n))),"b293450.txt");

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

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

Original entry on oeis.org

1, 2, 2, 6, 2, 36, 2, 30, 12, 60, 2, 2700, 2, 180, 120, 210, 2, 7560, 2, 6300, 360, 252, 2, 661500, 20, 420, 168, 94500, 2, 23814000, 2, 2310, 504, 132, 600, 43659000, 2, 396, 840, 2425500, 2, 187110000, 2, 207900, 352800, 1980, 2, 560290500, 60, 194040, 264, 485100, 2, 115259760, 840, 254677500, 792, 4620, 2, 264737261250000, 2, 13860

Offset: 1

Antti Karttunen, Oct 03 2017

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

Cf. A001065, A002110, A019565, A048675, A091954, A292257, A293215 (restricted growth sequence transform).

Cf. also A293216, A293221, A293222, A293225, A293231, A293442, A300830, A300831, A300832, A300833, A300834.

A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A293214(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(d))); m; };

a(n) = Product_{d|n, dA019565(d).
a(n) = A300830(n) * A300831(n) * A300832(n). - Antti Karttunen, Mar 16 2018
Other identities.
For n >= 0, a(2^n) = A002110(n).
For n >= 1:
A048675(a(n)) = A001065(n).
A001222(a(n)) = A292257(n).
A007814(a(n)) = A091954(n).
A087207(a(n)) = A218403(n).
A248663(a(n)) = A227320(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)).

A319991 a(n) = Product_{d|n, dA019565(d)^[1 == d mod 3].

Original entry on oeis.org

1, 2, 2, 2, 2, 2, 2, 10, 2, 2, 2, 10, 2, 60, 2, 10, 2, 2, 2, 210, 60, 2, 2, 10, 2, 140, 2, 300, 2, 42, 2, 110, 2, 2, 60, 10, 2, 132, 140, 210, 2, 60, 2, 1650, 2, 2, 2, 110, 60, 6468, 2, 700, 2, 2, 2, 115500, 132, 2, 2, 210, 2, 4620, 60, 110, 140, 330, 2, 390, 2, 1260, 2, 10, 2, 260, 308, 660, 60, 140, 2, 210210, 2, 2, 2, 115500, 2, 1092, 2

Offset: 1

Antti Karttunen, Oct 03 2018

nonn

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

Cf. A019565, A293214, A293895, A293897, A319990, A319992, A320001, A320011 (rgs-transform).

Cf. also A293221.

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
A319991(n) = { my(m=1); fordiv(n,d,if((dA019565(d))); m; };

a(n) = Product_{d|n, dA019565(d)^[1 == d mod 3].
a(n) = A293214(n) / (A319990(n)*A319992(n)).
For all n >= 1:
A007814(a(n)) = A320001(n).
A048675(a(n)) = A293897(n).
A195017(a(n)) = A293895(n) mod 3.

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)).

A300834 a(n) = Product_{d|n, dA019565(A003714(d)), where A003714(n) is the n-th Fibbinary number.

Original entry on oeis.org

1, 2, 2, 6, 2, 30, 2, 60, 10, 42, 2, 4200, 2, 126, 70, 660, 2, 9240, 2, 13860, 210, 330, 2, 5082000, 14, 78, 220, 32760, 2, 3783780, 2, 42900, 550, 780, 294, 924924000, 2, 1092, 130, 41621580, 2, 3898440, 2, 112200, 60060, 306, 2, 28078050000, 42, 235620, 1300, 92820, 2, 200119920, 770, 128648520, 1820, 1122, 2, 424964656116000, 2, 3366

Offset: 1

Antti Karttunen, Mar 18 2018

nonn

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

Cf. A003714, A019565, A300835 (rgs-transform of this sequence), A300836.

Cf. also A293214, A293216, A293221, A293222, A293231.

A072649(n) = { my(m); if(n<1, 0, m=0; until(fibonacci(m)>n, m++); m-2); }; \\ From A072649
A003714(n) = { my(s=0,w); while(n>2, w = A072649(n); s += 2^(w-1); n -= fibonacci(w+1)); (s+n); }
A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A300834(n) = { my(m=1); fordiv(n,d,if(d < n,m *= A019565(A003714(d)))); m; };

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

For n >= 1, A001222(a(n)) = A300836(n).

A296071 a(n) = Product_{d|n, dA019565(A289813(A295882(d))); a product obtained from the 1's present in balanced ternary representation of the deficiencies of the proper divisors of n.

Original entry on oeis.org

1, 2, 2, 4, 2, 12, 2, 8, 6, 24, 2, 24, 2, 20, 36, 16, 2, 60, 2, 144, 30, 40, 2, 48, 12, 60, 30, 240, 2, 1080, 2, 32, 60, 56, 60, 120, 2, 28, 90, 576, 2, 3600, 2, 400, 900, 168, 2, 96, 10, 1008, 84, 1200, 2, 420, 120, 480, 42, 56, 2, 4320, 2, 84, 1500, 64, 180, 4200, 2, 784, 252, 90720, 2, 1200, 2, 140, 2520, 784, 100, 75600, 2, 1152, 210, 840, 2

Offset: 1

Antti Karttunen, Dec 04 2017

nonn

Used as a part of filter A296073.

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

Cf. A033879, A019565, A289813, A295882, A296072, A296073.

Cf. also A293221.

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

(define (A296071 n) (let loop ((m 1) (props (proper-divisors n))) (cond ((null? props) m) (else (loop (* m (A019565 (A289813 (A295882 (car props))))) (cdr props))))))
(define (proper-divisors n) (reverse (cdr (reverse (divisors n)))))
(define (divisors n) (let loop ((k n) (divs (list))) (cond ((zero? k) divs) ((zero? (modulo n k)) (loop (- k 1) (cons k divs))) (else (loop (- k 1) divs)))))

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

A319708 a(n) = Product_{d|n, dA276086(d).

Original entry on oeis.org

1, 2, 2, 6, 2, 36, 2, 54, 12, 108, 2, 1620, 2, 60, 216, 810, 2, 5400, 2, 43740, 120, 540, 2, 607500, 36, 300, 360, 40500, 2, 21870000, 2, 182250, 1080, 2700, 360, 151875000, 2, 1500, 600, 246037500, 2, 101250000, 2, 5467500, 972000, 13500, 2, 85429687500, 20, 6075000, 5400, 5062500, 2, 2531250000, 3240, 3417187500, 3000, 67500, 2

Offset: 1

Antti Karttunen, Oct 03 2018

nonn

Antti Karttunen, Table of n, a(n) for n = 1..2310
Index entries for sequences related to primorial base

Cf. A276085, A276086, A319709 (rgs-transform).

Cf. A293214, A293221, A293222, A300834 for similar constructions for other bases.

A276086(n) = { my(i=0,m=1,pr=1,nextpr); while((n>0),i=i+1; nextpr = prime(i)*pr; if((n%nextpr),m*=(prime(i)^((n%nextpr)/pr));n-=(n%nextpr));pr=nextpr); m; };
A319708(n) = { my(m=1); fordiv(n, d, if(dA276086(d))); (m); };

a(n) = Product_{d|n, dA276086(d).
For all n >= 1:
A276085(a(n)) = A001065(n).
A001222(a(n)) = A319713(n).

cp's OEIS Frontend

A293226 Restricted growth sequence transform of A293225, a filter combining two products, the other formed from the 1-digits (A293221) and the other from the 2-digits (A293222) present in the ternary expansions of proper divisors of n.

Views

Author

Keywords

Comments

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

A293450 Restricted growth sequence transform of (3*A293225(n) + A010872(n)), a filter combining (n mod 3) with two products, the other formed from the 1-digits (A293221) and the other from the 2-digits (A293222) present in the ternary expansions of proper divisors of n.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

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

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

A319991 a(n) = Product_{d|n, dA019565(d)^[1 == d mod 3].

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

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

A300834 a(n) = Product_{d|n, dA019565(A003714(d)), where A003714(n) is the n-th Fibbinary number.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A296071 a(n) = Product_{d|n, dA019565(A289813(A295882(d))); a product obtained from the 1's present in balanced ternary representation of the deficiencies of the proper divisors of n.

Views

Author

Keywords

Comments

Links

Crossrefs