cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A016017 Smallest k such that 1/k can be written as a sum of exactly 2 unit fractions in n ways.

This page as a plain text file.
%I A016017 #45 Nov 08 2024 05:51:09
%S A016017 1,2,4,8,6,32,64,12,256,512,24,2048,36,30,16384,32768,96,72,262144,
%T A016017 192,1048576,2097152,60,8388608,216,768,67108864,288,1536,536870912,
%U A016017 1073741824,120,576,8589934592,6144,34359738368,68719476736,180,864
%N A016017 Smallest k such that 1/k can be written as a sum of exactly 2 unit fractions in n ways.
%C A016017 From _Jianing Song_, Aug 30 2021: (Start)
%C A016017 a(n) is the smallest number whose square has exactly 2n-1 divisors.
%C A016017 a(n) is the earliest occurrence of 2n-1 in A048691. (End)
%H A016017 David W. Wilson, <a href="/A016017/b016017.txt">Table of n, a(n) for n = 1..1000</a> (terms 122, 365 and 608 corrected by _Amiram Eldar_, Nov 08 2024)
%F A016017 a(n+1) <= 2^n.
%F A016017 From _Labos Elemer_, May 22 2001: (Start)
%F A016017 a(n) = sqrt(A061283(n)).
%F A016017 a(n) = sqrt(Min{k| A000005(k)=2n-1}).
%F A016017 a((p+1)/2) = 2^((p-1)/2) = 2^A005097(i) if p is the i-th odd prime. [Corrected by _Jianing Song_, Aug 30 2021] (End)
%F A016017 a(n) is the least k such that (tau(k^2) + 1)/2 = n. - _Vladeta Jovovic_, Aug 01 2001
%e A016017 a(1)=1 and a(2)=2 because 1/2 = 1/3 + 1/6 = 1/4 + 1/4.
%e A016017 a(3)=4 because 1/4 = 1/5 + 1/20 = 1/6 + 1/12 = 1/8 + 1/8.
%e A016017 a(4)=8 because 1/8 = 1/9 + 1/72 = 1/10 + 1/40 = 1/12 + 1/24 = 1/16 + 1/16.
%e A016017 a(5)=6 because 1/6 = 1/7 + 1/42 = 1/8 + 1/24 = 1/9 + 1/18 = 1/10 + 1/15 = 1/12 + 1/12.
%t A016017 f[j_, n_] := (Times @@ (j(Last /@ FactorInteger[n]) + 1) + j - 1)/j; t = Table[0, {50}]; Do[a = f[2, n]; If[a < 51 && t[[a]] == 0, t[[a]] = n; Print[{a, n}]], {n, 2^30}] (* _Robert G. Wilson v_, Aug 03 2005 *)
%o A016017 (PARI) a(n) = {k = 1; while (numdiv(k^2) != (2*n-1), k++); return (k); }; \\ _Amiram Eldar_, Jan 07 2019 after _Michel Marcus_ at A071571
%Y A016017 Identical to A071571 shifted right.
%Y A016017 Cf. A000005, A000290, A005408, A005179, A003680, A037992, A055079, A048691.
%K A016017 nonn
%O A016017 1,2
%A A016017 _Robert G. Wilson v_
%E A016017 Entry revised by _N. J. A. Sloane_, Aug 14 2005
%E A016017 Offset corrected by _David W. Wilson_, Dec 27 2018