6

四分位数を見つけるための次のコードがあります。

#include <stdio.h>
#include <stdlib.h>

typedef struct {
    double qrt[3];
    double *value;
    int count;
} t_data;

static void set_qrt(t_data *data, int qrt)
{
    int n, e;
    double d;

    d = qrt * 0.25 * data->count + 0.5;
    n = (int)d;
    e = n != d;
    data->qrt[qrt - 1] = data->value[n - 1];
    if (e) {
        data->qrt[qrt - 1] += data->value[n];
        data->qrt[qrt - 1] *= 0.5;
    }
}

static void set_qrts(t_data *data)
{
    set_qrt(data, 2);
    if (data->count > 1) {
        set_qrt(data, 1);
        set_qrt(data, 3);
    } else {
        data->qrt[0] = 0.0;
        data->qrt[2] = 0.0;
    }
}

static int comp(const void *pa, const void *pb)
{
    const double a = *(const double *)pa;
    const double b = *(const double *)pb;

    return (a > b) ? 1 : (a < b) ? -1 : 0;
}

int main(void)
{
    double values[] = {3.7, 8.9, 7.1, 5.4, 1.2, 6.8, 4.3, 2.7};
    t_data data;

    data.value = values;
    data.count = (int)(sizeof(values) / sizeof(double));
    qsort(data.value, data.count, sizeof(double), comp);
    set_qrts(&data);
    printf("Q1 = %.1f\nQ2 = %.1f\nQ3 = %.1f\n", data.qrt[0], data.qrt[1], data.qrt[2]);
}

は

d = qrt * 0.25 * data->count + 0.5;
n = (int)d;
e = n != d;

期待どおりに動作することが保証されていますか？(e == isinteger(d))

4

4 に答える 4

1

倍精度浮動小数点形式のマニティッサには 53 ビットがあり、そのうちの 1 つは暗黙的です。これは、2^0 から 2^53-1 の範囲のすべての正負の整数を表すことができることを意味します。

0 (ゼロ) は、独自の形式を持つ特殊なケースです。

0.25 の間隔になると、範囲は単純に 2^-2 から 2^51-0.25 と計算されます。これは、0.25 のすべての倍数が倍精度形式で正確に表現できるわけではありませんが、かなりの数の整数が正確に表現できるわけではないことを意味します。

したがって、正確に表現可能な間隔が 2^x の場合、表現可能な範囲は 2^x から 2^(53+x)-2^x です。

于 2013-07-30T13:24:24.103 に答える
1

dasblinkenlight は完全に正しいです。double/float 型と整数型は、IEEE754 に従って異なる方法で格納されます。興味のある方は、こちらの簡単なチュートリアルをご覧ください。

于 2013-07-30T10:47:28.760 に答える