python的浮點數精度

[Arton0306 (Ar藤)]

下方有python和c版本的兩種code
做的事是一樣的
但算出來 兩者在小數點後幾位卻不太一樣
而且c的才是較準確的
這是google code jam的題目 c的可以過 但python不行
題目 http://code.google.com/codejam/contest/dashboard?c=188266#s=p2

請問為什麼會這樣??

python版==================
import psyco
psyco.full()

floor=[1,1]
for i in range(2,11):
    floor.append(i*floor[-1])

def C(n,m):
    #print n,m
    return ((floor[n]+0.0)/floor[m]/floor[n-m])

def K(t,c,r,dp):
    if r in dp:
        return dp[r]

    sq=0.0
    for i in range(1,r+1):
        sq+=C(r,i)*C(t-r,c-i)/C(t,c)*(1+K(t,c,r-i,dp))
    cir=C(r,0)*C(t-r,c)/C(t,c)
    ret=(sq+cir)/(1-cir)
    dp[r]=ret
    return ret


casenu=int(raw_input())
for i in xrange(casenu):
    c,n=[int(e) for e in raw_input().split()]
    dp={}
    dp[0]=0
    print "Case #%d: %.7f" % (i+1, 1+K(c,n,c-n,dp))


c版本==============================================

#include <iostream>

using namespace std;

double floor(int x){
    if (x==0) return 1;
    double ret=1;
    for (int i=x; i>=1; i--){
        ret*=i;
    }
        return ret;
}

double com(int n, int m){
    return ((floor(n))/floor(m)/floor(n-m));
}

double C[11][11];

double K(int t, int c, int r){
    if (r==0) return 0;
    double sq=0.0;
    for (int i=1; i<=r; i++){
        sq+=C[r][i]*(C[t-r][c-i]/C[t][c])*(1+K(t,c,r-i));
    }
    double cir=C[r][0]*(C[t-r][c]/C[t][c]);
    double ret=(sq+cir)/(1-cir);
    return ret;
}

int main(){

        for (int i=0; i<11; i++){
                for (int j=0; j<11; j++){
                        if (j>i){
                                C[i][j]=0;
                        }
                        else {
                                C[i][j]=com(i,j);
                        }
                }
        }

    int casenum,c,n;
    cin>>casenum;
    for (int caseid=1; caseid<=casenum; caseid++){
                cin>>c>>n;
        printf("Case #%d: %.8f\n", caseid, 1+K(c,n,c-n));
    }

    return 0;
}

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 220.135.140.149
※ 編輯: Arton0306       來自: 220.135.140.149      (09/13 15:02)
→ Arton0306:就算c的不用double改用float就好 還是比python的準
→ Lucemia:python 的float 在C層是用double implement的
→ AndCycle:老闆你兩邊輸出的浮點位數不一樣多呀 _A_
推 AndCycle:http://www.andcycle.idv.tw/~andcycle/tmp/20090914.py
→ AndCycle:人家抓不出你的code問題在哪, 所以照C寫了一個py, 沒問題
→ AndCycle:所以你的py一定有哪裡做的事跟C不一樣, 請加油 _A_
→ Arton0306:浮點位數沒差 不過你的code確實是ok的 看來是我哪寫錯了
→ Arton0306:感謝你還花時間改寫!!!
→ Arton0306:我找到錯誤了 python版的C(n,m) 如果m是負的 就錯了
→ Arton0306:真是汗顏...囧>