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

@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;
}

Новые блоги и статьи Все статьи Все блоги /
Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Programma_Boinc 26.12.2025 Кто-нибудь знает, где можно бесплатно получить настольный компьютер или ноутбук? США. Нашел на реддите интересную статью под названием «Кто-нибудь знает, где получить бесплатный компьютер или. . .	Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Рецензия / Мнение/ Перевод Нашел на реддите интересную статью под названием The Thinkpad X220 Tablet is the best budget school laptop period . Ниже её машинный перевод. Thinkpad X220 Tablet —. . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro
Музыка, написанная Искусственным Интеллектом volvo 05.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .

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

Решение