A328129 Finite cardinalities of equivalence classes of real intervals with respect to the symmetric transitive closure of R(x,y) = "x is an integer multiple of y".
1, 2, 4, 7, 12, 24, 41, 58, 75, 92, 109, 214, 319, 424, 529, 634, 1176, 1718, 2260, 2802, 3344, 3886, 4428, 4970, 9411
Offset: 1
Examples
The real interval [1.0, 2.0] has equivalence classes including {1.5} and {1.0, 2.0}, so |{1.5}| = 1 and |{1.0, 2.0}| = 2 are in the sequence. To establish whether 3 is in the sequence, let x be the smallest member (by absolute value) of an equivalence class E on an interval I. E has cardinality greater than 1 only if 2*x is in I. If 3*x is not in I, there are no other multiples or submultiples of {x, 2*x} in I, so E has cardinality less than 3. Otherwise 3*x is in I, 3/2*x is in I, so {x, 2*x, 3*x, 3/2*x} is a subset of E and E has cardinality greater than 3. So 3 is not in the sequence. Equivalence classes of cardinality 7 take the form {x, 2*x, 3*x, 3/2*x, 4*x, 4/3*x, 8/3*x} for real x <> 0. Equivalence classes of cardinality 12 take the form {x, 2*x, 3*x, 3/2*x, 4*x, 4/3*x, 8/3*x, 9/2*x, 9/4*x, 9/8*x, 27/8*x, 27/16*x} for real x <> 0. In the table below the columns are as follows. |E|: cardinality of equivalence class; I_min: smallest interval with lower bound 1 that has such an equivalence class; I_max: largest interval with lower bound 1 that has such an equivalence class; u_max: upper bound of I_max in decimal notation. |E| I_min I_max u_max 1 [1, 1] (1, 4) 4.0 2 [1, 2] (1, 9/2) 4.5 4 [1, 3] (1, 16/3) 5.3333333333... 7 [1, 4] (1, 81/16) 5.0625 12 [1, 9/2] (1, 50/9) 5.5555555555... 24 [1, 5] (1, 6561/1280) 5.12578125 41 [1, 81/16] (1, 531441/102400) 5.189853515625 58 [1, 6561/1280] (1, 3^16/(2^16*5^3)) 5.2547266845... ... 214 [1, 16/3] (1, 2^52*5^8/3^43) 5.3592727015... ... oo (1, 27/5) not defined not defined
Links
- Peter Munn, The disappearance of finite classes
- Eric Weisstein's World of Mathematics, Closure
- Eric Weisstein's World of Mathematics, Equivalence Class
- Wikipedia, Transversal (combinatorics)
Programs
-
PARI
{ my (points = Set([1]), newmax); print (1," ",matsize(points)[2]); for (n=2,25, newmax = vecmin(setminus(setbinop((x,y)->x*y,[2,3,5],points),points)); points = setunion(setbinop((x,y)->x/y,[newmax],points),points); print (n," ",matsize(points)[2]); )} \\ Peter Munn, Jun 29 2022
Formula
a(n) = |S_n| where S_n = S_(n-1) U {min( {2*x, 3*x, 5*x : x in S_(n-1)} \ S_(n-1)) / y : y in S_(n-1)} n >= 2, S_1 = {1.0}.
Extensions
a(21)-a(25) from Peter Munn, Jun 29 2022
Comments