{"id":3087,"date":"2017-03-14T06:00:11","date_gmt":"2017-03-14T04:00:11","guid":{"rendered":"http:\/\/www.glc.us.es\/~jalonso\/exercitium\/?p=3087"},"modified":"2021-04-25T16:13:31","modified_gmt":"2021-04-25T14:13:31","slug":"distribucion-de-diferencias-de-digitos-consecutivos-de-pi-2017","status":"publish","type":"post","link":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/distribucion-de-diferencias-de-digitos-consecutivos-de-pi-2017\/","title":{"rendered":"Distribuci\u00f3n de diferencias de d\u00edgitos consecutivos de pi"},"content":{"rendered":"<p>La distribuci\u00f3n de las diferencias de los d\u00edgitos consecutivos para los 18 primeros d\u00edgitos de pi se calcula como sigue: los primeros 18 d\u00edgitos de pi son<\/p>\n<pre lang=\"text\">\n   3, 1, 4, 1, 5, 9, 2, 6, 5, 3, 5, 8, 9, 7, 9, 3, 2, 3\n<\/pre>\n<p>Las diferencias de sus elementos consecutivos es<\/p>\n<pre lang=\"text\">\n   2, -3, 3, -4, -4, 7, -4, 1, 2, -2, -3, -1, 2, -2, 6, 1, -1\n<\/pre>\n<p>y la distribuci\u00f3n de sus frecuencias en el intervalo [-9,9] es<\/p>\n<pre lang=\"text\">\n   0, 0, 0, 0, 0, 3, 2, 2, 2, 0, 2, 3, 1, 0, 0, 1, 1, 0, 0\n<\/pre>\n<p>es decir, el desde el -9 a -5 no aparecen, el -4 aparece 3 veces, el -2 aparece 2 veces y as\u00ed sucesivamente.<\/p>\n<p>Definir las funciones<\/p>\n<pre lang=\"text\">\n   distribucionDDCpi :: Int -> [Int]\n   graficas :: [Int] -> FilePath -> IO ()\n<\/pre>\n<p>tales que<\/p>\n<ul>\n<li>(distribucionDDCpi n) es la distribuci\u00f3n de las diferencias de los d\u00edgitos consecutivos para los primeros n d\u00edgitos de pi. Por ejemplo,<\/li>\n<\/ul>\n<pre lang=\"text\">\n   \u03bb> distribucionDDCpi 18\n   [0,0,0,0,0,3,2,2,2,0,2,3,1,0,0,1,1,0,0]\n   \u03bb> distribucionDDCpi 100\n   [1,2,1,7,7,7,6,5,8,6,7,14,4,9,3,6,4,1,0]\n   \u03bb> distribucionDDCpi 200\n   [3,6,2,13,14,12,11,12,15,17,15,19,11,17,8,13,9,2,0]\n   \u03bb> distribucionDDCpi 1000\n   [16,25,23,44,57,61,55,75,92,98,80,88,64,65,42,54,39,14,8]\n   \u03bb> distribucionDDCpi 5000\n   [67,99,130,196,245,314,361,391,453,468,447,407,377,304,242,221,134,97,47]\n<\/pre>\n<ul>\n<li>(graficas ns f) dibuja en el fichero f las gr\u00e1ficas de las distribuciones de las diferencias de los d\u00edgitos consecutivos para los primeros n d\u00edgitos de pi, para n en ns. Por ejemplo, al evaluar (graficas [100,250..4000] \u00abdistribucionDDCpi.png\u00bb se escribe en el fichero \u00abdistribucionDDCpi.png\u00bb la siguiente gr\u00e1fica<br \/>\n<a href=\"https:\/\/i0.wp.com\/www.glc.us.es\/~jalonso\/exercitium\/wp-content\/uploads\/2017\/03\/Distribucion_de_diferencias_de_digitos_consecutivos_de_pi.png\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/i0.wp.com\/www.glc.us.es\/~jalonso\/exercitium\/wp-content\/uploads\/2017\/03\/Distribucion_de_diferencias_de_digitos_consecutivos_de_pi.png?resize=640%2C480\" alt=\"Distribucion_de_diferencias_de_digitos_consecutivos_de_pi\" width=\"640\" height=\"480\" class=\"aligncenter size-full wp-image-3088\" srcset=\"https:\/\/i0.wp.com\/www.glc.us.es\/~jalonso\/exercitium\/wp-content\/uploads\/2017\/03\/Distribucion_de_diferencias_de_digitos_consecutivos_de_pi.png?w=640&amp;ssl=1 640w, https:\/\/i0.wp.com\/www.glc.us.es\/~jalonso\/exercitium\/wp-content\/uploads\/2017\/03\/Distribucion_de_diferencias_de_digitos_consecutivos_de_pi.png?resize=300%2C225&amp;ssl=1 300w, https:\/\/i0.wp.com\/www.glc.us.es\/~jalonso\/exercitium\/wp-content\/uploads\/2017\/03\/Distribucion_de_diferencias_de_digitos_consecutivos_de_pi.png?resize=100%2C75&amp;ssl=1 100w, https:\/\/i0.wp.com\/www.glc.us.es\/~jalonso\/exercitium\/wp-content\/uploads\/2017\/03\/Distribucion_de_diferencias_de_digitos_consecutivos_de_pi.png?resize=150%2C112&amp;ssl=1 150w\" sizes=\"(max-width: 640px) 100vw, 640px\" data-recalc-dims=\"1\" \/><\/a><\/li>\n<\/ul>\n<p><strong>Nota:<\/strong> Se puede usar la librer\u00eda <a href=\"https:\/\/hackage.haskell.org\/package\/numbers-3000.2.0.1\/docs\/Data-Number-CReal.html\">Data.Number.CReal<\/a>.<\/p>\n<h4>Soluciones<\/h4>\n<pre lang=\"haskell\">\nimport Data.Number.CReal\nimport Graphics.Gnuplot.Simple\nimport Data.Array\n\n--    \u03bb> digitosPi 18\n--    [3,1,4,1,5,9,2,6,5,3,5,8,9,7,9,3,2,3]\ndigitosPi :: Int -> [Int]\ndigitosPi n = init [read [c] | c <- (x:xs)]\n  where (x:_:xs) = showCReal n pi\n\n--    \u03bb> diferenciasConsecutivos (digitosPi 18)\n--    [2,-3,3,-4,-4,7,-4,1,2,-2,-3,-1,2,-2,6,1,-1]\ndiferenciasConsecutivos :: Num a => [a] -> [a]\ndiferenciasConsecutivos xs =\n  zipWith (-) xs (tail xs)\n\ndistribucionDDCpi :: Int -> [Int]\ndistribucionDDCpi =\n  distribucion . diferenciasConsecutivos . digitosPi\n  where distribucion xs =\n          elems (accumArray (+) 0 (-9,9) (zip xs (repeat 1)))\n\ngraficas :: [Int] -> FilePath -> IO ()\ngraficas ns f = \n  plotLists [Key Nothing, PNG f]\n            [puntos n | n <- ns]\n  where puntos :: Int -> [(Int,Int)]\n        puntos n = zip [-9..9] (distribucionDDCpi n)\n<\/pre>\n","protected":false},"excerpt":{"rendered":"<p>La distribuci\u00f3n de las diferencias de los d\u00edgitos consecutivos para los 18 primeros d\u00edgitos de pi se calcula como sigue: los primeros 18 d\u00edgitos de pi son 3, 1, 4, 1, 5, 9, 2, 6, 5, 3, 5, 8, 9, 7, 9, 3, 2, 3 Las diferencias de sus elementos consecutivos es 2, -3, 3,&#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":[7],"tags":[249,8,245,376,285,42,11,375,95,49,378,45,9,76],"jetpack_featured_media_url":"","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts\/3087"}],"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=3087"}],"version-history":[{"count":6,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts\/3087\/revisions"}],"predecessor-version":[{"id":3120,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/posts\/3087\/revisions\/3120"}],"wp:attachment":[{"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/media?parent=3087"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/categories?post=3087"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.glc.us.es\/~jalonso\/exercitium\/wp-json\/wp\/v2\/tags?post=3087"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}