A217584
Numbers k such that d(k^2)/d(k) is an integer, where d(k) is the number of divisors of k.
Original entry on oeis.org
1, 144, 324, 400, 784, 1936, 2025, 2500, 2704, 3600, 3969, 4624, 5625, 5776, 7056, 8100, 8464, 9604, 9801, 13456, 13689, 15376, 15876, 17424, 19600, 21609, 21904, 22500, 23409, 24336, 26896, 29241, 29584, 30625, 35344, 39204, 41616, 42849, 44944, 48400, 51984
Offset: 1
d(1^2)/d(1) = d(1)/d(1) = 1 an integer, so 1 belongs to the sequence.
144^2 has 45 divisors: 1, 2, 3, 4, 6, 8, 9, 12, ..., 20736, while 144 has 15 divisors: 1, 2, 3, 4, 6, 8, 9, 12, ..., 144; 45/15 = 3 and so 144 is in the sequence.
- Marcin E. Kuczma, International Mathematical Olympiads, 1986-1999, The Mathematical Association of America, 2003, pages 134-135.
Cf.
A339055 (values taken by d(a(n)^2)/d(a(n))),
A339056 (smallest k such that d(k^2)/d(k) = n-th odd).
-
Select[Range[1000], IntegerQ[DivisorSigma[0, #^2]/DivisorSigma[0, #]] &] (* Alonso del Arte, Oct 07 2012 *)
Select[Range[228]^2, Divisible[DivisorSigma[0, #^2], DivisorSigma[0, #]] &] (* Amiram Eldar, May 23 2020 *)
-
dn2dn(n)= {for (i=1, n, if (denominator(numdiv(i^2)/numdiv(i))==1, print1(i,", ");););}
A190106
Numbers with prime factorization p^2*q^3*r^3 where p, q, and r are distinct primes.
Original entry on oeis.org
5400, 9000, 10584, 13500, 24696, 26136, 36504, 37044, 49000, 62424, 68600, 77976, 95832, 114264, 121000, 143748, 158184, 165375, 169000, 171500, 181656, 207576, 231525, 237276, 266200, 289000, 295704, 332024, 353736, 361000, 363096
Offset: 1
-
f[n_]:=Sort[Last/@FactorInteger[n]]=={2,3,3};Select[Range[500000],f]
-
list(lim)=my(v=List(),t1,t2);forprime(p=2, (lim\4)^(1/6), t1=p^3;forprime(q=p+1, (lim\t1)^(1/3), t2=t1*q^3;forprime(r=2, sqrt(lim\t2), if(p==r||q==r, next);listput(v,t2*r^2)))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011
-
from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A190106(n):
def bisection(f,kmin=0,kmax=1):
while f(kmax) > kmax: kmax <<= 1
kmin = kmax >> 1
while kmax-kmin > 1:
kmid = kmax+kmin>>1
if f(kmid) <= kmid:
kmax = kmid
else:
kmin = kmid
return kmax
def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=integer_nthroot(x//r**2,3)[0])))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(isqrt(x)+1))+sum(primepi(integer_nthroot(x//p**5,3)[0]) for p in primerange(integer_nthroot(x,5)[0]+1))-primepi(integer_nthroot(x,8)[0])
return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025
A175752
Numbers with 45 divisors.
Original entry on oeis.org
3600, 7056, 8100, 15876, 17424, 19600, 20736, 22500, 24336, 39204, 41616, 48400, 51984, 54756, 67600, 76176, 86436, 93636, 94864, 99225, 104976, 115600, 116964, 121104, 122500, 132496, 138384, 144400, 147456, 160000, 171396, 197136
Offset: 1
Showing 1-3 of 3 results.
Comments