Факторизация методом Ферма

@AlexProg · Регистрация: 03.11.2012

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

не пойму в чем проблема , реализую факторизацию методом Ферма , не пойму почему не работает

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
#include<iostream>
#include<clocale>
#include<cmath>
using namespace std;
void faktor(double n);
int POW(int k,int d);
int main()
{   
    double n=19691;
    faktor(n);  
    getchar();
    return 0;
}
void faktor(double n)
{
    double A;
    int m;
    int p;
    int q;
    int B;
    m = sqrt(n);
    for(int x=0;x<m;++x)
        A=m+x;
        B=sqrt(POW(A,2)-n);
        p=A+B;
        q=A-B;
        if(p*q==n)
        {
            cout<<"простые делитеои числа "<<n<<"числа "<<p<<" и "<<q;        
        }
}
 
int POW(int k,int d)
{
    int e=k;
    for(int i=0;i<d-1;++i)
        k=k*e;
    return k;
}

@Andrey1302 · 26.03.2018, 01:05

Там же все очень просто. Вот на си и на асме (вторая работает раза в 2-3 быстрее):

C
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
unsigned __int128 fer_c(unsigned __int128 n) {//Ферма__128 на С++
unsigned __int128 q=sqr128(n); __int128 a=q*q-n;
for(n=0;a!=0;) a>0 ? a-=(n++<<1)+1 : a+=(q++<<1)+1;
q==n ? q<<=1 : q-=n; if(q==1) return 0; return q;}
 
unsigned __int128 fer_asm(unsigned __int128 n) {unsigned __int64 ql=sqr128(n), qh;// Ферма__128 на asm
asm("xor %%rcx,%%rcx\n" "mov %%rbx,%%rax\n" "mul %%rbx\n" "sub %%rsi,%%rax\n" "sbb %%rdi,%%rdx\n" "xor %%rsi,%%rsi\n" "xor %%rdi,%%rdi\n"
".p2align 4\n" "MF128:\n" "test %%rdx,%%rdx\n" "jl MFS128\n" "jnz MFS0128\n" "test %%rax,%%rax\n" "jz RF128\n"
".p2align 4\n" "MFS0128:\n" "mov %%rsi,%%r8\n" "mov %%rdi,%%r9\n" "shl %%r9\n" "shl %%r8\n" "adc $0,%%r9\n" "add $1,%%r8\n" "adc $0,%%r9\n"
"sub %%r8,%%rax\n" "sbb %%r9,%%rdx\n" "add $1,%%rsi\n" "adc $0,%%rdi\n" "jmp MF128\n"
".p2align 4\n" "MFS128:\n" "mov %%rbx,%%r8\n" "mov %%rcx,%%r9\n" "shl %%r9\n" "shl %%r8\n" "adc $0,%%r9\n" "add $1,%%r8\n" "adc $0,%%r9\n"
"add %%r8,%%rax\n" "adc %%r9,%%rdx\n" "add $1,%%rbx\n" "adc $0,%%rcx\n" "jmp MF128\n"
"RF128:\n" "cmp %%rdi,%%rcx\n" "jne MRF128\n" "cmp %%rsi,%%rbx\n" "jne MRF128\n" "shl %%rcx\n" "shl %%rbx\n" "adc $0,%%rcx\n" "jmp EF128\n"
"MRF128:\n" "sub %%rsi,%%rbx\n" "sbb %%rdi,%%rcx\n"
"EF128:\n" :"=c"(qh),"=b"(ql): "S"((unsigned __int64)n),"D"((unsigned __int64)(n>>64)),"b"(ql));
n=qh; n=(n<<64)+ql; if(n==1) return 0; return n;}

@perevertysh · 26.03.2018, 01:22

вы тут всей группой перед семинаром отрываетесь?

Факторизация больших чисел

Факторизация числа

Новые блоги и статьи Все статьи Все блоги /
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Рецензия / Мнение/ Перевод https:/ / **********/ gallery/ thinkpad-x220-tablet-porn-gzoEAjs . . .	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 04.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 - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .

@perevertysh 5 / 5 / 6 Регистрация: 23.03.2018 Сообщений: 98
	26.03.2018, 01:22
	вы тут всей группой перед семинаром отрываетесь? Факторизация больших чисел Факторизация числа 0