A Book on C 1장 연습문제 17번...

lsj0713의 이미지

/* Compute a better average. */

#include <stdio.h>

int main(void)
{
    int i;
    double x;
    double avg = 0.0;    /* a better average */
    double navg;    /* a naive average */
    double sum = 0.0;

    printf("%5s%17s%17s%17s\n%5s%17s%17s%17s\n\n",
        "Count", "Item", "Average", "Naive avg",
        "-----", "----", "-------", "---------");
    for (i = 1; scanf("%lf", &x) == 1; ++i )
    {
        avg += (x - avg) / i;
        sum += x;
        navg = sum / i;
        printf("%5d%17e%17e%17e\n", i, x, avg, navg);
    }
    return 0;
}

17. 이 연습 문제는 연습 문제 16번에 계속되는 것이다. better_average 프로그램의 입력으로 위와 같은 큰 숫자가 아닌 평범한 숫자를 하용한다면 Average 열과 Naive Avg열은 같을 것처럼 보인다. 두 열의 값이 달라지는 경우를 찾아보아라. 즉, sum이 오버플로 되지 않을 때에도 better_average가 실제로 더 좋은 방법이라는 것을 실험적으로 증명하여라.

A Book on C 1장 연습문제 17번입니다. 도저히 오버플로 이외의 이 알고리즘의 장점을 찾을 수가 없습니다. 해답 아시는 분은 도움 주시길 부탁드립니다. 다른 문제는 다 풀겠는데 이 문제는 지금에 이르러서도 도저히 해답이 나오질 않는군요-_-;

답변해주실 분들께 미리 감사드립니다.

winner의 이미지

double 인 것을 보니 유효숫자문제인 것 같은데요.

#include <stdio.h>

int main(void)
{
    int i;
    double x;
    double avg = 0.0;    /* a better average */
    double navg;    /* a naive average */
    double sum = 0.0;

    printf("%5s%24s%24s%24s\n%5s%24s%24s%24s\n\n",
        "Count", "Item", "Average", "Naive avg",
        "-----", "----", "-------", "---------");
    for (i = 1; scanf("%lf", &x) == 1; ++i )
    {
        avg += (x - avg) / i;
        sum += x;
        navg = sum / i;
        printf("%5d%25.15e%25.15e%25.15e\n", i, x, avg, navg);
    }
    return 0;
}


Count                    Item                 Average               Naive avg
-----                    ----                 -------               ---------

0.1
    1    1.000000000000000e-01    1.000000000000000e-01    1.000000000000000e-01
0.2
    2    2.000000000000000e-01    1.500000000000000e-01    1.500000000000000e-01
0.3
    3    3.000000000000000e-01    2.000000000000000e-01    2.000000000000000e-01
0.4
    4    4.000000000000000e-01    2.500000000000000e-01    2.500000000000000e-01
0.5
    5    5.000000000000000e-01    3.000000000000000e-01    3.000000000000000e-01
0.6
    6    6.000000000000000e-01    3.500000000000000e-01    3.500000000000000e-01
0.7
    7    7.000000000000000e-01    4.000000000000000e-01    4.000000000000000e-01
0.8
    8    8.000000000000000e-01    4.500000000000000e-01    4.500000000000000e-01
0.9
    9    9.000000000000000e-01    4.999999999999999e-01    5.000000000000000e-01
pynoos의 이미지

float, double 형은 잘 안써봐서 수치해석책 좀 찾아보니..

http://www.google.co.kr/search?hl=ko&ie=UTF-8&oe=UTF-8&newwindow=1&q=catastrophic+cancellation+float+point&lr=

이것과 관련한 문제가 아닐까요?

lampo의 이미지

a book on C 연습문제 사이트입니다. ^.^

http://cosmos.ssu.ac.kr/abc/4ed/

흠. 그곳에서도 17 번 정답을 달아논 분은 없군요...

지리즈의 이미지

그렇군요... 부동소수점 연산임으로 비슷한 크기의 숫자간의 계산이 더 정확해 진다는 뜻이네요...

그렇다면... 달라지는 답을 구하기 위해서는 sum이
부동소수점에서 허용하는 자리수 보다가 길어지게 되는 순간 값이 차이가
발생하겠군요..

pynoos wrote:
float, double 형은 잘 안써봐서 수치해석책 좀 찾아보니..

http://www.google.co.kr/search?hl=ko&ie=UTF-8&oe=UTF-8&newwindow=1&q=catastrophic+cancellation+float+point&lr=

이것과 관련한 문제가 아닐까요?

There is no spoon. Neo from the Matrix 1999.

댓글 달기

Filtered HTML

  • 텍스트에 BBCode 태그를 사용할 수 있습니다. URL은 자동으로 링크 됩니다.
  • 사용할 수 있는 HTML 태그: <p><div><span><br><a><em><strong><del><ins><b><i><u><s><pre><code><cite><blockquote><ul><ol><li><dl><dt><dd><table><tr><td><th><thead><tbody><h1><h2><h3><h4><h5><h6><img><embed><object><param><hr>
  • 다음 태그를 이용하여 소스 코드 구문 강조를 할 수 있습니다: <code>, <blockcode>, <apache>, <applescript>, <autoconf>, <awk>, <bash>, <c>, <cpp>, <css>, <diff>, <drupal5>, <drupal6>, <gdb>, <html>, <html5>, <java>, <javascript>, <ldif>, <lua>, <make>, <mysql>, <perl>, <perl6>, <php>, <pgsql>, <proftpd>, <python>, <reg>, <spec>, <ruby>. 지원하는 태그 형식: <foo>, [foo].
  • web 주소와/이메일 주소를 클릭할 수 있는 링크로 자동으로 바꿉니다.

BBCode

  • 텍스트에 BBCode 태그를 사용할 수 있습니다. URL은 자동으로 링크 됩니다.
  • 다음 태그를 이용하여 소스 코드 구문 강조를 할 수 있습니다: <code>, <blockcode>, <apache>, <applescript>, <autoconf>, <awk>, <bash>, <c>, <cpp>, <css>, <diff>, <drupal5>, <drupal6>, <gdb>, <html>, <html5>, <java>, <javascript>, <ldif>, <lua>, <make>, <mysql>, <perl>, <perl6>, <php>, <pgsql>, <proftpd>, <python>, <reg>, <spec>, <ruby>. 지원하는 태그 형식: <foo>, [foo].
  • 사용할 수 있는 HTML 태그: <p><div><span><br><a><em><strong><del><ins><b><i><u><s><pre><code><cite><blockquote><ul><ol><li><dl><dt><dd><table><tr><td><th><thead><tbody><h1><h2><h3><h4><h5><h6><img><embed><object><param>
  • web 주소와/이메일 주소를 클릭할 수 있는 링크로 자동으로 바꿉니다.

Textile

  • 다음 태그를 이용하여 소스 코드 구문 강조를 할 수 있습니다: <code>, <blockcode>, <apache>, <applescript>, <autoconf>, <awk>, <bash>, <c>, <cpp>, <css>, <diff>, <drupal5>, <drupal6>, <gdb>, <html>, <html5>, <java>, <javascript>, <ldif>, <lua>, <make>, <mysql>, <perl>, <perl6>, <php>, <pgsql>, <proftpd>, <python>, <reg>, <spec>, <ruby>. 지원하는 태그 형식: <foo>, [foo].
  • You can use Textile markup to format text.
  • 사용할 수 있는 HTML 태그: <p><div><span><br><a><em><strong><del><ins><b><i><u><s><pre><code><cite><blockquote><ul><ol><li><dl><dt><dd><table><tr><td><th><thead><tbody><h1><h2><h3><h4><h5><h6><img><embed><object><param><hr>

Markdown

  • 다음 태그를 이용하여 소스 코드 구문 강조를 할 수 있습니다: <code>, <blockcode>, <apache>, <applescript>, <autoconf>, <awk>, <bash>, <c>, <cpp>, <css>, <diff>, <drupal5>, <drupal6>, <gdb>, <html>, <html5>, <java>, <javascript>, <ldif>, <lua>, <make>, <mysql>, <perl>, <perl6>, <php>, <pgsql>, <proftpd>, <python>, <reg>, <spec>, <ruby>. 지원하는 태그 형식: <foo>, [foo].
  • Quick Tips:
    • Two or more spaces at a line's end = Line break
    • Double returns = Paragraph
    • *Single asterisks* or _single underscores_ = Emphasis
    • **Double** or __double__ = Strong
    • This is [a link](http://the.link.example.com "The optional title text")
    For complete details on the Markdown syntax, see the Markdown documentation and Markdown Extra documentation for tables, footnotes, and more.
  • web 주소와/이메일 주소를 클릭할 수 있는 링크로 자동으로 바꿉니다.
  • 사용할 수 있는 HTML 태그: <p><div><span><br><a><em><strong><del><ins><b><i><u><s><pre><code><cite><blockquote><ul><ol><li><dl><dt><dd><table><tr><td><th><thead><tbody><h1><h2><h3><h4><h5><h6><img><embed><object><param><hr>

Plain text

  • HTML 태그를 사용할 수 없습니다.
  • web 주소와/이메일 주소를 클릭할 수 있는 링크로 자동으로 바꿉니다.
  • 줄과 단락은 자동으로 분리됩니다.
댓글 첨부 파일
이 댓글에 이미지나 파일을 업로드 합니다.
파일 크기는 8 MB보다 작아야 합니다.
허용할 파일 형식: txt pdf doc xls gif jpg jpeg mp3 png rar zip.
CAPTCHA
이것은 자동으로 스팸을 올리는 것을 막기 위해서 제공됩니다.