Matematyka dyskretna 1/Wykład 7: Funkcje tworzące: Różnice pomiędzy wersjami

Z Studia Informatyczne
Przejdź do nawigacjiPrzejdź do wyszukiwania
Pitab (dyskusja | edycje)
Nie podano opisu zmian
Pitab (dyskusja | edycje)
Linia 170: Linia 170:


Przyjrzyjmy się teraz rozwinięciu w szereg funkcji <math>\left(1+x\right)^y</math>, gdzie <math>y\in\mathbb{R}</math> jest parametrem. Rozwinięcie takie okaże się bardzo przydatne w rozwiązywaniu wielu przykładów. Aby poznać ciąg odpowiadający tej funkcji wprowadźmy definicję.
Przyjrzyjmy się teraz rozwinięciu w szereg funkcji <math>\left(1+x\right)^y</math>, gdzie <math>y\in\mathbb{R}</math> jest parametrem. Rozwinięcie takie okaże się bardzo przydatne w rozwiązywaniu wielu przykładów. Aby poznać ciąg odpowiadający tej funkcji wprowadźmy definicję.
{{kotwica|usymdwu|'''Uogólniony symbol dwumianowy'''}} <math>{ y \choose n }</math>, gdzie <math>y\in\mathbb{R}</math> oraz <math>n\in\mathbb{N}</math> jest oznaczeniem na
<center><math>{ y \choose n }\ =\ \frac{y^{\underline{n}}}{n!}\ =\
\frac{y\cdot\left(y-1\right)\cdot\ldots\cdot\left(y-\left(n-1\right)\right)}{1\cdot2\cdot\ldots\cdot\left(n-1\right)\cdot n}.
</math></center>
{{uwaga|||
Oczywiście dla <math>y\in\mathbb{N}</math> spełniającego dodatkowo <math>y\geq n</math>, uogólniony symbol dwumianowy <math>{ y \choose n }</math> jest liczbą <math>n</math>-elementowych podzbiorów zbioru <math>y</math>-elementowego.
}}
{{twierdzenie|7.4|tw 7.4|
Dla liczby rzeczywistej <math>y</math> oraz liczby naturalnej <math>n</math> zachodzi
<center><math>\left(1+x\right)^y=\sum_{n=0}^{\infty}{ y \choose n }x^n.
</math></center>
}}
{{wniosek|7.5|wn 7.5|
Dla liczby naturalnej <math>m</math> zachodzi
<center><math>\frac{1}{\left(1-x\right)^{m+1}}=\sum_{n=0}^{\infty}{ m+n \choose n }x^n.
</math></center>}}
{{dowod|||
Dowód zostawiony jest jako ćwiczenie ['''ex''']['''ex newton for integer'''].
}}
{{przyklad|||
Policzmy sumę
<center><math>\sum_{k=0}^nk^2=1+4+9+\ldots+n^2.
</math></center>
Zacznijmy od znalezienia zwartej postaci funkcji tworzącej
<math>\fGen{G}{x}=\sum_{n=0}^{\infty}n^2x^n</math>. Korzystając z Wniosku [[#wn_7.5|7.5] otrzymujemy:
{{wzor|wzor_8|8|
<math>\frac{1}{1-x}&=&\sum_{n=0}^{\infty}{n \choose n}x^n=\sum_{n=0}^{\infty}x^n,</math>}}
{{wzor|wzor_9|9|
\frac{1}{\left(1-x\right)^2}&=&\sum_{n=0}^{\infty}{n+1 \choose n
}x^n\ =\ \sum_{n=0}^{\infty}nx^n+\sum_{n=0}^{\infty}x^n.
</math>}}
Po przekształceniu równości ([[#wzor_9|9]]) uzyskuje się
{{wzor|wzor_10|10|
<math>
\sum_{n=0}^{\infty}nx^n= \frac{1}{\left(1-x\right)^2}
-\frac{1}{1-x}.
</math>}}
Powołując się ponownie na Wniosek [[#wn_7.5|7.5]] otrzymujemy
<center><math>\frac{1}{\left(1-x\right)^3}
=\sum_{n=0}^{\infty}{ n+2 \choose n}x^n
=\frac{1}{2}\sum_{n=0}^{\infty}n^2x^n+\frac{3}{2}\sum_{n=0}^{\infty}nx^n+\sum_{n=0}^{\infty}x^n,
</math></center>
co w połączeniu z równościami ([[#wzor_9|9]]) oraz ([[#wzor_10|10]])
daje zwartą postać funkcji tworzącej <math>\fGen{G}{x}</math> dla ciągu <math>1,4,9,\ldots,n^2,\ldots</math>:
<center><math>\fGen{G}{x}=\sum_{n=0}^{\infty}n^2x^n
=\frac{2}{\left(1-x\right)^3}-\frac{3}{\left(1-x\right)^2}+\frac{1}{1-x}.
</math></center>
Naszym zadaniem było jednakże policzenie funkcji tworzącej <math>H(x)</math> dla ciągu <math>1,1+4,1+4+9,\ldots,1+4+9+\ldots+n^2,\ldots</math>, tzn. ciągu sum początkowych wyrazów ciągu <math>1,4,9,\ldots,n^2,\ldots</math>. Aby uzyskać <math>\fGen{H}{x}</math> wystarczy więc skorzystać ze wzoru ([[#wzor_7|7]]) i podzielić <math>\fGen{G}{x}</math> przez <math>1-x</math>.
Tak więc poszukiwanym rozwiązaniem są współczynniki funkcji tworzącej
<center><math>\fGen{H}{x}=\frac{\fGen{G}{x}}{1-x}
=\frac{2}{\left(1-x\right)^4}-\frac{3}{\left(1-x\right)^3}+\frac{1}{\left(1-x\right)^2}.
</math></center>
Korzystając po raz kolejny z Wniosku [[#wn_7.5|7.5]] otrzymujemy
<center><math>\aligned\fGen{H}{x}
&=&2\sum_{n=0}^{\infty}{n+3 \choose n}x^n-3\sum_{n=0}^{\infty}{n+2 \choose n}x^n+\sum_{n=0}^{\infty}{n+1 \choose n}x^n\\
&=&\sum_{n=0}^{\infty}\left(\frac{1}{3}n^3+\frac{1}{2}n^2+\frac{1}{6}n\right)x^n.
\endaligned</math></center>
W konsekwencji zachodzi równość
<center><math>\sum_{k=1}^nk^2=\vect{x^n}\fGen{H}{x}=\frac{2n^3+3n+n}{6}.
</math></center>
{{przyklad|||
Wracamy  do przykładu z monetami. Występowały tam funkcje tworzące postaci
<center><math>\fGen{A_k}{x} = 1+x^k+x^{2k}+x^{3k}+\ldots,
</math></center>
dla <math>k=1,5,10,25</math> i <math>50</math>.
Z równości ([[#wzor_7|7]]) wiemy, że
<center><math>1+x^k+x^{2k}+x^{3k}+\ldots, =\frac{1}{1-x^k}
</math></center>
tak więc:
<center><math>\aligned\fGen{A}{x}= \fGen{A_1}{x}&=& \frac{1}{1-x},\\
\fGen{B}{x}= \fGen{A}{x}\cdot \fGen{A_5}{x} &=&\frac{\fGen{A}{x}}{1-x^5},\\
\fGen{C}{x}= \fGen{B}{x}\cdot \fGen{A_{10}}{x} &=&\frac{\fGen{B}{x}}{1-x^{10}},\\
\fGen{D}{x}= \fGen{C}{x}\cdot \fGen{A_{25}}{x} &=&\frac{\fGen{C}{x}}{1-x^{25}},\\
\fGen{E}{x}= \fGen{D}{x}\cdot \fGen{A_{50}}{x} &=&\frac{\fGen{D}{x}}{1-x^{50}},
\endaligned</math></center>
skąd natychmiast:
<center><math>\aligned\fGen{A}{x}&=&1+x\fGen{A}{x},\\
\fGen{B}{x}&=&\fGen{A}{x}+x^5\fGen{B}{x},\\
\fGen{C}{x}&=&\fGen{B}{x}+x^{10}\fGen{C}{x},\\
\fGen{C}{x}&=&\fGen{D}{x}+x^{25}\fGen{C}{x},\\
\fGen{D}{x}&=&\fGen{E}{x}+x^{50}\fGen{D}{x}.
\endaligned</math></center>
Równości te dają zależności między współczynnikami:
<center><math>a_n=1,\quad b_n=a_n+b_{n-5},\quad c_n=b_n+c_{n-10},\quad
d_n=c_n+d_{n-25},\quad e_n=d_n+e_{n-50}.
</math></center>
Wykorzystując te zależności rekurencyjne możemy wypełnić następującą tabelę:
{| border=1
|+ <span style="font-variant:small-caps">Uzupelnij tytul</span>
|-
|
<math>n</math> || 0 || 5 || 10 || 15 || 2 || 25 || 30 || 35 || 40 || 45 || 50 || 55 || 60 || 65 || 70 || 75 || 80 || 85 || 90 || 95 || 100
|-
|
<math>a_n</math> || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1 || 1
|-
|
<math>b_n</math> || 1 || 2 || 3 || 4 || 5 || 6 || 7 || 8 || 9 || 10 || 11 || 12 || 13 || 14 || 15 || 16 || 17 || 18 || 19 || 20 || 21
|-
|
<math>c_n</math> || 1 || 2 || 4 || 6 || 9 || 12 || 16 || 10 || 25 || 30 || 36 || 42 || 49 || 56 || 64 || 72 || 81 ||  || 100 ||  || 121
|-
|
<math>d_n</math> || 1 ||  ||  ||  ||  || 13 ||    ||    ||    ||    || 49 ||    ||    ||    ||  || 121 ||    ||  ||    ||  || 242
|-
|
<math>e_n</math> || 1 ||  ||  ||  ||  ||    ||    ||    ||    ||    || 50 ||    ||    ||    ||  ||    ||    ||  ||    ||  || 292
|-
|
|}
Pół dolara można rozmienić na <math>50</math> sposobów.
Z kolei rozmieniać jednego dolara można na aż <math>292</math> sposoby.
Do problemu tego wrócimy jeszcze w następnym wykładzie.}}
==Funkcje tworzące w rozwiązywaniu zależności rekurencyjnych==

Wersja z 09:51, 21 sie 2006

Przykład

Słynny matematyk Georg Pólya rozważał problem polegający na policzeniu wszystkich możliwych sposobów, na które można rozmienić 50 centów używając jednocentówek (1), pięciocentówek (5), dziesięciocentówek (10), ćwierćdolarówek (25), oraz półdolarówki (50). Rozważania te doprowadziły go do użycia analitycznych metod funkcji tworzących w zaproponowanym przez niego rozwiązaniu. W tym i następnym wykładzie poznamy te metody i zobaczymy jak mogą być pomocne w zliczaniu rożnych obiektów kombinatorycznych.

Wracając do problemu rozmieniania monet, wygodnie nam będzie posiadać jeszcze monetę [0], którą możemy interpretować jako brak monet. Wypiszmy teraz (nadużywając trochę notacji) nieskończoną sumę wszystkich możliwości rozmiany dowolnej kwoty za pomocą jednocentówek


Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle A_1=[0]+\moneta(1)+\moneta(1)\moneta(1)+\moneta(1)\moneta(1)\moneta(1)+\moneta(1)\moneta(1)\moneta(1)\moneta(1)+\ldots }


i analogicznie przeanalizujmy sumę dla pieciocentówek


Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle A_5=[0]+\moneta(5)+\moneta(5)\moneta(5)+\moneta(5)\moneta(5)\moneta(5)+\moneta(5)\moneta(5)\moneta(5)\moneta(5)+\ldots }


Wtedy zbiór par A1×A5 jest zbiorem wszystkich możliwości rozmiany kwoty mając do dyspozycji dowolnie wiele jednocentówek oraz pięciocentówek.


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned B= A_1 \times A_5 &=&\left([0]+\moneta(1)+\moneta(1)\moneta(1)+\moneta(1)\moneta(1)\moneta(1)+\moneta(1)\moneta(1)\moneta(1)\moneta(1)+\ldots\right)\\ &&\times\left([0]+\moneta(5)+\moneta(5)\moneta(5)+\moneta(5)\moneta(5)\moneta(5)+\moneta(5)\moneta(5)\moneta(5)\moneta(5)+\ldots\right)\\ &=&[0]+\moneta(1)+\moneta(5)+\moneta(1)\moneta(1)+\moneta(1)\moneta(5)+\moneta(5)\moneta(5)+\moneta(1)\moneta(1)\moneta(1)+\ldots \endaligned}


Sumy wszystkich możliwości rozmiany za pomocą dziesięciocentówek Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(10)} , ćwierćdolarówek Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(25)} , oraz półdolarówek Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(50)} wyglądają następująco:


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned A_{10} &=& [0]+\moneta(10)+\moneta(10)\moneta(10)+\moneta(10)\moneta(10)\moneta(10)+\moneta(10)\moneta(10)\moneta(10)\moneta(10)+\ldots\\ A_{25} &=& [0]+\moneta(25)+\moneta(25)\moneta(25)+\moneta(25)\moneta(25)\moneta(25)+\moneta(25)\moneta(25)\moneta(25)\moneta(25)+\ldots\\ A_{50} &=& [0]+\moneta(50)+\moneta(50)\moneta(50)+\moneta(50)\moneta(50)\moneta(50)+\moneta(50)\moneta(50)\moneta(50)\moneta(50)+\ldots. \endaligned}


Dodając kolejno monety Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(10)} , Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(25)} , i na końcu Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(50)} do możliwych rozmian uzyskujemy odpowiednio:


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned C&=&B\times\left([0]+\moneta(10)+\moneta(10)\moneta(10)+\moneta(10)\moneta(10)\moneta(10)+\moneta(10)\moneta(10)\moneta(10)\moneta(10)+\ldots\right)\\ D&=&C\times\left([0]+\moneta(25)+\moneta(25)\moneta(25)+\moneta(25)\moneta(25)\moneta(25)+\moneta(25)\moneta(25)\moneta(25)\moneta(25)+\ldots\right)\\ E&=&D\times\left([0]+\moneta(50)+\moneta(50)\moneta(50)+\moneta(50)\moneta(50)\moneta(50)+\moneta(50)\moneta(50)\moneta(50)\moneta(50)+\ldots\right)\\ &=&[0]+\moneta(1)+\moneta(5)+\moneta(10)+\moneta(25)+\moneta(50)+\moneta(1)\moneta(1)+\moneta(1)\moneta(5)+\moneta(1)\moneta(10)+\ldots \endaligned}


Grupując teraz składniki sumy E w podsumy o tych samych wartościach dostajemy wyrażenie:


Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \begin{array} {rcl} E&=&\big(\moneta(1)\big)+\big(\moneta(1)\moneta(1)\big)+\big(\moneta(1)\moneta(1)\moneta(1)\big)+\big(\moneta(1)\moneta(1)\moneta(1)\moneta(1)\big)\\ &&+\big(\moneta(1)\moneta(1)\moneta(1)\moneta(1)\moneta(1)+\moneta(5)\big)\\ &&+\big(\moneta(1)\moneta(1)\moneta(1)\moneta(1)\moneta(1)\moneta(1)+\moneta(5)\moneta(1)\big)+\ldots \end{array} }      (1)


Zliczając zaś tylko składniki w podsumie odpowiadającej wartości n centów, otrzymujemy liczbę sposobów, na które można rozmienić n centów przy użyciu monet Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(1)} , Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(5)} , Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(10)} , Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(25)} , oraz Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(50)} . Pomysłem pochodzącym od Pólya, było zastąpienie monety Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(1)} przez zmienną x, monety Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(5)} przez xxxxx=x5 i analogicznie Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(10)} przez x10, Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(25)} przez x25, oraz Parser nie mógł rozpoznać (nieznana funkcja „\moneta”): {\displaystyle \moneta(50)} przez x50. Uzyskujemy w ten sposób nieskończony szereg zmiennej x:


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned\fGen{E}{x}&=&\left(1+x+x^2+x^3\ldots\right)\cdot\left(1+x^5+x^{10}+x^{15}\ldots\right)\cdot\left(1+x^{10}+x^{20}+x^{30}\ldots\right)\\ &&\cdot\left(1+x^{25}+x^{50}+x^{75}\ldots\right)\cdot\left(1+x^{50}+x^{100}+x^{150}\ldots\right)\\ &=&1+x+x^2+x^3+x^4+2x^5+2x^6+2x^7+2x^8+2x^9+4x^{10}+\ldots \endaligned}


Godne zauważenia jest, że liczba różnych możliwych sposobów rozmiany n centów (równa liczbie grup monet w odpowiednim nawiasie we wzorze (1)) jest równa współczynnikowi stojącemu przy jednomianie xn.

Funkcja tworząca Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} dla ciągu liczb rzeczywistych (lub zespolonych) (g0,g1,g2,g3,) to szereg funkcyjny zmiennej rzeczywistej (lub zespolonej) x postaci


Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=\sum_{n=0}^{\infty}{g_nx^n}=g_0+g_1x+g_2x^2+g_3x^3+g_4x^4+\ldots. }


Na oznaczenie współczynnika n-tego wyrazu szeregu Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} używać będziemy oznaczenia Parser nie mógł rozpoznać (nieznana funkcja „\vect”): {\displaystyle \vect{x^n}\fGen{G}{x}=g_n} .

Uwaga Jak traktowac funkcje tworzące

Na funkcje tworzące można spojrzeć dwoiście. Pierwszym sposobem jest potraktowanie Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} jako szeregu liczb rzeczywistych (lub ogólniej zespolonych). Oczywistym pytaniem jest tu kwestia zbieżności szeregu Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=\sum_{n=0}^{\infty}{g_nx^n}} . Z wykładu Analiza Matematyczna wiemy, że szereg Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} jest zbieżny wtedy i tylko wtedy, gdy istnieje stała M0 ograniczająca wszystkie skończone początkowe sumy, tzn.


|g0|+|g1x|+|g2x2|++|gnxn|M


zachodzi dla dowolnego n0. Ponadto jeśli dla pewnej liczby x0 szereg Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x_0}=g_0+g_1x_0+g_2x_0^2+\ldots} jest zbieżny, to i także szereg Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x_1}=g_0+g_1x_1+g_2x_1^2+\ldots} jest zbieżny dla dowolnego x1 spełniającego |x1||x0|. Możemy więc określić promień zbieżności szeregu jako taką liczbę Parser nie mógł rozpoznać (nieznana funkcja „\vect”): {\displaystyle r\in\mathbb{R}_*\cup{\left\{ {\infty} \right\}\ }=\vect{0,+\infty}} , że

  • jeśli x<r, to Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} jest zbieżny;
  • jeśli x>r, to Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} jest rozbieżny.

Szereg Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=g_0+g_1x+g_2x^2+\ldots} można więc potraktować jako funkcję


G:(r,r),


o wartościach Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=\lim_{n\rightarrow\infty}{\left(g_0+g_1x+g_2x^2+\ldots+g_nx^n\right)}.} Oczywiście Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{0}=g_0} , więc dla x=0 szereg Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} jest zbieżny.

Drugim podejściem, bardziej użytecznym w praktycznych obliczeniach i przekształceniach jest spojrzenie na szereg Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=g_0+g_1x+g_2x^2+\ldots} jako formę zapisu ciągu (g0,g1,g2,), czyli jedynie jako ciąg symboli. Równości pomiędzy odpowiednimi wzorami służą rozwiązaniu problemów kombinatorycznych, tak więc traktujemy je jako równości dwu wyrażeń, a nie jako równość dwu funkcji rzeczywistych, pomimo że mają one uzasadnienia w języku analizy matematycznej.

Jak zobaczymy na wielu przykładach, funkcje tworzące są bardzo użytecznym narzędziem przy wyznaczaniu wartości elementów ciągu. Jeśli bowiem Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=g_0+g_1x+g_2x^2+\ldots} jest funkcją tworzącą ciągu (g0,g1,g2,g3,), oraz w jakiś sposób będziemy w stanie poznać postać zwartą funkcji G(x), to rozwijając tę postać zwartą w szereg Taylora, poznamy kolejne współczynniki tego rozwinięcia. A współczynniki te, to właśnie kolejne wyrazy naszego ciągu.

Będziemy się zajmowali jedynie tymi funkcjami, dla których promień zbieżności r>0. Ponadto będziemy pomijać problem zbieżności oraz wartość r promienia zbieżności, skupiając się jedynie na przekształceniach wzorów. Poniżej zebrane zostały te własności, które często wykorzystywane są w takich przekształceniach.

Obserwacja 7.1

Dla dwu funkcji tworzących Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{F}{x}=f_0+f_1x+f_2x^2+\ldots} oraz Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=g_0+g_1x+g_2x^2+\ldots} mamy:


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned\fGen{F}{x}=\fGen{G}{x}&\Leftrightarrow& f_0=g_0,\ f_1=g_1,\ f_2=g_2,\ \ldots\\ &&\\ \alpha\cdot\fGen{F}{x}+\beta\cdot\fGen{G}{x}&=& \sum_{n=0}^{\infty}{\left(\alpha\cdot f_n+\beta\cdot g_n\right)x^n}\\ &=&\left(\alpha\cdot f_0+\beta\cdot g_0\right) + \left(\alpha\cdot f_1+\beta\cdot g_1\right)x + \left(\alpha\cdot f_2+\beta\cdot g_2\right)x^2 + \ldots\\ &&\\ \fGen{F}{x}\cdot\fGen{G}{x}&=&\sum_{n=0}^{\infty}\left(\sum_{k=0}^n f_k g_{n-k}\right) x^n\\ &=& f_0g_0 + \left(f_0g_1+f_1g_0\right)x\\ && + \left(f_0g_2+f_1g_1+f_2g_0\right)x^2\\ && + \left(f_0g_3+f_1g_2+f_2g_1+f_3g_0\right)x^3+\ldots\\ \endaligned}


Wyrażenie Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{F}{x}\cdot\fGen{G}{x}} nazywać będziemy splotem szeregów Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{F}{x}} oraz Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} .

Twierdzenie 7.2

Funkcja tworząca postaci


Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=g_0+g_1x+g_2x^2+g_3x^3+\ldots }


ma odwrotną względem mnożenia (splotu), tzn. istnieje funkcja tworząca Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{U}{x}} taka, że Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{U}{x}\fGen{G}{x}=1} , wtedy i tylko wtedy, gdy g00.

Następne własności są bardzo pomocne w dokonywanych przekształceniach funkcji tworzących.

Obserwacja 7.3

Dla dwu funkcji tworzących Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{F}{x}=f_0+f_1x+f_2x^2+\ldots} oraz Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=g_0+g_1x+g_2x^2+\ldots} mamy:


Parser nie mógł rozpoznać (nieznana funkcja „\alignedx”): {\displaystyle \alignedx^m\fGen{G}{x}&=&0+\ldots+0x^{m-1}+g_0x^m+g_1x^{m+1}+g_2x^{m+2}+\ldots\\ \frac{\fGen{G}{x}-\sum_{i=0}^{m-1}{g_ix^i}}{x^{m}}&=&g_m+g_{m+1}x+g_{m+2}x^{2}+g_{m+3}x^{3}+g_{m+4}x^{4}+\ldots\\ \fGen{G}{\alpha x}&=&g_0+g_1\alpha x+g_2\alpha^2x^2+g_3\alpha^3x^3+g_4\alpha^4x^4+\ldots\\ \fGen{G'}{x}&=&g_1+2g_2x+3g_3x^2+4g_4x^3+5g_5x^4+\ldots\\ \int \fGen{G}{x}dx &=& 0+g_0x+\frac{1}{2}g_1x^2+\frac{1}{3}g_2x^3+\frac{1}{4}g_3x^4+\ldots\\ \frac{\fGen{G}{x}}{1-x}&=&g_0+\left(g_0+g_1\right)x+\left(g_0+g_1+g_2\right)x^2+\ldots \endaligned}


Funkcje tworzące w zliczaniu

Widzieliśmy już, że dla n


(1+x)m=(m0)x0+(m1)x+(m2)x2++(mm1)xm1+(mm)xm=n=0(mn)xn.


Przyjrzyjmy się teraz rozwinięciu w szereg funkcji (1+x)y, gdzie y jest parametrem. Rozwinięcie takie okaże się bardzo przydatne w rozwiązywaniu wielu przykładów. Aby poznać ciąg odpowiadający tej funkcji wprowadźmy definicję.

Uogólniony symbol dwumianowy (yn), gdzie y oraz n jest oznaczeniem na


Parser nie mógł rozpoznać (błąd składni): {\displaystyle { y \choose n }\ =\ \frac{y^{\underline{n}}}{n!}\ =\ \frac{y\cdot\left(y-1\right)\cdot\ldots\cdot\left(y-\left(n-1\right)\right)}{1\cdot2\cdot\ldots\cdot\left(n-1\right)\cdot n}. }


Uwaga

Oczywiście dla y spełniającego dodatkowo yn, uogólniony symbol dwumianowy (yn) jest liczbą n-elementowych podzbiorów zbioru y-elementowego.

Twierdzenie 7.4

Dla liczby rzeczywistej y oraz liczby naturalnej n zachodzi


(1+x)y=n=0(yn)xn.


Wniosek 7.5

Dla liczby naturalnej m zachodzi


1(1x)m+1=n=0(m+nn)xn.


Dowód

Dowód zostawiony jest jako ćwiczenie [ex][ex newton for integer].

{{przyklad||| Policzmy sumę


k=0nk2=1+4+9++n2.


Zacznijmy od znalezienia zwartej postaci funkcji tworzącej Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=\sum_{n=0}^{\infty}n^2x^n} . Korzystając z Wniosku [[#wn_7.5|7.5] otrzymujemy:


Parser nie mógł rozpoznać (błąd składni): {\displaystyle \frac{1}{1-x}&=&\sum_{n=0}^{\infty}{n \choose n}x^n=\sum_{n=0}^{\infty}x^n,}      (8)

{{{3}}}      (9)


Po przekształceniu równości (9) uzyskuje się


n=0nxn=1(1x)211x.      (10)


Powołując się ponownie na Wniosek 7.5 otrzymujemy


1(1x)3=n=0(n+2n)xn=12n=0n2xn+32n=0nxn+n=0xn,


co w połączeniu z równościami (9) oraz (10) daje zwartą postać funkcji tworzącej Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} dla ciągu 1,4,9,,n2,:


Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}=\sum_{n=0}^{\infty}n^2x^n =\frac{2}{\left(1-x\right)^3}-\frac{3}{\left(1-x\right)^2}+\frac{1}{1-x}. }


Naszym zadaniem było jednakże policzenie funkcji tworzącej H(x) dla ciągu 1,1+4,1+4+9,,1+4+9++n2,, tzn. ciągu sum początkowych wyrazów ciągu 1,4,9,,n2,. Aby uzyskać Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{H}{x}} wystarczy więc skorzystać ze wzoru (7) i podzielić Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{G}{x}} przez 1x. Tak więc poszukiwanym rozwiązaniem są współczynniki funkcji tworzącej


Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{H}{x}=\frac{\fGen{G}{x}}{1-x} =\frac{2}{\left(1-x\right)^4}-\frac{3}{\left(1-x\right)^3}+\frac{1}{\left(1-x\right)^2}. }


Korzystając po raz kolejny z Wniosku 7.5 otrzymujemy


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned\fGen{H}{x} &=&2\sum_{n=0}^{\infty}{n+3 \choose n}x^n-3\sum_{n=0}^{\infty}{n+2 \choose n}x^n+\sum_{n=0}^{\infty}{n+1 \choose n}x^n\\ &=&\sum_{n=0}^{\infty}\left(\frac{1}{3}n^3+\frac{1}{2}n^2+\frac{1}{6}n\right)x^n. \endaligned}


W konsekwencji zachodzi równość


Parser nie mógł rozpoznać (nieznana funkcja „\vect”): {\displaystyle \sum_{k=1}^nk^2=\vect{x^n}\fGen{H}{x}=\frac{2n^3+3n+n}{6}. }


Przykład

Wracamy do przykładu z monetami. Występowały tam funkcje tworzące postaci


Parser nie mógł rozpoznać (nieznana funkcja „\fGen”): {\displaystyle \fGen{A_k}{x} = 1+x^k+x^{2k}+x^{3k}+\ldots, }


dla k=1,5,10,25 i 50. Z równości (7) wiemy, że


1+xk+x2k+x3k+,=11xk


tak więc:


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned\fGen{A}{x}= \fGen{A_1}{x}&=& \frac{1}{1-x},\\ \fGen{B}{x}= \fGen{A}{x}\cdot \fGen{A_5}{x} &=&\frac{\fGen{A}{x}}{1-x^5},\\ \fGen{C}{x}= \fGen{B}{x}\cdot \fGen{A_{10}}{x} &=&\frac{\fGen{B}{x}}{1-x^{10}},\\ \fGen{D}{x}= \fGen{C}{x}\cdot \fGen{A_{25}}{x} &=&\frac{\fGen{C}{x}}{1-x^{25}},\\ \fGen{E}{x}= \fGen{D}{x}\cdot \fGen{A_{50}}{x} &=&\frac{\fGen{D}{x}}{1-x^{50}}, \endaligned}


skąd natychmiast:


Parser nie mógł rozpoznać (nieznana funkcja „\aligned”): {\displaystyle \aligned\fGen{A}{x}&=&1+x\fGen{A}{x},\\ \fGen{B}{x}&=&\fGen{A}{x}+x^5\fGen{B}{x},\\ \fGen{C}{x}&=&\fGen{B}{x}+x^{10}\fGen{C}{x},\\ \fGen{C}{x}&=&\fGen{D}{x}+x^{25}\fGen{C}{x},\\ \fGen{D}{x}&=&\fGen{E}{x}+x^{50}\fGen{D}{x}. \endaligned}


Równości te dają zależności między współczynnikami:


an=1,bn=an+bn5,cn=bn+cn10,dn=cn+dn25,en=dn+en50.


Wykorzystując te zależności rekurencyjne możemy wypełnić następującą tabelę:


{

Funkcje tworzące w rozwiązywaniu zależności rekurencyjnych