русский
Germany.ruForen → Архив Досок→ Programmierung

Поиск простого числа

262  1 2 3 alle
taksos завсегдатай08.11.06 07:54
taksos
08.11.06 07:54 
Задача: найти простое число по его порядковому номеру в последовательности простых чисел, вот такая прога. Не хочет работать. Пожалуйста, тыкните в ошибку!
program prostoe;
const n=30000;
var
a:array[1..n] of integer;
chs:integer;
sum,num,i,b:integer;
function prov (c:integer):boolean;
var s:byte;
del:integer;
begin
s:=0;
for del:=2 to c do
begin
if c mod del=0 then inc(s);
if s<>0 then break;
end;
if s=0 then prov:=true;
end;
begin
readln(chs);
num:=1;
for b:=1 to chs do
begin
if prov(num)=true then a:=num;
inc(num);
end;
readln;
end.
Проблема в том, что массив при пошаговой проверке не заполняется, только первый элемент становится равным единице и все.
#1 
Murr коренной житель08.11.06 09:57
Murr
NEW 08.11.06 09:57 
in Antwort taksos 08.11.06 07:54
if prov(num)=true then a:=num;
-----
??? - что такое 'a' в этом контексте?
#2 
  scorpi_ скептик08.11.06 09:59
NEW 08.11.06 09:59 
in Antwort taksos 08.11.06 07:54, Zuletzt geändert 08.11.06 10:45 (scorpi_)
Это что за бредовый алгоритм? Посмотри в инете решето Эратоcфена, и перепиши программу соответствующим образом.
#3 
  scorpi_ скептик08.11.06 10:01
NEW 08.11.06 10:01 
in Antwort Murr 08.11.06 09:57
Переменная из четвёртой строчки
#4 
Murr коренной житель08.11.06 10:07
Murr
NEW 08.11.06 10:07 
in Antwort scorpi_ 08.11.06 10:01
Эээ... все никак не привыкнешь, что Я, при подобных вопросах, обычно не даю прямого ответа, но вынуждаю подумать...
Странно, что компайлер на этой строке не матюгнулся...
#5 
barmaglot местный житель08.11.06 10:16
barmaglot
NEW 08.11.06 10:16 
in Antwort taksos 08.11.06 07:54, Zuletzt geändert 08.11.06 10:20 (barmaglot)
В ответ на:
Пожалуйста, тыкните в ошибку!

[Недоброжелатель сказал бы, что ошибка -- в днк прокладки:). Шутка:). Ничего личного:).] Пардон. Вычеркнуто после просмотра профиля спрашивавшего
А по сути -- алгоритм выбран странный, далёкий от оптимального. Перефразируя условие задачи -- нужно искать простые числа одно за другим, 1, 2, 3, 5, ... до тех пор, пока счётчик текущего найденного простого числа не сравняется с наперёд заданным значением, и это самое текущее простое число и будет ответом. "Шагать" по ряду натуральных чисел нужно с шагом 2 (чётные числа простыми не бывают), а для проверки на простоту -- делить только на уже найденные простые числа. Вот и вся любовь:).
#6 
  scorpi_ скептик08.11.06 10:18
NEW 08.11.06 10:18 
in Antwort barmaglot 08.11.06 10:16
Делить здесь абсолютно ничего не надо. Достаточно прибавлять и сравнивать.
#7 
  scorpi_ скептик08.11.06 10:22
NEW 08.11.06 10:22 
in Antwort Murr 08.11.06 10:07
Думать здесь надо в первую очередь над алгоритмом.
#8 
barmaglot местный житель08.11.06 10:28
barmaglot
NEW 08.11.06 10:28 
in Antwort scorpi_ 08.11.06 10:18
В ответ на:
Достаточно прибавлять и сравнивать.

Я, наверное, дикий. Тему прошу развернуть. В особенности "сравнивать".
#9 
Murr коренной житель08.11.06 10:29
Murr
NEW 08.11.06 10:29 
in Antwort scorpi_ 08.11.06 10:22
Думать здесь надо в первую очередь над алгоритмом.
------
Думать там надо ой как много над чем - над использоваными индентификаторами, над стилем написания, над алгоритмом, над реализацией алгоритма с точки зрения использования памяти, над... тем, об чем писал бармаглот...
Но все эти думания не дадут ответа на заданный вопрос.
#10 
Murr коренной житель08.11.06 10:30
Murr
NEW 08.11.06 10:30 
in Antwort barmaglot 08.11.06 10:28
Почитай Решето Эратосфена - там только сложение и сравнение...
#11 
barmaglot местный житель08.11.06 10:36
barmaglot
NEW 08.11.06 10:36 
in Antwort Murr 08.11.06 10:30
О да! Только задолбаться можно, следя за текущим счётчиком-вычёркивателем:).
Кроме того, мальчегу нужны "первые сто простых чисел", а не "простые числа из интервала от одного до ста", а значит априори сходу и не скажешь, до какого числа решетить. ИМХО "мой" алгоритм оптимален с точки зрения затрат на "думание". Ах да, делить конечно можно не на все найденные, а только на те, которые .LE. SQRT (iCurrCounter).
#12 
Murr коренной житель08.11.06 10:40
Murr
NEW 08.11.06 10:40 
in Antwort barmaglot 08.11.06 10:36
ИМХО "мой" алгоритм оптимален с точки зрения затрат на "думание".
------
Ты будешь несколько удивлен, но Решето "работает" примерно на три-четыре порядка быстрее... Так что затраты на думание - не существенны...
#13 
barmaglot местный житель08.11.06 10:47
barmaglot
NEW 08.11.06 10:47 
in Antwort Murr 08.11.06 10:40
Хи. Спорю на ящик пива, что для чисел, которые влезаютЪ в интеджер, 4 порядка -- гон. Это первоЕ.
Также прошу представить достоверный метод априорной оценки "длины" решета для наперёд заданного номера простого числа. Это второЕ.
Ну что, "Гиннес" против "Францисканера" :) ?
#14 
Murr коренной житель08.11.06 10:58
Murr
NEW 08.11.06 10:58 
in Antwort barmaglot 08.11.06 10:47
Спорю
------
Сравни сначала набор операций для обоих случаев.
оценки "длины" решета
------
Верхняя оценка устроит? Для чисел в пределах Int32 - не более 4294967296...
Эээ... зазипованный набор простых чисел в этом диаппазоне весит порядка 70Мб...
#15 
  scorpi_ скептик08.11.06 10:58
NEW 08.11.06 10:58 
in Antwort barmaglot 08.11.06 10:47, Zuletzt geändert 08.11.06 11:04 (scorpi_)
Также прошу представить достоверный метод априорной оценки "длины" решета для наперёд заданного номера простого числа.
http://de.wikipedia.org/wiki/Primzahlsatz/
#16 
Simple Nothing is f*cked08.11.06 10:59
Simple
NEW 08.11.06 10:59 
in Antwort barmaglot 08.11.06 10:47
В ответ на:
Ну что, "Гиннес" против "Францисканера" :) ?

Хиииитренький.
SCNR.
#17 
barmaglot местный житель08.11.06 11:00
barmaglot
NEW 08.11.06 11:00 
in Antwort Murr 08.11.06 10:58
В ответ на:
Верхняя оценка устроит? Для чисел в пределах Int32 - не более 4294967296...

Такая -- не устроит. Даю вторую попытку .
#18 
Murr коренной житель08.11.06 11:02
Murr
NEW 08.11.06 11:02 
in Antwort barmaglot 08.11.06 11:00
Даю вторую попытку
------
Эээ, нет... теперь твоя очередь оценивать... скажем, количество операций деления для нахождения простого числа с заданным номером...
#19 
barmaglot местный житель08.11.06 12:04
barmaglot
NEW 08.11.06 12:04 
in Antwort Murr 08.11.06 11:02
Ну что ж, из уважения к Вашим сединам сыграю один раз по чужим правилам:).
Я писал о числах, влезающих в интеджер, т. е. до порядка 32 тысяч. В интервале 1..32000 лежит около 2500 простых. Таким образом, если оценивать сверху размером типа, при реализации решета вычеркнуть прийдётся где-то 29500 чисел.
Если, как я предложил, каждое нечётное делить на все простые, которые меньше квадратного корня из него, то делений прийдётся сделать не более чем 16000 умножить на количество простых, которые не превосходят корень из 32000 (это количество -- 41, но на самом деле результат будет меньше, это оценка сверху).
IDL> print, 1.6e4*41.
656000.
Теперь положим, что одно деление жрёт ресурса как 20 сложений, и поделим на заявленные четыре порядка::
IDL> print, 20.*1.6e4*41./1e4
1312.00
Итого:: если бы песня про четыре порядка была правильной, то решето в этом случае потребовало бы меньше, чем 1312 вычёркиваний.
Вердикт:: для типа интеджер четыре порядка -- гон, что и требовалось доказать. Гиннес прошу высылать в Мюнхен . С моей стороны спор окончен.
#20 
Murr коренной житель08.11.06 12:21
Murr
NEW 08.11.06 12:21 
in Antwort barmaglot 08.11.06 12:04
Теперь положим, что
------
Там, при делении, примерно 350 тактов, против 16 при сложении. Но, там не только деление...
Ну а вычеркивание (рассылка) - одна машинная команда, плюс еще три на подготовку следующей...
Видишь ли, во времена первого пня я прогнал оба варианта, предварительно оптимизировав их до предела...
#21 
barmaglot местный житель08.11.06 12:27
barmaglot
NEW 08.11.06 12:27 
in Antwort Murr 08.11.06 12:21
В ответ на:
... одна машинная команда, плюс еще три на подготовку следующей

Фактор четыре не спасёт гиганта мысли. Где мой Гиннес ?
#22 
Murr коренной житель08.11.06 12:49
Murr
NEW 08.11.06 12:49 
in Antwort barmaglot 08.11.06 12:27
Фактор четыре...
------
Ну хорошо.
Требование - найти простые числа в пределах... ну скажем... хммм... 2^82
Пиши оба варианта. Оптимизируй. Выставляй на общее обозрение для убивания неэффективностей... результаты - сравним, оценим, так сказать, фактор...
P.S. 2^82 дано с учем того, что какой-то из компайлеров поддерживает тип 2^80
#23 
taksos завсегдатай08.11.06 12:50
taksos
NEW 08.11.06 12:50 
in Antwort Murr 08.11.06 09:57
Отправитель: Murr
Заголовок: Re: Поиск простого числа
if prov(num)=true then a:=num;
-----
??? - что такое 'a' в этом контексте?
опечатка, там должно стоять a
мне пока не до идеальности алгоритма, хочу добиться чтобы он работал
я программированием занимаюсь с гулькин нос, так что не ругайте сильно
#24 
barmaglot местный житель08.11.06 12:53
barmaglot
NEW 08.11.06 12:53 
in Antwort Murr 08.11.06 12:49
Батюшки, Любезный, до Вас не ещё дошло? Спор был о том, что в рамках типа интеджер четыре порядка -- гон. Убедительно доказано. Спор окончен, трёп прродолжается. А пиво -- где?
Прямой вопрос:: пиво -- будет? Или словоблудство одно?
#25 
Murr коренной житель08.11.06 12:56
Murr
NEW 08.11.06 12:56 
in Antwort taksos 08.11.06 12:50

if prov(num)=true then a:=num;
-----
??? - что такое 'a' в этом контексте?
опечатка, там должно стоять a
=====
???
Думай снова над этой же строкой. И посмотри на замечание Скорпи - по поводу 4-й строки.
#26 
Murr коренной житель08.11.06 12:59
Murr
NEW 08.11.06 12:59 
in Antwort barmaglot 08.11.06 12:53
Так оно, любезный, ОДИНАКОВО, что для 32-х, что для 82-х... в последнем случае только лишь эффективно подчеркивается разница во времени вычислений...
Прямой вопрос:: пиво -- будет?
------
А код?
#27 
taksos завсегдатай08.11.06 13:00
taksos
NEW 08.11.06 13:00 
in Antwort barmaglot 08.11.06 12:53
ясно, уменьшим порядок числа, можно тогда использовать extended
п.с. это была не лпечатка форум не отбражает квадратные скобки и что в них (а "квадратная скобка" b "квадратная скобка" )
#28 
barmaglot местный житель08.11.06 13:04
barmaglot
NEW 08.11.06 13:04 
in Antwort Murr 08.11.06 12:59
Код не требуется, была приведена оценка количества операций сверху.
Считаю Вас неблагородным человеком, не умеющим держать слово и от дальнейшей полемики отказываюсь. Ящик пива "Гиннес", который Вы мне проиграли, дарю любому из участников этого форума, который сможет его из Вас выдоить .
#29 
Simple Nothing is f*cked08.11.06 13:05
Simple
NEW 08.11.06 13:05 
in Antwort taksos 08.11.06 13:00
Пользуйся тэгом pre.
#30 
Murr коренной житель08.11.06 13:11
Murr
NEW 08.11.06 13:11 
in Antwort taksos 08.11.06 13:00
В том же месте тебе нужно использовать не b, а отдельную переменную, которая увеличивается после добавления числа в массив.
#31 
Murr коренной житель08.11.06 13:17
Murr
NEW 08.11.06 13:17 
in Antwort barmaglot 08.11.06 13:04
Код не требуется, была приведена оценка количества операций сверху.
-----
на первом пне, для 32-битных целых, де факто, время счета при помощи деления - более двух месяцев, а Решетом - менее пяти дней... разница таки более трех порядков... хотя и двоичных...
P.S. Если не сложно - где именно было принято пари? Бо, в заведомо выигрышных(проигрышных) ситуациях Я не спорю...
#32 
taksos завсегдатай08.11.06 15:44
taksos
NEW 08.11.06 15:44 
in Antwort Murr 08.11.06 13:11
Murr, спасибо доделал, теперь нормально работает, сейчас буду делать с решетом Эратосфена
#33 
  scorpi_ скептик08.11.06 17:24
NEW 08.11.06 17:24 
in Antwort Murr 08.11.06 13:17
Результат щас будет. Я тут решил миллионное простое число посчитать. Решетом - 0.92с (на Athlon64 3000+), а вот когда будет готова бармаглотова версия... минут 5 она уже считает. Пиво пополам?
#34 
Murr коренной житель08.11.06 17:26
Murr
NEW 08.11.06 17:26 
in Antwort scorpi_ 08.11.06 17:24, Zuletzt geändert 08.11.06 17:28 (Murr)
Забирай целиком - я пиво не пользую, а Хеннесси он не предлагал...
Когда? Думаю, что не ранее послезавтра...
#35 
  scorpi_ скептик08.11.06 17:33
NEW 08.11.06 17:33 
in Antwort Murr 08.11.06 17:26
Ну, если за час не кончит, поставлю число поменьше... Или на ночь запущу. Три порядка - 16 минут, четыре - 3 часа. Но судя по расходу памяти она только к 400 тысячам подбирается Это за 14 минут. А дальше ведь только хуже будет...
#36 
Murr коренной житель08.11.06 17:43
Murr
NEW 08.11.06 17:43 
in Antwort scorpi_ 08.11.06 17:33
Я ведь реальные цифирьки по времени приводил.
#37 
Murr коренной житель08.11.06 17:50
Murr
NEW 08.11.06 17:50 
in Antwort scorpi_ 08.11.06 17:24
Решетом - 0.92с
----
Эээ... если ты шел "снизу", т.е. от 0 к 10^6, то можно еще примерно 0.10-0.15 снять идя"сверху", от миллиона... чисто за счет отсутствия сравнения с верхним пределом... и нулем...
#38 
  scorpi_ скептик08.11.06 17:55
NEW 08.11.06 17:55 
in Antwort Murr 08.11.06 17:50, Zuletzt geändert 08.11.06 18:13 (scorpi_)
Каким образом? Мы генерируем простые числа начиная с трёх...
Вот собственно -
В ответ на:
#include <iostream>
#include <vector>
#include <ctime>
const unsigned long MAX_PRIME = 0x19999998;
unsigned long get_prime_by_sieve( unsigned long number )
{
if ( number == 1 )
return 2;
std::vector<bool> sieve = std::vector<bool>( number * 10, 0 );
unsigned long i = 0;
size_t sieve_size = 2 * sieve.size();
for( unsigned long counter = 1; counter < number; ++counter )
{
while( sieve[ ++i ] );
unsigned int np = 2 * i + 1;
for( unsigned long j = np * np; j < sieve_size; j += np )
if ( j & 1 )
sieve[ j >> 1 ] = true;
}
return 2 * i + 1;
}
unsigned long get_prime_by_barmaglot( unsigned long number )
{
if ( number == 1 )
return 2;
std::vector<unsigned long> found_primes( 1, 2 );
found_primes.reserve( number * 10 );
unsigned long i = 3;
for( unsigned long counter = 1; counter < number; i += 2 )
{
unsigned long j = 0;
for( ; j < found_primes.size(); ++j )
if ( i % found_primes[j] == 0 )
break;
if ( j == found_primes.size() )
{
found_primes.push_back( i );
++counter;
}
}
return i - 2;
}
int main()
{
std::cout << "Berechnung der millionsten Primzahl" << std::endl;
clock_t start = clock();
unsigned long prime = get_prime_by_sieve( 1000000 );
double elapsed = (double)( clock() - start ) / CLOCKS_PER_SEC;
std::cout << "Sieve, Primzahl: " << prime << " time elapsed: " << elapsed << std::endl;
start = clock();
prime = get_prime_by_barmaglot( 1000000 );
elapsed = (double)( clock() - start ) / CLOCKS_PER_SEC;
std::cout << "Barmaglot, Primzahl: " << prime << " time elapsed: " << elapsed << std::endl;
}

#39 
  scorpi_ скептик08.11.06 18:04
NEW 08.11.06 18:04 
in Antwort barmaglot 08.11.06 12:04, Zuletzt geändert 08.11.06 19:54 (scorpi_)
Итого:: если бы песня про четыре порядка была правильной, то решето в этом случае потребовало бы меньше, чем 1312 вычёркиваний.
Гыыыыы, извини, ты о теории сложности хоть какое-нибудь представление имеешь? Решето это O(n*log n), а твой алгоритм O(n^2), поэтому коэффициенты (которые весьма не в твою пользу, ибо в решете сплошь однотактовые операции, а у тебя дорогостоящее деление) считать излишне.
PS Сорри, не заметил ограничения по корням. Тогда O(n^1.5) - всё равно хуже чем O(n*log n).
PPS Я писал о числах, влезающих в интеджер, т. е. до порядка 32 тысяч. - всё ещё на 8086 сидишь? Бедненьий... У меня так интеджер 64-битовый
#40 
Murr коренной житель08.11.06 18:22
Murr
NEW 08.11.06 18:22 
in Antwort scorpi_ 08.11.06 17:55
Каким образом?
-----
За счет префиксного декремента, устанавливающего флаг равенства нулю. Правда для того, чтобы оценка (0.10) была верной там должен быть только развернутый код.
Вместо:
for( unsigned long counter = 1; counter < number; ++counter )
for( unsigned long counter = number; --counter; /* no */ )
И так же со вторым циклом.
#41 
Murr коренной житель08.11.06 18:27
Murr
NEW 08.11.06 18:27 
in Antwort scorpi_ 08.11.06 18:04
а твой алгоритм O(n^3)
------
он предлагал понизить до О(n*sqrt(n)), если я не ошибаюсь...
#42 
  scorpi_ скептик08.11.06 18:51
NEW 08.11.06 18:51 
in Antwort Murr 08.11.06 18:22
Да, можно. Только измеримых результатов не даёт, так как это одна однотактовая команда.
#43 
  scorpi_ скептик08.11.06 19:40
NEW 08.11.06 19:40 
in Antwort Murr 08.11.06 18:27, Zuletzt geändert 13.11.06 12:54 (scorpi_)
А, я его постинга сразу не заметил... Тогда да, О(n^1,5), первоначальный вариант кстати был конечно же О(n^2), а не в куб.
Короче вот -
В ответ на:
unsigned long get_prime_by_barmaglot( unsigned long number )
{
if ( number == 1 )
return 2;
std::vector<unsigned long> found_primes( 1, 2 );
found_primes.reserve( number );
unsigned long i = 3;
for( unsigned long counter = number - 1; counter; i += 2 )
{
for( unsigned long j = 0; j < found_primes.size(); ++j )
if ( i < found_primes[j] * found_primes[j] || i % found_primes[j] == 0 )
break;
found_primes.push_back( i );
--counter;
}
return i - 2;
}


Даёт ответ за 15сек, против 0.85 для решета (number =1000000)
#44 
Murr коренной житель08.11.06 20:06
Murr
NEW 08.11.06 20:06 
in Antwort scorpi_ 08.11.06 18:51, Zuletzt geändert 08.11.06 20:09 (Murr)
так как это одна однотактовая команда.
-----
На базовом i8086, DEC и ИБМ360 в одну команду укладываться не удавалось...
#45 
  scorpi_ скептик08.11.06 20:12
08.11.06 20:12 
in Antwort Murr 08.11.06 20:06
Такого антиквариата я не застал.
#46 
barmaglot местный житель09.11.06 09:54
barmaglot
NEW 09.11.06 09:54 
in Antwort scorpi_ 08.11.06 18:04
Лень мне вникать во всё написанное. Ты опроверг утверждение про четыре порядка для простых до 32 тысяч? Если вдруг да (весьма сомнительно ) -- вечером как будет время напишу и проверю сам.
Что касается
В ответ на:
всё ещё на 8086 сидишь? Бедненьий...

то "сижу" я на таком (http://www.rzg.mpg.de/computing/IBM_P5/hardware.html), от чего знающему человеку впору захлебнуться слюнями .
#47 
Murr коренной житель09.11.06 10:38
Murr
NEW 09.11.06 10:38 
in Antwort barmaglot 09.11.06 09:54
вечером как будет время напишу и проверю сам.
------
Тебе предлагалось это сделать еще вчера...
#48 
1 2 3 alle