Arch Linux

[C] - prüfen ob variable kubikzahl ist

SiD

Hi

ich möchte gerne prüfen ob eine zahl x ein kubikzahl ist. Zum Beispiel für x=56623104

x^(1/3)=384

ich habe bis jetzt das hier

#include <stdio.h>
#include <math.h>

#define EIN_DRITTEL 0.33333333333333333

int main(void)
{
    double a=56623104;

    printf("%.17lf\n",pow(a,EIN_DRITTEL));
    return 0;
}

als Ergebnis erhalte ich 383.99999999999988631, was ja an der internen verarbeitung von double liegt.
Wie kann ich jetzt z.B bis zur (sagen wir mal 10ten Stelle) prüfen und wenn alles neunen sind sagen, das ergebnis ist 384, bzw. eine ganzzahl und somit a eine Kubikzahl?

Kinch

So spontan: Du rundest dein Ergebnis auf eine Ganzzahl auf oder ab. Dann testest du, ob die Differenz zwischen deiner Rundung und deinem Ergebnis unterhalb eines gewissen Schwellenwertes ist. Also in Pseudocode:

kubik = x ^ (1/3);
if(abs(round(kubik) -kubik) < 0.0000000001)
  then Kubikzahl
  else keine Kubikzahl

for_mat_C

hallo!

folgendes sollte funktionieren:

#include <stdio.h>
#include <math.h>

int main()
{
double zahl1, zahl2 = 56623104;
zahl1 =(double) 1/3;
printf ("Ergebnis: %f\n", pow(zahl2, zahl1));

return 0;
}

dann bekommst du auch "Ergebnis: 384.000000"

aber nicht den -lm schalter für gcc vergessen!

gruß
for_mat_C

piet

Hi,
ich kann dein Problem nicht nachvollziehen:

int main ( int argc, char* argv[] ){
        double base  = 27.0;
        double exp   = 1.0/3;
        double result= 0.0;

        printf("%f^%f=%f\n",base,exp,result);

        if( (argc>2) ) {
                base=atof(argv[1]);
                exp =atof(argv[2]);
        }

        result=pow(base,exp);

        printf("%f^%f=%f\n",base,exp,result);

        return 0;
}

ergibt bei mir:

./testPowDouble 56623104 0.33333333333
27.000000^0.333333=3.000000
56623104.000000^0.333333=384.000000

@ Alle
Wenn ihr Variablen hartcodiert mit Inhalten füllt, ist die Chance sehr hoch dass der Compiler optimiert und das Ergebnis nicht wirklich repräsentativ ist.

@for_mat_C
Vermutlich compiliert dein Quellcode auch ohne -lm, oder?
Dein zahl1 =(double) 1/3; funktioniert bei mir so auch nicht, ein zahl1 =1.0/3; sollte helfen. (Mach mal printf-Ausgaben, dann merkst du es...)

@SiD
Gewöhne dir bitte an, gleich vernünftige Namen für Variablen, Funktionen usw. zu vergeben.

Nicht vergessen double sind bis 16 Stellen genau, ab dann wird gelogen dass es nur so kracht.
SiD könnte auch die Zahlen nach der Berechnung in integer umwandeln und die Rechnung rückwärts durchführen und so eine Kontrolle realisieren. Natürlich mit einer eigenen pow-Funktion: int sidsPow (int base, int exp)

cu

for_mat_C

naja, es wurde ja auch nicht die problemstellung mit kommandozeilenübergabe geschildert
😉 aber du hast recht, deine zeilen/proc haben mehr stil

nein, ohne einbinden von libm bekomme ich immer folgende fehlermeldung

undefined reference to `pow'

auch mit der folgenden variante:

#include <stdio.h>
#include <math.h>

int main()
{
double zahl1 = 1.0/3;
double zahl2 = 56623104;
printf ("Ergebnis: %f\n", pow(zahl2, zahl1));
return 0;
}

vielleicht liegt das an der älteren version von gcc . hab hier gerade nen debian mit "gcc version 4.3.2" zur hand.

gruß
for_mat_C