C言語で3乗根を計算する方法：cbrt関数の使い方とpow関数との違いを解説

2026年5月17日

数学・アルゴリズム

C言語を用いて数値計算プログラムを作成する際、ある数値の3乗根（立方根）を求めたい場面は多々あります。

幾何学的な計算や物理シミュレーション、あるいは統計処理など、3乗根の計算は幅広い分野で必要とされる基本的な操作の一つです。

C言語の標準ライブラリには、これらの計算を効率的に行うための関数が用意されていますが、選択する手法によって精度や挙動が異なる場合があります。

この記事では、主にmath.hヘッダで提供されているcbrt関数を中心に、pow関数との違いや具体的な実装方法について詳しく解説します。

目次 [ close ]

math.hの活用と3乗根計算の基本

C言語で数学的な計算を行うためには、まず標準ライブラリであるmath.hをインクルードする必要があります。

このヘッダファイルには、三角関数や対数関数、累乗計算などの高度な数学関数が定義されています。

3乗根の計算においても、自力で複雑なアルゴリズムを構築する前に、これらの標準関数を正しく理解して活用することが推奨されます。

標準関数を利用することで、コードの可読性が高まるだけでなく、コンパイラによる最適化の恩恵も受けやすくなります。

特に浮動小数点数の扱いに長けた標準関数は、数値計算における精度の維持とエラー処理が適切に設計されている点が大きなメリットです。

cbrt関数の基本的な使い方

C99規格以降のC言語では、3乗根を直接計算するためのcbrt関数が導入されました。

cbrtという名称は「cube root（立方根）」の略称からきています。

この関数は非常にシンプルでありながら、数学的に正確な処理を目的として設計されています。

cbrt関数の定義と戻り値

cbrt関数は、引数として与えられた数値の3乗根を計算し、その結果を返します。

引数の型に合わせて、いくつかのバリエーションが存在します。

double cbrt(double x);：基本となるdouble型の計算。
float cbrtf(float x);：メモリ節約や速度優先のfloat型の計算。
long double cbrtl(long double x);：より高精度な計算が必要な場合に使用。

この関数の最大の特徴は、負の数に対しても正しい3乗根を返せるという点にあります。

例えば、-8の3乗根は-2ですが、cbrt(-8.0)は正確に-2.0を返却します。

cbrt関数の具体的なソースコード例

実際にcbrt関数を使用して、数値の3乗根を表示するプログラムを作成してみましょう。

C言語

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

int main(void) {
    double value1 = 27.0;
    double value2 = -64.0;
    double value3 = 10.0;

    // cbrt関数を使用して3乗根を算出
    double result1 = cbrt(value1);
    double result2 = cbrt(value2);
    double result3 = cbrt(value3);

    printf("%f の3乗根は %f です\n", value1, result1);
    printf("%f の3乗根は %f です\n", value2, result2);
    printf("%f の3乗根は %f です\n", value3, result3);

    return 0;
}

実行結果

27.000000 の3乗根は 3.000000 です
-64.000000 の3乗根は -4.000000 です
10.000000 の3乗根は 2.154435 です

このプログラムでは、正の整数、負の整数、そして割り切れない数値の3パターンで計算を行っています。

実行結果からわかる通り、負の値に対しても数学的に妥当な結果が得られています。

pow関数を使った3乗根の算出

C言語で累乗を計算する際によく使われるのがpow関数です。

pow(x, y)は「xのy乗」を計算する関数であり、これを利用して3乗根を求めることも可能です。

数学的には「xの3乗根」は「xの1/3乗」と等しいため、第2引数に1.0 / 3.0を指定することで代用できます。

pow関数の仕組みと引数の指定方法

pow関数を利用する場合は、指数部分の指定に注意が必要です。

1 / 3と記述してしまうと、整数同士の除算となり結果が0になってしまいます。

そのため、必ず浮動小数点数として1.0 / 3.0と記述しなければなりません。

C言語

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

int main(void) {
    double x = 125.0;
    // 1/3乗として計算
    double result = pow(x, 1.0 / 3.0);

    printf("%f の3乗根（pow経由）: %f\n", x, result);

    return 0;
}

実行結果

125.000000 の3乗根（pow経由）: 5.000000

pow関数における負の数の取り扱い

pow関数を使用して3乗根を求める際には、重大な制約が存在します。

多くの環境において、pow関数の第1引数に負の値を渡し、第2引数に整数でない値を渡すとドメインエラー（数学的定義域エラー）が発生します。

これは、pow(x, y)が内部的に指数関数と対数関数を用いてexp(y * log(x))として計算されることが多いためです。

負の数の対数は実数の範囲で定義されないため、プログラムがNaN（Not a Number）を返す原因となります。

したがって、負の数の3乗根を求める可能性がある場合は、pow関数の使用は避けるべきです。

cbrt関数とpow関数の違いと比較

3乗根を求めるにあたり、cbrt関数とpow関数のどちらを採用すべきかは重要です。

基本的には専用関数であるcbrtを使用するのがベストプラクティスですが、その理由を比較表で確認してみましょう。

比較項目	cbrt関数	pow関数 (x, 1.0/3.0)
負の値の処理	可能（正しい結果を返す）	不可能（NaNになる場合が多い）
計算精度	3乗根に最適化されている	一般的な累乗計算のため誤差が出やすい
処理速度	高速に動作するよう設計されている	汎用的なため、相対的に遅い場合がある
可読性	意図が明確（3乗根を求めている）	指数計算の意味を知っている必要がある

計算精度とパフォーマンスの差

cbrt関数は3乗根専用のアルゴリズムを内部で持っているため、極めて高い精度を維持します。

一方、pow関数は任意の累乗を計算する汎用関数であるため、1.0 / 3.0という割り切れない値を指数に指定した時点で、微小な誤差が蓄積されやすくなります。

また、計算コストの面でも、専用関数であるcbrtの方がCPUサイクルを効率的に利用できる傾向にあります。

近年のコンパイラは非常に優秀ですが、プログラムの意図を明確に伝えるためにもcbrtを選択するのが賢明です。

どちらの関数を選択すべきか

結論として、C99以降の環境であれば、迷わずcbrt関数を選択してください。

もし、非常に古い規格（C89/C90など）に従わなければならない制約がある場合のみ、pow関数を検討することになります。

ただし、その際も負の値が入力される可能性があるならば、事前に入力値を絶対値に変換し、後から符号を戻すといった処理を自前で実装する必要があります。

ニュートン法による3乗根の独自実装

ライブラリ関数を使用せずに3乗根を求める方法を知ることは、アルゴリズムの理解を深める上で非常に役立ちます。

その代表的な手法が「ニュートン法（ニュートン・ラフソン法）」です。

アルゴリズムの原理

ニュートン法は、関数 $f(x) = 0$ となる解を近似的に求めるための反復計算手法です。

3乗根を求める場合、方程式 $x^3 – a = 0$ を解くことに相当します。

この関数の接線を利用して次の近似値を求める更新式は、以下のように導出されます。

x_next = (2 * x + a / (x * x)) / 3

この計算を繰り返すことで、急速に真の値へと収束していきます。

ニュートン法を用いたプログラム例

標準ライブラリに頼らず、ニュートン法を用いて3乗根を計算する関数を作成してみましょう。

C言語

#include <stdio.h>

double my_cbrt(double a) {
    if (a == 0.0) return 0.0;
    
    // 負の数への対応
    double sign = (a < 0) ? -1.0 : 1.0;
    if (a < 0) a = -a;

    double x = a; // 初期値
    double prev;

    // 近似値が収束するまで繰り返す
    do {
        prev = x;
        x = (2.0 * x + a / (x * x)) / 3.0;
    } while (x != prev); // 完全に値が変わらなくなるまで（簡易的な判定）

    return x * sign;
}

int main(void) {
    double n = 27.0;
    printf("%f の独自3乗根計算: %f\n", n, my_cbrt(n));
    return 0;
}

実行結果

27.000000 の独自3乗根計算: 3.000000

このコードでは、浮動小数点数の特性を利用して値が変化しなくなるまでループを回していますが、実務では特定の誤差範囲（イプシロン）を下回った時点で終了させるのが一般的です。

自作の関数を持つことで、数学ライブラリが制限されている組み込み環境などでも3乗根計算が可能になります。

浮動小数点数型と精度の選択

C言語における数値計算では、使用するデータ型によって得られる結果の精度が変わります。

3乗根の計算においても、どの程度の桁数が必要なのかを考慮して型を選択する必要があります。

float型、double型、long double型の使い分け

それぞれの型の特性を理解し、適切に使い分けることが重要です。

float型：有効桁数は約7桁です。メモリ消費を抑えたい場合や、グラフィックス処理などで速度を優先する場合にcbrtfとともに使用します。
double型：有効桁数は約15桁です。C言語の数値計算における標準的な型であり、cbrt関数はこの型を対象としています。
long double型：有効桁数は環境に依存しますが、多くのシステムで18桁以上の精度を持ちます。科学技術計算など極めて高い精度が求められる場合にcbrtlを使用します。

現代のコンピュータではdouble型の計算がハードウェアレベルで高速化されているため、特別な理由がない限りはdouble型を選択するのが最もバランスが良いとされています。

精度の低いfloat型を使用すると、3乗根のような累乗計算を繰り返す中で誤差が無視できないレベルに膨らむ可能性があるため注意してください。

まとめ

C言語で3乗根を計算する最も確実で効率的な方法は、math.hヘッダで定義されているcbrt関数を使用することです。

cbrt関数は、負の数に対する正しい処理、高い計算精度、そして優れた実行速度を兼ね備えています。

一方でpow関数を用いた方法は、負の数の入力によるエラーのリスクや精度の面で劣るため、代用として使用する際には細心の注意が必要です。

また、ニュートン法のようなアルゴリズムを理解しておくことは、ライブラリが使えない特殊な環境下での開発において強力な武器となります。

それぞれの関数の特性を正しく把握し、プロジェクトの要件に合わせて最適な計算手法を選択するようにしましょう。