急求！！！float数乘法溢出问题！

[Civiano]

如图所示

A_ij取0或1,theta是个(0,1)区间上的小数,i,j的区间为(0,N),s的区间为(0,C)。在做累乘时,如果N很大,常常会导致分子分母都为0。
另一方面，由于分母的存在，似乎不能采用取对数再运算的方法。
求各位大神指点！

[iasky]

这个不是很懂，建议参考numpy这个库。 http://www.numpy.org/

[panghuhu250]

首先, 把记号简化一下, 令b_s = pi_r*product of(theta...), 你的公式可以表示为: b_r/(sum of b_s), 分子分母同除以b_r,得到1/sum(b_s/b_r). 只要能计算b_s/b_r就行了, b_s/b_r完全是乘除法, 可以通过log计算: 先算

log(pi_r)+sum of A_ij*log(theta_rj) - log(pi_s) - sum of A_ij*log(theta_sj)

再取指数. 分母中有一个是1, 其它的b_s/b_r也一定大于0, 所以分母一定不会是0.

[panghuhu250]

log(pi_r)+sum of A_ij*log(theta_rj) - log(pi_s) - sum of A_ij*log(theta_sj)

这个反了,应该是:

-log(pi_r) - sum of A_ij*log(theta_rj) + log(pi_s) + sum of A_ij*log(theta_sj)

[maniachhz]

数学运算太神秘.....

[Bird_1989]

完全看不懂的运算