Дискуссионный математический форумМатематический форум
Математический форум Math Help Planet

Обсуждение и решение задач по математике, физике, химии, экономике

Теоретический раздел
Часовой пояс: UTC + 3 часа [ Летнее время ]
новый онлайн-сервис
число, сумма и дата прописью

Часовой пояс: UTC + 3 часа [ Летнее время ]




Начать новую тему Ответить на тему  [ Сообщений: 22 ]  На страницу 1, 2, 3  След.
Автор Сообщение
 Заголовок сообщения: Массив простых чисел
СообщениеДобавлено: 30 май 2019, 19:36 
Не в сети
Начинающий
Зарегистрирован:
27 апр 2017, 21:12
Сообщений: 38
Cпасибо сказано: 19
Спасибо получено:
1 раз в 1 сообщении
Очков репутации: 2

Добавить очки репутацииУменьшить очки репутации
Как сформировать массив простых чисел? Какие есть алгоритмы? Длина массива задана. Есть мысль сравнить длину прохода от двух до случайно полученного числа. У простого она максимальна.

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 30 май 2019, 20:36 
Не в сети
Последняя инстанция
Зарегистрирован:
06 дек 2014, 09:11
Сообщений: 7070
Cпасибо сказано: 115
Спасибо получено:
1662 раз в 1508 сообщениях
Очков репутации: 283

Добавить очки репутацииУменьшить очки репутации
Какого прохода?

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 31 май 2019, 08:25 
Не в сети
Light & Truth
Аватара пользователя
Зарегистрирован:
28 апр 2016, 13:40
Сообщений: 3550
Cпасибо сказано: 5
Спасибо получено:
624 раз в 591 сообщениях
Очков репутации: 98

Добавить очки репутацииУменьшить очки репутации
Провести последовательные целочисленные деления от 2 до n (проверяемое число) и проверять остаток.

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 31 май 2019, 15:40 
Не в сети
Гений
Зарегистрирован:
02 июн 2018, 08:50
Сообщений: 659
Cпасибо сказано: 21
Спасибо получено:
105 раз в 103 сообщениях
Очков репутации: 16

Добавить очки репутацииУменьшить очки репутации
Если относительно просто то можно так

program massiv
implicit none

!Объявим переменные,
!в том числе целочисленный массив из 50-ти элементов

integer :: a(50), x, i
logical :: l, f

i=1; x=2

!Гоним числа вверх по единичке, пока не наберем 50 простых чисел

do while (i <= 50)
l = f(x)
if (l) then
a(i) = x
i=i+1
end if
x=x+1
end do

!Напечатаем полученные числа

print 100, a

100 format (10i5)
end program massiv !Конец программы

!Логическая функция определения числа на простоту

logical function f(x)
implicit none
integer, intent (in) :: x
integer :: temp, i

f = .true.
temp = x**(0.5)

do i = 2, temp
if (mod (x, i)==0) then
f =.false.
exit
end if
end do

end function f

результат
Код:
    2    3    5    7   11   13   17   19   23   29
   31   37   41   43   47   53   59   61   67   71
   73   79   83   89   97  101  103  107  109  113
  127  131  137  139  149  151  157  163  167  173
  179  181  191  193  197  199  211  223  227  229


Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 31 май 2019, 16:35 
Не в сети
Light & Truth
Зарегистрирован:
12 окт 2017, 13:50
Сообщений: 2358
Cпасибо сказано: 94
Спасибо получено:
709 раз в 684 сообщениях
Очков репутации: 200

Добавить очки репутацииУменьшить очки репутации
pacha писал(а):
Какие есть алгоритмы? Длина массива задана.

Решето Ератостена Вам известно? Если да, то :
1)задайте какого то большое [math]N >[/math] длину массива;
2) Через решето Ератостена найдите все простые чисель [math]\leqslant N[/math] и на кажды шаг следите
достигли ли длину массива;
3.1) если достигли, то идите к концу алгоритма;
3.2) если не достигли длину массива, но достигли [math]N-1[/math], то увеличите горная граница цикла решето
Ератостена на [math]= 2N[/math] и продолжайте от 2) дальше - [math]N+1,N+2,...,2N[/math]

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 31 май 2019, 18:12 
Не в сети
Последняя инстанция
Зарегистрирован:
14 июн 2011, 08:15
Сообщений: 3565
Cпасибо сказано: 50
Спасибо получено:
502 раз в 465 сообщениях
Очков репутации: 23

Добавить очки репутацииУменьшить очки репутации
Чтобы создать массив простых чисел до N, надо иметь
массив простых чисел до [math]\sqrt N[/math]

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 31 май 2019, 18:50 
Не в сети
Beautiful Mind
Аватара пользователя
Зарегистрирован:
09 авг 2018, 23:20
Сообщений: 1011
Cпасибо сказано: 32
Спасибо получено:
121 раз в 116 сообщениях
Очков репутации: 8

Добавить очки репутацииУменьшить очки репутации
pacha писал(а):
сравнить длину прохода от двух до случайно полученного числа. У простого она максимальна.

Простые числа статистикой не поборешь. В любом вероятностном алгоритме будут пропуски.
Обычный алгоритм, видимо, и самый быстрый:
1. Есть два простых числа: 2 и 3 (первоначальный Х) .
2. Прибавляем к Х двойку.
3. Вычисляем квадратный корень из Х. Это ограничитель на последнее простое число, делимость на которое мы проверяем.
4. Делим Х на простые, начиная с 3 до корня из шага 3.
5. Если нулевой остаток получился, то немедленное прерывание: число составное. Переход к шагу 2.
6. Если нулевого остатка не случилось, то добавляем Х в наш массив. Переход к шагу 2.

Вернуться к началу
 Профиль  
Cпасибо сказано 
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 01 июн 2019, 03:03 
Не в сети
Гений
Зарегистрирован:
02 июн 2018, 08:50
Сообщений: 659
Cпасибо сказано: 21
Спасибо получено:
105 раз в 103 сообщениях
Очков репутации: 16

Добавить очки репутацииУменьшить очки репутации
Если логическую функцию немного усложнить (как описал atlakatl), то получим приличное ускорение работы. В этом варианте определяем число на четность. Если оно четное, то сразу же говорим что оно не простое и остаток от деления перебором уже не ищем, для нечетных чисел цикл выполняется с шагом 2 начиная с попытки деления на 3.

logical function f(x)
implicit none
integer, intent (in) :: x
integer :: temp, i

f = .true.

if (x /= 2) then
if (mod(x, 2) == 0) then
f = .false.
return
end if
end if

temp = x**(0.5)

do i = 3, temp, 2
if (mod(x, i) == 0) then
f = .false.
exit
end if
end do

end function f


массив первых 1000 простых чисел

▼ массив
Код:
     2     3     5     7    11    13    17    19    23    29    31    37    41    43    47
    53    59    61    67    71    73    79    83    89    97   101   103   107   109   113
   127   131   137   139   149   151   157   163   167   173   179   181   191   193   197
   199   211   223   227   229   233   239   241   251   257   263   269   271   277   281
   283   293   307   311   313   317   331   337   347   349   353   359   367   373   379
   383   389   397   401   409   419   421   431   433   439   443   449   457   461   463
   467   479   487   491   499   503   509   521   523   541   547   557   563   569   571
   577   587   593   599   601   607   613   617   619   631   641   643   647   653   659
   661   673   677   683   691   701   709   719   727   733   739   743   751   757   761
   769   773   787   797   809   811   821   823   827   829   839   853   857   859   863
   877   881   883   887   907   911   919   929   937   941   947   953   967   971   977
   983   991   997  1009  1013  1019  1021  1031  1033  1039  1049  1051  1061  1063  1069
  1087  1091  1093  1097  1103  1109  1117  1123  1129  1151  1153  1163  1171  1181  1187
  1193  1201  1213  1217  1223  1229  1231  1237  1249  1259  1277  1279  1283  1289  1291
  1297  1301  1303  1307  1319  1321  1327  1361  1367  1373  1381  1399  1409  1423  1427
  1429  1433  1439  1447  1451  1453  1459  1471  1481  1483  1487  1489  1493  1499  1511
  1523  1531  1543  1549  1553  1559  1567  1571  1579  1583  1597  1601  1607  1609  1613
  1619  1621  1627  1637  1657  1663  1667  1669  1693  1697  1699  1709  1721  1723  1733
  1741  1747  1753  1759  1777  1783  1787  1789  1801  1811  1823  1831  1847  1861  1867
  1871  1873  1877  1879  1889  1901  1907  1913  1931  1933  1949  1951  1973  1979  1987
  1993  1997  1999  2003  2011  2017  2027  2029  2039  2053  2063  2069  2081  2083  2087
  2089  2099  2111  2113  2129  2131  2137  2141  2143  2153  2161  2179  2203  2207  2213
  2221  2237  2239  2243  2251  2267  2269  2273  2281  2287  2293  2297  2309  2311  2333
  2339  2341  2347  2351  2357  2371  2377  2381  2383  2389  2393  2399  2411  2417  2423
  2437  2441  2447  2459  2467  2473  2477  2503  2521  2531  2539  2543  2549  2551  2557
  2579  2591  2593  2609  2617  2621  2633  2647  2657  2659  2663  2671  2677  2683  2687
  2689  2693  2699  2707  2711  2713  2719  2729  2731  2741  2749  2753  2767  2777  2789
  2791  2797  2801  2803  2819  2833  2837  2843  2851  2857  2861  2879  2887  2897  2903
  2909  2917  2927  2939  2953  2957  2963  2969  2971  2999  3001  3011  3019  3023  3037
  3041  3049  3061  3067  3079  3083  3089  3109  3119  3121  3137  3163  3167  3169  3181
  3187  3191  3203  3209  3217  3221  3229  3251  3253  3257  3259  3271  3299  3301  3307
  3313  3319  3323  3329  3331  3343  3347  3359  3361  3371  3373  3389  3391  3407  3413
  3433  3449  3457  3461  3463  3467  3469  3491  3499  3511  3517  3527  3529  3533  3539
  3541  3547  3557  3559  3571  3581  3583  3593  3607  3613  3617  3623  3631  3637  3643
  3659  3671  3673  3677  3691  3697  3701  3709  3719  3727  3733  3739  3761  3767  3769
  3779  3793  3797  3803  3821  3823  3833  3847  3851  3853  3863  3877  3881  3889  3907
  3911  3917  3919  3923  3929  3931  3943  3947  3967  3989  4001  4003  4007  4013  4019
  4021  4027  4049  4051  4057  4073  4079  4091  4093  4099  4111  4127  4129  4133  4139
  4153  4157  4159  4177  4201  4211  4217  4219  4229  4231  4241  4243  4253  4259  4261
  4271  4273  4283  4289  4297  4327  4337  4339  4349  4357  4363  4373  4391  4397  4409
  4421  4423  4441  4447  4451  4457  4463  4481  4483  4493  4507  4513  4517  4519  4523
  4547  4549  4561  4567  4583  4591  4597  4603  4621  4637  4639  4643  4649  4651  4657
  4663  4673  4679  4691  4703  4721  4723  4729  4733  4751  4759  4783  4787  4789  4793
  4799  4801  4813  4817  4831  4861  4871  4877  4889  4903  4909  4919  4931  4933  4937
  4943  4951  4957  4967  4969  4973  4987  4993  4999  5003  5009  5011  5021  5023  5039
  5051  5059  5077  5081  5087  5099  5101  5107  5113  5119  5147  5153  5167  5171  5179
  5189  5197  5209  5227  5231  5233  5237  5261  5273  5279  5281  5297  5303  5309  5323
  5333  5347  5351  5381  5387  5393  5399  5407  5413  5417  5419  5431  5437  5441  5443
  5449  5471  5477  5479  5483  5501  5503  5507  5519  5521  5527  5531  5557  5563  5569
  5573  5581  5591  5623  5639  5641  5647  5651  5653  5657  5659  5669  5683  5689  5693
  5701  5711  5717  5737  5741  5743  5749  5779  5783  5791  5801  5807  5813  5821  5827
  5839  5843  5849  5851  5857  5861  5867  5869  5879  5881  5897  5903  5923  5927  5939
  5953  5981  5987  6007  6011  6029  6037  6043  6047  6053  6067  6073  6079  6089  6091
  6101  6113  6121  6131  6133  6143  6151  6163  6173  6197  6199  6203  6211  6217  6221
  6229  6247  6257  6263  6269  6271  6277  6287  6299  6301  6311  6317  6323  6329  6337
  6343  6353  6359  6361  6367  6373  6379  6389  6397  6421  6427  6449  6451  6469  6473
  6481  6491  6521  6529  6547  6551  6553  6563  6569  6571  6577  6581  6599  6607  6619
  6637  6653  6659  6661  6673  6679  6689  6691  6701  6703  6709  6719  6733  6737  6761
  6763  6779  6781  6791  6793  6803  6823  6827  6829  6833  6841  6857  6863  6869  6871
  6883  6899  6907  6911  6917  6947  6949  6959  6961  6967  6971  6977  6983  6991  6997
  7001  7013  7019  7027  7039  7043  7057  7069  7079  7103  7109  7121  7127  7129  7151
  7159  7177  7187  7193  7207  7211  7213  7219  7229  7237  7243  7247  7253  7283  7297
  7307  7309  7321  7331  7333  7349  7351  7369  7393  7411  7417  7433  7451  7457  7459
  7477  7481  7487  7489  7499  7507  7517  7523  7529  7537  7541  7547  7549  7559  7561
  7573  7577  7583  7589  7591  7603  7607  7621  7639  7643  7649  7669  7673  7681  7687
  7691  7699  7703  7717  7723  7727  7741  7753  7757  7759  7789  7793  7817  7823  7829
  7841  7853  7867  7873  7877  7879  7883  7901  7907  7919



atlakatl писал(а):
4. Делим Х на простые, начиная с 3 до корня из шага 3.

Хотел только уточнить, именно на простые предлагаете делить? У меня просто с шагом 2.

Вернуться к началу
 Профиль  
Cпасибо сказано 
За это сообщение пользователю Emphatic18 "Спасибо" сказали:
atlakatl
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 01 июн 2019, 03:55 
Не в сети
Beautiful Mind
Аватара пользователя
Зарегистрирован:
09 авг 2018, 23:20
Сообщений: 1011
Cпасибо сказано: 32
Спасибо получено:
121 раз в 116 сообщениях
Очков репутации: 8

Добавить очки репутацииУменьшить очки репутации
Emphatic18 писал(а):
именно на простые предлагаете делить?

Цэ ж основная теорема арифметики. Если число делится на 26, то оно делится и на 2, и на 13.
Обратно, если число не делится ни на 2, ни на 13, то оно не делится на 26.
Ну и проверять чётные числа незачем. Просто увеличиваем очередного кандидата на 2.

Вернуться к началу
 Профиль  
Cпасибо сказано 
За это сообщение пользователю atlakatl "Спасибо" сказали:
Emphatic18
 Заголовок сообщения: Re: Массив простых чисел
СообщениеДобавлено: 01 июн 2019, 04:40 
Не в сети
Гений
Зарегистрирован:
02 июн 2018, 08:50
Сообщений: 659
Cпасибо сказано: 21
Спасибо получено:
105 раз в 103 сообщениях
Очков репутации: 16

Добавить очки репутацииУменьшить очки репутации
Тогда вся программа получается такой:

▼ код
program massiv
implicit none

integer :: a(50), x, i
logical :: l, f

i=1; x=3; a(1) = 2; a(2) = 3

!Гоним числа вверх с шагом 2, пока не наберем 50 простых чисел

do while (i<=50)
l = f(x, a)
if (l) then
a(i+1) = x
i=i+1
end if
x=x+2
end do

!Напечатаем полученные числа

print 100, a
100 format (10i5)
end program massiv !Конец программы

!Логическая функция определения числа на простоту

logical function f(x, a)
implicit none
integer, intent (in) :: x, a(50)
integer :: temp, i

f = .true.
temp = x**(0.5)

do i = 1, temp
if (mod (x, a(i)) == 0) then
f = .false.
exit
end if
end do

end function f


результат работы тот же
но в этом случае функция применима только в этой основной программе

Вернуться к началу
 Профиль  
Cпасибо сказано 
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему    На страницу 1, 2, 3  След.  Страница 1 из 3 [ Сообщений: 22 ]

 Похожие темы   Автор   Ответы   Просмотры   Последнее сообщение 
Множество простых чисел и пар простых чисел-близнецов бескон

в форуме Размышления по поводу и без

korolchukvasily

2

257

28 июн 2023, 11:23

Дан массив целых чисел. Найти максимальное число. Сделать в

в форуме Информатика и Компьютерные науки

Nesquik60

0

641

29 июн 2014, 14:02

Сумма простых чисел

в форуме Интересные задачи участников форума MHP

MurChik

15

222

28 янв 2024, 10:52

Группы простых чисел

в форуме Теория чисел

vorvalm

5

1063

03 дек 2014, 15:00

Формула простых чисел?

в форуме Теория чисел

Ferma

18

1087

05 дек 2018, 21:11

Формула для простых чисел

в форуме Размышления по поводу и без

Tirpa

30

1343

22 авг 2019, 23:30

Тройки простых чисел

в форуме Теория чисел

Claudia

5

591

18 июн 2018, 13:13

Свойства простых чисел

в форуме Палата №6

Galina Alexandrovna

13

1569

21 июл 2016, 07:14

Формула для простых чисел

в форуме Размышления по поводу и без

ammo77

1

241

31 янв 2020, 12:22

Последовательность простых чисел

в форуме Теория чисел

DeD

2

654

28 мар 2017, 01:43


Часовой пояс: UTC + 3 часа [ Летнее время ]



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 9


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Перейти:  

Яндекс.Метрика

Copyright © 2010-2023 MathHelpPlanet.com. All rights reserved