[问题] 大一程式设计作业请教

楼主: mpyh12345 (嘉义金城武)   2018-11-22 01:29:04
不好意思,作业应该要自己做的,
但是有笔测资输出了我怎么想都想不出来的结果。
https://i.imgur.com/vYuDM5e.png
怕格式乱掉 贴截图。
题目是要算组合 C n取k。
我一开始是先把分子跟分母分别算出来之后在相除,但这题有限制不能overflow。
于是图片上的做法我的想法就是假设C5取2,就是1*(5/2*4/1),但是因为只能够改函式部分
,cin的n,k,m一开始就是int,所以我在函式计算里面强制把n以及k转换成double。
问题来了,输入了一堆测资大部分都正确,结果C 8取3出错,正确应该是56,但是输出结果
跑出了55这样的奇妙结果,百思不得其解这个数字到底怎么跑出来的,所以想请各位帮我看
哪里出了问题。
作者: alan23273850   2018-11-22 01:40:00
真神奇,我自己试出来也是 55,很肯定是从 double转成 long long 的问题,但是解释不出理由用 float 的话就没错
作者: TMDTMD2487 (ㄚ冰)   2018-11-22 01:45:00
浮点数的运算都是有误差的,不然你把每一轮的c用double型态印出来看
作者: alan23273850   2018-11-22 01:47:00
https://goo.gl/5vCq8X 这边有人有尝试过
作者: oToToT (屁孩)   2018-11-22 01:47:00
return c+0.5应该就好了
作者: ckc1ark (伪物)   2018-11-22 01:49:00
c会是55.99999999999999289457
作者: alan23273850   2018-11-22 01:49:00
rounding error,所以有没有误差跟数字大小未必有关那 rounding error 会愈运算差愈多,所谓的误差传递
作者: ckc1ark (伪物)   2018-11-22 01:50:00
建议用8/1*7/2*6/3 这样整数除法不会出问题 只怕太大而已
作者: alan23273850   2018-11-22 01:50:00
这也就是为什么你直接把 56.0 转过去会对,但用运算的方式会错就是。clark大指的是 8*7*6/1/2/3 吧
作者: TMDTMD2487 (ㄚ冰)   2018-11-22 01:53:00
请问一下为什么不用递回做top down, 怕太慢就从下build up起来就好了C(n, m) = C(n – 1, m – 1) + C(n – 1, m)
作者: alan23273850   2018-11-22 01:55:00
他马的大大说的就是 dynamic programming 解
作者: TMDTMD2487 (ㄚ冰)   2018-11-22 01:56:00
这公式高中就交过了吧 这样不会超过ㄅDP难的地方是从问题找递回
作者: ckc1ark (伪物)   2018-11-22 01:57:00
8/1*7/2*6/3没错 过程中都会是整数
作者: TMDTMD2487 (ㄚ冰)   2018-11-22 01:59:00
n*m的阵列从0,0开始算到n,m 作业试试这样写 不行的话递回慢慢做也行
作者: alan23273850   2018-11-22 02:54:00
喔喔 是我误会clark大的意思ㄌ顺便提醒一下这种题目测资范围会影响解法的一定要给,不然就不是好题目
作者: IcecreamHsu (冰淇淋)   2018-11-23 02:08:00
for i=1; i<=n; p *= (m-i+1)/i; Cm取n的话这样就好噢抱歉推文有写了 请忽略

Links booklink

Contact Us: admin [ a t ] ucptt.com