A130279 Smallest number having exactly n square divisors.
1, 4, 16, 36, 256, 144, 4096, 576, 1296, 2304, 1048576, 3600, 16777216, 36864, 20736, 14400, 4294967296, 32400, 68719476736, 57600, 331776, 9437184, 17592186044416, 129600, 1679616, 150994944, 810000, 921600, 72057594037927936
Offset: 1
Keywords
A298735 Number of odd squares dividing n.
1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 1
Offset: 1
Comments
The smallest integer with exactly m odd square divisors is A357450(m). - Bernard Schott, Oct 03 2022
Examples
a(81) = 3 because 81 has 5 divisors {1, 3, 9, 27, 81} among which 3 are odd squares {1, 9, 81}.
Links
- Antti Karttunen, Table of n, a(n) for n = 1..65537
Crossrefs
Programs
-
Mathematica
nmax = 105; Rest[CoefficientList[Series[Sum[x^(2 k - 1)^2/(1 - x^(2 k - 1)^2), {k, 1, nmax}], {x, 0, nmax}], x]] a[n_] := Length[Select[Divisors[n], IntegerQ[Sqrt[#]] && OddQ[#] &]]; Table[a[n], {n, 1, 105}] f[2, e_] := 1; f[p_, e_] := Floor[e/2] + 1; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 11 2020 *)
-
PARI
a(n)=factorback(apply(e->e\2+1, factor(n/2^valuation(n,2))[, 2])) \\ Rémy Sigrist, Jan 26 2018
Formula
G.f.: Sum_{k>=1} x^((2*k-1)^2)/(1 - x^((2*k-1)^2)).
Multiplicative with a(2^e) = 1 and a(p^e) = floor(e/2) + 1 for p > 2. - Amiram Eldar, Sep 11 2020
Asymptotic mean: lim_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi^2/8 (A111003). - Amiram Eldar, Sep 25 2022
Extensions
Keyword mult added by Rémy Sigrist, Jan 26 2018
A358252 a(n) is the least number with exactly n non-unitary square divisors.
1, 8, 32, 128, 288, 864, 1152, 2592, 4608, 13824, 10368, 20736, 28800, 41472, 64800, 279936, 115200, 331776, 345600, 663552, 259200, 1679616, 518400, 1620000, 1166400, 4860000, 1036800, 17915904, 2073600, 15552000, 6998400, 26873856, 4147200, 53747712, 8294400
Offset: 0
Keywords
Comments
Examples
a(1) = 8 since 8 is the least number that has exactly one non-unitary square divisor, 4.
Links
- Amiram Eldar, Table of n, a(n) for n = 0..176
Crossrefs
Programs
-
Mathematica
f1[p_, e_] := 1 + Floor[e/2]; f2[p_, e_] := 2^(1 - Mod[e, 2]); f[1] = 0; f[n_] := Times @@ f1 @@@ (fct = FactorInteger[n]) - Times @@ f2 @@@ fct; seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s]; seq[21, 10^6]
-
PARI
s(n) = {my(f = factor(n)); prod(i = 1, #f~, 1 + floor(f[i,2]/2)) - 2^sum(i = 1, #f~, 1 - f[i,2]%2);} lista(len, nmax) = {my(v = vector(len), c = 0, n = 1, i); while(c < len && n < nmax, i = s(n) + 1; if(i <= len && v[i] == 0, c++; v[i] = n); n++); v};
A358262 a(n) is the least number with exactly n noninfinitary square divisors.
1, 16, 144, 256, 3600, 1296, 2304, 65536, 129600, 16777216, 32400, 20736, 57600, 331776, 589824, 4294967296, 6350400, 1099511627776, 150994944, 810000, 1587600, 1679616, 518400, 5308416, 2822400, 84934656, 8294400, 26873856, 14745600, 21743271936, 38654705664
Offset: 0
Keywords
Comments
Examples
a(1) = 16 since 16 is the least number with exactly one noninfinitary divisor, 4.
Links
- Amiram Eldar, Table of n, a(n) for n = 0..50
Crossrefs
Programs
-
Mathematica
f1[p_, e_] := 1 + Floor[e/2]; f2[p_, e_] := 2^DigitCount[If[OddQ[e], e - 1, e], 2, 1]; f[1] = 0; f[n_] := Times @@ f1 @@@ (fct = FactorInteger[n]) - Times @@ f2 @@@ fct; seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s]; seq[15, 2*10^7]
-
PARI
s(n) = {my(f = factor(n)); prod(i=1, #f~, 1+f[i,2]\2) - prod(i=1, #f~, 2^hammingweight(if(f[i,2]%2, f[i,2]-1, f[i,2])))}; lista(len, nmax) = {my(v = vector(len), c = 0, n = 1, i); while(c < len && n < nmax, i = s(n) + 1; if(i <= len && v[i] == 0, c++; v[i] = n); n++); v};
A361418 a(n) is the least number with exactly n noninfinitary divisors.
1, 4, 12, 16, 60, 36, 48, 256, 360, 4096, 180, 144, 240, 576, 768, 65536, 2520, 1048576, 12288, 900, 1260, 1296, 720, 2304, 1680, 9216, 2880, 5184, 3840, 147456, 196608, 36864, 27720, 46656, 3145728, 4398046511104, 61440, 3600, 6300, 18014398509481984, 10080, 20736
Offset: 0
Keywords
Comments
Examples
a(1) = 4 since 4 is the least number with exactly one noninfinitary divisor, 2.
Links
- Amiram Eldar, Table of n, a(n) for n = 0..50
Crossrefs
Programs
-
Mathematica
f[1] = 0; f[n_] := DivisorSigma[0, n] - Times @@ Flatten[2^DigitCount[#, 2, 1] & /@ FactorInteger[n][[;; , 2]]]; seq[len_, nmax_] := Module[{s = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = f[n] + 1; If[i <= len && s[[i]] == 0, c++; s[[i]] = n]; n++]; s]; seq[35, 10^7]
-
PARI
s(n) = {my(f = factor(n)); numdiv(f) - prod(i = 1, #f~, 2^hammingweight(f[i,2]));} lista(len, nmax) = {my(v = vector(len), c = 0, n = 1, i); while(c < len && n < nmax, i = s(n) + 1; if(i <= len && v[i] == 0, c++; v[i] = n); n++); v};
Comments
Links
Crossrefs
Programs
PARI
Formula