LINUX.ORG.RU

Ответ на: комментарий от Selecter

C рулит :))
#include <gmp.h>
#include <stdarg.h>

mpz_t n;

int main()
{
int i;
mpz_init (n);
mpz_set_ui (n,1);
for(i=1; i<=1003; i++)
{
mpz_mul_si(n,n,i);
}
gmp_printf("%Zd\n",n);
mpz_clear(n);
}

Selecter ★★★★
() автор топика
Ответ на: комментарий от Selecter

Особенно прикольно факториал миллиона вычислять =))) гы

Selecter ★★★★
() автор топика

BTW, gmp - не самая лучшая библиотека на этот счет...

Есть очень хорошая библиотека, пИсаная на Фортране (автор - Бейли) http://crd.lbl.gov/~dhbailey/

Есть (глюкавые) врапперы к ней на ЦеПП

Касательно ЦеПП можно стартануть отсюда: http://www.ginac.de/

Die-Hard ★★★★★
()
Ответ на: комментарий от human0id

human0id (23.05.2004 0:49:35):

> А зачем таки понадобилась такая точность?

Часто надо бывает! Например: Криптография, аналитические вычисления, восстановление точности при эстраполяции, и т.п.

Например, алгоритм PSLQ (http://mathworld.wolfram.com/PSLQAlgorithm.html) позволяет иногда восстановить аналитическое выражение через трансцендентные функции, но для его применения требуется знать 50-100 знаков после запятой...

Die-Hard ★★★★★
()
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.