作者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