OFFSET
1,2
COMMENTS
This sequence arises in the analytically obtained strict upper bound of the set of equivalent resistances formed by any conceivable network (series/parallel or bridge, or non-planar) of n equal resistors. Consequently it provides an strict upper bound of the sequences: A048211, A153588, A174283, A174284, A174285 and A174286. This sequence provides a better strict upper bound than A176499 but is harder to compute. [Corrected by Antoine Mathys, May 07 2019]
From Hugo Pfoertner, Jan 24 2021: (Start)
The claim that this sequence is a strict upper bound for the number of representable resistance values of any conceivable network is wrong. It only applies to purely serial-parallel networks, but it already fails when bridges are allowed, as described in A174283. Even more so if arbitrary nonplanar networks are allowed as in A337517. See the linked illustrations of the respective quotients.
LINKS
Antoine Mathys, Table of n, a(n) for n = 1..40
Antoni Amengual, The intriguing properties of the equivalent resistances of n equal resistors combined in series and in parallel, American Journal of Physics, 68(2), 175-179 (February 2000).
Sameen Ahmed Khan, The bounds of the set of equivalent resistances of n equal resistors combined in series and in parallel, arXiv:1004.3346v1 [physics.gen-ph], (20 April 2010).
Sameen Ahmed Khan, Mathematica notebook
Hugo Pfoertner, Ratio for series-parallel networks, Plot2 of A048211(n)/a(n).
Hugo Pfoertner, Ratio for networks with bridges, Plot2 of A174283(n)/a(n).
Hugo Pfoertner, Ratio for arbitrary networks, Plot2 of A337517(n)/a(n).
EXAMPLE
n = 5, I = [1/5, 1], m = Fibonacci(5 + 1) = 8, Farey(8) = 23, Farey(8; I) = 19
MATHEMATICA
a[n_ /; n<4] := 2^(n-1); a[n_] := Module[{m = Fibonacci[n+1], v}, v = Reap[ Do[Sow[j/i], {i, n+1, m}, {j, 1, (i-1)/n}]][[2, 1]]; Total[ EulerPhi[ Range[m]]] - Length[v // Union]];
Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 1, 23}] (* Jean-François Alcover, Aug 30 2018, after Antoine Mathys *)
PROG
(PARI) farey(n) = sum(i=1, n, eulerphi(i)) + 1;
a(n) = my(m=fibonacci(n + 1), count=0); for(b=n+1, m, for(a=1, (b-1)/n, if(gcd(a, b)==1, count++))); farey(m) - 1 - count; \\ Antoine Mathys, May 07 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Sameen Ahmed Khan, Apr 21 2010
EXTENSIONS
a(19)-a(27) from Antoine Mathys, Aug 10 2018
a(28)-a(31) from Antoine Mathys, May 07 2019
STATUS
approved