作者chrisjon (随机数能吃吗？)
看板C_and_CPP
标题[问题] 哪里出问题了..=.=a?
时间Tue Apr 28 22:51:42 2009
#include <stdio.h>
#include <math.h>
int main()
{
float f;
int lcg_p,lcg_a;
int Lcg_seed;
lcg_p = 2147483647;
lcg_a = 16807;
scanf("%d", &Lcg_seed);
printf("%d\n", Lcg_seed);
Lcg_seed = Lcg_seed*lcg_a % lcg_p;
printf("%d\n", Lcg_seed);
f = Lcg_seed / lcg_p;
printf("%f\n", f);
return Lcg_seed;
}
==================================
简单的数学式表示就是
X[i] = a*X[i-1] mod p
mod 为取余数
==================================
我重新修改过
compile成功
但是
f完全没有计算，print出来是0.00000
另外，return是整个程式计算完之後，要回传的值嘛?
所以我return後，应该会回传Lcg_seed的值不是吗?
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 123.195.128.31
1F：推 BSpowerx:取余数限定用在整数。一定要将取余数的数改成int型态 04/28 22:54
2F：→ sunneo:浮点数的余数要用 math.h的 fmod 04/28 22:56
3F：→ sunneo:另外main的回传值必须是整数 04/28 22:56
4F：→ chrisjon:可是，如果取int而不用double的话，p算不会爆量吗? 04/28 23:17
5F：→ chrisjon:例如int上限100，而a*b % 3 = 101 % 3 = 2 不会爆掉吗? 04/28 23:18
6F：推 littleshan: (a * b) mod c = (a mod c)*(b mod c) mod c 04/28 23:24
7F：→ chrisjon:a<c,b<c，但有可能a*b>c，而c是2^31 - 1，谢谢 04/28 23:36
※ 编辑: chrisjon 来自: 123.195.128.31 (04/29 00:07)
8F：→ MOONRAKER:f = Lcg_seed / (float) lcg_p; 新手级问题 04/29 02:07
谢谢，只学过其他语言，临时因用到的套装软体需要而来了解c语言
在图书馆借几书要在几个月内就要弄好，真的有点赶，好在有学过其他语言
能比较快一点进入状况。
真不好意思问这些新手问题，但周围又没其他人会c语言，所以只能
麻烦板上各位先进帮我这临阵磨枪的新手一个忙了^^"
感恩
※ 编辑: chrisjon 来自: 123.195.128.31 (04/29 04:41)
9F：→ MOONRAKER:整数/整数会用整数除法算，所以你至少要把其中一者强制 04/29 09:53
10F：→ MOONRAKER:转为浮点，他才会用浮点除法计算 04/29 09:54
11F：→ chrisjon:唔...这就麻烦了..Lcg_seed和lcg_p前部份要设成int 04/29 22:15
12F：→ chrisjon:这样才能计算%，现在又要设成float...囧" 04/29 22:15
13F：→ sunneo:推过了 fmod ... 04/30 00:32
14F：→ sunneo:会用内建函式就省了麻烦了不是 ? 04/30 00:32
15F：→ sunneo:还是你压根不想使用内建函式？ 04/30 00:33
16F：推 VictorTom:也搞不好是作业, 老师规定不准用内建函数....XD 04/30 00:38
17F：→ VictorTom:遥想以前计组, 老师作业是要我们自己写把ASCII 123.456 04/30 00:38
18F：→ VictorTom:这样的user input, 转换并写入float f去, IEEE754啊Orz 04/30 00:39
19F：→ sunneo:不过那可真令人疑惑了呢 因为他连math.h都include了... 04/30 01:11
20F：→ VictorTom:我看他好像没真的call到math.h的lib就是了, 不过小弟也 04/30 01:23
21F：→ VictorTom:只是乱猜啦; 真要的话还是得请原po自己说明清楚了:) 04/30 01:24
22F：→ chrisjon:fmod那个是我漏看了...^^"拍洗 04/30 02:18
如果要说是作业，其实也不为过...因为是论文...~.~
http://www.iro.umontreal.ca/~simardr/testu01/tu01.html
pdf载点
http://www.iro.umontreal.ca/~simardr/testu01/guideshorttestu01.pdf
不过这是guide short testu01
我是看 guide long testu01
我现在在看第二章的部份
苦於英文及组合能力太差....单字查完看得懂，全段意思组不起来...囧"
只知道里面有提到 "非常简单"、"使用外来的产生器做检定很简单"
.....= ="
然後大概也看得懂部份程式，但是就是...不会使用....~.~
我现在先试着把我已经用Maple语言写好的程式码转成C语言来跑
然後再运用cygwin去测转好的程式有没有问题
说明书嘛...再慢慢看...总会看得通的....
※ 编辑: chrisjon 来自: 123.195.128.31 (04/30 02:31)
23F：推 sunneo:但是不管怎麽说 X[i] = a*X[i-1] mod p 04/30 02:38
24F：→ sunneo:其中X[ i ]如果是浮点数 那麽这个式子应该是 04/30 02:39
25F：→ sunneo:X[i] = a * fmod( X[i-1], p ) 04/30 02:39
26F：→ sunneo:其中fmod就是浮点数该使用的余数函式 04/30 02:39
27F：→ sunneo:然後 f = Lcg_seed / lcg_p 这个式子你会永远得到0 04/30 02:41
28F：→ chrisjon:找到了fmod的说明，谢谢 用法了解中 04/30 02:57
29F：→ chrisjon:原来 /(float)lcg_p 也是一种用法= =" 04/30 03:13