两个大阶乘相除溢出。急救。

[阿博]

我要求一个组合问题。公式为 N! / ( k!*(n-k)! )
当Ｎ＝２９２时，并且仅当 k = 132 ...160 时，我的算法走不下去了。
我的算法已经作如下优化处理
（１）分子分母已经约去。如：假设 n-k > k ，则公式变为：
n*(n-1)*(n-2)*...*(n-k+1)
-------------------------
k*(k-1)*(k-2)*...*1
（２）分子分母都把偶数除以２，即： n,n-1,n-2 ... , k,k-1,k-2，... 凡是偶数都先除以２再参与计算。

　发现：当 k=140时，还能有结果。当 k=145，ＪＶＭ提示：Infinity

　－－－－　求求兄弟们帮忙，建议一个好方法　。以下是源码：

public static double zhuhe(int n , int k ){
double N1 = 1d ; // 存放公式中的被除数 N! 约去 max(K!,(N-k)!)后的值。
double N2 = 1d ; // 存放公式中 min(K!,(N-K)!)中的值。
int i_max = Math.max(k,n-k);
int i_min = Math.min(k,n-k);
int i_N1_oushu = 0; // 记录分子(N1)最小值是否偶数。1表示偶数，０表示奇数。
if( (float)i_max/2.0 > Math.round((float)i_max/2.0) ){
i_N1_oushu = 1 ; // i_max是奇数。分子(N1)最小值(i_max+1)是偶数
}
int i_N2_oushu = 0; // 分母肯定是奇数开始，因为从１开始。

int m = 1;
while ( m <= i_min ){
if( i_N1_oushu == 1 ){
N1*=((m+i_max)/2);
i_N1_oushu = 0 ;
}
else {
N1 *= (m + i_max);
i_N1_oushu = 1 ;
}
if( i_N2_oushu == 1) {
N2*=m/2;
i_N2_oushu = 0 ;
}
else{
N2 *= m;
i_N2_oushu = 1 ;
}
m++;
}
return N1/N2;
}

[qiyongjun2003]

用大整型java.math.BigInteger.

[追求自由]

UP

[lee_nan]

n (n-1) (n-2) (n-k+1)
-* ---- * ---- *..... * ------
k (k-1) (k-2) 1
可否？

[阿博]

谢谢。。。正在采用　java.math.BigDecimal　。

---　　我马上改用 BigInteger

[flyforlove]

你不应该先求N！然后再计算。应该每项都计算完后再相乘，像lee_nan() ( )所说的那样。

否则本来结果不会溢出，但是在计算过程中却出现了溢出的值，这很不合理。

[阿博]

昨晚已经改用BigInteger完成。溢出问题解决了，但是误差问题出来了。

但是结合其它部分的功能求出的概率的总和不为1（或者不接近1）了。
原来采用double时，在100以内的基数进行计算时，概率总和是非常接近1的。

郁闷！！！！！！

谢谢 flyforlove(为情飞)、lee_nan() ，double相乘，也经常出现相当大的误差的，所以，我一直没有采用这种做法。
不过，今天我得试一试。

[过客猫2022]

数值太大，就不用这个，因为精度有限，网上有大把的例子，介绍高精度计算!

[阿博]

我就是在google.com 搜了，找不到有关 浮点。
记得大学时，关于数学的计算，有 MathLab 软件。

－－－　现在一搜“高精度计算”，确实有很多。如：超大整数高精度快速计算器

谢谢各位，欢迎继续提供意见。

等我把问题完全解决了，我再来结帐。