2 / 2 / 0
Регистрация: 10.09.2016
Сообщений: 235
1

Код Алгоритма Полларда как устранить большие константа?

10.01.2019, 02:00. Показов 1950. Ответов 1
Метки нет (Все метки)

Приветствую Всех! Очень интересуюсь фа́кторингом в сети нашел статью "Алгоритм Ро Полларда для простой факторизации"

При целом положительном число n мы находим его делитель.


C++ (Qt)
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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
/* C++ program to find a prime factor of composite using 
   Pollard's Rho algorithm */
#include<bits/stdc++.h> 
using namespace std; 
  
/* Function to calculate (base^exponent)%modulus */
long long int modular_pow(long long int base, int exponent, 
                          long long int modulus) 
{ 
    /* initialize result */
    long long int result = 1; 
  
    while (exponent > 0) 
    { 
        /* if y is odd, multiply base with result */
        if (exponent & 1) 
            result = (result * base) % modulus; 
  
        /* exponent = exponent/2 */
        exponent = exponent >> 1; 
  
        /* base = base * base */
        base = (base * base) % modulus; 
    } 
    return result; 
} 
  
/* method to return prime divisor for n */
long long int PollardRho(long long int n) 
{ 
    /* initialize random seed */
    srand (time(NULL)); 
  
    /* no prime divisor for 1 */
    if (n==1) return n; 
  
    /* even number means one of the divisors is 2 */
    if (n % 2 == 0) return 2; 
  
    /* we will pick from the range [2, N) */
    long long int x = (rand()%(n-2))+2; 
    long long int y = x; 
  
    /* the constant in f(x). 
     * Algorithm can be re-run with a different c 
     * if it throws failure for a composite. */
    long long int c = (rand()%(n-1))+1; 
  
    /* Initialize candidate divisor (or result) */
    long long int d = 1;   
  
    /* until the prime factor isn't obtained. 
       If n is prime, return n */
    while (d==1) 
    { 
        /* Tortoise Move: x(i+1) = f(x(i)) */
        x = (modular_pow(x, 2, n) + c + n)%n; 
  
        /* Hare Move: y(i+1) = f(f(y(i))) */
        y = (modular_pow(y, 2, n) + c + n)%n; 
        y = (modular_pow(y, 2, n) + c + n)%n; 
  
        /* check gcd of |x-y| and n */
        d = __gcd(abs(x-y), n); 
  
        /* retry if the algorithm fails to find prime factor 
         * with chosen x and c */
        if (d==n) return PollardRho(n); 
    } 
  
    return d; 
} 
  
/* driver function */
int main() 
{ 
    long long int n = 10967535067; 
    printf("One of the divisors for %lld is %lld.", 
          n, PollardRho(n)); 
    return 0; 
}
Получается что один из делителей для 10967535067 - 104729

C++ (Qt)
1
2
3
$g++ -o main *.cpp
$main
One of the divisors for 10967535067 is 104723.
Но вопрос когда я задаю очень большое значение
Допустим:

C++ (Qt)
1
2
3
....
    long long int n = 56152467775622139884518581101516654419856441754822443617277137044268690521212; 
....

Я получаю что целочисленная константа слишком велика для своего типа


C++ (Qt)
1
2
3
4
5
6
$g++ -o main *.cpp
main.cpp:77:23: warning: integer constant is too large for its type
     long long int n = 56152467775622139884518581101516654419856441754822443617277137044268690521212;
                       ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
$main
One of the divisors for -376121745022567300 is 2.
Я новичок в программирование что нужно исправить в коде чтобы реализовать это?
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
10.01.2019, 02:00
Ответы с готовыми решениями:

Даны функциональная константа, предикатная константа и определённое количество аксиом. Как построить модель
Даны функциональная константа, предикатная константа и определённое количество аксиом. Как построит...

Как устранить ошибку «Система Windows остановила это устройство код 43»?
Всем привет ! Устанавливал на ноутбук HP драйвера.До установки драйверов у меня были драйвера они...

C# Есть код алгоритма Дейкстры на C++, как его можно преобразовать на язык C#?
using System; using System.Collections.Generic; using System.Linq; using System.Text; ...

Как преобразовать в строке все маленькие буквы в большие а большие в маленькие?
Дана строка .Преобразовать в ней все маленькие буквы в большие а большие в маленькие. Вот что я...

1
1389 / 1019 / 323
Регистрация: 28.07.2012
Сообщений: 2,805
10.01.2019, 10:00 2
Цитата Сообщение от DewCooper Посмотреть сообщение
что нужно исправить в коде чтобы реализовать это?
Тебе нужна длинная арифметика. Готовых реализаций в Интернете довольно много: boost, gmp и т.п.
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
10.01.2019, 10:00
Помогаю со студенческими работами здесь

Не получается запустить длинный код Алгоритма Гомори, код правильный.
Собственно как запустить код через С++Builder 6 #include&lt;ctype.h&gt; #include&lt;string.h&gt;...

ро-метод Полларда
Здравствуйте! Задание такое: Реализовать ро-метод Полларда факторизации челых чисел на примере 32...

Ошибка константа. как исправить?
помогите пожалуйста, как убрать эту ошибку?

Метод факторизации Полларда (p-1)
Весь форум облазил но не нашёл, пришлось зарегаться. Очень нужно реализовать в Pascal ABC метод...


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

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2023, CyberForum.ru