mapFolding 0.15.4__py3-none-any.whl → 0.16.0__py3-none-any.whl

mapFolding/algorithms/oeisIDbyFormula.py

@@ -1,113 +1,347 @@
+"""Compute a(n) for an OEIS ID by computing other OEIS IDs.
+
+TODO Implement A178961 for unknown values of A001010
+TODO A223094 For n >= 3: a(n) = n! - Sum_{k=3..n-1} (a(k)*n!/k!) - A000682(n+1). - _Roger Ford_, Aug 24 2024
+TODO A301620 a(n) = Sum_{k=3..floor((n+3)/2)} (A259689(n+1,k)*(k-2)). - _Roger Ford_, Dec 10 2018
+"""
 from functools import cache
-from mapFolding import countFolds
-from mapFolding.algorithms.matrixMeanders import doTheNeedful
+from mapFolding import countFolds, dictionaryOEISMeanders
+from mapFolding.basecamp import A000682, A005316
 
 @cache
 def A000136(n: int) -> int:
-	return n * A000682(n)
+    """
+    Compute A000136(n) as a function of A000682.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A000136 is: "Number of ways of folding a strip of n labeled stamps."
+
+    The domain of A000136 starts at 1, therefore for values of `n` < 1, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 46.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Number of ways of folding a strip of n labeled stamps.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A000136
+    """
+    return n * A000682(n)
 
 def A000560(n: int) -> int:
-	return A000682(n + 1) // 2
+    """
+    Compute A000560(n) as a function of A000682.
 
-def A000682getCurveLocations(n: int) -> dict[int, int]:
-	curveLocationsMAXIMUM: int = 1 << (2 * n + 4)
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A000560 is: "Number of ways of folding a strip of n labeled stamps."
 
-	curveStart: int = 5 - (n & 0b1) * 4
-	listCurveLocations: list[int] = [(curveStart << 1) | curveStart]
+    The domain of A000560 starts at 2, therefore for values of `n` < 2, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 45.
 
-	while listCurveLocations[-1] < curveLocationsMAXIMUM:
-		curveStart = (curveStart << 4) | 0b101
-		listCurveLocations.append((curveStart << 1) | curveStart)
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
 
-	return dict.fromkeys(listCurveLocations, 1)
+    Returns
+    -------
+    a(n) : int
+        Number of ways of folding a strip of n labeled stamps.
 
-@cache
-def A000682(n: int) -> int:
-	return doTheNeedful(n - 1, A000682getCurveLocations(n - 1))
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A000560
+    """
+    return A000682(n + 1) // 2
 
 def A001010(n: int) -> int:
-	"""Formulas.
-
-	a(2n-1) = 2*A007822(n)
-	OddQ[n], 2*A007822[[(n - 1)/2 + 1]]]
-
-	a(2n) = 2*A000682(n+1)
-	EvenQ[n], 2*A000682[[n/2 + 1]]
-	"""
-	if n == 1:
-		foldsTotal = 1
-	elif n & 0b1:
-		foldsTotal = 2 * countFolds(oeisID='A007822', oeis_n=(n - 1)//2 + 1, flow='theorem2Numba')
-	else:
-		foldsTotal = 2 * A000682(n // 2 + 1)
-	return foldsTotal
+    """
+    Compute A001010(n) as a function of A000682 or A007822.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A001010 is: "Number of symmetric foldings of a strip of n stamps."
+
+    The domain of A001010 starts at 1, therefore for values of `n` < 1, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 53.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Number of symmetric foldings of a strip of n stamps.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A001010
+    """
+    if n == 1:
+        foldsTotal = 1
+    elif n & 1:
+        foldsTotal = 2 * countFolds(oeisID='A007822', oeis_n=(n - 1) // 2 + 1, flow='theorem2Numba')
+    else:
+        foldsTotal = 2 * A000682(n // 2 + 1)
+    return foldsTotal
 
 def A001011(n: int) -> int:
-	if n == 0:
-		foldsTotal = 1
-	else:
-		foldsTotal = (A001010(n) + A000136(n)) // 4
-	return foldsTotal
-
-def A005316getCurveLocations(n: int) -> dict[int, int]:
-	if n & 0b1:
-		return {22: 1}
-	else:
-		return {15: 1}
+    """
+    Compute A001011(n) as a function of A000136 and A001010.
 
-@cache
-def A005316(n: int) -> int:
-	return doTheNeedful(n-1, A005316getCurveLocations(n-1))
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A001011 is: "Number of ways to fold a strip of n blank stamps."
+
+    The domain of A001011 starts at 1, therefore for values of `n` < 1, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 46.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Number of ways to fold a strip of n blank stamps.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A001011
+    """
+    if n == 0:
+        foldsTotal = 1
+    else:
+        foldsTotal = (A001010(n) + A000136(n)) // 4
+    return foldsTotal
 
 @cache
 def A005315(n: int) -> int:
-	if n == 1:
-		foldsTotal = 1
-	else:
-		foldsTotal = A005316(2 * n - 1)
-	return foldsTotal
+    """
+    Compute A005315(n) as a function of A005316.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A005315 is: "Closed meandric numbers (or meanders): number of ways a loop can cross a road 2n times."
+
+    The domain of A005315 starts at 0, therefore for values of `n` < 0, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 29.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Closed meandric numbers (or meanders): number of ways a loop can cross a road 2n times.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A005315
+    """
+    if n == 1:
+        foldsTotal = 1
+    else:
+        foldsTotal = A005316(2 * n - 1)
+    return foldsTotal
 
 def A060206(n: int) -> int:
-	return A000682(2 * n + 1)
+    """
+    Compute A060206(n) as a function of A000682.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A060206 is: "Number of rotationally symmetric closed meanders of length 4n+2."
+
+    The domain of A060206 starts at 0, therefore for values of `n` < 0, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 21.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Number of rotationally symmetric closed meanders of length 4n+2.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A060206
+    """
+    return A000682(2 * n + 1)
 
 def A077460(n: int) -> int:
-	"""Formulas.
+    """
+    Compute A077460(n) as a function of A005315, A005316, and A060206.
 
-	a[0] = a[1] = 1;
-	a[n_] := If[OddQ[n], (A005316[[n + 1]] + A005316[[2n]] + A000682[[n]])/4
-	a(2n+1) = (A005315(2n+1) + A005316(2n+1) + A060206(n)) / 4.
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A077460 is: "Number of nonisomorphic ways a loop can cross a road (running East-West) 2n times."
 
-	a(2n) = (A005315(2n) + 2 * A005316(2n)) / 4.
-	(A005316[[2n]] + 2 A005316[[n + 1]])/4];
+    The domain of A077460 starts at 0, therefore for values of `n` < 0, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 21.
 
-	"""
-	if n in {0, 1}:
-		foldsTotal = 1
-	elif n & 0b1:
-		foldsTotal = (A005315(n) + A005316(n) + A060206((n - 1) // 2)) // 4
-	else:
-		foldsTotal = (A005315(n) + 2 * A005316(n)) // 4
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
 
-	return foldsTotal
+    Returns
+    -------
+    a(n) : int
+        Number of nonisomorphic ways a loop can cross a road (running East-West) 2n times.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A077460
+    """
+    if n in {0, 1}:
+        foldsTotal = 1
+    elif n & 1:
+        foldsTotal = (A005315(n) + A005316(n) + A060206((n - 1) // 2)) // 4
+    else:
+        foldsTotal = (A005315(n) + 2 * A005316(n)) // 4
+    return foldsTotal
 
 def A078591(n: int) -> int:
-	return A005315(n) // 2
+    """
+    Compute A078591(n) as a function of A005315.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A078591 is: "Number of nonisomorphic ways a loop can cross a road (running East-West) 2n times."
+
+    The domain of A078591 starts at 0, therefore for values of `n` < 0, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 29.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Number of nonisomorphic ways a loop can cross a road (running East-West) 2n times.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A078591
+    """
+    return A005315(n) // 2
 
 def A178961(n: int) -> int:
-	from mapFolding.oeis import dictionaryOEISMeanders  # noqa: PLC0415
-	A001010valuesKnown: dict[int, int] = dictionaryOEISMeanders['A001010']['valuesKnown']
-	foldsTotal: int = 0
-	for n下i in range(1, n+1):
-		foldsTotal += A001010valuesKnown[n下i]
-	return foldsTotal
+    """
+    Compute A178961(n) as a function of A001010.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A178961 is: "Partial sums of A001010."
+
+    The domain of A178961 starts at 1, therefore for values of `n` < 1, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 53.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Partial sums of A001010.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A178961
+    """
+    A001010valuesKnown: dict[int, int] = dictionaryOEISMeanders['A001010']['valuesKnown']
+    foldsTotal: int = 0
+    for n下i in range(1, n + 1):
+        foldsTotal += A001010valuesKnown[n下i]
+    return foldsTotal
 
 def A223094(n: int) -> int:
-	return A000136(n) - A000682(n + 1)
-# TODO A223094 For n >= 3: a(n) = n! - Sum_{k=3..n-1} (a(k)*n!/k!) - A000682(n+1). - _Roger Ford_, Aug 24 2024
+    """
+    Compute A223094(n) as a function of A000136 and A000682.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A223094 is: "Number of foldings of n labeled stamps in which leaf n is inwards."
+
+    The domain of A223094 starts at 1, therefore for values of `n` < 1, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 44.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Number of foldings of n labeled stamps in which leaf n is inwards.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A223094
+    """
+    return A000136(n) - A000682(n + 1)
 
 def A259702(n: int) -> int:
-	return A000682(n) // 2 - A000682(n - 1)
+    """
+    Compute A259702(n) as a function of A000682.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A259702 is: "Row sums of A259701 except first column."
+
+    The domain of A259702 starts at 2, therefore for values of `n` < 2, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 33.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        Row sums of A259701 except first column.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A259702
+    """
+    return A000682(n) // 2 - A000682(n - 1)
 
 def A301620(n: int) -> int:
-	return A000682(n + 2) - 2 * A000682(n + 1)
-# TODO A301620 a(n) = Sum_{k=3..floor((n+3)/2)} (A259689(n+1,k)*(k-2)). - _Roger Ford_, Dec 10 2018
+    """
+    Compute A301620(n) as a function of A000682.
+
+    *The On-Line Encyclopedia of Integer Sequences* (OEIS) description of A301620 is: "a(n) is the total number of top arches with exactly one covering arch for semi-meanders with n top arches."
+
+    The domain of A301620 starts at 1, therefore for values of `n` < 1, a(n) is undefined. The smallest value of n for which a(n)
+    has not yet been computed is 44.
+
+    Parameters
+    ----------
+    n : int
+        Index (n-dex) for a(n) in the sequence of values. "n" (lower case) and "a(n)" are conventions in mathematics.
+
+    Returns
+    -------
+    a(n) : int
+        a(n) is the total number of top arches with exactly one covering arch for semi-meanders with n top arches.
+
+    Would You Like to Know More?
+    ----------------------------
+    OEIS : webpage
+        https://oeis.org/A301620
+    """
+    return A000682(n + 2) - 2 * A000682(n + 1)

mapFolding/algorithms/zCuzDocStoopidoeisIDbyFormula.py

@@ -0,0 +1,84 @@
+"""Compute a(n) for an OEIS ID by computing other OEIS IDs.
+
+TODO Implement A178961 for unknown values of A001010
+TODO A223094 For n >= 3: a(n) = n! - Sum_{k=3..n-1} (a(k)*n!/k!) - A000682(n+1). - _Roger Ford_, Aug 24 2024
+TODO A301620 a(n) = Sum_{k=3..floor((n+3)/2)} (A259689(n+1,k)*(k-2)). - _Roger Ford_, Dec 10 2018
+"""
+from functools import cache
+from mapFolding import countFolds, dictionaryOEISMeanders
+from mapFolding.basecamp import A000682, A005316
+
+@cache
+def A000136(n: int) -> int:
+	"""A000682"""
+	return n * A000682(n)
+
+def A000560(n: int) -> int:
+	"""A000682"""
+	return A000682(n + 1) // 2
+
+def A001010(n: int) -> int:
+	"""A000682 or A007822"""
+	if n == 1:
+		foldsTotal = 1
+	elif n & 0b1:
+		foldsTotal = 2 * countFolds(oeisID='A007822', oeis_n=(n - 1)//2 + 1, flow='theorem2Numba')
+	else:
+		foldsTotal = 2 * A000682(n // 2 + 1)
+	return foldsTotal
+
+def A001011(n: int) -> int:
+	"""A000136 and A001010"""
+	if n == 0:
+		foldsTotal = 1
+	else:
+		foldsTotal = (A001010(n) + A000136(n)) // 4
+	return foldsTotal
+
+@cache
+def A005315(n: int) -> int:
+	"""A005316"""
+	if n == 1:
+		foldsTotal = 1
+	else:
+		foldsTotal = A005316(2 * n - 1)
+	return foldsTotal
+
+def A060206(n: int) -> int:
+	"""A000682"""
+	return A000682(2 * n + 1)
+
+def A077460(n: int) -> int:
+	"""A005315, A005316, and A060206"""
+	if n in {0, 1}:
+		foldsTotal = 1
+	elif n & 0b1:
+		foldsTotal = (A005315(n) + A005316(n) + A060206((n - 1) // 2)) // 4
+	else:
+		foldsTotal = (A005315(n) + 2 * A005316(n)) // 4
+
+	return foldsTotal
+
+def A078591(n: int) -> int:
+	"""A005315"""
+	return A005315(n) // 2
+
+def A178961(n: int) -> int:
+	"""A001010"""
+	A001010valuesKnown: dict[int, int] = dictionaryOEISMeanders['A001010']['valuesKnown']
+	foldsTotal: int = 0
+	for n下i in range(1, n+1):
+		foldsTotal += A001010valuesKnown[n下i]
+	return foldsTotal
+
+def A223094(n: int) -> int:
+	"""A000136 and A000682"""
+	return A000136(n) - A000682(n + 1)
+
+def A259702(n: int) -> int:
+	"""A000682"""
+	return A000682(n) // 2 - A000682(n - 1)
+
+def A301620(n: int) -> int:
+	"""A000682"""
+	return A000682(n + 2) - 2 * A000682(n + 1)