Алгоритм Ленстры

@theseventhline · Регистрация: 16.03.2015

Студворк — интернет-сервис помощи студентам

Нужна помощь в реализации алгоритма, не совсем все ясно в псевдокоде

Комментарий модератора

П.5.18.Правил
Запрещено размещать задания и решения в виде картинок и других файлов с их текстом.
Перепечатывайте тексты в тему.

@theseventhline · 13.04.2017, 13:03 **[ТС]**

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
57
58
59
60
61
62
63
64
#include <iostream>
 
using namespace std;
 
int gcd(int a, int b)
{
    int c;
    while (a != 0) {
        c = a; a = b % a;  b = c;
    }
    return b;
}
 
int main(){
    int r, s, n, con, r1, r2;
    int a[100], b[100], c[100], q[100];
    cin >> r >> " " >> s >> " " >> n;
    cout << endl;
    gcd(r, s); 
    r1 = (1 % s)/r;
    gcd(n, s);
    if (r2 >= 0 && r2 < s){
        r2 = r1 * n % s;
    }
    for (int i = 0; i < n; i++){
        a[0] = s;
        b[0] = c[i];
        c[0] = 0;
        if (a[1] > 0 && a[1] <= s){
            a[1] = (r1*r2) % s;
            b[1] = 1;
            c[1] = ((n - r * r2) / s) * r1 % s;
        }
        if (i >= 2){
            a[i] = a[i - 2] - q[i] * a[i - 1];
            b[i] = b[i - 2] - q[i] * b[i - 1];
            c[i] = c[i - 2] - q[i] * c[i - 1] % s;
        }
 
        con = c[i] % s;
 
        for (int i = 1; i < n; i + 2){
            if (abs(con) < s)
                con = c[i] % s;
        }
        for (int i = 2; i < n; i + 2){
            if (con >= 2 * a[i] * b[i])
                con = c[i] % s;
        }
        int x;
        x = (con - )
            if (a[i] == 0) break;
    }
    for (int i = 0; i < n; i++){
        
 
        
    }
    
 
 
    system("pause");
    return 0;
}

вот недописанный код

@theseventhline · 24.04.2017, 20:48 **[ТС]**

Готовый код, вдруг кому-то пригодится

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
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
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <algorithm>
#include <vector>
#include <string>
#include <cmath>
 
using namespace std;
 
 
int gcd(int a, int b)
{
    int c;
    while (a != 0) {
        c = a; a = b % a;  b = c;
    }
    return b;
}
 
int obr(int r, int s) {
    for (int i = 0; i < 100000; i++)
        if (((i * r) % s) == 1)
            return i;
    return -1;
}
 
int a[1000000], b[1000000], c[1000000];
 
int main() {
#ifdef _DEBUG
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
#endif
    ios_base::sync_with_stdio(false);
 
    int r, s, n;
 
    cin >> r >> s >> n;
 
    if (gcd(r, s) != 1)
        return !(cout << "ERROR: r, s not simple!");
 
    int r_star = obr(r, s);
    int r_sh = r_star * n % s;
 
    // инициализация
    a[0] = s;
    b[0] = 0;
    c[0] = 0;
 
    a[1] = r_star * r_sh % s;
    b[1] = 1;
    c[1] = (n - r * r_sh) / s * r_star % s;
 
 
 
    for (int i = 0; i < n; i++) {
 
        if (i >= 2) {
            int q = a[i - 2] / a[i - 1];
            if (i % 2 != 0)
                q = 1;
            a[i] = a[i - 2] - q * a[i - 1];
            b[i] = b[i - 2] - q * b[i - 1];
            c[i] = c[i - 2] - q * c[i - 1] % s;
        }
 
 
        vector<int> con(2, -1);
 
        for (int j = 1; j < s - 1; j += 2) {
            if (j == c[i] % s && abs(j) < s)
                if (con[0] == -1)
                    con[0] = j;
                else
                    con[1] = j;
        }
        for (int j = 2; j < s; j += 2) {
            if (j == c[i] % s
                && 2 * a[i] * b[i] <= j
                && j <= n / (s * s) + a[i] * b[i]) {
                if (con[0] == -1)
                    con[0] = j;
                else
                    con[1] = j;
            }
        }
 
 
        for (int j = 0; j < 2; j++) {
            if (con[j] == -1) continue;
            cout << con[j] << endl;
 
            int d, x1, x2, x = 0, y = 0;
            d = (con[j] * s + a[i] * r + b[i] * r_sh) * (con[j] * s + a[i] * r + b[i] * r_sh) - 4 * a[i] * b[i] * n;
            if (d < 0)
                break;
            x1 = (con[j] * s + a[i] * r + b[i] * r_sh + sqrt((double)d)) / 2;
            x2 = (con[j] * s + a[i] * r + b[i] * r_sh - sqrt((double)d)) / 2;
 
            if (a[i] == 0)
                break;
                
            if (a[i] != 0 && b[i] != 0 && ((x1 % a[i] == 0 && x2 % b[i] == 0) || (x1 % b[i] == 0 && x2 % a[i] == 0))) {
                x = (x1 - r * a[i]) / s * a[i];
                y = (x2 - r_sh * b[i]) / s * b[i];
            }
            if (x > 0 && y > 0) {
                cout << "dopolnenie: " << x*s + r << endl;
            }
            
        }
 
        if (a[i] == 0)
            break;
    }
    return 0;
}

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .

Алгоритм Ленстры

Решение