作者Arton0306 (Ar藤)
看板Python
標題python的浮點數精度
時間Sun Sep 13 15:00:56 2009
下方有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)
1F：→ Arton0306:就算c的不用double改用float就好 還是比python的準 09/13 15:05
2F：→ Lucemia:python 的float 在C層是用double implement的 09/14 01:06
3F：→ AndCycle:老闆你兩邊輸出的浮點位數不一樣多呀 _A_ 09/14 03:42
5F：→ AndCycle:人家抓不出你的code問題在哪, 所以照C寫了一個py, 沒問題 09/14 04:45
6F：→ AndCycle:所以你的py一定有哪裡做的事跟C不一樣, 請加油 _A_ 09/14 04:52
7F：→ Arton0306:浮點位數沒差 不過你的code確實是ok的 看來是我哪寫錯了 09/14 05:41
8F：→ Arton0306:感謝你還花時間改寫!!! 09/14 05:42
9F：→ Arton0306:我找到錯誤了 python版的C(n,m) 如果m是負的 就錯了 09/14 06:17
10F：→ Arton0306:真是汗顏...囧> 09/14 06:18