[討論] randi 超過上限 該如何處理

[JACK19920102 ()]

我需要在 0到20!
或0到30! 甚至0到50!

之間隨機生成一個正整數


但是randi()函數 括弧裡面 最多不可超過2^53
階層最大只能求到

randi(factorial(18)) =  5.4416e+14

randi(factorial(19)) 就超出randi的範圍了( 19! > 2^53 )

請問這有辦法解決嗎??

感謝大家回答 前五個推文的id各送100p幣 謝謝


--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.114.54.67
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1474531511.A.F9A.html
※ 編輯: JACK19920102 (140.114.54.67), 09/22/2016 16:06:11
推 sunev: randi(10,1,ceil(gammaln(51)/log(10)))-1
→ JACK19920102: s大 我試著用這個語法 但是寫出來的好像不是耶
→ JACK19920102: 還是我用錯了? 不太知道該怎麼用這個語法
→ JACK19920102: 如果是20!階層的話呢?
→ sunev: 好像?
→ celestialgod: 不能在0到log(20!)之間取一個隨機數 在轉回去取
→ celestialgod: ceiling嗎
→ celestialgod: 或是floor, round之類都好，保證是正數就好
→ celestialgod: 如果下界是1話，取ceiling，下界是0取floor
→ celestialgod: 這樣做好像機率會有問題Orz
推 sunev: 搞不好人家要研究benford's law啊