Форум программистов «Весельчак У»
  *
Добро пожаловать, Гость. Пожалуйста, войдите или зарегистрируйтесь.
Вам не пришло письмо с кодом активации?

  • Рекомендуем проверить настройки временной зоны в вашем профиле (страница "Внешний вид форума", пункт "Часовой пояс:").
  • У нас больше нет рассылок. Если вам приходят письма от наших бывших рассылок mail.ru и subscribe.ru, то знайте, что это не мы рассылаем.
   Начало  
Наши сайты
Помощь Поиск Календарь Почта Войти Регистрация  
 
Страниц: [1]   Вниз
  Печать  
Автор Тема: Что то не то  (Прочитано 10776 раз)
0 Пользователей и 1 Гость смотрят эту тему.
vika
Постоялец

ru
Offline Offline
Пол: Женский

« : 17-05-2012 21:40 » 

Вот код набросала, вроде для небольших степеней и чисел все считает нормально, но почему -то при больших значения степеней, например 2^150 выдает ересь... Самое обидное, что проверить даже нечем, тупо все числа что ли вбивать в онлайн калькулятор.... Переполнение что ли в нем... Взгляните пожалуйста на код и ткине носом, что не так.

Код: (Perl)
#!/usr/bin/perl -w
$| = 1;

$a = 2;
$p = 131;
$d = 2;
$d1 = 5;
$d2 = 13;
open (M, ">> rezultat.txt");

while ($a < $p) {
  $one = $a ** (($p - 1) / $d);
  $otvet = $one % $p;
  $two = $a ** (($p - 1) / $d1);
  $otvet1 = $two % $p;
  $three = $a ** (($p - 1) / $d2);
  $otvet2 = $three % $p;
  print(M "Проверка числа $a дала результаты $otvet ; $otvet1 ; $otvet2\n\n");
  ++$a;
}

close (M);

Это поиск примитивного элемента в поле Галуа, там надо в степень в возводить число и искать остаток от деления...
« Последнее редактирование: 18-05-2012 06:37 от RXL » Записан
RXL
Технический
Администратор

ru
Offline Offline
Пол: Мужской

WWW
« Ответ #1 : 18-05-2012 06:45 » 

Внутренне оператор ** использует функцию pow из стандартной библиотеки Си. Она оперирует типом double.

Почитай мануал по числам в Perl: man perlnumber
Почитай о формате IEEE754 64-бита. У тебя происходит потеря точности.

Попробуй использовать пакет Math::BigFloat или Math::BigInt. По моему тут логичнее Math::BigInt — вычисления у тебя целочисленные. Будь готова к падению производительности на несколько порядков.
« Последнее редактирование: 18-05-2012 12:09 от RXL » Записан

... мы преодолеваем эту трудность без синтеза распределенных прототипов. (с) Жуков М.С.
Страниц: [1]   Вверх
  Печать  
 

Powered by SMF 1.1.21 | SMF © 2015, Simple Machines