login
For odd numbers 2n - 1, half the difference between the largest divisor not exceeding the square root, and the least divisor not less than the square root.
6

%I #99 Nov 20 2024 17:51:01

%S 0,1,2,3,0,5,6,1,8,9,2,11,0,3,14,15,4,1,18,5,20,21,2,23,0,7,26,3,8,29,

%T 30,1,4,33,10,35,36,5,2,39,0,41,6,13,44,3,14,7,48,1,50,51,4,53,54,17,

%U 56,9,2,5,0,19,10,63,20,65,6,3,68,69,22,1,12,7,74,75,4,13,78,25,8,81,2,83,0,5,86,9,28,89

%N For odd numbers 2n - 1, half the difference between the largest divisor not exceeding the square root, and the least divisor not less than the square root.

%C We consider 2n-1 which has only odd divisors, so any difference among them is always even.

%C a(n) = 0 if and only if 2n - 1 is a square.

%C From _Charles Kusniec_, Oct 31 2024: (Start)

%C a(n) = 1 iff 2n is a square.

%C a(n) = 2 iff 2n-1 = (2k-1)^2 - 4.

%C a(n) = 3 iff 2n-1 = (2k)^2 - 9.

%C Since 2n-1 is odd, 2n-1 = d*D can be written 2n-1 = (x-y)*(x+y) = x^2 - y^2 with integers x and y, from which a(n) is the smallest possible y, which also means smallest possible x and that being x = A377499(n).

%C Record values occur when 2n-1 is an odd prime (n in A006254), and then a(n) = n-1 (and A377499(n) = n). (End)

%H Seiichi Manyama, <a href="/A219695/b219695.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = (A056737(2n-1))/2 = (A033677(2n-1) - A033676(2n-1))/2.

%F a(n) = sqrt(A377499(n)^2 - (2n-1)). - _Charles Kusniec_, Oct 31 2024

%e For n = 2, consider divisors of 2n - 1 = 3 which are {1, 3}. The least one greater than or equal to sqrt(3) is 3, the largest one less than or equal to sqrt(3) is 1; whence a(2) = (3 - 1)/2 = 1.

%e For n = 14, consider divisors of 2n - 1 = 27 which are {1, 3, 9, 27}. The least one greater than or equal to sqrt(27) is 9, the largest one less than or equal to sqrt(27) is 3; whence a(14) = (9 - 3)/2 = 3.

%e For n = 1, 5, 13, 25, ..., the number 2n - 1 equals the square 1, 9, 25, 49, ...; so the two beforementioned "median divisors" coincide with the square root, and a(n) = 0/2 = 0.

%t Table[(Divisors[n][[(Length[Divisors[n]] - Boole[IntegerQ[Sqrt[n]]])/2 + 1]] - Divisors[n][[(Length[Divisors[n]] + Boole[IntegerQ[Sqrt[n]]])/2]])/2, {n, 1, 199, 2}] (* _Alonso del Arte_, Nov 25 2012, corrected March 21 2024, with help from _Giorgos Kalogeropoulos_ *)

%t A219695[n_] := (d = Divisors[2n - 1]; l = Floor[Length@d/2 + 1]; (d[[l]] - d[[-l]])/2); Array[A219695, 100] (* _Giorgos Kalogeropoulos_, Mar 15 2024 *)

%o (PARI)

%o A056737(n)={n=divisors(n); n[(2+#n)\2]-n[(1+#n)\2]}

%o A219695(n)=A056737(2*n-1)/2 \\ _M. F. Hasler_, Nov 25 2012

%Y Cf. A033676, A033677, A056737, A001844, A016754, A377499.

%Y Cf. A006254 (indices of record highs).

%K nonn,changed

%O 1,3

%A _M. F. Hasler_, Nov 25 2012