Умножение длинного числа на короткое

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

Разобьем этот алгоритм на две части:

1) Умножим каждую цифру длинного числа a на короткое число k и сохраним результаты в массиве c, не заботясь о переносах:

c[i] := a[i] * k;

2) Пройдем по полученному массиву и “сделаем переносы”:

for i:=1 to 1000 do begin //следующие строки важно выполнить именно в таком порядке

c[i+1]:=c[i+1]+(c[i] div 10);

c[i]:=c[i] mod 10;

end;

Умножение длинного числа на длинное

Алгоритм будет похоже на предыдущий, с той лишь разницей, что сначала мы умножим каждую цифру первого числа на каждую цифру второго числа:

for i:=1 to 1000 do

for j:=1 to 1000 do

c[…] := с[…] + a[i] * b[j];

(Именно здесь мы используем то, что объявлен массив длины 200 – это упростит код.)

Подумаем, что нужно поставить вместо многоточия.

При умножении 1-й цифры a на первую цифру b результат нужно положить в первую цифру c (разряд единиц).

При умножении 1-й цифры a на вторую цифру b результат нужно положить во вторую цифру c (в разряд десятков).

При умножении 2-й цифры a на вторую цифру b результат нужно положить в третью цифру с (разряд сотен).

Можно заметить, что при умножении i-й цифры a на j-ю цифру b результат нужно положить в (i+j-1)-ю цифру с, то есть на место каждого многоточия надо вписать i+j-1.

После этого мы точно так же, как и в предыдущем разделе, проходим по массиву с и делаем переносы.

Несколько слов о реализации

Рассмотрев (в самом простом варианте) основные алгоритмы, поговорим немного об их реализации. Ясно, что каждый алгоритм удобно оформить в виде соответствующей функции или процедуры. Для того, чтобы передавать им массив, объявим новый тип:

type

MyNumber = array [1..1 000] of integer;

После этого будем передавать аргументы наших арифметических операций как константы (они не будут меняться), а переменную, в которую будет попадать результат – с модификатором var:

Procedure Add (const a,b : MyNumber; var c : MyNumber);

Для удобства можно ввести константу – максимальную длину числа:

Const

Maxlen = 1000;

В этом случае нам не придется вносить исправления по всей программе, если мы захотим изменить это число.

Заключение

Мы рассмотрели самый простой (и потому самый неэффективный) способ работы с длинными числами. Тем не менее, во многих приложениях его вполне хватает. Если же требуется увеличить скорость работы, то можно двигаться в нескольких направлениях:

- хранить не по одной, а по несколько последовательных цифр числа в одном элементе массива;

- обрабатывать массив не до конца, а вместо этого хранить реальную длину каждого числа;

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

ЗАДАЧИ НА ДЛИННУЮ АРИФМЕТИКУ

Шестаков А.П.

http://comp-science.narod.ru/DL-AR/dl-ar.html

Рассмотрим достаточно популярную в программировании задачу на работу с "длинными" числами. Реально с "астрономическими" или "микроскопическими" числами приходится сталкиваться не так уж и часто. Тем не менее, упражнения, рассматриваемые в этой публикации, могут послужить хорошей тренировкой в области программирования и занять достойное место в классах с углубленным изучением информатики или на кружках по программированию. Алгоритмы, представленные ниже, записаны на Turbo Pascal, версия 7.0. При желании или необходимости они могут легко быть адаптированы к любой другой программной среде.

Диапазон представления целых чисел (Integer, Word, LongInt) ограничен, о чем не раз уже говорилось (впрочем, для действительных величин это замечание тоже актуально). Поэтому при решении задач всегда приходится действовать с оглядкой, — как бы не допустить возникновения ошибки выхода за диапазон или переполнения. Например, вычисляя факториал (n! = 1 * 2 * 3 * … * n), в диапазоне представления величин типа Integer удастся правильно получить только 7! = 5040, а в диапазоне представления типа LongInt — 12! = 479001600. Для больших значений, конечно, можно использовать действительные типы данных, но это уже не гарантирует точного результата. Поэтому полезно для получения точных значений при действиях с многозначными числами разработать другие способы представления таких чисел, алгоритмы выполнения арифметических и других операций, процедуры ввода и вывода результатов и т.д.

Покажем реализацию решения такого рода задач на примере умножения одного многозначного числа на другое. Именно эта арифметическая операция наиболее часто используется при решении других задач.

Наиболее естественным способом представления многозначного числа является запись каждого его разряда в виде отдельного элемента линейного массива (или списка, где память под цифру будет отводиться по мере надобности, в то время как в массиве приходится заранее задавать максимальное количество элементов в нем). Пусть (для удобства дальнейших действий) разряды "длинного" числа при записи в массив нумеруются с единицы, начиная с разряда единиц, т.е. цифра из разряда единиц — элемент массива с номером один, цифра из разряда десятков — элемент массива с номером два и т.д. Определим для работы с "длинными" неотрицательными числами тип данных:

Const MNax = 2000; Type Digit = 0..9; DlChislo = Array[1..Nmax] Of Digit;

Для решения поставленной задачи необходимо уметь выполнять следующие действия:

1) ввод "длинного" числа;

2) собственно умножение двух "длинных" чисел;

3) вывод "длинного" числа;

4) определение количества цифр в записи числа.

Каждую из подзадач реализуем в виде отдельной подпрограммы. Начнем с ввода. Ввести большое число целесообразно в виде строки, а в дальнейшем преобразовать в массив цифр. В процедуре учтен указанный выше способ размещения "длинного" числа в массиве, т.е. с точки зрения пользователя число записывается как бы в обратном порядке.

{Процедура преобразования длинного числа, записанного в виде строки, в массив цифр; переменная OK принимает значение True, если в записи числа нет посторонних символов, отличных от десятичных цифр, иначе — false} Procedure Translate(S : String; Var A : DlChislo; Var OK : Boolean); Var I : Word; Begin Zero(A); I := Length(S); OK := True; While (I >= 1) And OK Do Begin If S[I] In ['0'..'9'] Then A[Length(S)- I + 1]:= Ord(S[I]) - 48 Else OK := False; I := I - 1 End End;

В процедуре вызывается подпрограмма Zero(A), назначение которой — запись нуля в каждый разряд длинного числа. Вот текст этой процедуры:

{Процедура обнуления длинного числа} Procedure Zero(Var A : DlChislo); Var I : Integer; Begin For I := 1 To NMax Do A[I] := 0; End;

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

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

{Функция определения количества цифр в записи длинного числа} Function Dlina(C : DlChislo) : Integer; Var I : Integer; Begin I := NMax; While (I > 1) And (C[I] = 0) Do I := I - 1; Dlina := I End;

При ее разработке было использовано следующее соображение: если число не равно нулю, то количество цифр в его записи равно номеру первой цифры, отличной от нуля, если просмотр числа осуществляется от старшего разряда к младшему. Если же длинное число равно нулю, то получается, что количество цифр в его записи равно одной, что и требовалось.

Ну и, наконец, главная процедура, ради которой и была проделана вся предшествующая работа. При составлении алгоритма используется идея умножения "столбиком", хотя в нашем варианте сложение выполняется не по окончанию умножения, а по ходу его, т.е. перемножив очередные цифры, сразу же добавляем результирующую цифру в нужный разряд и формируем перенос в следующий разряд.

{Процедура умножения длинных чисел. A, B — множители, C — произведение} Procedure Multiplication(A, B : DlChislo; Var C : DlChislo); Var I, J : Integer; P : Digit; VspRez : 0..99; Begin Zero(C); For I := 1 To Dlina(A) Do {Цикл по количеству цифр в первом числе} Begin P := 0; {Первоначально перенос равен нулю} For J := 1 To Dlina(B) Do {Цикл по количеству цифр во втором числе} Begin VspRez := A[I] * B[J] + P + C[I + J - 1]; C[I + J - 1] := VspRez Mod 10; {Очередное значение цифры в разряде I + J - 1} P := VspRez Div 10 {Перенос в следующий разряд} End; C[I + J] := P {последний перенос может быть отличен от нуля, запишем его в пока ещё свободный разряд} End End;

Сейчас приведем листинг программы целиком.

Program DlUmn;Const NMax = 2000;Type Digit = 0..9; DlChislo = Array[1..Nmax] Of Digit;Var S : String; M, N, R, F : DlChislo; I, MaxF : Word; Logic : Boolean;{Процедура обнуления длинного числа}Procedure Zero(Var A : DlChislo);Var I : Integer; Begin For I := 1 To NMax Do A[I] := 0; End;{Функция определения количества цифр в записи длинного числа}Function Dlina(C : DlChislo) : Integer;Var I : Integer; Begin I := NMax; While (I > 1) And (C[I] = 0) Do I := I - 1; Dlina := I End;{Процедура печати длинного числа}Procedure Print(A : DlChislo);Var I : Integer; Begin For I := Dlina(A) DownTo 1 Do Write(A[I] : 1); WriteLn End;{Процедура преобразования длинного числа в массив цифр}Procedure Translate(S : String; Var A : DlChislo; Var OK : Boolean);Var I : Word; Begin Zero(A); I := Length(S); OK := True; While (I >= 1) And OK Do Begin If S[I] In ['0'..'9'] Then A[Length(S) - I+ 1] := Ord(S[I]) - 48 Else OK := False; I := I - 1 End End;Procedure Multiplication(A, B : DlChislo; Var C : DlChislo);Var I, J : Integer; P : Digit; VspRez : 0..99; Begin Zero(C); For I := 1 To Dlina(A) Do Begin P := 0; For J := 1 To Dlina(B) Do Begin VspRez := A[I] * B[J] + P + C[I + J - 1]; C[I + J - 1] := VspRez Mod 10; P := VspRez Div 10 End; C[I + J] := P End End;{Основная программа}Begin Repeat {повторяем ввод, пока число не будет введено правильно} Write('Введите первый множитель: '); ReadLn(S); Translate(S, M, Logic) Until Logic; Repeat Write('Введите второй множитель: '); ReadLn(S); Translate(S, N, Logic) Until Logic; Multiplication(M, N, R); Print(R)End.

В приведенном листинге Print — процедура вывода длинного числа. Предоставим читателю самостоятельно разобраться в алгоритме ее работы.

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

Вот измененный фрагмент основной программы, решающий поставленную задачу.

Begin MaxF := 810; Zero(F); F[1] := 1; For I := 1 To MaxF Do Begin Str(I, S); {преобразование числа I к строковому типу S} Translate(S, M, Logic); Multiplication(F, M, F); Print(F); WriteLn('Факториал числа ', I : 4, ' содержит ', Dlina(F), ' цифр.') End End.

Расчеты показали, что можно вычислять факториалы до значения 810! включительно, в записи которого 1999 цифр. Далее вновь возникает переполнение.

Ниже будет предложен список задач для самостоятельного выполнения. Из них, по мнению автора, наибольшую сложность представляют реализации алгоритмов деления одного длинного числа на другое и извлечение квадратного корня. Алгоритм извлечения квадратного корня подробно описан в справочнике В.А. Гусева и А.Г. Мордковича [7]. В некоторых случаях составленные программы могут выступать как подпрограммы при разработке алгоритмов решения других, более сложных (как в примере с факториалом), задач. Кроме авторских задач и задач из списка литературы здесь приведены задания из олимпиад школьников по программированию, проводившихся в Пермской области в 1989-99 гг.

Задачи для самостоятельного решения

1. Составить программу сравнения двух многозначных чисел (количество знаков в записи чисел более 20).
2. Составить программу, суммирующую два натуральных многозначных числа с количеством знаков более 20.
3. Составить программу вычисления степени aⁿ, если a > MaxInt, n > 10.
4. Составить программу вычисления числа 2⁶⁴ – 1, в результате сохранить все цифры.
5. Составить программу вычисления 100!.
6. Составить программу извлечения точного квадратного корня из n-разрядного числа (n > 40).
7. Составить программу вычисления точного значения n!, где n > 12.
8. Составить программу вычисления точного значения nⁿ, где n > 10.
9. Составить программу деления числа a на число b, если a, b — многозначные числа.
10. Вычислить 100! + 2¹⁰⁰.
11. Вычислить 100! – 2¹⁰⁰.
12. Вычислить 7¹²³.
13. Встречаются ли среди цифр числа 2¹¹²¹³ – 1 две подряд идущие девятки?
14. Вычислить 2^–200.
15. Составить программу нахождения частного и остатка от деления m-значного числа на n–значное (m, n > 20).
16. Выяснить, какое из чисел a^m, bⁿ больше и на сколько (a, b <= 40000; m, n <= 10).
17. Найти n знаков в десятичной записи квадратного корня из целого числа m (n >= 50).
18. Найти количество делителей n-значного натурального числа (n > 20).
19. Вычислить точное значение (n!)! (n >= 3).
20. Составить программу вычисления точного значения суммы 1! + 2! + 3! + ... + n! при n > 10.
21. Составить программу вычисления точного значения суммы дробей

при n > 10. Ответ должен быть представлен в виде несократимой дроби P / Q, где P, Q — натуральные числа.

1. Вычислить точное значение (nⁿ)! при n >= 3.
2. Составить программу вычисления точного значения суммы первых n членов последовательности 1, k, k², k³, ..., kⁿ (n > MaxInt). Указание: используйте формулу суммы n членов геометрической прогрессии.
3. Составить программу вычисления точного значения суммы первых n членов последовательности чисел, кратных данному натуральному числу k (n > MaxInt). Указание: используйте формулу суммы n членов арифметической прогрессии.
4. Вычислить точное значение суммы 1² + 2² + 3² + … + n² (n >= 20000).
5. Вычислить точное значение суммы 1ⁿ + 2ⁿ + 3ⁿ + ... + nⁿ (n >= 10).
6. Найти первое простое число, которое больше 10¹¹.
7. Составить программу вычисления точного значения многочлена a_nxⁿ + a_n - 1xⁿ - 1 + ... + a₁x + a₀, где a_i и x — целые числа больше 10¹¹.
8. Найти наибольший общий делитель и наименьшее общее кратное чисел m и n (m, n >= 10¹¹).
9. Проверить, являются ли числа m и n (m, n >= 10¹¹) взаимно простыми.
10. Докажите, что число 2¹⁹⁹³⁶ * (2¹⁹⁹³⁷ – 1) является совершенным, т.е. равно сумме всех своих делителей, кроме самого себя.
11. "Вращающееся число". Написать программу, которая находит число, обладающее следующими свойствами:
    1) число оканчивается на 5;
    2) при умножении его на 5 образуется новое число, которое может быть получено из исходного вычеркиванием цифры 5 на конце и переписыванием ее в начало числа.
12. Дана последовательность, заданная рекуррентной формулой

a_{n + 1} = 7a_n mod 2023, a₁ = 1,

где x mod y означает остаток от деления x на y. Написать программу, вычисляющую a_n при 1 <= n <= 1000000000000000000000.

1. Дана последовательность

Написать программу, находящую точное значение a_n при 1 <= n <= 150.

Пример. При n = 58 получаем a_n = 10359022039470231387111424.

1. Напишите программу перевода многозначного числа (с количеством знаков больше 20) в системы счисления с основанием два, восемь, шестнадцать.
2. Разложить на простые множители натуральное число с количеством знаков более 11.
3. Умножение периодической дроби. Задана некоторая положительная правильная периодическая дробь Q и натуральное число N. Числа Q и N таковы, что количество цифр, используемых для их описания, не превосходит 100. При изображении дроби Q периодическая часть заключается в круглые скобки.

Требуется написать программу, которая определяет результат умножения Q на N, то есть непериодическую часть и минимальный период числа Q * N.

В случае получения результата умножения в виде конечной дроби скобки опускаются.

Пример работы правильной программы
Введите периодическую дробь: 0.1(6)
Введите натуральное число: 2
Ответ: 0.(3)

Список использованной литературы

1. Абрамов С.А., Гнездилова Г.Г., Капустина Е.Н., Селюн М.И. Задачи по программированию. М.: Наука, 1988.
2. Олимпиады по информатике. Задачи и решения. Методические рекомендации для учителей и учащихся школ. Красноярск, 1991.
3. Пильщиков В.Н. Сборник упражнений по языку Паскаль. М.: Наука, 1989.
4. Касаткин В.Н. Информация. Алгоритмы. ЭВМ. М.: Просвещение, 1991.
5. Хонсбергер Р. Математические изюминки. М.: Наука, 1992.
6. Семакин И.Г., Шестаков А.П. Лекции по программированию. Пермь: изд-во ПГУ, 1998.
7. Гусев В.А., Мордкович А.Г. Математика. Справочные материалы. М.: Просвещение, 1990.
8. Гладков В.П. Курс лабораторных работ по программированию. Пермь: изд-во ПГТУ, 1998.

Проверочная работа

1 вариант

1. Представить число 5012 в виде суммы наименьшего количества степеней двойки. 5012 = 2¹² + 2⁹ +2⁸ +2⁷ +2⁴+2²
2. Сформулируйте алгоритм перевода десятичной дроби в двоичную дробь.
3. Число 1394₁₆ записать в восьмеричной системе счисления. 11624₈
4. Байт 11001100 представляет число, записанное в обратном дополнительном коде. Чему равно это число в десятичной системе счисления? -52
5. Для кодирования букв А, Б, В, Г решили использовать двухразрядные последовательные двоичные числа (от 00 до 11, соответственно). Запишите шестнадцатеричным кодом последовательность ГВАБА. 384₁₆
6. Игрок в карты шифрует свои выигрыши/проигрыши в системе счисления с цифрами {♠, ♣, ♥, ♦}. Шифр простой: выигрыш кодируется ♥♥, проигрыш ♠♠. Как следует декодировать фразу «♦♣.♦.♥♣ ♠♠ ♣♦♦♣♠». 13.03.09 проиграл 500

2 вариант

1. Число 0,525 записать в виде двоичной периодической дроби. 0.100(0011)
2. Сумму чисел 43₈ + 56₁₆ представьте в двоичной системе счисления. 1111001₂
3. Сформулируйте алгоритм перевода восьмеричного целого числа в двоичное.
4. Запишите 8-ми разрядное представление числа -52 в обратном дополнительном коде. 11001100
5. Для кодирования букв А, Б, В, Г решили использовать двухразрядные последовательные двоичные числа (от 00 до 11, соответственно). Какое слово закодировано шестнадцатеричным кодом 384₁₆? (последовательность ГВАБА. )
6. Игрок в карты шифрует свои выигрыши/проигрыши в системе счисления с цифрами {♠, ♣, ♥, ♦}. Шифр простой: выигрыш кодируется ♥♥, проигрыш ♠♠. Как будет кодироваться фраза «13.03.09 выиграл 500». ♦♣.♦.♥♣ ♥♥ ♣♦♦♣♠

Дополнительные задания

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

a. Из 16-й системы счисления 2-ю. Число А8₁₆.

b. Из 2-й системы счисления в 8-ю. Число 1001001₂.

Даны два числа в десятичной системе счисления. Перевести каждое число в двоичную систему счисления и вычислить их сумму, разность, произведение, частное и остаток в двоичной системе счисления.

a. 13 и 2.

b. 15 и 4.

a. Число, записанное в десятичной системе счисления, оканчивается цифрой 5. Будет ли оно делиться на 5, если записать его в троичной системе счисления? Привести объяснения.

b. В некоторой позиционной системе счисления имеются цифры {J, K, L}. Чему равно число KKJLJK в десятичной системе счисления?

a. В каких системах счисления число 10 – нечетное. Привести объяснения и пример.

Найти основание p позиционной системы счисления и цифру n, если верно равенство 33m5n + 2n443 = 55424. Здесь m – максимальная цифра в этой систем