реклама на сайте
подробности

 
 
> Простые числа, где можно найти ряд
ig_z
сообщение Oct 2 2006, 12:21
Сообщение #1


Местный
***

Группа: Свой
Сообщений: 437
Регистрация: 27-08-04
Пользователь №: 551



Подскажите, где можно найти ряд простых чисел. От 1000 и дальше. Поиском найти не смог.
Go to the top of the page
 
+Quote Post
 
Start new topic
Ответов
DRUID3
сообщение Oct 2 2006, 13:21
Сообщение #2


山伏
*****

Группа: Свой
Сообщений: 1 827
Регистрация: 3-08-06
Из: Kyyiv
Пользователь №: 19 294



Цитата(ig_z @ Oct 2 2006, 15:21) *
Подскажите, где можно найти ряд простых чисел. От 1000 и дальше. Поиском найти не смог.

biggrin.gif Вы бы еще формулу попросили...

Код
//////////////////////////////////////////////////////////////////////////////
//
//  Finding prime numbers
//  (c) Johna Smith, 1996
//
//  Method description:
//   We take a number and try to divide it. If we can divide it
//   without remainder - this is not prime number.
//   We can take into account only odd numbers, because we can
//   divide all even number by 2. Also we can store all prime
//   numbers that are already found in an array and try to divide
//   all new numbers only by numbers from this array.
//   If we want to find all prime numbers less than N the size of
//   the array should be sqrt(N)/2
//
//////////////////////////////////////////////////////////////////////////////

#include <stdio.h>

#define N   160  // so we can find all prime numbers that are less than 100000
#define M   25   // check all numbers less than 250

int Simple[N];
int k=0;
enum {yes,no} simple;

void main(void)
{
// it's easy: 2 and 3 are prime
if (M>=2) printf("2\n");                    // 2 is simple 'cause we can divide it only by itself and 1
Simple[k++]=2;
if (M>=3) printf("3\n");
Simple[k++]=3;

// but what we can say about other numbers:
for(int i=5; i<=M; i+=2)
{
   simple=yes;
   for(int j=0; j<k; j++)
   {
     if (Simple[j]*Simple[j]>i) break;     // other Simple[j] is too big for i
     if ((i%Simple[j])==0) simple=no;    // there's no remainder - not prime
   }
   if (simple==yes)
   {
     printf("%d\n",i);
     Simple[k++]=i;
   }
}
}

Это я копирнул с "винграда", раздел "алгоритмы".

P.S.: Кстати, о формуле простых чисел ( biggrin.gif laugh.gif ), вот человек получил кое какие аппроксимации ряда (я сам не проверял, т.к. довольно холоден к фундаментальной математике). Но если надумаете, как их применить в прикладных задачах - Вам и карты в руки... wink.gif


--------------------
Нас помнят пока мы мешаем другим...
//--------------------------------------------------------
Хороший блатной - мертвый...
//--------------------------------------------------------
Нет старик, это те дроиды которых я ищу...
Go to the top of the page
 
+Quote Post
ig_z
сообщение Oct 2 2006, 13:53
Сообщение #3


Местный
***

Группа: Свой
Сообщений: 437
Регистрация: 27-08-04
Пользователь №: 551



Цитата(DRUID3 @ Oct 2 2006, 16:21) *
Цитата(ig_z @ Oct 2 2006, 15:21) *

Подскажите, где можно найти ряд простых чисел. От 1000 и дальше. Поиском найти не смог.

biggrin.gif Вы бы еще формулу попросили...



Спасибо за ответ. Я еще более холоден к данной теме. Просто просьба со стороны.

Что касается формульного расчета, то действительно гарантировано простое число на выходе. Но, если я не ошибаюсь, нет гарантии, что будет рассчитан весь ряд без пропусков. Вроде я слышал, что раньше публиковали полные ряды, получаемые тупым перебором. Вот такое бы найти.
Go to the top of the page
 
+Quote Post

Сообщений в этой теме
- ig_z   Простые числа   Oct 2 2006, 12:21
- - Andrew10   В Mathematica-5.2 функция Prime[n] дает n-ое прост...   Oct 2 2006, 14:35
|- - DRUID3   Цитата(Andrew10 @ Oct 2 2006, 17:35) В Ma...   Oct 2 2006, 15:06
- - Stanislav   В матлабе: P=primes(n); Здесь n - верхний предел п...   Oct 2 2006, 15:11
- - Andrew10   Цитата дык кто Вам мешает копирнуть ряд сюда и пом...   Oct 2 2006, 18:17
- - Krys   Извиняюсь, если офтопик, но ещё в школе на олимпиа...   Oct 27 2006, 08:08
|- - exSSerge   Цитата(Krys @ Oct 27 2006, 15:08) Извиняю...   Oct 27 2006, 18:44
- - net   да ладно вам в известной книге приведена очень ...   Oct 27 2006, 19:11
|- - BEZU   Цитата(net @ Oct 27 2006, 23:11) 2*n+1 Эт...   Oct 31 2006, 09:43
|- - net   Цитата(BEZU @ Oct 31 2006, 12:43) Цитата(...   Oct 31 2006, 16:43
- - Krys   exSSerge, не очень Вас понял. Нельзя ли прямо на п...   Oct 28 2006, 07:19
|- - pdk   всегото одна строчка в Maple: writedata("c:...   Oct 31 2006, 07:55
|- - SKov   Цитата(Krys @ Oct 28 2006, 10:19) exSSerg...   Nov 1 2006, 14:22
|- - Oldring   Цитата(SKov @ Nov 1 2006, 17:22) Цитата(K...   Nov 1 2006, 14:51
||- - SKov   Цитата(Oldring @ Nov 1 2006, 17:51) Цитат...   Nov 2 2006, 07:49
||- - exSSerge   Цитата(SKov @ Nov 2 2006, 13:49) В чуть б...   Nov 2 2006, 09:24
||- - Oldring   Цитата(exSSerge @ Nov 2 2006, 12:24) Цита...   Nov 2 2006, 09:49
|- - -=ВН=-   Цитата(SKov @ Nov 1 2006, 17:22) Цитата(K...   Nov 3 2006, 09:05
|- - -=ВН=-   Цитата(-=ВН=- @ Nov 3 2006, 12...   Nov 6 2006, 23:16
- - BEZU   Не читали... Каюсь! судя по всему, сглупил... ...   Oct 31 2006, 23:03
|- - UMP   В книге Дьяконов В.П. Энциклопедия Mathcad 200i, 1...   Nov 2 2006, 20:09
- - -=ВН=-   Цитата(ig_z @ Oct 2 2006, 15:21) Подскажи...   Nov 3 2006, 08:17


Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 22nd July 2025 - 00:49
Рейтинг@Mail.ru


Страница сгенерированна за 0.01389 секунд с 7
ELECTRONIX ©2004-2016