With Ada.Text_IO; SELECT `desciption` FROM `db_blogs` WHERE `author` = 'SkG';

Criba de Eratostenes en Ada (ACTUALIZADO)

Viendo el panorama que hay en clase con respecto a la practica de Metodología de la Programación de la Criba de eratostenes me decido a crear esta entrada explicando como se realiza dicha practica.

Este es el código, en los comentarios explico cada parte:

----------------------------------------------------------------
-- Nombre: Criba de Eratostenes
-- Fecha: 14/12/2008
-- Descripción: Creación de una lista que contenga todos los
-- primos menores que el numero que el usuario introduce
-- Requiere: Modulo de manejo de listas (lista.ads y lista.adb)
----------------------------------------------------------------
-- Autor: [SkG]
----------------------------------------------------------------
with Ada.Text_IO; use Ada.Text_IO;
with Ada.Integer_Text_IO; use Ada.Integer_Text_IO;
with Ada.Numerics.Elementary_Functions; use Ada.Numerics.Elementary_Functions;
with Lista; -- Llamada al modulo de manejo de listas que el profesor nos dio

procedure criba is
   ----------------------------------------------------------------
   -- Declaración de tipos y paquetes
   ----------------------------------------------------------------
   -- Declaramos el paquete ListaCriba como una lista de tipo positivo y
   -- la 'usamos' para poder interactuar con ella
   package ListaCriba is new Lista(positive); use ListaCriba;
   -- Creamos un Subtipo 'LCriba' que sea del tipo de las listas
   subtype LCriba is ListaCriba.Tipo;

   ----------------------------------------------------------------
   -- Declaración de funciones y procedimientos
   ----------------------------------------------------------------
   -- Esta funcion la tenemos de practicas anteriores y devuelve si un
   -- numero es primo o no mediante True o False (Booleano)
   function Primo(N,Div:positive) return boolean is
   -- POST: Res = True si not Div(X,N) para todo valor de X|2 < = X < sqrt(N)
   begin
      if Div > integer(sqrt(float(N))) then
         return True;
      else
         return (N mod Div/=0) and then Primo(N,Div+1);
      end if;
   end Primo;

   -- Esta función devuelve una lista que va construyendo según si el
   -- numero es primo o no mediante la función 'Cons' que está en el paquete 'Lista'
   -- El parámetro Inc es para saltar los números en tandas de 2, 4, 2, 4 ...
   -- y no tener que comprobar los pares superiores al 2 que no son primos
   function Criba_Aux(N,Cand,Inc:positive) return LCriba is
   begin
      if Cand > N then
         return Vacia;
      elsif Primo(Cand,2) then
         return Cons(Cand, Criba_Aux(N, Cand+Inc, 6-Inc));
      else
         return Criba_Aux(N, Cand+Inc, 6-Inc);
      end if;
   end Criba_Aux;

   -- Incluimos 2 y 3 en la criba y empezamos desde el 5
   function Criba(N:positive) return LCriba is
   -- POST: Devuelve todos los primos anteriores al numero especificado
   begin
      ----------------------------------------------------------------
      -- NOTA: Si usamos el código que dio el profesor no funciona.
      --           Como no es posible introducir una lista vacía por teclado
      --           no es necesario comprobar si es vacía, si lo intentamos
      --           da error y no podemos compilar.
      --           También hay que añadir las opciones de si es 2 o 3
      --           ya que si pones un numero anterior a 3 devuelve siempre
      --           una lista con 2 y 3.
      ----------------------------------------------------------------
      if N = 2 then
         return Cons(2,Vacia);
      elsif N = 3 then
         return Cons(2,Cons(3,Vacia));
      else
         return Cons(2,Cons(3,Criba_Aux(N,5,2)));
      end if;
   end Criba;

   -- Estas dos funciones SIEMPRE son necesarias si usamos el modulo lista
   -- Esta función nos permite sacar por pantalla la lista. Pondríamos como
   -- queremos ver los valores como si de la parte ejecutiva se tratase
   procedure ImprPrimo(N:positive) is
      begin
      Put(N,0);
   end ImprPrimo;

   -- Esta otra es llamada desde la parte ejecutiva y hace que se 'Escriba' la lista,
   -- para ello llama a la anterior que es la que devuelve los valores.
   procedure ImprLista is new Escribir(ImprPrimo);

   ----------------------------------------------------------------
   -- Declaración de variables y/o constantes
   ----------------------------------------------------------------
   N: natural;

   ----------------------------------------------------------------
   -- Parte ejecutiva
   ----------------------------------------------------------------
begin
   -- Aquí supongo que si sabréis que hacer
   Put("Introduce un numero: "); Get(N); New_Line;
   -- Pedimos que nos Imprima una lista con la Criba hasta el número 'N'
   if N = 0 then
      Put("No hay numeros naturales antes del 0");
   elsif N < 2 then
      Put("No hay numeros primos antes del "); Put(N,0);
   else
      Put("Primos anteriores a ese numero:"); New_Line;
      ImprCriba(Criba(N));
   end if;
end criba;

Espero que os sea útil para comprender la práctica.

Saludos.