cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A329605 Number of divisors of A108951(n), where A108951 is fully multiplicative with a(prime(i)) = prime(i)# = Product_{i=1..i} A000040(i).

Original entry on oeis.org

1, 2, 4, 3, 8, 6, 16, 4, 9, 12, 32, 8, 64, 24, 18, 5, 128, 12, 256, 16, 36, 48, 512, 10, 27, 96, 16, 32, 1024, 24, 2048, 6, 72, 192, 54, 15, 4096, 384, 144, 20, 8192, 48, 16384, 64, 32, 768, 32768, 12, 81, 36, 288, 128, 65536, 20, 108, 40, 576, 1536, 131072, 30, 262144, 3072, 64, 7, 216, 96, 524288, 256, 1152, 72, 1048576, 18, 2097152, 6144, 48
Offset: 1

Views

Author

Antti Karttunen, Nov 18 2019

Keywords

Links

Crossrefs

Cf. A000005, A000040, A000142, A002110, A010052, A034386, A052126, A056239, A108951, A329902, A329378, A329382, A329614, A329617, A331283, A331286 (odd part).
Cf. A329606 (rgs-transform), A329608, A331284 (ordinal transform).
Cf. A331285 (the position where for the first time some term has occurred n times in this sequence).

Programs

  • Mathematica
    Block[{a}, a[n_] := a[n] = Module[{f = FactorInteger[n], p, e}, If[Length[f] > 1, Times @@ a /@ Power @@@ f, {{p, e}} = f; Times @@ (Prime[Range[PrimePi[p]]]^e)]]; a[1] = 1; Array[DivisorSigma[0, a@ #] &, 75]] (* Michael De Vlieger, Jan 24 2020, after Jean-François Alcover at A108951 *)
  • PARI
    A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A329605(n) = numdiv(A108951(n));
  • PARI
    A329605(n) = if(1==n,1,my(f=factor(n),e=1,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= e^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); (m)); \\ Antti Karttunen, Jan 14 2020
  • PARI
    A329605(n) = if(1==n,1,my(f=factor(n),e=0,d); forstep(i=#f~,1,-1, e += f[i,2]; d = (primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1]))); f[i,1] = (e+1); f[i,2] = d); factorback(f)); \\ Antti Karttunen, Jan 14 2020

Formula

a(n) = A000005(A108951(n)).
a(n) >= A329382(n) >= A329617(n) >= A329378(n).
A020639(a(n)) = A329614(n).
From Antti Karttunen, Jan 14 2020: (Start)
a(A052126(n)) = A329382(n).
a(A002110(n)) = A000142(1+n), for all n >= 0.
a(n) > A056239(n).
a(A329902(n)) = A002183(n).
A000265(a(n)) = A331286(n).
gcd(n,a(n)) = A331283(n).
If n = p(k1)^e(k1) * p(k2)^e(k2) * p(k3)^e(k3) * ... * p(kx)^e(kx), with p(n) = A000040(n) and k1 > k2 > ... > kx, then a(n) = (1+e(k1))^(k1-k2) * (1+e(k1)+e(k2))^(k2-k3) * ... * (1+e(k1)+e(k2)+...+e(kx))^kx.
A000035(a(n)) = A000035(A000005(n)) = A010052(n).
(End)

A329608 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = A329605(n) for all other n, except for odd primes p, f(p) = 0.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Nov 18 2019

Keywords

Comments

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

Links

Crossrefs

Cf. A000005, A108951, A305801, A329605, A329606, A329614.

Programs

  • PARI
    up_to = 100000;
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; };
A034386(n) = prod(i=1, primepi(n), prime(i));
A108951(n) = { my(f=factor(n)); prod(i=1, #f~, A034386(f[i, 1])^f[i, 2]) };  \\ From A108951
A329605(n) = numdiv(A108951(n));
Aux329608(n) = if((n%2)&&isprime(n),0,A329605(n));
v329608 = rgs_transform(vector(up_to, n, Aux329608(n)));
A329608(n) = v329608[n];

A331284 Number of values of k, 1 <= k <= n, with A329605(k) = A329605(n), where A329605 is the number of divisors of primorial inflation of n (A108951).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 3, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 3, 1, 1, 3, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 3, 2, 1, 2, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 3, 1, 1, 2, 1, 1, 3
Offset: 1

Views

Author

Antti Karttunen, Jan 14 2020

Keywords

Comments

Ordinal transform of A329605, or equally, of A329606.

Links

Crossrefs

Cf. A000005, A108951, A329605, A329606, A331285 (positions of the first occurrences of each n, also positions of records).
Cf. also A067004.

Programs

  • Mathematica
    c[n_] := c[n] = If[n == 1, 1, Module[{f = FactorInteger[n], p, e}, If[Length[f] > 1, Times @@ c /@ Power @@@ f, {{p, e}} = f; Times @@ (Prime[Range[PrimePi[p]]]^e)]]];
A329605[n_] := DivisorSigma[0, c[n]];
Module[{b}, b[_] = 0;
a[n_] := With[{t = A329605[n]}, b[t] = b[t] + 1]];
Array[a, 105] (* Jean-François Alcover, Jan 12 2022 *)
  • PARI
    up_to = 65537;
ordinal_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), pt); for(i=1, length(invec), if(mapisdefined(om,invec[i]), pt = mapget(om, invec[i]), pt = 0); outvec[i] = (1+pt); mapput(om,invec[i],(1+pt))); outvec; };
A329605(n) = if(1==n,1,my(f=factor(n),e=1,m=1); forstep(i=#f~,1,-1, e += f[i,2]; m *= e^(primepi(f[i,1])-if(1==i,0,primepi(f[i-1,1])))); (m));
v331284 = ordinal_transform(vector(up_to, n, A329605(n)));
A331284(n) = v331284[n];

Formula

a(A331285(n)) = n for all n.

A331285 a(n) is the index of the first occurrence of n in A331284.

Original entry on oeis.org

1, 8, 27, 108, 180, 396, 672, 1056, 1372, 1760, 2352, 3087, 3696, 3744, 4896, 5733, 7497, 9724, 11907, 13600, 15200, 18513, 19773, 23940, 24752, 28917, 32319, 33534, 42282, 45472, 47500, 52668, 55890, 59976, 66048, 74240, 77792, 81144, 86944, 100035, 105248, 109368, 122825, 127908, 134368, 144648, 156325, 168948, 175770
Offset: 1

Views

Author

Antti Karttunen, Jan 14 2020

Keywords

Comments

Also positions of records in A331284.
a(n) is the least k such that in range 1 .. k of A329605 (equally: in A329606[1..k]) there can be found exactly n occurrences of some term. In A329605 these "champion terms" are A329605(a(n)): 1, 4, 16, 24, 48, 192, 96, 192, 384, 288, 288, 576, 576, 576, 1152, 2304, 4608, 9216, 576, 2304, ..., that appear all to be 3-smooth numbers (in A003586).
For example, a(5)=180 and 48 is the first term in A329605 to occur for five times, as A329605(22) = A329605(42) = A329605(75) = A329605(112) = A329605(180) = 48.

Crossrefs

Cf. A331284 (a left inverse).
Cf. A003586, A329605, A329606.

Formula

A331284(a(n)) = n for all n >= 1.
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.