C++ - Сумма бесконечного ряда с факториалом

[Angre]

Найти сумму ряда с заданной точностью e = 0,001, общий член которого:
a с индексем n= n! / pow(n,n)
Проблема в том, что вычислить надо не с помощью функции, а с помощью for, while или do-while.

#include <iostream.h>
#include <math.h>
int main()
{
int n, f= 1, i= 1;
double a;
do
{
f= f * i;
i++;
do
{
a= a + f/ pow(n,n);
n++;}
while(a > 0.001);}
while( i <= n );
cout« a;

return 0;
}

И такой:
#include <iostream.h>
#include <math.h>
int main()
{
int n= 1, i, f;
double a= 1.00;
do
{for(i=1; i <= n; i++ )
{n++;
f= f * i;}
a= a + f/pow(n,n);
n++;
}
while( a > 0.001 );
cout« a;
return 0;
}

[min@y™]

и? где проблема?

[challengerr]

В первом нужно инициализировать a перед началом суммирования и инициализировать n.Заданная точность это разница между двумя последовательно вычисленными значениями a, а не само а, поэтому сравнение а с 0.001 некорректно. В классе cout нет переопределенного оператора «, а есть оператор <<. Факториал растет очень быстро, в общем случае нужна длинная арифметика.

[8Observer8]

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

#include <iostream>
#include <stdexcept>

unsigned int factorial( unsigned int n ) throw (std::out_of_range)
{
    static const unsigned int table[] = { 1, 1, 2, 6, 24, 120, 720,
                                          5040, 40320, 362880, 3628800,
                                          39916800, 479001600 };
    if ( n >= sizeof table / sizeof *table )
        throw std::out_of_range( "Range error: the number must be from the range [0, 12]" );
    return table[n];
}

int main( )
{
    int f = 0;
    int n = 12;

    try {
        f = factorial( n );
        std::cout << f << std::endl;
    } catch ( const std::out_of_range& e ) {
        std::cout << e.what( ) << std::endl;
    }

    return 0;
}

[evg_m]

По поводу факториала могу подсказать.

И нарваться на переполнение! (вот такой он факториал вредный больше 10! вряд ли кто посчитает без ухищрений (длинная арифметика).

n:=????
s:=0;
repeat
r:=1;  for k:=1 to n do r:=r*k/n;  /// r:=(1/n)*(2/n)*...*(n/n);
s:=s+r;
until r<0.001;

Это конечно же не С.

[8Observer8]

И нарваться на переполнение!

В коде из моего сообщения нет переполнения. Там факториал уже подсчитан. Функция factorial() выдаёт значение из таблицы. Если пользователь введёт 13, то ему программа покажет сообщение о некорректных данных. Такой способ быстрее, чем расчёт

[evg_m]

В коде из моего сообщения нет переполнения.

Если пользователь введёт 13, то ему программа покажет сообщение о некорректных данных.

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

Такой способ быстрее, чем расчёт

Ощутимая разница будет заметна при миллионах повторных расчетов 13!
В нашем случае удастся провести только тринадцать.(1!... 13!)
А дальше "...некорректные данные".

Мой код позволит рассчитать и за пределами 13! БЕЗ каких либо выкрутасов (вроде длинной арифметики).