OFFSET
1,2
COMMENTS
a(n) is least term of A054994 with exactly n divisors. - Ray Chandler, Jan 05 2012
From Jianing Song, Apr 24 2021: (Start)
a(n) is the smallest k such that A004018(k) = 4n.
Also a(n) is the smallest index of n in A002654.
a(n) is the smallest term in A004613 that has exactly n divisors.
This is a subsequence of A054994. (End)
LINKS
Ray Chandler, Table of n, a(n) for n = 1..1432 (a(1433) exceeds 1000 digits).
FORMULA
A000446(n) = min(a(2n-1), a(2n)) for n > 1.
A124980(n) = min(a(2n-1), a(2n)).
A016032(n) = min(2*a(2n-1), a(2n), a(2n+1)).
A093195(n) = min(a(2n), a(2n+1)).
From Jianing Song, Apr 24 2021: (Start)
If the factorization of n into primes is n = Product_{i=1..r} p_i with p_1 >= p_2 >= ... >= p_r, then a(n) <= (q_1)^((p_1)-1) * (q_2)^((p_2)-1) * ... * (q_r)^((p_r)-1), where q_1 < q_2 < ... < q_r are the first r primes congruent to 1 modulo 4. The smallest n such that the equality does not hold is n = 16.
a(n) <= 5^(n-1) for all n, where the equality holds if and only if n = 1 or n is a prime.
a(p*q) = 5^(p-1) * 13^(q-1) for primes p >= q. (End)
EXAMPLE
4225 = 5^2 * 13^2 is the smallest number all of whose prime factors are congruent to 1 modulo 4 with exactly 9 divisors, so a(9) = 4225. - Jianing Song, Apr 24 2021
MATHEMATICA
(* This program is not convenient to compute huge terms - A054994 is assumed to be computed with maxTerm = 10^16 *) a[n_] := Catch[ For[k = 1, k <= Length[A054994], k++, If[DivisorSigma[0, A054994[[k]]] == n, Throw[A054994[[k]]]]]]; Table[a[n], {n, 1, 28}] (* Jean-François Alcover, Jan 21 2013, after Ray Chandler *)
PROG
(PARI) primelist(d, r, l) = my(v=vector(l), i=0); if(l>0, forprime(p=2, oo, if(Mod(p, d)==r, i++; v[i]=p; if(i==l, break())))); v
prodR(n, maxf)=my(dfs=divisors(n), a=[], r); for(i=2, #dfs, if( dfs[i]<=maxf, if(dfs[i]==n, a=concat(a, [[n]]), r=prodR(n/dfs[i], min(dfs[i], maxf)); for(j=1, #r, a=concat(a, [concat(dfs[i], r[j])]))))); a
A018782(n)=my(pf=prodR(n, n), a=1, b, v=primelist(4, 1, bigomega(n))); for(i=1, #pf, b=prod(j=1, length(pf[i]), v[j]^(pf[i][j]-1)); if(b<a || i==1, a=b)); a \\ Jianing Song, Apr 25 2021, following R. J. Mathar's program for A005179.
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved