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.

Previous Showing 91-100 of 158 results. Next

A336389 The least positive integer k for which A336835(k) >= n, where A336835(k) is the number of iterations of x -> A003961(x) needed before the result is deficient (sigma(x) < 2x), when starting from x=k.

Original entry on oeis.org

1, 6, 120, 19399380, 195534950863140268380, 538938984694949877040715541221415046162838700, 216487559804430601784907786655491617909711008142914104790481010259258659171900
Offset: 0

Views

Author

Antti Karttunen, Aug 07 2020

Keywords

Comments

For n > 0, the least k such that for at least n-1 iterations of map x -> A003961(x), starting from x=k, x stays nondeficient. In other words, from each a(n) starts a chain of at least n nondeficient numbers (A023196) obtained by successive prime shifts, e.g, for a(3) we have: 19399380 -> 334639305 -> 5391411025, where -> stands for applying A003961, the prime shift towards larger primes.
After 1 all other terms here are even, because if an odd number k is nondeficient, then A064989(k) is nondeficient also, where A064989 is the prime shift towards smaller primes. Moreover, because A047802 is defined for every n >= 0, also this sequence is.
From Peter Munn, Aug 13 2020 (Start)
Upper bounds for a(4) and a(5) are:
a(4) <= 195534950863140268380 = A064989(A064989(A064989(20169691981106018776756331))) = A337202(3).
a(5) <= 538938984694949877040715541221415046162838700 = A064989^4((A047802(4)*17*19)/137).
(End)
From David A. Corneth, Aug 21 2020: (Start)
Subsequence of A025487.
Let prime(n)# be the n-th primorial number, A002110(n) = A034386(prime(n)). Then:
a(6) <= 191# * 7#;
a(7) <= 311# * 5#;
a(8) <= 457# * 5#.
(End)
That each term occurs in A025487 follows because (1), the abundancy index of prime(i)^e is larger than that of prime(i+1)^e, that is, sigma(prime(i)^e)/prime(i)^e > sigma(prime(i+1)^e)/prime(i+1)^e, and (2) because the abundancy index of p^(e+d) * q^e is larger than that of p^e * q^(e+d), where p and q are distinct primes, p < q, and e, d > 0. Thus, for any n, we can first find a "prime-factorization compressed version" of it, A071364(n), and then sort the exponents to the non-ascending order with A046523 (and actually, A046523(A071364(n)) = A046523(n), so we need to apply just A046523), to get a term x of A025487, that certainly have the abundancy index >= n [and this inequivalence stays same for their successive prime shifts as well, the abundancy index of A003961(x) being at least that of A003961(n), etc.], and as A046523(n) <= n for all n, it is guaranteed that the least k for which A336835(k) >= n are found from A025487, which is the range of A046523.

Crossrefs

From term a(2) = 120 onward a subsequence of A337386.

Programs

  • PARI
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
    A336835(n) = { my(i=0); while(sigma(n) >= (n+n), i++; n = A003961(n)); (i); };
    A336389(n) = for(i=1,oo,if(A336835(i)>=n,return(i)));

Formula

For n >= 0, A336835(a(n)) >= n.
For all n >= 1, a(n) <= A337202(n-1) [= 2*A246277(A047802(n-1))].
a(n) = A025487(A337477(n)).
a(n) = A108951(A337478(n)).

Extensions

a(4) - a(6) from combined work of David A. Corneth and Peter Munn Aug 13-26 2020

A284001 a(n) = A005361(A283477(n)).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 4, 6, 1, 2, 4, 6, 8, 12, 18, 24, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480, 600, 720, 1, 2, 4, 6, 8, 12, 18, 24, 16, 24, 36, 48, 54, 72, 96, 120, 32, 48, 72, 96, 108, 144, 192, 240, 162, 216, 288, 360, 384, 480
Offset: 0

Views

Author

Antti Karttunen, Mar 18 2017

Keywords

Comments

a(n) is the product of elements of the multiset that covers an initial interval of positive integers with multiplicities equal to the parts of the n-th composition in standard order (graded reverse-lexicographic, A066099). This composition is obtained by taking the set of positions of 1's in the reversed binary expansion of n, prepending 0, taking first differences, and reversing again. For example, the 13th composition is (1,2,1) giving the multiset {1,2,2,3} with product 12, so a(13) = 12. - Gus Wiseman, Apr 26 2020

Crossrefs

Row products of A095684.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Weighted sum is A029931.
- Necklaces are A065609.
- Sum is A070939.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon words are A275692.
- Distinct parts are counted by A334028.

Programs

  • Mathematica
    Table[Times @@ FactorInteger[#][[All, -1]] &[Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e == 1 :> {Times @@ Prime@ Range@ PrimePi@ p, e}] &[Times @@ Prime@ Flatten@ Position[#, 1] &@ Reverse@ IntegerDigits[n, 2]]], {n, 0, 93}] (* Michael De Vlieger, Mar 18 2017 *)
  • PARI
    A005361(n) = factorback(factor(n)[, 2]); \\ From A005361
    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
    A019565(n) = {my(j,v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
    A283477(n) = A108951(A019565(n));
    A284001(n) = A005361(A283477(n));
    
  • Scheme
    (define (A284001 n) (A005361 (A283477 n)))

Formula

a(n) = A005361(A283477(n)).
a(n) = A003963(A057335(n)). - Gus Wiseman, Apr 26 2020
a(n) = A284005(A053645(n)) for n > 0 with a(0) = 1. - Mikhail Kurkov, Jun 05 2021 [verification needed]

A329045 Lexicographically earliest infinite sequence such that a(i) = a(j) => A046523(A329044(i)) = A046523(A329044(j)) for all i, j.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Nov 08 2019

Keywords

Comments

Restricted growth sequence transform of function f(n) = A046523(A329044(n)).
For all i, j:
A305800(i) = A305800(j) => a(i) = a(j),
a(i) = a(j) => A324888(i) = A324888(j),
a(i) = a(j) => A329046(i) = A329046(j).

Crossrefs

Programs

  • PARI
    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; };
    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
    A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
    A324886(n) = A276086(A108951(n));
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A329044(n) = A064989(A324886(n));
    A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); }; \\ From A046523
    v329045 = rgs_transform(vector(up_to, n, A046523(A329044(n))));
    A329045(n) = v329045[n];

A329345 Lexicographically earliest infinite sequence such that a(i) = a(j) => A246277(A329044(i)) = A246277(A329044(j)) for all i, j.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Nov 11 2019

Keywords

Comments

Restricted growth sequence transform of function f(n) = A246277(A329044(n)).
For all i, j:
A305800(i) = A305800(j) => a(i) = a(j),
a(i) = a(j) => A329045(i) = A329045(j),
a(i) = a(j) => A329343(i) = A329343(j),
a(i) = a(j) => A329348(i) = A329348(j),
a(i) = a(j) => A329349(i) = A329349(j).

Crossrefs

Programs

  • PARI
    up_to = 1024;
    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
    A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
    A324886(n) = A276086(A108951(n));
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A329044(n) = A064989(A324886(n));
    A246277(n) = if(1==n, 0, my(f = factor(n), k = primepi(f[1,1])-1); for (i=1, #f~, f[i,1] = prime(primepi(f[i,1])-k)); factorback(f)/2);
    v329345 = rgs_transform(vector(up_to, n, A246277(A329044(n))));
    A329345(n) = v329345[n];

A324895 Largest proper divisor of A276086(n); a(0) = 1.

Original entry on oeis.org

1, 1, 1, 3, 3, 9, 1, 5, 5, 15, 15, 45, 5, 25, 25, 75, 75, 225, 25, 125, 125, 375, 375, 1125, 125, 625, 625, 1875, 1875, 5625, 1, 7, 7, 21, 21, 63, 7, 35, 35, 105, 105, 315, 35, 175, 175, 525, 525, 1575, 175, 875, 875, 2625, 2625, 7875, 875, 4375, 4375, 13125, 13125, 39375, 7, 49, 49, 147, 147, 441, 49, 245, 245, 735, 735, 2205, 245, 1225, 1225
Offset: 0

Views

Author

Antti Karttunen, Mar 30 2019

Keywords

Crossrefs

Programs

  • PARI
    A032742(n) = if(1==n,n,n/vecmin(factor(n)[,1]));
    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; };
    A324895(n) = A032742(A276086(n));

Formula

a(n) = A032742(A276086(n)).
For n >= 0, a(A002110(n)) = 1.
For n >= 1, a(n) = A276086(n) / A053669(n) = A276086(A276151(n)).
For n >= 1, A276085(a(n)) = A276151(n).
For n >= 1, a(A108951(n)) = A324896(n).

A329888 a(n) = A329900(A329602(n)); Heinz number of the even bisection (even-indexed parts) of the integer partition with Heinz number n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 3, 2, 1, 2, 1, 2, 3, 4, 1, 3, 1, 2, 3, 2, 1, 4, 5, 2, 3, 2, 1, 3, 1, 4, 3, 2, 5, 6, 1, 2, 3, 4, 1, 3, 1, 2, 3, 2, 1, 4, 7, 5, 3, 2, 1, 6, 5, 4, 3, 2, 1, 6, 1, 2, 3, 8, 5, 3, 1, 2, 3, 5, 1, 6, 1, 2, 5, 2, 7, 3, 1, 4, 9, 2, 1, 6, 5, 2, 3, 4, 1, 6, 7, 2, 3, 2, 5, 8, 1, 7, 3, 10, 1, 3, 1, 4, 5
Offset: 1

Views

Author

Antti Karttunen, Dec 22 2019

Keywords

Comments

From Gus Wiseman, Aug 05 2021 and Antti Karttunen, Oct 13 2021: (Start)
Also the product of primes at even positions in the weakly decreasing list (with multiplicity) of prime factors of n. For example, the prime factors of 108 are (3,3,3,2,2), with even bisection (3,2), with product 6, so a(108) = 6.
Proof: A108951(n) gives a number with the same largest prime factor (A006530) and its exponent (A071178) as in n, and with each smaller prime p = 2, 3, 5, 7, ... < A006530(n) having as its exponent the partial sum of the exponents of all prime factors >= p present in n (with primes not present in n having the exponent 0). Then applying A000188 replaces each such "partial sum exponent" k with floor(k/2). Finally, A319626 replaces those halved exponents with their first differences (here the exponent of the largest prime present stays intact, because the next larger prime's exponent is 0 in n). It should be easy to see that if prime q is not present in n (i.e., does not divide it), then neither it is present in a(n). Moreover, if the partial sum exponent of q is odd and only one larger than the partial sum exponent of the next larger prime factor of n, then q will not be present in a(n), while in all other cases q is present in a(n). See also the last example.
(End)

Examples

			From _Gus Wiseman_, Aug 15 2021: (Start)
The list of all numbers with image 12 and their corresponding prime factors begins:
  144: (3,3,2,2,2,2)
  216: (3,3,3,2,2,2)
  240: (5,3,2,2,2,2)
  288: (3,3,2,2,2,2,2)
  336: (7,3,2,2,2,2)
  360: (5,3,3,2,2,2)
(End)
The positions from the left are indexed as 1, 2, 3, ..., etc, so e.g., for 240 we pick the second, the fourth and the sixth prime factor, 3, 2 and 2, to obtain a(240) = 3*2*2 = 12. For 288, we similarly pick the second (3), the fourth (2) and the sixth (2) to obtain a(288) = 3*2*2 = 12. - _Antti Karttunen_, Oct 13 2021
Consider n = 11945934 = 2*3*3*3*7*11*13*13*17. Its primorial inflation is A108951(11945934) = 96478365991115908800000 = 2^9 * 3^8 * 5^5 * 7^5 * 11^4 * 13^3 * 17^1. Applying A000188 to this halves each exponent (floored down if the exponent is odd), leaving the factors 2^4 * 3^4 * 5^2 * 7^2 * 11^2 * 13^1 = 2497294800. Then applying A319626 to this number retains the largest prime factor (and its exponent), and subtracts from the exponent of each of the rest of primes the exponent of the next larger prime, so from 2^4 * 3^4 * 5^2 * 7^2 * 11^2 * 13^1 we get 2^(4-4) * 3^(4-2) * 5^(2-2) * 7^(2-2) * 11^(2-1) * 13^1 = 3^2 * 11^1 * 13^1 = 1287 = a(11945934), which is obtained also by selecting every second prime from the list [17, 13, 13, 11, 7, 3, 3, 3, 2] and taking their product. - _Antti Karttunen_, Oct 15 2021
		

Crossrefs

A left inverse of A000290.
Positions of 1's are A008578.
Positions of primes are A168645.
The sum of prime indices of a(n) is A346700(n).
The odd version is A346701.
The odd non-reverse version is A346703.
The non-reverse version is A346704.
The version for standard compositions is A346705, odd A346702.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A001414 adds up prime factors, row sums of A027746.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A346633 adds up the even bisection of standard compositions.
A346698 adds up the even bisection of prime indices.

Programs

  • Mathematica
    Table[Times@@Last/@Partition[Reverse[Flatten[Apply[ConstantArray,FactorInteger[n],{1}]]],2],{n,100}] (* Gus Wiseman, Oct 13 2021 *)
  • PARI
    A329888(n) = A329900(A329602(n));
    
  • PARI
    A329888(n) = if(1==n,n,my(f=factor(n),m=1,p=0); forstep(k=#f~,1,-1,while(f[k,2], m *= f[k,1]^(p%2); f[k,2]--; p++)); (m)); \\ (After Wiseman's new interpretation) - Antti Karttunen, Sep 21 2021

Formula

A108951(a(n)) = A329602(n).
a(n^2) = n for all n >= 1.
a(n) * A346701(n) = n. - Gus Wiseman, Aug 07 2021
A056239(a(n)) = A346700(n). - Gus Wiseman, Aug 07 2021
Antti Karttunen, Sep 21 2021
From Antti Karttunen, Oct 13 2021: (Start)
For all x in A102750, a(x) = a(A253553(x)). (End)

Extensions

Name amended with Gus Wiseman's new interpretation - Antti Karttunen, Oct 13 2021

A331188 Primorial inflation of A052126(n), where A052126(n) = n/(largest prime dividing n).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 4, 6, 2, 1, 4, 1, 2, 6, 8, 1, 12, 1, 4, 6, 2, 1, 8, 30, 2, 36, 4, 1, 12, 1, 16, 6, 2, 30, 24, 1, 2, 6, 8, 1, 12, 1, 4, 36, 2, 1, 16, 210, 60, 6, 4, 1, 72, 30, 8, 6, 2, 1, 24, 1, 2, 36, 32, 30, 12, 1, 4, 6, 60, 1, 48, 1, 2, 180, 4, 210, 12, 1, 16, 216, 2, 1, 24, 30, 2, 6, 8, 1, 72, 210, 4, 6, 2, 30, 32
Offset: 1

Views

Author

Antti Karttunen, Jan 14 2020

Keywords

Comments

The primorial inflation of n, A108951(n), divided by its largest squarefree divisor, which is also its largest primorial divisor.

Crossrefs

Programs

  • PARI
    A002110(n) = prod(i=1,n,prime(i));
    A331188(n) = if(1==n, n, my(f=factor(n)); prod(i=1, #f~, A002110(primepi(f[i, 1]))^(f[i, 2]-(#f~==i))));

Formula

a(n) = A108951(A052126(n)).
a(n) = A003557(A108951(n)).
a(n) = A111701(A108951(n)) = A108951(n) / A002110(A061395(n)).
Other identities. For all >= 1:
A000005(a(n)) = A329382(n) = A005361(A108951(n)).
a(n) mod A117366(n) = A329348(n).

A181817 a(n) is the smallest integer that, when divided by any divisor of A025487(n), yields a member of A025487.

Original entry on oeis.org

1, 2, 4, 12, 8, 24, 16, 48, 360, 32, 144, 96, 720, 64, 288, 192, 1440, 128, 576, 4320, 384, 75600, 1728, 2880, 256, 1152, 8640, 768, 151200, 3456, 5760, 512, 2304, 17280, 1536, 302400, 6912, 129600, 11520, 1024, 51840, 4608, 907200, 20736, 34560, 3072, 604800, 13824, 259200, 23040, 2048
Offset: 1

Views

Author

Matthew Vandermast, Nov 30 2010

Keywords

Comments

A permutation of A181818.

Examples

			For any divisor d of 6 (d = 1, 2, 3, 6), 12/d (12, 6, 4, 2) is always a member of A025487. 12 is the smallest number with this relationship to 6; therefore, since 6 = A025487(4), a(4) = 12.
		

Crossrefs

Formula

If A025487(n) = Product prime(i)^e(i), then a(n) = Product A002110(i)^e(i). I.e., a(n) = A108951(A025487(n)).
If A025487(n) = Product A002110(i)^e(i), then a(n) = Product A006939(i)^e(i).
a(n) = A025487(n) * A181816(n).

A284003 a(n) = A007913(A283477(n)) = A019565(A006068(n)).

Original entry on oeis.org

1, 2, 6, 3, 30, 15, 5, 10, 210, 105, 35, 70, 7, 14, 42, 21, 2310, 1155, 385, 770, 77, 154, 462, 231, 11, 22, 66, 33, 330, 165, 55, 110, 30030, 15015, 5005, 10010, 1001, 2002, 6006, 3003, 143, 286, 858, 429, 4290, 2145, 715, 1430, 13, 26, 78, 39, 390, 195, 65, 130, 2730, 1365, 455, 910, 91, 182, 546, 273, 510510, 255255, 85085, 170170, 17017
Offset: 0

Views

Author

Antti Karttunen, Mar 18 2017

Keywords

Comments

A squarefree analog of A302783. Each term is either a divisor or a multiple of the next one. In contrast to A302033 at each step the previous term can be multiplied (or divided), not just by a single prime, but possibly by a product of several distinct ones, A019565(A000975(k)). E.g., a(3) = 3, a(4) = 2*5*a(3) = 30. - Antti Karttunen, Apr 17 2018

Crossrefs

Programs

Formula

a(n) = A007913(A283477(n)).
Other identities. For all n >= 0:
A048675(a(n)) = A006068(n).
A046523(a(n)) = A284004(n).
It seems that A001222(a(n)) = A209281(n).
a(n) = A019565(A006068(n)) = A302033(A064707(n)). - Antti Karttunen, Apr 16 2018

Extensions

Name amended with a second formula by Antti Karttunen, Apr 16 2018

A071187 Smallest prime factor of number of divisors of n; a(1) = 1.

Original entry on oeis.org

1, 2, 2, 3, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 5, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 7, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 5, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2
Offset: 1

Views

Author

Reinhard Zumkeller, May 15 2002

Keywords

Comments

a(n) = 2 for nonsquare n. - David A. Corneth, Jul 24 2017

Examples

			324 = 18^2 = 2^2 * 3^4 has (2+1)*(4+1) = 3 * 5 = 15 divisors, thus a(324) = A020639(15) = 3. - _Antti Karttunen_, Nov 18 2019
		

Crossrefs

Differs from A329614 for the first time at n=324, where a(324) = 3, while A329614(324) = 5. A329613 gives the positions of differences.

Programs

  • Mathematica
    a[n_] := FactorInteger[DivisorSigma[0, n]][[1, 1]]; Array[a, 90] (* Jean-François Alcover, Oct 01 2016 *)
  • PARI
    A071187(n) = if(1==n, n, my(f = factor(numdiv(n))); vecmin(f[, 1])); \\ Antti Karttunen, Jul 24 2017
    
  • PARI
    first(n) = my(v = vector(n, i, 2)); for(i=1,sqrtint(n), v[i^2] = numdiv(i^2)); v

Formula

a(n) = A020639(A000005(n)).
a(A108951(n)) = A329614(n). - Antti Karttunen, Nov 17 2019
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2. - Amiram Eldar, Jan 15 2024

Extensions

Data section extended up to term a(105) by Antti Karttunen, Nov 17 2019
Previous Showing 91-100 of 158 results. Next