Выводит неточный результат

@c++\noob · Регистрация: 13.11.2010

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

Где-то есть ошибка. Выводит неточный результат. Ввожу : 1 0.001 Выводит: 2.71806 , а надо 2.71825

Исправьте пожалуйста.

C++

#include <stdio.h>
#include <math.h>
int factorial(int a)// funcion factorial
{
int g=1;
for(int i=1; i<=a; i++)
{
g=g*i;
}
return(g);
}
 
int main()
{
FILE *f;
double x,eps;
f = fopen("input.txt","r");
fscanf(f,"%lf %lf", &x, &eps);
fclose(f);
double b=1.0;
 
int p;
for(p=1; pow(x,p)/factorial(p)>eps; p++)
{
b=b+pow(x,p)/factorial(p);
}
 
f = fopen("output.txt","w");
fprintf(f, "%.5lf", b);
fclose(f);
return 0;
}

valeriikozlov · 20.12.2010, 12:09

c++\noob, Все правильно выводит Ваша программа.
Вот если ввести: 1 0.0001, тогда выдаст тоже правильный результат: 2.71825

@c++\noob · 20.12.2010, 13:59 **[ТС]**

ну тогда извините, недосмотрел =)

Добавлено через 14 минут
нет именно при вводе 1 0.001 должно быть 2.71825. Отправлял задачу на econtest - не приняли((

valeriikozlov · 20.12.2010, 14:03

Условие задачи напишите.

@c++\noob · 20.12.2010, 14:10 **[ТС]**

valeriikozlov, нет именно при вводе 1 0.001 должно быть 2.71825. Отправлял задачу на econtest - не приняли((

Добавлено через 22 секунды
щас секунду

Добавлено через 33 секунды
извините за второе сообщение нечаянно отправил)

Добавлено через 3 минуты
Не используя стандартные функции (за исключением ABS), вычислить с точностью Eps:

y=e^x=1 + x/1! + x^2/2! + ... + x^n/n! + ...

Считать, что требуемая точность достигнута, если очередное слагаемое по модулю меньше Eps, все последующие слагаемые можно уже не учитывать

valeriikozlov · 20.12.2010, 14:22

Сообщение от c++\noob

Считать, что требуемая точность достигнута, если очередное слагаемое по модулю меньше Eps, все последующие слагаемые можно уже не учитывать

Вот ключевая фраза. Если очередное слагаемое по модулю меньше Eps (его нужно складывать с b),
а вот все последующие слагаемые можно уже не учитывать.

@c++\noob · 20.12.2010, 14:36 **[ТС]**

напишите ,пожалуйста, как надо

@kazak · 20.12.2010, 14:49

Сообщение от c++\noob

for(p=1; pow(x,p)/factorial(p)>eps; p++)

Тут лучше делать с do{}while

@c++\noob · 20.12.2010, 15:24 **[ТС]**

Переписал через do{}while - результат всё равно тот же. =(

valeriikozlov · 20.12.2010, 16:19

C

#include <stdio.h>
#include <math.h>
int factorial(int a)// funcion factorial
{
int g=1;
for(int i=1; i<=a; i++)
{
g=g*i;
}
return(g);
}
 
int main()
{
FILE *f;
double x,eps;
f = fopen("input.txt","r");
fscanf(f,"%lf %lf", &x, &eps);
fclose(f);
double b=1.0;
 
int p, fl=1;
for(p=1; fl==1; p++)
{
b=b+pow(x,p)/factorial(p);
if(pow(x,p)/factorial(p)<eps)
    fl=0;
}
 
f = fopen("output.txt","w");
fprintf(f, "%.5lf", b);
fclose(f);
return 0;
}

@c++\noob · 20.12.2010, 19:23 **[ТС]**

Извините пожалуйста. Я наверно вас уже достал, но опять где-то есть недочёт (econtest не принимает). Ввожу 3.14 0.000001, выбивает 22.19869, а надо 23.10387.

Добавлено через 6 минут

Сообщение от valeriikozlov

C

#include <stdio.h>
#include <math.h>
int factorial(int a)// funcion factorial
{
int g=1;
for(int i=1; i<=a; i++)
{
g=g*i;
}
return(g);
}
 
int main()
{
FILE *f;
double x,eps;
f = fopen("input.txt","r");
fscanf(f,"%lf %lf", &x, &eps);
fclose(f);
double b=1.0;
 
int p, fl=1;
for(p=1; fl==1; p++)
{
b=b+pow(x,p)/factorial(p);
if(pow(x,p)/factorial(p)<eps)
    fl=0;
}
 
f = fopen("output.txt","w");
fprintf(f, "%.5lf", b);
fclose(f);
return 0;
}

valeriikozlov · 20.12.2010, 19:47

Для данного набора входных данных: 3.14 0.000001 типа int для вычисления факториала числа не хватает (функция factorial() возвращает отрицательное значение), пробуйте так:

C

#include <stdio.h>
#include <math.h>
double factorial(double a)// funcion factorial
{
double g=1;
for(int i=1; i<=a; i++)
{
g=g*i;
}
return(g);
}
 
int main()
{
FILE *f;
double x,eps;
f = fopen("input.txt","r");
fscanf(f,"%lf %lf", &x, &eps);
fclose(f);
double b=1.0;
 
double p;
int fl=1;
for(p=1; fl==1; p++)
{
b=b+pow(x,p)/factorial(p);
if(pow(x,p)/factorial(p)<eps)
        fl=0;
}
 
f = fopen("output.txt","w");
fprintf(f, "%.5lf", b);
fclose(f);
return 0;
}

@Напильнег · 20.12.2010, 20:07

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

Добавлено через 5 минут
В данном случае рекурентное соотношение такое:

$https://www.cyberforum.ru/cgi-bin/latex.cgi?{y}_{0}=1$ , $https://www.cyberforum.ru/cgi-bin/latex.cgi?{y}_{i}={y}_{i-1}*x/i$

И никаких сусликов!

@c++\noob · 20.12.2010, 21:22 **[ТС]**

Сообщение от valeriikozlov

Для данного набора входных данных: 3.14 0.000001 типа int для вычисления факториала числа не хватает (функция factorial() возвращает отрицательное значение), пробуйте так:

C

#include <stdio.h>
#include <math.h>
double factorial(double a)// funcion factorial
{
double g=1;
for(int i=1; i<=a; i++)
{
g=g*i;
}
return(g);
}
 
int main()
{
FILE *f;
double x,eps;
f = fopen("input.txt","r");
fscanf(f,"%lf %lf", &x, &eps);
fclose(f);
double b=1.0;
 
double p;
int fl=1;
for(p=1; fl==1; p++)
{
b=b+pow(x,p)/factorial(p);
if(pow(x,p)/factorial(p)<eps)
        fl=0;
}
 
f = fopen("output.txt","w");
fprintf(f, "%.5lf", b);
fclose(f);
return 0;
}

econtest опять выдаёт ошибку. ввод 2 300 , вывод 3.00000, надо 1.00000 . ((( хз что сней делать

valeriikozlov · 21.12.2010, 07:28

c++\noob,

Сообщение от c++\noob

ввод 2 300 , вывод 3.00000, надо 1.00000

Вы не ошибаетесь?
x=2, eps=300
b=1+2/1=3, цикл останавливается после вычисления первого слагаемого 2/1 (т.к. оно меньше 300).
Кстати Вам хороший алгоритм написал Напильнег (но даже в этом случае его алгоритм выдаст тоже 3.00000).
А ссылку на задачу можете написать?

@Напильнег · 21.12.2010, 11:23

Сообщение от valeriikozlov

x=2, eps=300
b=1+2/1=3, цикл останавливается после вычисления первого слагаемого 2/1 (т.к. оно меньше 300).

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

Сообщение от valeriikozlov

Кстати Вам хороший алгоритм написал Напильнег (но даже в этом случае его алгоритм выдаст тоже 3.00000).

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

Например, так:

C++

  y=0; 
  i=0; yi=1;
  while (1) {
    y=y+yi; 
    if (yi<eps) break;
    i++; yi=yi*x/i;
  }

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

valeriikozlov · 21.12.2010, 11:46

Да я не учел что первое слагаемое 1 тоже учитывается:

C

#include <stdio.h>
 
int main()
{
FILE *f;
double x,eps;
f = fopen("input.txt","r");
fscanf(f,"%lf %lf", &x, &eps);
fclose(f);
double b=1.0, y=0.0;
double p;
int fl=1;
for(p=1; fl==1; p++)
{
y=y+b;
if(b<eps)
    fl=0;
b=b*x/p;
}
 
f = fopen("output.txt","w");
fprintf(f, "%.5lf", y);
fclose(f);
return 0;
}

@c++\noob · 22.12.2010, 10:42 **[ТС]**

valeriikozlov, спасибоо вам огромное всё заработало!!!

valeriikozlov 4728 / 2549 / 757 Регистрация: 18.08.2009 Сообщений: 4,568
	20.12.2010, 12:09	2
	c++\noob, Все правильно выводит Ваша программа. Вот если ввести: 1 0.0001, тогда выдаст тоже правильный результат: 2.71825 1

@c++\noob -2 / 2 / 1 Регистрация: 13.11.2010 Сообщений: 52
	20.12.2010, 13:59 [ТС]	3
	ну тогда извините, недосмотрел =) Добавлено через 14 минут нет именно при вводе 1 0.001 должно быть 2.71825. Отправлял задачу на econtest - не приняли(( 0

valeriikozlov 4728 / 2549 / 757 Регистрация: 18.08.2009 Сообщений: 4,568
	20.12.2010, 14:03	4
	Условие задачи напишите. 1

@c++\noob -2 / 2 / 1 Регистрация: 13.11.2010 Сообщений: 52
	20.12.2010, 14:10 [ТС]	5
	valeriikozlov, нет именно при вводе 1 0.001 должно быть 2.71825. Отправлял задачу на econtest - не приняли(( Добавлено через 22 секунды щас секунду Добавлено через 33 секунды извините за второе сообщение нечаянно отправил) Добавлено через 3 минуты Не используя стандартные функции (за исключением ABS), вычислить с точностью Eps: y=e^x=1 + x/1! + x^2/2! + ... + x^n/n! + ... Считать, что требуемая точность достигнута, если очередное слагаемое по модулю меньше Eps, все последующие слагаемые можно уже не учитывать 0

valeriikozlov 4728 / 2549 / 757 Регистрация: 18.08.2009 Сообщений: 4,568
	20.12.2010, 14:22	6
	Сообщение от c++\noob Считать, что требуемая точность достигнута, если очередное слагаемое по модулю меньше Eps, все последующие слагаемые можно уже не учитывать Вот ключевая фраза. Если очередное слагаемое по модулю меньше Eps (его нужно складывать с b), а вот все последующие слагаемые можно уже не учитывать. 1

@c++\noob -2 / 2 / 1 Регистрация: 13.11.2010 Сообщений: 52
	20.12.2010, 14:36 [ТС]	7
	напишите ,пожалуйста, как надо 0

@kazak 3565 / 2712 / 347 Регистрация: 11.03.2009 Сообщений: 6,246
	20.12.2010, 14:49	8
	Сообщение от c++\noob for(p=1; pow(x,p)/factorial(p)>eps; p++) Тут лучше делать с do{}while 0

@c++\noob -2 / 2 / 1 Регистрация: 13.11.2010 Сообщений: 52
	20.12.2010, 15:24 [ТС]	9
	Переписал через do{}while - результат всё равно тот же. =( 0

@Напильнег 481 / 119 / 17 Регистрация: 30.09.2010 Сообщений: 473
	20.12.2010, 20:07	13
	Double тоже надолго не хватит. Правильное решение таких задач - не вычислять факториал отдельно вообще, и pow по условию задачи тоже нельзя использовать кстати, а вычислять следующий член через предыдущий по рекурентной формуле. Добавлено через 5 минут В данном случае рекурентное соотношение такое: $https://www.cyberforum.ru/cgi-bin/latex.cgi?{y}_{0}=1$ , $https://www.cyberforum.ru/cgi-bin/latex.cgi?{y}_{i}={y}_{i-1}*x/i$ И никаких сусликов! 1

valeriikozlov 4728 / 2549 / 757 Регистрация: 18.08.2009 Сообщений: 4,568
	21.12.2010, 07:28	15
	c++\noob, Сообщение от c++\noob ввод 2 300 , вывод 3.00000, надо 1.00000 Вы не ошибаетесь? x=2, eps=300 b=1+2/1=3, цикл останавливается после вычисления первого слагаемого 2/1 (т.к. оно меньше 300). Кстати Вам хороший алгоритм написал Напильнег (но даже в этом случае его алгоритм выдаст тоже 3.00000). А ссылку на задачу можете написать? 0

@c++\noob -2 / 2 / 1 Регистрация: 13.11.2010 Сообщений: 52
	22.12.2010, 10:42 [ТС]	18
	valeriikozlov, спасибоо вам огромное всё заработало!!! 0