SORU
18 EYLÜL 2008, PERŞEMBE


Nasıl okunabilir kesirler için yüzer dönüştürmek için?

Hadi çıktı etmeliyiz 0.33, mesela "1/3".< / ^ br . Eğer varsa "" ihtiyacımız olan çıkış "2/5". 0.4

Fikri okunabilir insan kullanıcı anlamasını sağlamak için "x" anlayış verilerin daha iyi bir yol gibi. y dışı parçalar

Yüzdeleri iyi bir yedek olur ama eğer bunu yapmak için basit bir yolu varsa merak ettiğimi biliyor musun?

CEVAP
18 EYLÜL 2008, PERŞEMBE


David 2 ** C kodu Eppstein için soruyorsun tam olarak bulduk. Sürekli kesirler teorisi üzerine dayalı ve çok hızlı ve oldukça kompakt.

Bu pay ve payda belirli sınırlar için özelleştirilmiş sürümleri kullandım.

/*
** find rational approximation to given real number
** David Eppstein / UC Irvine / 8 Aug 1993
**
** With corrections from Arno Formella, May 2008
**
** usage: a.out r d
**   r is real number to approx
**   d is the maximum denominator allowed
**
** based on the theory of continued fractions
** if x = a1   1/(a2   1/(a3   1/(a4   ...)))
** then best approximation is found by truncating this series
** (with some adjustments in the last term).
**
** Note the fraction can be recovered as the first column of the matrix
**  ( a1 1 ) ( a2 1 ) ( a3 1 ) ...
**  ( 1  0 ) ( 1  0 ) ( 1  0 )
** Instead of keeping the sequence of continued fraction terms,
** we just keep the last partial product of these matrices.
*/

#include <stdio.h>

main(ac, av)
int ac;
char ** av;
{
    double atof();
    int atoi();
    void exit();

    long m[2][2];
    double x, startx;
    long maxden;
    long ai;

    /* read command line arguments */
    if (ac != 3) {
        fprintf(stderr, "usage: %s r d\n",av[0]);  // AF: argument missing
        exit(1);
    }
    startx = x = atof(av[1]);
    maxden = atoi(av[2]);

    /* initialize matrix */
    m[0][0] = m[1][1] = 1;
    m[0][1] = m[1][0] = 0;

    /* loop finding terms until denom gets too big */
    while (m[1][0] *  ( ai = (long)x )   m[1][1] <= maxden) {
        long t;
        t = m[0][0] * ai   m[0][1];
        m[0][1] = m[0][0];
        m[0][0] = t;
        t = m[1][0] * ai   m[1][1];
        m[1][1] = m[1][0];
        m[1][0] = t;
        if(x==(double)ai) break;     // AF: division by zero
        x = 1/(x - (double) ai);
        if(x>(double)0x7FFFFFFF) break;  // AF: representation failure
    } 

    /* now remaining x is between 0 and 1/ai */
    /* approx as either 0 or 1/m where m is max that will fit in maxden */
    /* first try zero */
    printf("%ld/%ld, error = %e\n", m[0][0], m[1][0],
           startx - ((double) m[0][0] / (double) m[1][0]));

    /* now try other possibility */
    ai = (maxden - m[1][1]) / m[1][0];
    m[0][0] = m[0][0] * ai   m[0][1];
    m[1][0] = m[1][0] * ai   m[1][1];
    printf("%ld/%ld, error = %e\n", m[0][0], m[1][0],
           startx - ((double) m[0][0] / (double) m[1][0]));
}

Bunu Paylaş:
  • Google+
  • E-Posta
Etiketler:

YORUMLAR

SPONSOR VİDEO

Rastgele Yazarlar

  • DeAdPiXel6667

    DeAdPiXel666

    2 Ocak 2010
  • Sergio Lafuente Rubio

    Sergio Lafue

    11 Aralık 2008
  • SPBedition

    SPBedition

    24 HAZİRAN 2013