Форум программистов, компьютерный форум CyberForum.ru

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
YouDoItWrong
46 / 46 / 7
Регистрация: 29.10.2011
Сообщений: 154
#1

Алгоритм Диксона - C++

14.05.2014, 13:37. Просмотров 824. Ответов 0
Метки нет (Все метки)

Что я делаю:
1) Генерируем базис
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
std::vector<std::size_t> &PrimeNumbersProcessor::genPrimesUpToMaxLimit(const std::size_t &limit)
{
    if (limit < 2)
    {
        return m_PrimeNumbers;
    }
    else if (limit == 2)
    {
        m_PrimeNumbers.push_back(2);
        return m_PrimeNumbers;
    }
 
    m_PrimeNumbers.push_back(2);
 
    for (std::size_t i = 3; i <= limit; i++)
    {
        if (isPrime(i))
        {
            m_PrimeNumbers.push_back(i);
        }
    }
    return m_PrimeNumbers;
}
2) Генерируем вектора B и A
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 std::size_t l = static_cast<std::size_t>(std::log(num));
    std::size_t limit = static_cast<std::size_t>(exp(std::sqrt(l * log(l))));
    limit = static_cast<std::size_t>(std::sqrt(limit));
 
    m_Primes.genPrimesUpToMaxLimit(limit);
 
    while (m_Avec.size() != (m_Primes.getGeneratedPrimesList().size() + 1))
    {
        std::size_t tmpB = genRandom();
        std::size_t tmpA = PrimeNumbersProcessor::modPow(tmpB, 2, m_Num);
        if (isSmooth(tmpA))
        {
            std::cout << "This B is ok" << tmpB << std::endl;
            m_Avec.push_back(tmpA);
            m_Bvec.push_back(tmpB);
        }
        else
        {
            std::cout << "This B is BAD" << tmpB << std::endl;
        }
    }
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
bool DixonFactorization::isSmooth(const std::size_t &num)
{
    CanonicalRepresentation tmp(num);
    for (std::size_t i = 0; i < tmp.getRepresentation().size(); i++)
    {
        if (!isInBasis(tmp.getRepresentation().at(i).num))
        {
            return false;
        }
    }
    return true;
 
}
3) Записываем степени из канонического представления чисел из А в матрицу по модулю 2, решаем СЛАУ ( http://eigen.tuxfamily.org/dox/class...1SparseQR.html )
4) На выходе получаем
Кликните здесь для просмотра всего текста
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
==================A====================
35574 2^1 3^1 7^2 11^2 
36504 2^3 3^3 13^2 
53760 2^9 3^1 5^1 7^1 
600 2^3 3^1 5^2 
43264 2^8 13^2 
2904 2^3 3^1 11^2 
66825 3^5 5^2 11^1 
==================B====================
40303
20953
10565
74030
72012
65683
72700
==================Basis================
2
3
5
7
11
13
=====MATRIX=====
1 1 1 1 0 1 0 
1 1 1 1 0 1 1 
0 0 1 0 0 0 0 
0 0 1 0 0 0 0 
0 0 0 0 0 0 1 
0 0 0 0 0 0 0 
==========X-Vector============
0 0 0 0 0 0 0

т.е. все вектора линейно независимы (все элементы в "right hand" матрице считаем равными 0). Как такое может быть? Может я что-то делаю не так?
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
14.05.2014, 13:37     Алгоритм Диксона
Посмотрите здесь:

Нужен алгоритм поиска пути в этом лабиринте (будь то волновой алгоритм или алгоритм правой/левой руки ) - C++
#include &quot;stdafx.h&quot; #include &lt;iostream&gt; #include &lt;conio.h&gt; using namespace std; void lab () { int s1 = 0; int s2 =...

Волновой алгоритм поиска (Алгоритм A* / Алгоритм А стар) - C++
Хочу разработать алгоритм для решения головоломки с подвижными дисками (перестановочная головоломка). Определение. Перестано́вочные...

Линейный алгоритм, Алгоритм с ветвлениями, Циклический алгоритм Линейный алгоритм - Pascal
Линейный алгоритм, Алгоритм с ветвлениями, Циклический алгоритм Линейный алгоритм 1. Объясни, что будет напечатано программой Program...

Построить алгоритм ДО и алгоритм ПОКА для вычислений значения функции на отрезке [a,b] с шагом h. - Free Pascal
Построить алгоритм ДО и алгоритм ПОКА для вычислений значения функции на отрезке с шагом h. Написать программу: F=3+tgx Мой...

Составить алгоритм-вычисление квадрата суммы двух чисел и алгоритм для вычисления функции - Pascal ABC
Здравствуйте!Мне нужно все с самого начала и точно,помогите пожалуйста! 1.составить алгоритм-вычисление квадрата суммы двух чисел.

Построить алгоритм ДО и алгоритм ПОКА дя вычислений значения функции на отрезке [а,b] с шагом h. Написать программу - Pascal
F=3+tg x

Написать алгоритм по блок схеме (Алгоритм метода Ньютона) - Pascal
Прошу помогите, очень срочно! Сижу на экзамене!

Алгоритм грубой силы или алгоритм Бойера-Мура - Turbo Pascal
Составить алгоритм поиска заданного слова в тексте. Слово и текст являются массивами символов заданной длины. Если заданное слово...

Алгоритм поиска образа в строке. Алгоритм Бойера-Мура - Python
# Лабораторная работа № 1 # Поиск образа в строке def forming_d(pattern): &quot;&quot;&quot; Формируем массив d.&quot;&quot;&quot; d = new_p...

Алгоритм удаления невидимых линий и алгоритм закраски Гуро - Графика и игры
Как применить алгоритм удаления невидимых линий и алгоритм закраски Гуро? в данном коде, заранее спасибо за помощь #include...


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

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru