Форум программистов, компьютерный форум, киберфорум
С++ для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.80/5: Рейтинг темы: голосов - 5, средняя оценка - 4.80
 Аватар для SerVal
37 / 36 / 9
Регистрация: 16.04.2015
Сообщений: 283

Решето Эратосфена (сегмент): медленно работает - как можно ускорить?

19.01.2017, 16:59. Показов 1048. Ответов 4
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Подсчёт числа простых чисел в диапазоне от "from" до "to"
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
typedef UINT64 Number;
 
__device__ Number eratosthenesSingleBlock(const Number from, const Number to)
{
    const Number memorySize = (to - from + 1) / 2;
 
    // initialize
    char* isPrime = new char[memorySize];
    memset(isPrime, 1, memorySize);
 
    for (Number i = 3; i*i <= to; i += 2)
    {
        // skip multiples of three: 9, 15, 21, 27, ...
        if (i >= 3 * 3 && i % 3 == 0)
            continue;
        // skip multiples of five
        if (i >= 5 * 5 && i % 5 == 0)
            continue;
        // skip multiples of seven
        if (i >= 7 * 7 && i % 7 == 0)
            continue;
        // skip multiples of eleven
        if (i >= 11 * 11 && i % 11 == 0)
            continue;
        // skip multiples of thirteen
        if (i >= 13 * 13 && i % 13 == 0)
            continue;
 
        // skip numbers before current slice
        Number minJ = ((from + i - 1) / i)*i;
        if (minJ < i*i)
            minJ = i*i;
 
        // start value must be odd
        if ((minJ & 1) == 0)
            minJ += i;
 
        // find all odd non-primes
        for (Number j = minJ; j <= to; j += 2 * i)
        {
            Number index = j - from;
            isPrime[index / 2] = 0;
        }
    }
 
    // count primes in this block
    UINT32 found = 0;
    for (UINT32 i = 0; i < memorySize; i++)
        found += isPrime[i];
    // 2 is not odd => include on demand
    if (from <= 2)
        found++;
 
    delete[] isPrime;
    return found;
}
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
19.01.2017, 16:59
Ответы с готовыми решениями:

Решето Эратосфена. Как ускорить?
Этот код не проходит задачу. доля секунды. как ускорить. или каким методом проидет #include &lt;bits/stdc++.h&gt; using namespace...

Прокомментируйте оставшиеся строчки, не понимаю как работает "Решето Эратосфена"
#include &lt;iostream&gt; //подключение стандартной библиотеки ввода-вывода using namespace std; //пространство имён. Чтобы писать &quot;cout...

Как разложить число на простые множители, используя решето Эратосфена?
Я только код для решета Эратосфена знаю(

4
Форумчанин
Эксперт CЭксперт С++
 Аватар для MrGluck
8216 / 5047 / 1437
Регистрация: 29.11.2010
Сообщений: 13,453
19.01.2017, 17:02
http://e-maxx.ru/algo/eratosthenes_sieve
Решето Эратосфена
Самый быстрый алгоритм Евклида вычисления НОД
0
 Аватар для SerVal
37 / 36 / 9
Регистрация: 16.04.2015
Сообщений: 283
19.01.2017, 17:14  [ТС]
Спасибо, это я уже смотрел. Нахождения простых чисел в диапазоне от Num1 до Num2 там нетути.
Все ищут от "печки" - с нуля или двойки.
0
Dimension
595 / 463 / 223
Регистрация: 08.04.2014
Сообщений: 1,710
19.01.2017, 19:16
SerVal, потому что решето по другому не будет работать ,вот вообще линейное решето http://e-maxx.ru/algo/prime_sieve_linear
думаю потом по массиву от нум1 до нум2 сами сможете пробежать
0
 Аватар для SerVal
37 / 36 / 9
Регистрация: 16.04.2015
Сообщений: 283
19.01.2017, 19:35  [ТС]
Цитата Сообщение от Dimension Посмотреть сообщение
потому что решето по другому не будет работать
Решето Эратосфена, приведённое в первом посте вполне себе работает.
*хотелось бы чтобы работало маленько побыстрее.
C++
1
2
3
4
5
6
7
8
9
10
11
D:\aProjects\CudaTest\x64\Release>PrimeSearchCuda.exe
 
Prime search:   4000000000000 .. 4001000000000
 
Accelerator: GeForce GTX 460
Found : 34462887
First prime :   4000000000039
Last prime  :   4000999999987
=============
Exec time : 110 seconds.
Primes/sec : 313663
Большой диапазон разбивается на много маленьких блоков (по 8 килобайт).
И для каждого блока решетом подсчитывается число простых. Потом суммируется для всех блоков.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
19.01.2017, 19:35
Помогаю со студенческими работами здесь

Решето Эратосфена
Кому надо - программа &quot;Решето Эратосфена&quot; на C++. Записывает в файл 1 000 000 первых простых чисел за 1/10 секунды (без вывода)!!! ...

Решето Эратосфена
Как можно реализовать? Подскажите плиз

Решето Эратосфена
Простое число — это любое целое число, которое точно делится без остатка только само на себя и на 1. Решето Эратосфена — это способ...

Решето Эратосфена
В общем задание посчитать количество простых чисел до заданного числа N. Написал такой алгоритм, работает только до 11 :cry: Уже час не...

Решето Эратосфена
Здравствуйте. Реализовал алгоритм &quot;Решето Эратосфена&quot; в виде класса. Взгляните, пожалуйста, и скажите, где я не прав. Спасибо. ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
5
Ответ Создать тему
Новые блоги и статьи
http://iceja.net/ математические сервисы
iceja 20.01.2026
Обновила свой сайт http:/ / iceja. net/ , приделала Fast Fourier Transform экстраполяцию сигналов. Однако предсказывает далеко не каждый сигнал (см ограничения http:/ / iceja. net/ fourier/ docs ). Также. . .
http://iceja.net/ сервер решения полиномов
iceja 18.01.2026
Выкатила http:/ / iceja. net/ сервер решения полиномов (находит действительные корни полиномов методом Штурма). На сайте документация по API, но скажу прямо VPS слабенький и 200 000 полиномов. . .
Расчёт переходных процессов в цепи постоянного тока
igorrr37 16.01.2026
/ * Дана цепь постоянного тока с R, L, C, k(ключ), U, E, J. Программа составляет систему уравнений по 1 и 2 законам Кирхгофа, решает её и находит переходные токи и напряжения на элементах схемы. . . .
Восстановить юзерскрипты Greasemonkey из бэкапа браузера
damix 15.01.2026
Если восстановить из бэкапа профиль Firefox после переустановки винды, то список юзерскриптов в Greasemonkey будет пустым. Но восстановить их можно так. Для этого понадобится консольная утилита. . .
Сукцессия микоризы: основная теория в виде двух уравнений.
anaschu 11.01.2026
https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/
WordPad для Windows 11
Jel 10.01.2026
WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .
Classic Notepad for Windows 11
Jel 10.01.2026
Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .
Почему дизайн решает?
Neotwalker 09.01.2026
В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru