#include <stdio.h>
#include <math.h>
#include <stdbool.h>
#include <stdint.h>

bool is_prime(uint64_t number)
{
        for (uint64_t factor = 2; factor <= number / factor; ++factor)
                if (number % factor == 0)
                        return false;

        return true;
}

int main()
{
        uint64_t sum = 0;
        uint64_t high = 0;

        printf("How high should I go? ");
        scanf("%lld", &high);

        for (uint64_t n = 0; ; ++n)
        {
                uint64_t current_power = (uint64_t)pow(2, n);

                sum += current_power;

                if (is_prime(sum) && sum >= 3)
                {
                        uint64_t perfect_number = sum * current_power;

                        if (perfect_number > high)
                                break;

                        printf("%lld\n", perfect_number);
                }
        }

        return 0;
}

with Ada.Text_IO, Ada.Integer_Text_IO;
use Ada.Text_IO, Ada.Integer_Text_IO;
procedure Perfect_Numbers is
        function Is_Prime(Number : in Integer) return Boolean is
                Factor : Integer := 2;
        begin
                while Factor <= Number / Factor
                loop
                        if Number mod Factor = 0 then
                                return False;
                        end if;
                       
                        Factor := Factor + 1;
                end loop;
               
                return True;
        end Is_Prime;
       
        Sum, Power, High : Integer := 0;
begin
        Put("How high would you like to go? ");
        Get(High);
       
        loop
                declare
                        Current_Power, Perfect_Number : Integer := 0;
                begin
                        Current_Power := 2 ** Power;
                        Sum := Sum + Current_Power;

                        if Is_Prime(Sum) then
                                Perfect_Number := Sum * Current_Power;

                                if Perfect_Number > High then
                                        exit;
                                elsif Sum >= 3 then
                                        Put(Perfect_Number);
                                        New_Line;
                                end if;
                        end if;

                        Power := Power + 1;
                end;
        end loop;
end Perfect_Numbers;

function isprime (n : int) : boolean
    for a : 2 .. floor (sqrt (n))
        if n mod a = 0 then
            result false
        end if
    end for
    result true
end isprime

for p : 2 .. 16
    var mersenne : int := 2 ** p - 1

    if isprime (mersenne) then
        put mersenne * 2 ** (p - 1)
    end if
end for