Processing math: 100%

Простые числа. Решето Эратосфена

Очень интересными с математической (и не только) точки зрения считаются простые числа. Для начала сформулируем несколько определений для дальнейшей работы.

Определение. Простое число — это натуральное число больше единицы и которое делится нацело только на единицу и на само себя. Таким образом, p считается простым, если pN,p>1,aN,a1,ap,p mod a0.

Определение. Натуральное число не являющиеся простым и больше 1 называется составным.

Примеры

  1. 3,5,7,23 — простые числа, что можно с легкостью проверить мысленно перебрав возможные делители для этих чисел. 177539 — тоже простое число, однако проверить это устным перебором делителей будет значительно сложнее.
  2. Любое четное число кроме 2 — составное, так как имеет как минимум один делитель помимо 1 и самого себя — 2.

Леммы

Сформулируем и докажем несколько лемм. Далее, если это потребуется, будем упоминать их как лемму и её номер в списке. Лемма (2), к примеру.

  1. Лемма. Пусть p и является наименьшим делителем (не считая 1) nN,n>1. Тогда p — простое число.
    Спойлер
  2. Лемма. Пусть p — наименьший (не считая 1) натуральный делитель составного числа n. Тогда pn.
    Спойлер

Решето Эратосфена

Алгоритм. Способ нахождения простых чисел до определенного n. Метод подразумевает фильтрацию чисел до n, отсеивая составные числа. Является псевдополиномиальным алгоритмом. Алгоритм заключается в следующем:

  1. Требуется выписать все числа от 2 до n.
  2. Изначально p=2.
  3. Далее вычеркнем все числа представимые в виде 2p,3p,4p, до n.
  4. Присвоим p следующее не вычеркнутое число. Будем повторять 3 и 4 шаги до тех пор, пока pn (по лемме (2)).
  5. Таким образом, все составные числа будут вычеркнуты и останутся только простые.

Замечание

Если внимательно взглянуть на алгоритм, можно заметить что мы начинаем вычеркивать с p2. Пусть kN,k>1 и k очередное простое (а значит не вычеркнутое) число в списке. А значит, что перед тем как p=k, мы вычеркнули (при условии что k>2) 2k, ведь на первом шаге мы вычеркнули все делящиеся на 2 числа. Если k>3, то и все делящиеся на 3 числа были уже вычеркнуты. То есть 3k уже вычеркнуто. Таким образом, все составные числа имеющие нетривиальные делители до k(k1) включительно уже вычеркнуты, поэтому искать число чтобы вычеркнуть стоит начиная от p2. Подробнее с модфикациями алгоритма можно ознакомится на википедии и e-max.

Пример

Найдем все простые числа до 20 с помощью решета Эратосфена. Для начала выпишем все числа. 2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20

Положим p=2 и уберем все числа от p2 до 20. Останется 2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20

Далее p=3, и мы снова убираем ненужные нам числа. 2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20

Брать следующее p не смысла, так как это будет 5, а 52>20. Таким образом мы нашли все простые числа до 20.

Тест на простые числа и решето Эратосфена

У вас есть возможность проверить то, как вы усвоили материал.

Литература

  1. Электронный конспект по алгебре. Автор Белозеров.Г.С.
  2. И.М.Виноградов. Основы теории чисел. 6-ое издание, 1952 год. стр.18-20.

Критерий совместности СЛАУ Кронекера-Капелли

Теорема Кронекера-Капелли. Критерий совместности системы линейных алгебраических уравнений. СЛАУ совместна тогда и только тогда, когда ранг матрицы системы равен рангу расширенной матрицы. То есть, если в СЛАУ r=rangA=rang˜A, где rangA — обозначает ранг матрицы системы, а rang˜A — ранг расширенной матрицы, тогда данная матрица совместна, причём система имеет единственное решение, если rangA=rang˜A=n, где n — число неизвестных, и бесконечное число решений, если rangA=rang˜A<n.

Необходимость. Пусть задана расширенная матрица ˜A:

˜A={a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2am1x1+am2x2++amnxn=bm

Скажем, что данная система совместна, в таком случае существуют числа (c1,c2,,cn), которые являются частным решением матрицы, при подстановке их в систему. Мы получим равенство:

b1b2bn=c1a11a21am1+c2a12a22am2++cna1na2namn

Следовательно, вектор-столбец свободных членов является линейной комбинацией столбцов (a1,a2,,an), матрицы A. Так же, мы можем заметить, что сколько бы мы раз не приписали или не вычеркнули строку(столбец), от этого не меняется ранг системы, из этого следует, что rangA=rang˜A.

Достаточность. Если rangA=rang˜A, то это означает, что у них один и тот же базисный минор. Тогда, согласно теореме о базисном миноре, последний столбец свободных членов – линейная комбинация столбцов базисного минора.

Следствие:

  1. rangA=rang˜A=n единственное решение.
  2. rangA=rang˜A<n бесконечное число решений.
  3. Количество главных переменных равно рангу системы.

Примеры решения задач

Рассмотрим примеры задач, в которых используеться критерий совместности rangA=rang˜A.

  1. {2x1x2+5x3=43x1x2+5x3=05x12x2+3x3=2

    Решение

    Сначала, приведем матрицу к треугольному виду.

    (215431505232)(125413502532)

    (115401040177)(115401040073)

    Элементарные преобразования не меняют ранга матриц, поэтому в результате выполненных действий, получены эквивалентные исходнной матрице системы A=(115010007) и расширенная матрица системы ˜A=(115401040073)

    rangA=rang˜A=3 значит, по теореме Кронекера-Капелли система совместна.

  2. {x1+x2x3=7x1+2x23x3=12x12x3=3

    Решение

    Приведем матрицу к ступенчистому виду:

    (111412302023)(111401240245)(1114012400013)

    ˜A=(1114012400013)=rang˜A=3

    A=(111012000)=rangA=2

    rangArang˜A. По теореме Кронекера-Капелли система линейных уравнений несовместна.

  3. {5x13x2+2x3+4x4=34x12x2+3x3+7x4=18x16x2x35x4=97x13x2+7x3+17x4=λ

    Решение

    Очевидно, что от значения λ зависит, будет ли матрица совместна или нет.

    Сначала приведем матрицу к треугольному ввиду:

    ˜A=(53243423718615973717λ)(111324237102719773717λ)

    (11132027197027197041438λ14)(11132027197000000000λ)

    При λ0: rang˜A=3, rangA=2. По теореме Кронекера-Капелли система линейных уравнений несовместна.

    При λ=0: rang˜A=2, rangA=2. По теореме Кронекера-Капелли система линейных уравнений совместна.

Критерий совместности СЛАУ Кронекера-Капелли

Тест на закрепление материала «Критерий совместности СЛАУ Кронекера-Капелли».

Литература

  1. Личный конспект, составленный на основе лекций Белозерова Г.С.
  2. Фадеев Д.К. Лекции по алгебре. М.: Наука, 1984.-416 с.  стр 119.
  3. Проскуряков И.В. Сборник задач по линейной алгебре. М.: Наука, 1984.-384 с.  стр 101-103.

Существование ортонормированного базиса

Определение. Ортонормированный базис (ОНБ) — это базисная система векторов, которая ортогональна и нормирована.

Определение. Ортогональная система векторов — это система состоящая либо из только одного ненулевого вектора, либо из нескольких ненулевых векторов, которые попарно ортогональны.

Определение. Любой вектор евклидова пространства, скалярный квадрат которого равен единице, называется нормированным. Причем любой ненулевой вектор можно нормировать. Если вектор a1=μa, при μ=(a,a)12, становится нормированным.

Определение. Система называется нормированной, если каждый вектор этой системы нормирован.

Теорема. (существование ОНБ в евклидовом пространстве) В любом конечномерном евклидовом пространстве можно найти ортонормированный базис.

Допустим, имеется система S=e1,e2,,en в евклидовом пространстве e0. Если мы возьмем произвольный вектор a из E и если бы ортонормированная система S=e1,e2,,en была бы базисом, то вектор a совпадал бы с вектором b. Тогда рассмотрим вектор ab при b=(a,e1)e1+(a,e2)e2++(a,en)en.

Тогда вектор ab: (ab,ek)=(ani=1(a,ei)ei,ek)=(a,ek)ni=1(a,ei)(ei,ek)=
=(a,ek)(a,ek)=0.
То есть вектор ab ортогонален ко всем векторам системы S=e1,e2,,en. Причем мы еще и доказали, что ab=0a=b.
Значит ЛНЗ система S=e1,e2,,en образует базис в евклидовом пространстве, т. к. векторы S линейно выражают векторы E. Таким образом, в любом конечномерном евклидовом пространстве мы можем найти ортонормированный базис, причем ортогонализировать его векторы можно процессом ортогонализации Грама-Шмидта, а нормировать по определению выше.

Смотрите также

  1. Личный конспект на основе лекций Белозерова Г.С. Евклидовы пространства
  2. Курош А.Г. Курс высшей алгебры. М.: Наука, 1968, Глава 8, §34, «Определение евклидова пространства. Ортонормированные базы», c.215
  3. Воеводин В.В. Линейная алгебра. М.: Наука, 1994, Глава 3, §28, «Ортогональность», стр. 93

Существование ортонормированного базиса

Тест на знание темы «Существование ортонормированного базиса»

Векторное произведение векторов, свойства, координатное представление

Векторное произведение векторов

Определение. Если наблюдатель, идя против часовой стрелке сначала встречает вектор c, затем встречает вектор a, затем вектор b, то тройка векторов {a,b,c} называется правой (рис. 1), если же наблюдатель шел по часовой стрелке и встретил вектора в той же последовательности, то тогда тройка векторов {a,b,c} называется левой.

Определение с использованием руки (мнемоническое правило). Если обозначить указательный палец как a, средний палец как b, а большой палец как их произведение, т.е. c, то расположение пальцев на правой руке является правой тройкой векторов, а на левой руке левой тройкой векторов.

На рисунке 1 показано как будет выглядеть правая тройка векторов.

рис. 1

Определение. Векторным произведением неколлинеарных векторов a и b называется вектор c, такой, что

  1. |c|=|a||b|sinφ, где φ — угол между векторами a и b;
  2. Вектор c ортогонален вектору a и вектору b;
  3. Тройка векторов {a,b,c} правая.

Векторное произведение a и b обозначается как [a,b].

Свойства векторного произведения

  • [a,b]=[b,a] (антикоммутативность).

    Смотря на определение видно, что произведения a×b и b×a имеют одинаковую длину. Так же они имеют противоположное направление из-за того, что sinφ нечетен.

  • [λa,b]=λ[a,b] (ассоциативность).

    Докажем данное св-во для случая λ>0, а для λ<0, доказательство проводится аналогично. Легко заметить, что при λ>0 вектор λ(a×b) имеет то же направление, что и a×b (обратное при λ<0). Теперь нам надо доказать равенство длин этих произведений. |(a×b)|=|λ||a×b|=λ|a||b|sin(a;b),

    |(λa)×b|=|λa||b|sin(a;b)=λ|a||b|sin(a;b).

  • a×(b+c)=a×b+a×c (дистрибутивность).
  • Условие коллинеарности векторов.

    Для того, чтобы два ненулевых вектора были коллинеарны, необходимо и достаточно, чтобы их векторное произведение равнялось нулевому вектору. ab,|a|0,|b|0a×b=0.

    Необходимость. Очевидно, что если вектора a и b коллинеарны, то синус угла между ними равен нулю, отсюда видим, что по определению, векторное произведение равно нулю.
    Достаточность. Теперь докажем в обратную сторону: если a×b=0, то |a||b|sin(a;b)=0 один из сомножителей равен нулю. Так как ни один из векторов не равен нулю, то sin(a;b)=0, т.е. либо ^(a;b)=0, либо ^(a;b)=π и значит ab.

    Следствие: векторный квадрат равен нулевому вектору.

  • Геометрический смысл векторного произведения.

    Модуль векторного произведения равен площади параллелограмма построенного на перемножаемых векторах (рис. 2).

    Если посмотреть векторного произведения |a×b|=|a||b|sin(a;b), то мы видим общеизвестную формулу площади параллелограмма со сторонами, длины которых равны |a| и |b|.

    рис. 2

Координатное представление векторного произведения

Для того, чтобы выразить результат векторного произведения векторов a=(ax,ay,az) и b=(bx,by,bz) в координатах надо сначала найти все парные векторные произведения единичных векторов i,j,k. i×i=j×j=k×k=0,

i×j=k,j×k=i,k×i=j,
j×i=k,k×j=i,i×k=j.
a×b=(axi+ayj+azk)×(bxi+byj+bzk)=
=axbyi×j+axbzi×k+aybxj×i+aybzj×k+azbxk×i+
+azbyk×j=axbykaxbzjaybxk+aybzi+azbxj
azbvi=(aybzazby)i(axbzazbx)j+(axbyaybx)k.

Легко заметить, что разности, стоящие в скобочках, равны определителям второго порядка. a×b=|ayazbybz|i|axazbxbz|j+|axaybxby|k.
Итак, видим, что справа от знака равно записано разложение определителя третьего порядка по первой строке. a×b=|ijkaxayazbxbybz|.
То есть c=(|ayazbybz|,|axazbxbz|,|axaybxby|).

Примеры решения задач

  1. Найти модуль векторного произведения векторов a=(0,3,4) и b=(5,12,0),φ=π3.
    Решение

  2. Найти координаты вектора c, который является результатом векторного произведения векторов a=(1,2,3) и b=(3,4,6).
    Решение
  3. Найти длины и координаты всех векторов получившихся в результате векторного умножения векторов a=(2,3,4),b=(1,3,7),c=(0,0,3) зная, что sin(a,b)=12,sin(a,c)=13,sin(b,c)=56.
    Решение
  4. Найти площадь треугольника, у которого заданы координаты его вершин. A=(1,2,3),B=(5,112),C=(3,6,4).
    Решение
  5. Найти площадь параллелограмма, построенного на векторах a=(1,3,4),AB, если A=(3,8,6),B=(2,4,7) и угол между ними равен φ=π6.
    Решение

Список литературы

  1. Ефимов Н.В.: Краткий курс аналитической геометрии, стр. 154-163
  2. Постников М.М. Аналитическая геометрия, стр 133-134
  3. Личный конспект на основе лекций Белозерова Г.С.

Векторное произведение векторов

Тест для проверки знаний по теме «Векторное произведение векторов»

Простейшие сведения о непрерывных дробях и их свойствах

Пусть число \usepackageamsfontsxR. Тогда обозначим через q1 наибольшее целое число, меньшее x. Если x не целое число, то мы получим равенство вида x=q1+1x1, так как дробь 1x1<1, то x1>1, и тогда аналогично для x1 находим такое целое q2<x1, получаем x1=q2+1x2, возвращаясь к первому равенству x=q1+1q2+1x2. Продолжая этот процесс будем получать представления последующих xk: x2=q3+1x3,

x3=q4+1x4,
xi=qi+1+1xi+1,

В итоге и получим непрерывную дробь: \usepackageamsmathx=q1+1q2+1q3++1qn1+1qn+1xn.

Далее, нам стоит рассмотреть два случая: первый — \usepackageamsfontsxQ, т. е. xрациональное число и второй — \usepackageamsfontsxRQ, т. е. xиррациональное число. Почему важны именно эти случаи?

По определению, рациональное число представимо в виде несократимой дроби mn, где \usepackageamsfontsmZ,nN, а, значит, и разложение, представленное сверху, должно быть конечным и, более того, может быть получено благодаря алгоритму Евклида.

С иррациональным числом получим ситуацию обратную — процесс можно будет продолжать неограниченно долго т. к. на каждом этапе xi будет иррационально.

Пусть xi — иррационально, тогда xi=qi+1+1xi+1,

сумма qi+1+1xi+1 — иррациональна, однако qi+1 является целым по определению, которое мы дали ему выше дробь 1xi+1 должна быть иррациональной. А это означает, что и xi+1 — иррациональное число.

Т. е. получаем, что иррациональность xi влечёт за собой иррациональность xi+1, а т. к. изначальное число x — иррационально, то и все xj, при j=1,2,3 — иррациональны.

Как было упомянуто ранее, если \usepackageamsfontsxQ, то его разложение в непрерывную дробь можно получить с помощью алгоритма Евклида.

Перед описанием алгоритма стоит ввести понятие неполного частного — это целые числа вида qi,i=¯1,n.

Опишем сам алгоритм:

Суть алгоритма заключается в том, что на каждом шаге мы будем непосредственно получать одно из неполных частных — qi, а также отношение riri1 (начиная со второго шага).

Пусть нам задано рациональное число, тогда его можно записать в виде несократимой дроби mn, где \usepackageamsfontsmZ,nN. Тогда, первый шаг: m=nq1+r1mn=q1+1mr1,

узнали значение q1, а так же получили возможность вычислить значение r1. Второй шаг: n=r1q2+r2nr1=q2+1r1r2,
узнали значение q2, а так же получили возможность вычислить значение r2. Продолжая алгоритм далее: r1=r2q3+r3r1r2=q3+1r2r3,r2=r3q4+r4r2r3=q4+1r3r4,rn2=rn1qn+rnrn2rn1=qn+1rn1rn,rn1=rnqn+1,rn1rn=qn+1.
Заканчиваем алгоритм тогда, когда получим, что очередная дробь ri1ri будет целым числом и, соответственно, qi+1 будет полным частным.

Так как найдены все неполные частные, то дробь mn можно представить в виде: \usepackageamsmathx=q1+1q2+1q3++1qn1+1qn+1qn+1.

С помощью алгоритма Евклида есть возможность найти разложение в непрерывную дробь, однако, иногда промежуточные результаты важнее конечного, а именно: δ1=q1,δ2=q1+1q2,δ3=q1+1q2+1q3,

δi называются подходящими дробями. Несложно заметить зависимость δi+1 от δi — если в записи δi число qi заменить на сумму qi+1qi+1, то мы получим δi+1.

Подходящие дроби будут нас интересовать тем, что они образуют последовательность, которая приближается к изначальному числу. Понятно, что зная все неполные частные (после применения алгоритма Евклида) можно вычислить значения всех подходящих дробей, однако, это не очень удобно и долго.

Введем специальные обозначения для нахождения значений подходящих дробей: δi=PiQi. При этом положим, что P0=1,P1=q1 и Q0=0,Q1=1. Так же стоит отметить, что в силу того, что для рационального числа x=mn непрерывная дробь конечна, то и количество подходящих дробей будет конечно, а это означает что существует равенство mn=PiQi. А так как подходящие дроби так же являются несократимыми, то равенство можно упростить до m=Pi и n=Qi. Тогда получим, что: δ1=q11=P1Q1,

δ2=q1+1q2=q2q1+1q21+0=q2P1+P0q2Q1+Q0=P2Q2,
δ3=q1+1q2+1q3=q1(q2+1q3)+1q2+1q3=q1(q2q3+1)+q3q2q3+1=
=q1q2q3+q1+q3q2q3=q3(q1q2+1)+q1q3q2+1=q3P2+P1q3Q2+Q1.

Несложно заметить рекуррентное выражение для PiQi: Pi=qiPi1+Pi2Qi=qiQi1+Qi2.

Докажем это с помощью математической индукции.

База индукции: δ3=P3Q3=q3P2+P1q3Q2+Q1.

Предположим, что δk=PkQk=qkPk1+Pk2qkQk1+Qk2.

Тогда: δk+1=(qk+1qk+1)Pk1+Pk2(qk+1qk+1)Qk1+Qk2=Pk1qk+1+qkPk1+Pk2Qk1qk+1+qkQk1+Qk2=

(выполним замену по предположению индукции) =Pk1qk+1+PkQk1qk+1+Qk=Pk1+Pkqk+1Qk1+Qkqk+1=Pk+1Qk+1,
что и требовалось доказать.

Имеет место следующее свойство подходящих дробей:

Лемма. При n>0 имеет место равенство PnQn1Pn1Qn=(1)n.

Проверим значение левой части при n=1, получим: P1Q0P0Q1=1,

далее вычислим значение левой части при увеличении индекса на 1, т. е. при n+1, получим: Pn+1QnPnQn+1=(qn+1Pn+Pn1)QnPn(qn+1Qn+Qn1)==Pn1QnPnQn1,
получили выражение противоположное заданному в условии. А, значит, при изменении индекса на единицу меняется и знак выражения, а т. к. первое значение 1, то и получаем требуемое.

Примеры решения задач

Рассмотрим примеры задач в которых могут быть использованы непрерывные дроби. Рекомендуется сначала решать примеры самому, а только затем сверить решение с представленным ниже.

  1. Разложить число x=8913 в непрерывную дробь.
    Решение

    Применяя алгоритм Евклида получим: 89=136+11,q1=6;

    13=111+2,q2=1;
    11=25+1,q3=5;
    2=12,q4=2.
    Нашли все qi можем записать x как непрерывную дробь: x=6+11+15+12

  2. Найти все подходящие дроби числа x=12719.
    Решение

    Для этого используем рекуррентные формулы подходящих дробей. Воспользуемся алгоритмом Евклида для поиска всех qi: 127=196+13,q1=6;

    19=131+6,q2=1;
    13=62+1,q3=2;
    6=16,q4=6.

    Далее будем выписывать подходящие дроби в порядке возрастания индекса: δ1=q11=61
    δ2=q2P1+P0q2Q1+Q0=16+111+0=71,
    продолжая расчеты получим: δ3=27+621+1=203
    и, наконец, δ4=620+763+1=12719=x,
    как и ожидалось четвертая подходящая дробь равна заданному числу, т. к. максимальный индекс qi был равен четырём.

  3. Разложить в непрерывную дробь иррациональное число 7.
    Решение

    Для этого воспользуемся разложением, которое было представлено в теме первым. А именно
    x0=7=q1+1x1=2+(72)

    1x1=72x1=172=7+23=1+713=1+1x2,
    x2=371=3(7+1)6=7+12=1+712=1+1x3,
    x3=271=2(7+1)6=7+13=1+723=1+1x4,
    x4=372=3(7+2)3=7+2=4+(72).

    Однако, слагаемое вида 72 у нас уже было, а значит мы пришли к циклу. Выпишем все неполные частные, они же — целые части дробей. Получим: 7=[2,¯1,1,1,4] — часть чисел находятся под чертой т. к. они находятся в цикле.

  4. Восстановить по заданным qi=[10,4,3,2,4] рациональное число x.
    Решение

    Для этого воспользуемся разложением, полученным в результате алгоритма Евклида:
    Получим дробь: 10+14+13+12+14,

    её значением и будет искомым x. Посчитав значение этой дроби получим, что x=1361133.

  5. Восстановить по заданным qi=[¯2,9] иррациональное число x.
    Решение

    Так как число x — иррациональное, то его непрерывная дробь будет бесконечной, поэтому воспользоваться методом из предыдущего примера не получится. Однако, так как мы видим по данным qi, что дробь зацикливается, то можем записать следующие выражение: x=2+19+1x

    приведем дробь к квадратному уравнению: 9x218x+2=0,
    D=324+72=396,x1,2=1±39618.
    Т. к. целая часть числа равна 2, то вариант x=139618 можно отбросить. И окончательный ответ: x=1+39618

Простейшие сведения о непрерывных дробях и их свойствах

Тест на знание темы «Простейшие сведения о непрерывных дробях и их свойствах».

Смотрите также

  1. Виноградов И.М. Основы теории чисел. cтр. 14-18
  2. Арнольд В.И. Цепные дроби.
  3. Теоретические материалы основанные на конспекте и лекциях Белозерова Г.С.