%I #38 Nov 11 2019 09:23:30
%S 6,8,8,15,21,11,13,27,16,35,16,27,16,27,55,27,16,16,16,65,27,16,77,16,
%T 85,16,29,91,31,16,95,16,37,115,16,119,16,41,43,133,16,47,16,143,125,
%U 16,125,16,53,161,16,59,16,61,125,187,16,67,16,203,125,16,209,71,16,125
%N The 5th divisor of numbers with 25 divisors.
%C n=1 means the first number that has 25 divisors (1296), 6 is the 5th divisor of 1296. The second number with 25 divisors is 10000 and its 5th divisor is 8
%C This is one example of such a sequence where the divisor index is the square root of the total number of divisors (self included).
%C Other examples would be the 6th divisor of numbers with 36 divisors, 7th divisor of numbers with 49 divisors, etc.
%C Choice of the square root is arbitrary.
%C All but 16 primes {2, 3, 5, 7, 17, 19, 23, 83, 89, 97, 101, 103, 107, 109, 113} are in this sequence; p^3 and p^4 are in this sequence for all prime p; pq is in this sequence for all prime p and q with p < q < p^2. No other terms are members. - _Charles R Greathouse IV_, Nov 28 2011
%H Reinhard Zumkeller, <a href="/A135581/b135581.txt">Table of n, a(n) for n = 1..10000</a>
%H Laurens Lapré, <a href="http://laurenslapre.nl/lapre_006.htm">Natural division</a>.
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Divisor_function">Divisor function</a>
%e a(1) = 6 because 6 is the 5th divisor of 1296 and 1296 is the first number with 25 divisors.
%e a(2) = 8 because 8 is the 5th divisor of 10000 and 10000 is the second number with 25 divisors.
%t upto=10^10;With[{max1=Ceiling[Power[upto, (4)^-1]],max2=Ceiling[ Power[ upto, (24)^-1]]},Take[Divisors[#][[5]]&/@Select[Union[Join[ Range[ max2]^24, Times@@@(Subsets[Range[max1],{2}]^4)]],DivisorSigma[0,#] == 25&], Ceiling[max1/4]]] (* _Harvey P. Dale_, Nov 25 2011 *)
%o (Haskell)
%o a135581 n = [d | d <- [1..], a137488 n `mod` d == 0] !! 4
%o -- _Reinhard Zumkeller_, Nov 29 2011
%Y Cf. A005179, A137488, A200722, A201266.
%K nonn,look,easy
%O 1,1
%A G. H. Ens (GerardEns(AT)gmail.com), Feb 24 2008
%E Corrected and extended by _R. J. Mathar_, Apr 21 2008. The original entries were wrong from the 16th term onwards.
|