{"id":7423,"date":"2022-10-07T06:00:23","date_gmt":"2022-10-07T04:00:23","guid":{"rendered":"http:\/\/www.glc.us.es\/~jalonso\/exercitium\/?p=7423"},"modified":"2022-12-14T14:15:29","modified_gmt":"2022-12-14T12:15:29","slug":"numeros-abundantes-impares","status":"publish","type":"post","link":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/numeros-abundantes-impares\/","title":{"rendered":"N\u00fameros abundantes impares"},"content":{"rendered":"<p>Definir la lista<\/p>\n<pre lang=\"text\">\n  abundantesImpares :: [Integer]\n<\/pre>\n<p>cuyos elementos son los n\u00fameros abundantes impares. Por ejemplo,<\/p>\n<pre lang=\"text\">\n   \u03bb> take 12 abundantesImpares\n   [945,1575,2205,2835,3465,4095,4725,5355,5775,5985,6435,6615]\n<\/pre>\n<p><b>Soluciones<\/b><\/p>\n<p>A continuaci\u00f3n se muestran las <a href=\"#haskell\">soluciones en Haskell<\/a> y las <a href=\"#python\">soluciones en Python<\/a>.<\/p>\n<p><a name=\"haskell\"><\/a><br \/>\n<b>Soluciones en Haskell<\/b><\/p>\n<pre lang=\"haskell\">\nimport Math.NumberTheory.ArithmeticFunctions (sigma)\nimport Test.QuickCheck\n\n\n-- 1\u00aa soluci\u00f3n\n-- ===========\n\nabundantesImpares1 :: [Integer]\nabundantesImpares1 = [x | x <- [1,3..], numeroAbundante1 x]\n\n-- (numeroAbundante n) se verifica si n es un n\u00famero abundante. Por\n-- ejemplo,\n--    numeroAbundante 5  == False\n--    numeroAbundante 12 == True\n--    numeroAbundante 28 == False\n--    numeroAbundante 30 == True\nnumeroAbundante1 :: Integer -> Bool\nnumeroAbundante1 x =\n  x < sumaDivisores1 x - x\n\n-- (sumaDivisores x) es la suma de los divisores de x. Por ejemplo,\n--    sumaDivisores 12                 ==  28\n--    sumaDivisores 25                 ==  31\nsumaDivisores1 :: Integer -> Integer\nsumaDivisores1 n = sum (divisores1 n)\n\n-- (divisores x) es la lista de los divisores de x. Por ejemplo,\n--    divisores 60  ==  [1,5,3,15,2,10,6,30,4,20,12,60]\ndivisores1 :: Integer -> [Integer]\ndivisores1 n = [x | x <- [1..n], n `rem` x == 0]\n\n-- 2\u00aa soluci\u00f3n\n-- ===========\n\nabundantesImpares2 :: [Integer]\nabundantesImpares2 = filter numeroAbundante1 [1,3..]\n\n-- 3\u00aa soluci\u00f3n\n-- ===========\n\n-- Sustituyendo la definici\u00f3n de numeroAbundante1 de las soluciones\n-- anteriores por cada una de las del ejercicio \"N\u00fameros abundantes\"\n-- https:\/\/bit.ly\/3xSlWDU se obtiene una nueva definici\u00f3n de abundantes\n-- impares. La usada en las definiciones anteriores es la menos\n-- eficiente y la que se usa en la siguiente definici\u00f3n es la m\u00e1s eficiente.\n\nabundantesImpares3 :: [Integer]\nabundantesImpares3 = filter numeroAbundante3 [1,3..]\n\nnumeroAbundante3 :: Integer -> Bool\nnumeroAbundante3 x =\n  x < sumaDivisores3 x - x\n\nsumaDivisores3 :: Integer -> Integer\nsumaDivisores3 = sigma 1\n\n-- Comprobaci\u00f3n de equivalencia\n-- ============================\n\n-- La propiedad es\nprop_abundantesImpares :: Positive Int -> Bool\nprop_abundantesImpares (Positive n) =\n  all (== take n abundantesImpares1)\n      [take n abundantesImpares2,\n       take n abundantesImpares3]\n\n-- La comprobaci\u00f3n es\n--    \u03bb> quickCheckWith (stdArgs {maxSize=10}) prop_abundantesImpares\n--    +++ OK, passed 100 tests.\n\n-- Comparaci\u00f3n de eficiencia\n-- =========================\n\n-- La comparaci\u00f3n es\n--    \u03bb> abundantesImpares1 !! 5\n--    4095\n--    (2.07 secs, 841,525,368 bytes)\n--    \u03bb> abundantesImpares2 !! 5\n--    4095\n--    (2.06 secs, 841,443,112 bytes)\n--    \u03bb> abundantesImpares3 !! 5\n--    4095\n--    (0.01 secs, 550,776 bytes)\n<\/pre>\n<p>El c\u00f3digo se encuentra en <a href=\"https:\/\/github.com\/jaalonso\/Exercitium\/blob\/main\/src\/Numeros_abundantes_impares.hs\">GitHub<\/a>.<\/p>\n<p><a name=\"python\"><\/a><br \/>\n<b>Soluciones en Python<\/b><\/p>\n<pre lang=\"python\">\nfrom timeit import Timer, default_timer\n\nfrom hypothesis import given\nfrom hypothesis import strategies as st\nfrom sympy import divisor_sigma\n\n# 1\u00aa soluci\u00f3n\n# ===========\n\ndef abundantesImpares1(n: int) -> list[int]:\n    return [x for x in range(1, n, 2) if numeroAbundante1(x)]\n\n# divisores(n) es la lista de los divisores del n\u00famero n. Por ejemplo,\n#    divisores(30)  ==  [1,2,3,5,6,10,15,30]\ndef divisores1(n: int) -> list[int]:\n    return [x for x in range(1, n + 1) if n % x == 0]\n\n# sumaDivisores(x) es la suma de los divisores de x. Por ejemplo,\n#    sumaDivisores(12)                ==  28\n#    sumaDivisores(25)                ==  31\ndef sumaDivisores1(n: int) -> int:\n    return sum(divisores1(n))\n\n# numeroAbundante(n) se verifica si n es un n\u00famero abundante. Por\n# ejemplo,\n#    numeroAbundante(5)  == False\n#    numeroAbundante(12) == True\n#    numeroAbundante(28) == False\n#    numeroAbundante(30) == True\ndef numeroAbundante1(x: int) -> bool:\n    return x < sumaDivisores1(x) - x\n\n# 2\u00aa soluci\u00f3n\n# ===========\n\ndef abundantesImpares2(n: int) -> list[int]:\n    return list(filter(numeroAbundante1, range(1, n, 2)))\n\n# 3\u00aa soluci\u00f3n\n# ===========\n#\n# Sustituyendo la definici\u00f3n de numeroAbundante1 de las soluciones\n# anteriores por cada una de las del ejercicio \"N\u00fameros abundantes\"\n# https:\/\/bit.ly\/3xSlWDU se obtiene una nueva definici\u00f3n de abundantes\n# impares. La usada en las definiciones anteriores es la menos\n# eficiente y la que se usa en la siguiente definici\u00f3n es la m\u00e1s eficiente.\n\ndef abundantesImpares3(n: int) -> list[int]:\n    return list(filter(numeroAbundante3, range(1, n, 2)))\n\ndef sumaDivisores3(n: int) -> int:\n    return divisor_sigma(n, 1)\n\ndef numeroAbundante3(x: int) -> bool:\n    return x < sumaDivisores3(x) - x\n\n# Comprobaci\u00f3n de equivalencia\n# ============================\n\n# La propiedad es\n@given(st.integers(min_value=1, max_value=1000))\ndef test_abundantesImpares(n: int) -> None:\n    r = abundantesImpares1(n)\n    assert abundantesImpares2(n) == r\n    assert abundantesImpares3(n) == r\n\n# La comprobaci\u00f3n es\n#    src> poetry run pytest -q numeros_abundantes_impares.py\n#    1 passed in 1.42s\n\n# Comparaci\u00f3n de eficiencia\n# =========================\n\ndef tiempo(e: str) -> None:\n    \"\"\"Tiempo (en segundos) de evaluar la expresi\u00f3n e.\"\"\"\n    t = Timer(e, \"\", default_timer, globals()).timeit(1)\n    print(f\"{t:0.2f} segundos\")\n\n# La comparaci\u00f3n es\n#    >>> tiempo('abundantesImpares1(10000)[5]')\n#    1.25 segundos\n#    >>> tiempo('abundantesImpares2(10000)[5]')\n#    1.22 segundos\n#    >>> tiempo('abundantesImpares3(10000)[5]')\n#    0.33 segundos\n<\/pre>\n<p>El c\u00f3digo se encuentra en <a href=\"https:\/\/github.com\/jaalonso\/Exercitium-Python\/blob\/main\/src\/numeros_abundantes_impares.py\">GitHub<\/a>.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Definir la lista abundantesImpares :: [Integer] cuyos elementos son los n\u00fameros abundantes impares. Por ejemplo, \u03bb> take 12 abundantesImpares [945,1575,2205,2835,3465,4095,4725,5355,5775,5985,6435,6615] Soluciones A continuaci\u00f3n se muestran las soluciones en Haskell y las soluciones en Python. Soluciones en Haskell import Math.NumberTheory.ArithmeticFunctions (sigma) import Test.QuickCheck &#8212; 1\u00aa soluci\u00f3n &#8212; =========== abundantesImpares1 :: [Integer] abundantesImpares1 = [x | x&#8230;<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"jetpack_post_was_ever_published":false,"_kad_post_transparent":"","_kad_post_title":"","_kad_post_layout":"","_kad_post_sidebar_id":"","_kad_post_content_style":"","_kad_post_vertical_padding":"","_kad_post_feature":"","_kad_post_feature_position":"","_kad_post_header":false,"_kad_post_footer":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"footnotes":"","_jetpack_memberships_contains_paid_content":false},"categories":[581],"tags":[],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts\/7423"}],"collection":[{"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/comments?post=7423"}],"version-history":[{"count":2,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts\/7423\/revisions"}],"predecessor-version":[{"id":7676,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts\/7423\/revisions\/7676"}],"wp:attachment":[{"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/media?parent=7423"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/categories?post=7423"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/tags?post=7423"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}