Формула Бернулли

@Алексей9912 · Регистрация: 18.11.2017

Студворк — интернет-сервис помощи студентам

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
var
k,n:integer;
p,c,ber,p1,q1:real;
k1,n1:longint;
function fac(x:integer):longint;
begin
if x=1 then fac:=1 else fac:=fac(x-1)*x;
end;
begin
Read(k,n,p);
k1:=fac(k);
n1:=fac(n);
p1:=Power(p,k);
q1:=Power(1-p,n-k);
c:=n1/fac(n-k)/k1;
ber:=c*p1*q1;
Writeln(ber);
end.

Ввод: 100 200 0.25 (или например 6 66 0.66)
Вывод: "не число"

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

@Joy · 15.03.2018, 18:37

100! = 9332621544394415268169923885626670049071 5968264381621468592963895217599993229915 6089414639761565182862536979208272237582 51185210916864000000000000000000000000

в longint не влезет

@Cyborg Drone · 15.03.2018, 22:36

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

$https://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> C_k^n=\frac{n!}{k!(n-k)!}=\frac{\prod_{i=k+1}^{n}i}{(n-k)!}<br />$

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

$https://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> s=\sum_{i=1}\lfloor m/p^i\rfloor<br />$

суммирование продолжается до тех пор, пока очередное слагаемое не станет равным 0.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11680&d=1772460536 Одним из. . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .
SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .	Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .

@Joy 2388 / 1300 / 1492 Регистрация: 29.08.2014 Сообщений: 4,665
	15.03.2018, 18:37
	100! = 9332621544394415268169923885626670049071 5968264381621468592963895217599993229915 6089414639761565182862536979208272237582 51185210916864000000000000000000000000 в longint не влезет 0

@Cyborg Drone Модератор 10434 / 5722 / 3405 Регистрация: 17.08.2012 Сообщений: 17,412
	15.03.2018, 22:36
	Чтобы не было переполнения, биномиальные коэффициенты считайте либо с помощью треугольника Паскаля, либо с помощью разложения числителя и знаменателя на простые множители. В последнем случае резонно сократить формулу биномиальных коэффициентов на k! $https://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> C_k^n=\frac{n!}{k!(n-k)!}=\frac{\prod_{i=k+1}^{n}i}{(n-k)!}<br />$ Степень каждого из простых чисел в разложении факториала в знаменателе можно выяснить с помощью факторизации факториала. Если число p - простое, и, если обозначить m=n-k то количество сомножителей p в факториале будет равно $https://www.cyberforum.ru/cgi-bin/latex.cgi?<br /> s=\sum_{i=1}\lfloor m/p^i\rfloor<br />$ суммирование продолжается до тех пор, пока очередное слагаемое не станет равным 0. 1