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.

A094342 Successive record-setters for tau(n+1)*tau(n-1)/tau(n)^2, where tau(n) is the number of divisors of n.

This page as a plain text file.
%I A094342 #6 Aug 28 2019 19:06:13
%S A094342 2,3,5,7,11,17,19,29,41,71,181,239,379,449,701,881,1429,1871,2729,
%T A094342 3079,4159,10529,11969,23561,40699,51679,90271,104651,146719,226799,
%U A094342 244529,252449,388961,403649,825551,906751,1276001,2408561,2648449,3807649,4058209,4406401
%N A094342 Successive record-setters for tau(n+1)*tau(n-1)/tau(n)^2, where tau(n) is the number of divisors of n.
%C A094342 Most terms are primes. These are numbers with few factors which are sandwiched between numbers with many factors. Terms <379 are same as those of A090481.
%H A094342 Amiram Eldar, <a href="/A094342/b094342.txt">Table of n, a(n) for n = 1..61</a>
%e A094342 tau(16)*tau(18)/tau(17)^2 = 5*6/2^2 = 15/2 and this is larger than for any n < 17, so 17 is in the sequence.
%p A094342 f := x -> tau(x-1)*tau(x+1)/tau(x)^2:?print m := 1: A := []: for k from 2 to 10^6 do if f(k) > m then m := f(k): A := [op(A), [k, f(k)]]: fi; od;
%t A094342 s = {}; d1 = 1; d2 = 2; rm = 0; Do[d3 = DivisorSigma[0, n]; r = d1*d3/d2^2; If[r > rm, rm = r; AppendTo[s, n - 1]]; d1 = d2; d2 = d3, {n, 3, 10000}]; s (* _Amiram Eldar_, Aug 28 2019 *)
%Y A094342 Cf. A090481.
%K A094342 easy,nonn
%O A094342 1,1
%A A094342 Isabel C. Lugo (isabel(AT)mit.edu), Jun 04 2004
%E A094342 a(1) = 2 and more terms added by _Amiram Eldar_, Aug 28 2019