- •Основная теорема теории чисел Теорема: 4
- •Алгоритм Ферма
- •Бинарный алгоритм
- •Разложение числа в цепную дробь
- •Вычисление подходящих дробей
- •Диафантовы уравнения
- •Сравнение
- •Арифметика сравнений
- •Функция Эйлера и ее свойства.
- •Теорема Эйлера-Ферма (Малая теорема Ферма)
- •Решение сравнений в первой степени.
- •Теорема: Вильсена
- •Китайская теорема об остатках
- •Многочлены
- •Наибольший общий делитель
- •Алгоритм Евклида для 2-х многочленов
- •Частный случай теоремы 23
- •Приближенная интерполяция
- •Метод наименьших квадратов
- •Полиномиальное кодирование.
- •Коды Шеннона-Фено
- •Алгоритм Хоффмана.
- •Лексикографический порядок
- •Антилексикографический порядок
- •Алгоритм лексикографического порядка
- •Бинарный код Грея
- •Числа Стирлинга первого и второго рода
- •Числа Стирлинга
- •Числа Белла
- •Разбиение чисел
- •Рекуррентное уравнение (Уравнение в конечных разностях)
- •Линейные рекуррентные уравнения порядка k
- •Нахождение частного решения неоднородного уравнения
- •Основы теории графов
- •Связность
- •Алгоритм построения Эйлеровой цепи.
- •Деревья
- •Коциклы
- •Элементы теории вероятности
Частный случай теоремы 23
X1 |
X2 |
… |
Xn |
y1 |
y2 |
… |
yn |
xiyj
Si(x) : Si(x)=yi i=1,2…n
S(x) – многочлен проходящий через все узлы интерполяции – (xi,yi)
<S(x)>X-Xi=yi - остаток от деления S(x) на (x-xi) (Схема Горнера)
S(x)=(x-xi)+c
C(x)=(x-x1) (x-x2)… (x-xn)
Ci(x)=( ) (x-xi-1) (x-xi+1)…
<Di(x)Ci(x)>X-Xi =1
<Di(x) (x-x1) (x-x2)… (x-xi-1) (x-xi+1)… (x-xn)>X-Xi =1
Di(x) (x-x1) (x-x2)… (x-xi-1) (x-xi+1)… (x-xn)=1
S(x)=i=1nDi(x) Ci(x) Ui(x)=i=1nyi - Многочлен Лагранжа
Через n точек можно провести многочлен deg<=n-1
Недостаток: При добавлении данных надо пересчитывать весь многочлен.
Интерполяционная формула Ньютона.
Интерполяционный метод Ньютона.
Подход:
X1 |
X2 |
… |
Xn |
y1 |
y2 |
… |
yn |
S1(x) S2(x) … Sn(x) (x-x1)
S1(x)=y1 (x-x2)
…. …
Sk(x)-Sk-1(x) (x-xk-1)
x1,x2….xk-1
* a||b1 и a||b2 (b1;b2)=1 => a||b1b2
по аналогии с *
Sk(x)-Sk-1(x)||(x-x1)(x-x2) … (x-xk-1)
Sk(x)-Sk-1(x)=Ak(x-x1)(x-x2) … (x-xk-1)
Sn(x)=A1+A2(x-x1)+A3(x-x1)(x-x2)+…+An(x-x1)(x-x2) … (x-xn-1)
y1=A1 | x=x1
y2=A1+A2(x2-x1) | x=x2
Пример:
x1 x2 x3
-1 0 1
y1 y2 y3 S3(x) =A1+A2(x-x1)+A3(x-x1)(x-x2)
3 -1 2 3=A1 | x=-1
-1=A1+A2(1) | x=0 | A2=-4
2=A1+A2(2)+A3(2)(1) | A3=7/2
S3(x)=3-4(x+1)+7/2(x+1)x
S4(x)=S3(x)+A4(x-x1)(x-x2)(x-x3)
Приближенная интерполяция (Метод наименьших квадратов).
Приближенная интерполяция
F(x1,x2…xn)
F(x,y)=sin x.cos x или F(x,y,z)=x2z+xy+xy2z3
;
;
Метод наименьших квадратов
S(xk)-yk=k
k=1nk min
Sm(x)=a0+a1x+a2x2+…+amxm
k=1nk2 =k=1n(a0+a1x+a2x2+…+amxm)2
Ф(a0,a1,…,am) min {a0,…,am}
(a,b)ax+b
k=17(axk+b-yk)2=Ф(a,b)
;
Ф(a,b)=(-3a+b+2)2+(-2a+b-4)2+(-a+b)2+(b-1)2+(a+b)2+(2a+b-3)2+(3a+b-2)2
(-3a+b+2)(-3)+(-2a+b-4)(-2)+(-a+b)(-1)+(b-1)(0)+(a+b)(1)+(2a+b-3)(2)+(3a+b-2)(3)=0
Размещения и сочетания без повторений. Биномиальные коэффициенты и их свойства. Треугольник Паскаля.
Определение: способ размещения в определенном порядке некоторого числа элементов. Суммы S называют размещения.
S={a,b,c}
{ab,ac,ba,ca,bc,cb}
{ab,ac,bc}
Akn=n(n-1)(n-2)…(n-k+1)
Ann=Pn-n!
Ckn=Akn/Pk=n!/(n-k)!k!
Свойства:
1)
k=0nCkn=2n
k=0n(-1)kCkn=0
Ckn=Cn-kn
Ckn=Ck-1n-1+Ckn-1
C00=Cnn=1
Треугольник Паскаля:
С0i, С1i, … Сii (x+y)i
1
1 1
1 2 1
1 3 3 1
1 4 6 4 1
1 5 10 10 5 1
Тождество Коши:
Сkm+n=s=0kCsmCk-sn
m=n=k
Ck2k=s=0k(Csn)2
Доказательство:
M,n,k – надо выбирать
s=0kCsmCk-sn
Размещения и сочетания с повторениями. Кодирование с исправлением ошибок. Граница Хемминга.
Размещение и сочетание с повторениями.
n-k
Akn=nk n, n, n … n - k раз
Сочетания
X={a1,a2,…an)
a1 x1 xi>=0
a2 x2 x1+x2+…+xn=k
.. ..
an xn
ai - xi+1 > 0
a1a1a1| a2a2a2|…|ananan|
x1+1 x2+1 xn+1
Cn-1n+k-1
Кодирование с исправлениями (кодирование Хэмминга)
a=(a1,…,an) ai{0,1}
P<2km r – не более, чем в r разрядах могут существовать ошибки.
P – число слов.
a=(a1,…,am) | d(a,b)={ai+bi} – количество разрядов, которые не совпадают друг с b=(b1,…,bm) | другом. – Расстояние от a до b.
Пример: a=(0,1,1,0) |
b=(1,1,0,1) | => d(a,b)=3. = функция расстояния между объектами
(A,B)M (A,B) M*M R+ (декартово произведение).
(A,B)>0
(A,B)=0 <=> A=B
(A,B)=(B,A)
(A,B)<=(A,C)+(C,B)
Докажем, что d(a,b) – расстояние т.е. обладает 1,2,3.
очевидно
очевидно
d(a,b)<=d(a,c)+d(c,b)
ai=bi; 0<= …
aibi; 1<= Пусть нет вклада ai=ci |
ci=bi | => ai=bi Значит наше предположение не верно aibi;
Теорема: 26
S={a(1), a(2), … , a(p)} –
D(a(i), a(j))>=2r+1
a(i) – передали i разрядное слово, не более, чем в r разрядах ошибка.
b – переданное слово.
Известно d(a(i),b)<=r
Ai={b|d(a(i),b)<=r - окрестность с центром в точке b.
b<i
Окрестности не пересекаются:
Пусть это не так, тогда:
d(a(i),w)<=r | => d(a(i), b(i))<=d(a(i),w)+d(a(j),w)<=2+2=22 Значит наше пред-
d(a(j),w)<=r | положение неверно.
. .
. . . слова
Окрестность, любое слово на
расстоянии <r стало P<rn слов.
2m, r P?
Ai ai Ai={a, d(a,ai)<=r}
|Ai|=C0m+ C1m+ C2m+…+Crm=rCkm – число слов в любой окрестности.
PrCkm<=2m => p<=2m/rCkm
Пример: Кодирования по Хэммингу.
m=3, r=1, p<=8/C03+C13=2
a(0)=(0,0,0) |
a(1)=(1,1,1) | d(a(0), a(1))=3=2k+1
A0={(0,0,0), (1,0,0), (0,1,0), (0,0,1)}
A1={(1,1,1), (1,1,0), (0,1,1), (1,0,1)}
(0,1,1) (1,1,1)
(0,1,0) (1,1,0)
(0,0,1) (1,0,1)
(0,0,0) (1,0,0)
Полиномиальное кодирование. Код Шеннона-Фэно и алгоритм Хаффмена.