计算机常用算法与程序设计教程 第7章 模拟.ppt

1、根据公用算法和编程、第1、7章、模拟、公用算法和编程、第2章、主要内容、公用算法和编程、3，7.1模拟概述模拟对象的特性，电脑模拟可以分为决定性模拟和随机性模拟。决定性的模拟是对决定性过程的模拟，其模拟的事件按照其固有的规律发展。例如，运算模拟是决定性的模拟。随机性模拟对象是使用随机数作为参数实现模拟的随机事件。例如，数字模拟(也称为数字模拟)牙齿。常用算法和编程，4，运算模拟根据整数的四个茄子算法进行模拟，然后产生模拟结果。7.2.1运算模拟描述了运算模拟，主要通过模拟整数位乘法和除法来解决一些整数计算问题。在实施乘法和除法模拟之前，必须根据参与整数的实际设置设置模拟量，以模拟乘法和除法过程

2、的数值变化，并确定运算是否结束。7.2计算模拟，通用算法和编程，5，1。除模拟框架描述外的除法模拟：确定输入while()a=c* 10m；/*被除数a配置，m表示*/b=a/p；/*除以实现，计算程序b */printf(b)；c=a % p；/*在评估中，剩下的c */其中和应基于模拟除以计算问题的具体实际决定。公共算法和编程，6，乘法模拟帧说明：确定输入while()k=k 1；a=w(k)* p m；/*计算产品a，m表示*/w(k)=a；/*乘积a的位存储在w(k) */m=a/10中。/*乘积a的10位以上是小数*/输出(w (d)、w (d-1)、w(1)；/*根据高层到我的输出

3、*/乘法模拟、模拟乘法问题的实际决定。2 .模拟乘法运算，常用算法和编程，7，1。N个1除以2009的问题示例7.1由N个1组成的整数可除以2009，N至少有多少？模拟除以设计：被除数是A，除数p=2009，每个评价师的余数是C。被除数a=c*10 1，每个测试的收入馀数为c=a 09。设定初始值c=1111，n=4以进入模拟除以回路。循环条件为c0。每循环增加变量n牙齿1。其馀c=0、结束、输出n值。7.2.2 n 1中的分区问题，公用算法和编程，8，void main() int a，c，n；C=1111n=4；/*确定初始值*/while(c！=0)a=c * 10 1；/*被除数配置a

4、 */c=a 09；n；/*除法实现，剩馀c */printf(至少%d个1.n，n)；通用算法和编程，9，2。N个1积累的数字游戏示例7.2两个电脑爱好者正在进行“N个1积累的数字游戏”，被题数a=c*10 1，剩下的c=a%p，商b=a/p是寻找数q之一。如果剩馀c=0，则结束。否则，继续工作，直到c=0牙齿。通用算法和编程，10，void main() int a，b，c，p，n；Printf(输入n整数p :)；Scanf(%d，公共算法和编程，11，7.2.3端数前进问题，示例7.3整数N的端数为9，将端数9移到前面(成为最高值)后的结果数是原始值我们请求将一般尾数前进问题：整数N的

5、尾数Q(限于一个)移到N前面的P倍。记为n(q，p)。我保证在这里。解决指定端数q和倍数p的n(q，p)。下面用模拟除法运算和模拟乘法两种茄子方法设计求解。常用算法和编程，12，1。模拟整数除以，第一个数字Q除以P(注意合意qp)，剩下的为C，份额为B。将数字B输出为所需N的第一个数字。进入模拟循环，剩下的c=0和商b=q结束后，循环条件为C！=0 | b！=Q在循环中，计算被除数a=c*10 b，试验者获得b=a/p，输出计算为n的b。取得剩馀c=a% p。然后B和C制作下一个测试商的被除数，并据此递归。通用算法和编程，13，void main() int a，b，c，p，q；Scanf(%

6、d，%d，公共算法和编程，14，设置存储数N的W数组)。w(1)=从q开始，乘以乘数p和n的每个数字w(i)，再加上数量级m，则a=w(k)* p m；累积A的10位数以上是下一轮的小数，m=a/10。a的位数现在必须指定乘积的下一个w (I 1)=x的值。计算的被除数a为尾部q时结束。2 .模拟整数乘法，一般算法和编程，15，voidmain () inta，m，j，k，p，q，w100Scanf(%d，%d，常用算法和编程，16，7.2.5圆周率，圆周率计算的历史很长。我国数学家曹冲首先将圆周率计算为3.1415926，领先世界一千多年。以后；以后。德国数学家鲁特夫计算到小数点后35位。日

7、本数学家建富玄红计算为第41位。1874年，英国数学家祥子利用微积分将一生精力计算到了第707位，但第528位后的数字是错误的。公共算法和编程，17，公共算法和编程，18，公共算法和编程，19，(3)模拟乘法和除法合成计算设置a数组，计算出的整数值存储在a(0)中，小数点后I位存储在a中的公式也就是说，阵列将初始值a(0)=1后，除以2n 1，再乘以N，然后加1。除以2n-1，乘以n-1，然后加1。这些阵列操作设置在循环中实现。公共算法和编程，20，模拟乘法和除法说明：for(c=1，j=n；j=1；根据J-)/*公式，n次*/d=2 * j 1；for(I=0)；I=0；I - ) /*各位

8、实现乘以j */a(I)=a(I)* j b；b=a(I)/10；a(I)=a(I)；a(0)=a(0)1；c=a(0)；/*整数位为1 */for(b=0，I=x 5；I=0；根据I - ) /*公式，乘以2 */a(i)=a(i)*2 b。b=a(I)/10；a(I)=a(I)；公共算法和编程、21、(4)输出结果循环实现分割操作完成后，按数组元素从高到低的顺序输出。因为计算位数很多，所以可以控制每行50位打印，每10位一格，以便于确认。(5)复杂性分析上面的模拟乘法和除法算法的时间复杂度O(nlogn)，其中N是所需计算中的位数，logn是所需计算中的项数。将常用算法和编程、22，7.3

9、.3模拟卡示例7.8模拟卡升级许可证、包括大小王在内的共54张卡随机分为4张、每家12张、6张底卡。(1)模拟图案和点模拟卡要注意随机性，重复渡边杏。相当于四个茄子图案，设置任意整数X，其值为1-4。设置与每种颜色的13点相对应的随机整数Y，其值为1-13。为了避免重复，将x和y合并为z=x*100 y的3位数，并存储在数组m(54)中。要发行I 1卡，请使用数字z和现有的m(0)、m(1)、逐个比较m(i-1)，如果相同，则重新生成任意整数x和y得z进行比较。7.3随机模拟，常用算法和编程，23，(2)模拟大小王注意到升级扑克中有大小王，大小王的出现也是随机的。如果将随机整数Y的值放宽到01

10、3，则Z可以牙齿为100，200，300，400。定义Z=200时对应于国王，定义z=100时对应于国王，如上所示进行打印和分配处理。如果Z=300或400，则返回并重新生成x和y。设置通用算法和编程，24，(3)比较周期，生成I卡并存储在m数组中，I 1卡：for(j=1；J=10000j)x=rand()% 4 1；y=rand()；z=x * 100y；if(z=300 | | z=400)continue；t=0；for(k=0)；k=I-1；k)if(z=MK)t=1；布列克；if(t=0)mi=z；布列克；公共算法和编程，25，(4)打印输出打印直接应用C语言中ASCII代码16的

11、字符显示大小王和每个图案。设定文字阵列D以列印点数时，将Y=1，13，12，11分别转换为A，K，Q，J。为了实现真正的随机性，请根据时间设置t=time () 000。Srand(t)初始化随机数生成器以实现实际随机目的。，经常使用的算法和编程，26，7.4操作过程模拟，例7.10法国数学家Poisson有12品脱(容量单位)酒，5品脱和8品脱的空杯子各有一瓶。两个空杯子牙齿，怎么能平分12品脱牙齿的酒？我们要解决一般的二等奔忙问题。分别使用bv和cv(单位整数)的两个空杯子，以最小倒计时平分总容量为偶数A的酒。其中正整数bv、cv和偶数A都是从键盘输入的。7.4.2泊松忙碌，公共算法和编程

12、，27，设置i=a/2(i是全局变量)，拆分过程中在：瓶A中设置的酒量是A，(0B-C顺序操作是当B杯空(b=0)时在B杯中至少设置一次装满Bcv-c，c杯子，操作b=cv-c，清空b杯子，c杯子满(c=cv)后，从c杯子倒回a瓶子。公共算法和编程，28，b=0和ACV-c)B-=(cv-c)；C=cv/*在b中填充c罩杯*/else c=b；b=0；/*从b倒c，清空b杯*/printf(mmmn，a，b，c)。公共算法和编程，29，(2)按A-C-B顺序工作的牙齿循环操作，基本上等于(1)通过交换C和B杯返回函数值探测(A，cv，bv)。同时，设计实现函数Practice(a，bv，cv)，并与实验函数进行比较，将N增量1操作更改为输出中间过