WikiDer > Рациональное сито

В математика, то рациональное сито генерал алгоритм за разложение целых чисел на простые множители. Это частный случай общее числовое поле сито. Пока меньше эффективный чем общий алгоритм, он концептуально проще. Он служит полезным первым шагом в понимании того, как работает сито общего числового поля.

Метод

Предположим, мы пытаемся учесть составное число п. Выбираем границу B, и определить факторная база (который мы будем называть п), множество всех простых чисел, меньших или равных B. Далее мы ищем положительные целые числа z так что оба z и г + п находятся B-гладкий - т.е. все их простые множители находятся в п. Поэтому мы можем написать для подходящих показателей ${ displaystyle a_ {k}}$,

${ displaystyle z = prod _ {p_ {i} in P} p_ {i} ^ {a_ {i}}}$

а также для подходящих ${ displaystyle b_ {i}}$, у нас есть

${ Displaystyle Z + N = prod _ {p_ {i} in P} p_ {i} ^ {b_ {i}}}$.

Но ${ displaystyle z}$ и ${ displaystyle z + n}$ конгруэнтны по модулю ${ displaystyle n}$, поэтому каждое такое целое число z что мы находим, дает мультипликативное соотношение (мод п) среди элементов п, т.е.

${ displaystyle prod _ {p_ {i} in P} p_ {i} ^ {a_ {i}} Equiv prod _ {p_ {i} in P} p_ {i} ^ {b_ {i} } { pmod {n}}}$

(где а_я и б_я являются неотрицательными целыми числами.)

Когда мы сгенерировали достаточно этих отношений (обычно достаточно, чтобы количество отношений было на несколько больше, чем размер п) можно использовать методы линейная алгебра умножить эти различные отношения таким образом, чтобы все показатели простых чисел были четными. Это даст нам соответствие квадратов формы²≡b² (мод п), который можно превратить в факторизацию п = gcd(а-б,п) × gcd (а+б,п). Эта факторизация может оказаться тривиальной (т.е. п=п× 1), и в этом случае мы должны повторить попытку с другой комбинацией отношений; но если повезет, мы получим нетривиальную пару множителей п, и алгоритм завершится.

Пример

Мы будем множить целое число п = 187 с использованием рационального сита. Мы произвольно попробуем значение B= 7, что дает факторную базу п = {2,3,5,7}. Первый шаг - проверить п на делимость каждым из членов п; Очевидно, что если п делится на одно из этих простых чисел, то мы уже закончили. Однако 187 не делится на 2, 3, 5 или 7. Затем мы ищем подходящие значения z; первые несколько - 2, 5, 9 и 56. Четыре подходящих значения z дают четыре мультипликативных отношения (mod 187):

• 2¹3⁰5⁰7⁰ = 2 ≡ 189 = 2⁰3³5⁰7¹.............(1)
• 2⁰3⁰5¹7⁰ = 5 ≡ 192 = 2⁶3¹5⁰7⁰.............(2)
• 2⁰3²5⁰7⁰ = 9 ≡ 196 = 2²3⁰5⁰7².............(3)
• 2³3⁰5⁰7¹ = 56 ≡ 243 = 2⁰3⁵5⁰7⁰.............(4)

Теперь есть несколько принципиально разных способов их комбинирования и получения четных показателей. Например,

• (1) × (4): после их умножения и сокращения общего множителя 7 (что мы можем сделать с 7, будучи членом п, уже было определено, что они взаимно просты с п^[1]), это сводится к 2⁴ ≡ 3⁸ (мод п) или 4² ≡ 81² (мод п). Результирующая факторизация: 187 = НОД (81-4 ​​187) × НОД (81 + 4 187) = 11 × 17.

В качестве альтернативы, уравнение (3) уже находится в правильной форме:

• (3): Это говорит 3² ≡ 14² (мод п), что дает факторизацию 187 = НОД (14-3,187) × НОД (14 + 3,187) = 11 × 17.

Ограничения алгоритма

Рациональное сито, как и сито общего числового поля, не может разложить на множители числа вида п^м, куда п это простое и м целое число. Однако это не является большой проблемой - такие числа статистически редки, и, кроме того, существует простой и быстрый процесс проверки того, имеет ли данное число такую ​​форму. Наверное, самый элегантный метод - проверить, ${ displaystyle lfloor n ^ {1 / b} rfloor ^ {b} = n}$ выполняется для любого 1 Метод Ньютона для удаления корня.^[2]

Самая большая проблема - найти достаточное количество z так что оба z и z+п находятся B-гладкий. Для любого данного B, доля чисел, которые B-smooth быстро уменьшается с увеличением числа. Так что если п большой (скажем, сто цифр), будет сложно или невозможно найти достаточно z чтобы алгоритм работал. Преимущество сита общего числового поля состоит в том, что нужно искать только гладкие числа порядка. п^1/d для некоторого положительного целого числа d (обычно 3 или 5), а не по порядку п как требуется здесь.

Рекомендации