Dlaczego komenda pow podaje zły wynik przy podstawie=5?

Question

Witam serdecznie 

Tym razem mam problem z komendą potęgowania - pow. Oto kod:

#include <iostream>
#include <math.h>

using namespace std;

int a, b;
long long wynik;

int main()
{
        cin>>a>>b;
        wynik=pow(a,b);
        cout<<wynik<<endl;

    return 0;
}

Dlaczego, gdy w takim programie wpiszę a=5, i na przykład b=3 to otrzymuję wynik 124? Z tego co zauważyłem, to błędny wynik otrzymuję tylko, gdy a=5. Dla 5^2 dostaję 24, dla 5^7 ->78124.

Dla 5^10 dostaję poprawny wynik, a dla 5^12 znów zły.

 

Jeśli wpiszę:

cout<<pow(5,3)<<endl;

to problemu nie ma, wynik jest poprawny. 

Answer 1

Bardzo znany i prosty do wyjaśnienia błąd:

Każdy egzemplarz funkcji pow() zwraca typ zmiennoprzecinkowy (http://en.cppreference.com/w/cpp/numeric/math/pow). Ta funkcja ma pewną ułomność - nie jest w 100% dokładna. Dla np. pow(5, 3) zobaczysz wynik 124, gdyż wypisałeś tą liczbę jako całkowitą. Funkcja pow zwróciła w rzeczywistości liczbę 124.(9). Jeżeli użyjesz typu double to wszystko będzie działać:

#include <iostream>
#include <cmath> // cmath a nie math.h !

using namespace std;

int a, b;
double wynik;

int main()
{
        cin>>a>>b;
        wynik=pow(a,b);
        cout<<wynik<<endl;
}

Comments

No rok bym siedział i bym tego nie wymyślił. Teraz działa. Dzięki!

---

Sam myślałem nad tym jakiś rok temu. Wykminiłem to po dwóch dniach :D

---

Jak tak niedokładna funkcja mogła się znaleźć w pod. bibliotece?   jest jakiś powód dlaczego ona zwraca w ogóle double ?

---

Tego pana się pytaj, nie mnie :v http://bit.ly/1MtQOu9

---

Prosta sprawa zawsze można napisać coś takiego

template < class klasa >
klasa potegowanie(klasa x, klasa y){
	klasa wynik = x;
	while (y > 1){
		wynik *= x;
		y--;
	}
	return wynik;
}

Ale taki szablon się nie sprawdzi jeśli y będzie liczbą niecałkowitą sprawdz taki kod:

#include<iostream>

using namespace std;

template < class klasa >
klasa potegowanie(klasa x, klasa y){
	klasa wynik = x;
	while (y > 1){
		wynik *= x;
		y--;
	}
	return wynik;
}

int main(){
	double x = 2;
	double y = 2.2;
	cout << "Potegowanie przy pomocy szablonu: " << potegowanie < double >(x,y) << "\nPotegowanie funkcja pow: " << pow(x, y);
	cin >> x;

}

 

---

Zawsze można użyć asercji do sprawdzenia typu:

template <typename Type>
Type pow(Type base, Type index){

    static_assert(std::is_integral<Type>::value, "Typ musi być całkowity");
    return (index) ? (pow(base, index - 1) * base) : (1);
}

Oraz funkcji rekurencyjnej ;)

Answer 2

Prawdopodobnie dlatego, że wynik zapisujesz do zmiennej typu long long. Spróbuj tak zadeklarować tę zmienną:

double wynik;