Язык Си в примерах/Треугольник Паскаля

Материал из Викиучебника

Перейти к: навигация, поиск
Язык Си в примерах

  1. Компиляция программ
  2. Простейшая программа «Hello World»
  3. Учимся складывать
  4. Максимум
  5. Таблица умножения
  6. ASCII коды символов
  7. Верхний регистр
  8. Скобочки
  9. Факториал
  10. Степень числа
  11. Треугольник Паскаля
  12. Корень уравнения
  13. Система счисления
  14. Сортировка
  15. Библиотека complex
  16. Сортировка на основе qsort
  17. RPN калькулятор
  18. RPN калькулятор на Bison
  19. Простая грамматика
  20. Задача «Расчёт сопротивления схемы»
  21. Простая реализация конечного автомата
  22. Использование аргументов командной строки
  23. Чтение и печать без использования stdio

Всем известны простые формулы

  • (a + b)2 = a2 + 2ab + b2
  • (a + b)3 = a3 + 3a2b + 3ab2 + b3

А как находить коэффициенты в разложении (a + b)n?

Рассмотрим следующую таблицу:

Pascal.png

То, что нарисовано справа, называется треугольником Паскаля — в n-й строчке этого треугольника находятся коэффициенты для разложения (a+b)^n\,\!.

Число номер k + 1 в n-й строчке называется биномиальным коэффициентом C_n^k\,\!.

Например, C_3^0 = 1, \quad C_5^2 = 10, \quad C_4^2 = 6.\,\!

Стороны треугольника Паскаля состоят из единичек. Каждое число внутри треугольника Паскаля равно сумме двух чисел, стоящих над ним справа и слева в предыдущей строчке:

C^{k+1}_{n+1} = C^{k+1}_{n}+C^k_{n}.

Эти числа возникают в задаче о числе сочетаний: C_n^k\,\! — это число способов выбрать k элементов из n различных. Например, чиcло байт, в которых ровно 3 единицы — это число C_8^3\,\! — число способов выбрать три бита, в которых будут стоять единицы, из воcьми бит байта.


Докажите, что
C_n^1=n,\quad C_n^k+C_n^{k+1} = C_{n+1}^{k+1}, \quad C_n^k = \frac{n!}{k! (n-k)!}.\,\!


Рассмотрим две программы, которые решают следующие задачи:

  1. Запрограммировать функцию C(n,k) = C_n^k\,\!.
  2. Вывести на экран n строчек треугольника Паскаля.


 /*
   Вычисление биномиальных коэффициентов.
 */
 #include <stdio.h>
 long C(long n,long k) {
    if(k == 0 || n == k) return 1;
    return C(n - 1, k - 1) + C(n - 1, k);
 }
 int main() {
    long n, k;
    scanf ("%ld%ld", &n, &k);
    printf ("%ld ", C(n, k));
    return 0;
 }
  • Сколько раз вызовется функция C(., .) при вычислении С(n, k)?
  • Докажите, что время вычисления C_n^k\,\! по приведенному алгоритму пропорционально C_n^k\,\!.
  • Оцените ассимптотику C_n^{n/2}\,\!, а именно, напишите программу, которая вычисляет \log C_n^{n/2}\,\! для n = 2,4,...,40 и нарисуте график

\log C_n^{n/2}\,\! от n.

 /*
     Вычисление <i>n</i>-й строки треугольника Паскаля.
 */
 #include <stdio.h>
 #define N 1000
 long c[N];
 int main () {
    long n, i, j;
    scanf ("%ld",&n);
    for(i = 1; i <= n ; i++) c[i] =0;
    c[0] = 1;
    for(j = 1 ; j <= n; j++)
       for(i = j; i >= 1 ; i--)
          c[i] = c[i-1] + c[i];
    for(i = 0; i <= n ; i++)
       printf ("%ld ", c[i]);
    return 0;
 }
  • Докажите, что указанный алгоритм вычисления n-й строчки треугольника Паскаля

работает быстрее, чем алгоритм вычисления C_n^k\,\! из предыдущей программы, а именно время работы пропорционально n^2\,\!.

  • Начиная с какого n самое большое число из n-й строчки треугольника Паскаля не умещается в тип long?