A351950 -id:A351950 - cp's OEIS frontend

A351952 a(n) = A351950(n) / A351577(n).

0, 1, 1, 5, 2, 7, 1, 7, 8, 31, 13, 41, 2, 9, 11, 37, 16, 47, 3, 11, 14, 43, 19, 53, 1, 9, 10, 41, 17, 55, 12, 59, 71, 247, 106, 317, 19, 73, 92, 289, 127, 359, 26, 87, 113, 331, 148, 401, 2, 11, 13, 47, 20, 61, 17, 69, 86, 277, 121, 347, 24, 83, 107, 319, 142, 389, 31, 97, 128, 361, 163, 431, 3, 13, 16, 53, 23, 67

Offset: 0

Antti Karttunen, Apr 01 2022

Compare how different the scatter plot is to that of A342002, albeit with a very similar definition.

Note: this is at least partly because the other uses linear and the other uses logarithmic scatter plot. - Antti Karttunen, Oct 23 2024

Antti Karttunen, Table of n, a(n) for n = 0..20160
Index entries for sequences related to factorial base representation

Cf. A000142 (positions of 1's), A003415, A003557, A083345, A225901, A276076, A342001, A342002, A351576, A351577, A351950, A351953, A351954.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A351952(n) = { my(u=A276076(n)); (A003415(u) / A003557(u)); };

a(n) = A351950(n) / A351577(n).
a(n) = A342001(A276076(n)) = A083345(A276076(n)).
a(n) = A342002(A351576(n)).
a(n) = A351953(A225901(n)).

A327860 Arithmetic derivative of the primorial base exp-function: a(n) = A003415(A276086(n)).

Original entry on oeis.org

0, 1, 1, 5, 6, 21, 1, 7, 8, 31, 39, 123, 10, 45, 55, 185, 240, 705, 75, 275, 350, 1075, 1425, 3975, 500, 1625, 2125, 6125, 8250, 22125, 1, 9, 10, 41, 51, 165, 12, 59, 71, 247, 318, 951, 95, 365, 460, 1445, 1905, 5385, 650, 2175, 2825, 8275, 11100, 30075, 4125, 12625, 16750, 46625, 63375, 166125, 14, 77, 91, 329, 420

Offset: 0

Antti Karttunen, Sep 30 2019

Are there any other fixed points after 0, 1, 7, 8 and 2556? (A328110, see also A351087 and A351088).
Out of the 30030 initial terms, 19220 are multiples of 5. (See A327865).
Proof that a(n) is even if and only if n is a multiple of 4: Consider Charlie Neder's Feb 25 2019 comment in A235992. As A276086 is never a multiple of 4, and as it toggles the parity, we only need to know when A001222(A276086(n)) = A276150(n) is even. The condition for that is given in the latter sequence by David A. Corneth's Feb 27 2019 comment. From this it also follows that A166486 gives similarly the parity of terms of A342002, A351083 and A345000. See also comment in A327858. - Antti Karttunen, May 01 2022

			2556 has primorial base expansion [1,1,1,1,0,0] as 1*A002110(5) + 1*A002110(4) + 1*A002110(3) + 1*A002110(2) = 2310 + 210 + 30 + 6 = 2556. That in turn is converted by A276086 to 13^1 * 11^1 * 7^1 * 5^1 = 5005, whose arithmetic derivative is 5' * 1001 + 1001' * 5 = 1*1001 + 311*5 = 2556, thus 2556 is one of the rare fixed points (A328110) of this sequence.

Antti Karttunen, Table of n, a(n) for n = 0..2310
Antti Karttunen, Data supplement: n, a(n) computed for n = 0..30030
Index entries for sequences related to primorial base

Cf. A002110 (positions of 1's), A003415, A048103, A276086, A327858, A327859, A327865, A328110 (fixed points), A328233 (positions of primes), A328242 (positions of squarefree terms), A328388, A328392, A328571, A328572, A329031, A329032, A329041, A342002.
Cf. A345000, A351074, A351075, A351076, A351077, A351080, A351083, A351084, A351087 (numbers k such that a(k) is a multiple of k), A351088.
Coincides with A329029 on positions given by A276156.
Cf. A166486 (a(n) mod 2), A353630 (a(n) mod 4).
Cf. A267263, A276150, A324650, A324653, A324655 for omega, bigomega, phi, sigma and tau applied to A276086(n).
Cf. also A351950 (analogous sequence).

Block[{b = MixedRadix[Reverse@ Prime@ Range@ 12]}, Array[Function[k, If[# < 2, 0, # Total[#2/#1 & @@@ FactorInteger[#]] ] &@ Abs[Times @@ Power @@@ # &@ Transpose@{Prime@ Range@ Length@ k, Reverse@ k}]]@ IntegerDigits[#, b] &, 65, 0]] (* Michael De Vlieger, Mar 12 2021 *)

A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415
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; };
A327860(n) = A003415(A276086(n));

A327860(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= (p^e); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); }; \\ (Standalone version) - Antti Karttunen, Nov 07 2019

a(n) = A003415(A276086(n)).
a(A002110(n)) = 1 for all n >= 0.
From Antti Karttunen, Nov 03 2019: (Start)
Whenever A329041(x,y) = 1, a(x + y) = A003415(A276086(x)*A276086(y)) = a(x)*A276086(y) + a(y)*A276086(x). For example, we have:
a(n) = a(A328841(n)+A328842(n)) = A329031(n)*A328572(n) + A329032(n)*A328571(n).
A051903(a(n)) = A328391(n).
A328114(a(n)) = A328392(n).
(End)
From Antti Karttunen, May 01 2022: (Start)
a(n) = A328572(n) * A342002(n).
For all n >= 0, A000035(a(n)) = A166486(n). [See comments]
(End)

Verbal description added to the definition by Antti Karttunen, May 01 2022

A276076 Factorial base exp-function: digits in factorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Original entry on oeis.org

1, 2, 3, 6, 9, 18, 5, 10, 15, 30, 45, 90, 25, 50, 75, 150, 225, 450, 125, 250, 375, 750, 1125, 2250, 7, 14, 21, 42, 63, 126, 35, 70, 105, 210, 315, 630, 175, 350, 525, 1050, 1575, 3150, 875, 1750, 2625, 5250, 7875, 15750, 49, 98, 147, 294, 441, 882, 245, 490, 735, 1470, 2205, 4410, 1225, 2450, 3675, 7350, 11025, 22050, 6125, 12250, 18375, 36750, 55125, 110250, 343

Offset: 0

Antti Karttunen, Aug 18 2016

These are prime-factorization representations of single-variable polynomials where the coefficient of term x^(k-1) (encoded as the exponent of prime(k) in the factorization of n) is equal to the digit in one-based position k of the factorial base representation of n. See the examples.

			   n  A007623   polynomial     encoded as             a(n)
   -------------------------------------------------------
   0       0    0-polynomial   (empty product)        = 1
   1       1    1*x^0          prime(1)^1             = 2
   2      10    1*x^1          prime(2)^1             = 3
   3      11    1*x^1 + 1*x^0  prime(2) * prime(1)    = 6
   4      20    2*x^1          prime(2)^2             = 9
   5      21    2*x^1 + 1*x^0  prime(2)^2 * prime(1)  = 18
   6     100    1*x^2          prime(3)^1             = 5
   7     101    1*x^2 + 1*x^0  prime(3) * prime(1)    = 10
and:
  23     321  3*x^2 + 2*x + 1  prime(3)^3 * prime(2)^2 * prime(1)
                                      = 5^3 * 3^2 * 2 = 2250.

Antti Karttunen, Table of n, a(n) for n = 0..5040
Indranil Ghosh, Python program for computing this sequence.
Index entries for sequences related to factorial base representation.

Cf. A000040, A007623, A084558, A099563, A257687, A276009.
Cf. A276075 (a left inverse).
Cf. A276078 (same terms in ascending order).
Cf. also A000142, A001221, A001222, A002110, A007489, A008836, A019565, A033312, A034968, A048675, A051903, A059590, A060130, A076954, A246359, A248663, A262725, A276073, A276074, A351576, A351577, A351950, A351951, A351952, A351954.
Cf. also A275733, A275734, A275735, A275725 for other such encodings of factorial base related polynomials, and A276086 for a primorial base analog.

a[n_] := Module[{k = n, m = 2, r, p = 2, q = 1}, While[{k, r} = QuotientRemainder[k, m]; k != 0|| r != 0, q *= p^r; p = NextPrime[p]; m++]; q]; Array[a, 100, 0] (* Amiram Eldar, Feb 07 2024 *)

a(0) = 1, for n >= 1, a(n) = A275733(n) * a(A276009(n)).
Or: for n >= 1, a(n) = a(A257687(n)) * A000040(A084558(n))^A099563(n).
Other identities.
For all n >= 0:
A276075(a(n)) = n.
A001221(a(n)) = A060130(n).
A001222(a(n)) = A034968(n).
A051903(a(n)) = A246359(n).
A048675(a(n)) = A276073(n).
A248663(a(n)) = A276074(n).
a(A007489(n)) = A002110(n).
a(A059590(n)) = A019565(n).
For all n >= 1:
a(A000142(n)) = A000040(n).
a(A033312(n)) = A076954(n-1).
From Antti Karttunen, Apr 18 2022: (Start)
a(n) = A276086(A351576(n)).
A276085(a(n)) = A351576(n)
A003557(a(n)) = A351577(n).
A003415(a(n)) = A351950(n).
A069359(a(n)) = A351951(n).
A083345(a(n)) = A342001(a(n)) = A351952(n).
A351945(a(n)) = A351954(n).
A181819(a(n)) = A275735(n).
(End)
lambda(a(n)) = A262725(n+1), where lambda is Liouville's function, A008836. - Antti Karttunen and Peter Munn, Aug 09 2024

Name changed by Antti Karttunen, Apr 18 2022

A351577 a(n) = A003557(A276076(n)).

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 5, 5, 5, 5, 15, 15, 25, 25, 25, 25, 75, 75, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 3, 3, 5, 5, 5, 5, 15, 15, 25, 25, 25, 25, 75, 75, 7, 7, 7, 7, 21, 21, 7, 7, 7, 7, 21, 21, 35, 35, 35, 35, 105, 105, 175, 175, 175, 175, 525, 525, 49, 49, 49, 49, 147, 147, 49, 49, 49, 49, 147, 147, 245

Offset: 0

Antti Karttunen, Apr 01 2022

Antti Karttunen, Table of n, a(n) for n = 0..10080
Index entries for sequences related to factorial base representation

Cf. A003557, A085731, A276009, A276076, A328572, A351576, A351950, A351952.

A003557(n) = (n/factorback(factorint(n)[, 1]));
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A351577(n) = A003557(A276076(n));

a(n) = A003557(A276076(n)).
a(n) = A276076(A276009(n)).
a(n) = A328572(A351576(n)).
a(n) = A085731(A276076(n)) = gcd(A276076(n), A351950(n)).

A351951 a(n) = A069359(A276076(n)).

Original entry on oeis.org

0, 1, 1, 5, 3, 15, 1, 7, 8, 31, 24, 93, 5, 35, 40, 155, 120, 465, 25, 175, 200, 775, 600, 2325, 1, 9, 10, 41, 30, 123, 12, 59, 71, 247, 213, 741, 60, 295, 355, 1235, 1065, 3705, 300, 1475, 1775, 6175, 5325, 18525, 7, 63, 70, 287, 210, 861, 84, 413, 497, 1729, 1491, 5187, 420, 2065, 2485, 8645, 7455, 25935, 2100

Offset: 0

Antti Karttunen, Apr 01 2022

Antti Karttunen, Table of n, a(n) for n = 0..10080
Index entries for sequences related to factorial base representation

Cf. A069359, A276076, A329029, A351576.

Coincides with A351950 on positions given by A059590.

A069359(n) = (n*sumdiv(n, d, isprime(d)/d)); \\ From A069359
A276076(n) = { my(i=0,m=1,f=1,nextf); while((n>0),i=i+1; nextf = (i+1)*f; if((n%nextf),m*=(prime(i)^((n%nextf)/f));n-=(n%nextf));f=nextf); m; };
A351951(n) = A069359(A276076(n));

a(n) = A069359(A276076(n)).

a(n) = A329029(A351576(n)).

cp's OEIS Frontend

A351952 a(n) = A351950(n) / A351577(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A327860 Arithmetic derivative of the primorial base exp-function: a(n) = A003415(A276086(n)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A276076 Factorial base exp-function: digits in factorial base representation of n become the exponents of successive prime factors whose product a(n) is.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A351577 a(n) = A003557(A276076(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A351951 a(n) = A069359(A276076(n)).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula