第九部分递归程序设计技术学习程序设计需要注意规律性的东西北京交通大学计算机与信息技术学院教师:林友芳本章内容递归与循环递归函数的执行过程递归函数效率北京交通大学计算机与信息技术学院教师:林友芳循环与递归循环程序用于描述需要重复进行计算高级语言里,也常见用递归来实现重复的计算。递归—recursion,recursivealgorithm函数或过程调用自身C语言允许递归,可以在函数内调用自身,常常使程序更简单清晰。北京交通大学计算机与信息技术学院教师:林友芳1.阶乘和乘幂例:定义计算整数阶乘的函数1×2×…×(n-1)×n本例中,乘的次数依赖于n,计算所需的次数定义时无法确定。这是一种典型循环情况计算“次数”依赖某些参数的值。北京交通大学计算机与信息技术学院教师:林友芳程序longfact1(longn){longfac,i;for(fac=1,i=1;i=n;i++)fac*=i;returnfac;}北京交通大学计算机与信息技术学院教师:林友芳阶乘函数的精确定义是一种递归定义的形式要解决规模为n的问题,要先解决规模为n-1的子问题,依此类推。如果高级语言允许递归定义函数,就可以直接翻译为程序。C允许递归定义在函数定义内直接或间接地调用被定义函数本身。nnnnn!()!1010北京交通大学计算机与信息技术学院教师:林友芳写成递归函数longfact(longn){returnn==0?1:n*fact(n-1);}longfact(longn){if(n=1)return1;returnn*fact(n-1);}北京交通大学计算机与信息技术学院教师:林友芳longfact(1){if(1=1)return1;return1*fact(1-1);}longfact(2){if(2=1)return1;return2*fact(2-1);}longfact(3){if(3=1)return1;return3*fact(3-1);}main(){printf(“%d”,fact(3));}蓝线:函数调用线路黄线:函数内部执行路线红线:函数执行结束返回主调函数的路线longfact(longn){if(n=1)return1;returnn*fact(n-1);}北京交通大学计算机与信息技术学院教师:林友芳递归与计算过程包含递归的程序产生的计算过程和性质更复杂,能完成很复杂的工作。递归调用只有一个调用表达式或语句,但是可能要许多步才能完成。实际参数的不同,会实际产生的递归调用次数(步数)也会有很大的不同。递归程序理解的理解比较困难递归的函数定义需要有条件语句去控制递归过程的最终结束直接给出结果的时候,递归结束;把对较复杂情况的计算归结为对更简单情况的计算,需要进行递归处理。北京交通大学计算机与信息技术学院教师:林友芳递归和循环基本运算、关系判断、条件表达式,加函数定义和递归定义构成了一个(理论上)“足够强的”的程序语言。循环程序可以改成递归实现递归程序也可以改成循环实现。2.Fibonacci序列计算与时间北京交通大学计算机与信息技术学院教师:林友芳定义Fibonacci(斐波那契)序列的递归定义F0=1,F1=1Fn=Fn-1+Fn-2(n1)1,1,2,3,5,8,13,21,34,65,99,…北京交通大学计算机与信息技术学院教师:林友芳用递归程序实现longfib(intn){returnn2?1:fib(n-1)+fib(n-2);}问题分析:这个程序好不好?一方面,很好!程序与数学定义的关系很清晰,正确性容易确认,定义易读易理解北京交通大学计算机与信息技术学院教师:林友芳例fib(5)调用过程fib(5)fib(4)fib(3)fib(3)fib(2)fib(2)fib(1)fib(1)fib(0)fib(1)fib(0)fib(2)fib(1)fib(1)fib(0)存在什么问题?北京交通大学计算机与信息技术学院教师:林友芳问题存在大量重复计算,参数越大重复计算越多。有关系吗?随着参数增大,计算中重复增长迅速,最快的微机上一分钟大约可以算出fib(45)参数加1,fib多用近一倍时间(指数增长)。最快的微机一小时算不出fib(55),算fib(100)要数万年计算需要时间,复杂计算需要很长时间。这是计算机的本质特征和弱点。说明它不是万能,有些事情“不能”做。北京交通大学计算机与信息技术学院教师:林友芳计算复杂度人们发现了许多实际问题,理论上说可用计算机解决(可写出计算它的程序),但对规模大的情况(“大的参数n”),人类永远等不到计算完成。这时能说问题解决了吗?计算中有一大类问题被称为的“难解问题”,其中有许多很实际的问题,如规划、调度、优化等。解决这些问题,需要理论和实际技术的研究。另外,对于许多问题的实用的有效算法,有极大的理论价值和实际价值。计算复杂性,难解问题,“P=NP?”问题。北京交通大学计算机与信息技术学院教师:林友芳阅读材料:NP问题AproblemisassignedtotheNP(nondeterministicpolynomialtime)classifitissolvableinpolynomialtimebyanondeterministicTuringmachine.AP-problem(whosesolutiontimeisboundedbyapolynomial)isalwaysalsoNP.IfaproblemisknowntobeNP,andasolutiontotheproblemissomehowknown,thendemonstratingthecorrectnessofthesolutioncanalwaysbereducedtoasingleP(polynomialtime)verification.IfPandNParenotequivalent,thenthesolutionofNP-problemsrequires(intheworstcase)anexhaustivesearch.Linearprogramming,longknowntobeNPandthoughtnottobeP,wasshowntobePbyL.Khachianin1979.ItisanimportantunsolvedproblemtodetermineifallapparentlyNPproblemsareactuallyP.AproblemissaidtobeNP-hardifanalgorithmforsolvingitcanbetranslatedintooneforsolvinganyotherNP-problem.ItismucheasiertoshowthataproblemisNPthantoshowthatitisNP-hard.AproblemwhichisbothNPandNP-hardiscalledanNP-completeproblem.北京交通大学计算机与信息技术学院教师:林友芳为计算过程计时统计程序或程序片段的计算时间有助于理解程序性质。许多语言或系统都提供了内部计时功能。有关函数在time.h,统计程序时间时程序头部应写#includetime.h在程序里计时,通常写表达式clock()/CLOCKS_PER_SEC得到从程序开始到表达式求值时所经历的秒数。北京交通大学计算机与信息技术学院教师:林友芳确定计算fib(45)所需要的时间的程序#includestdio.h#includetime.hlongfib(intn){returnn=1?1:fib(n-1)+fib(n-2);}intmain(){doublex;x=clock()/CLOCKS_PER_SEC;fib(45);x=clock()/CLOCKS_PER_SEC-x;printf(Timingfib(45):%f.\n,x);return0;}北京交通大学计算机与信息技术学院教师:林友芳Fibonacci数的迭代计算Fibonacci数的递推计算,易见1)f1和f2是12)知道fn-2和fn-1连续两个Fibonacci数,就可算出下一个fn递推计算方式逐个往后推,可用循环实现递推方案longfib1(intn){longf1=1,f2=1,f3,i;if(n=1)return1;for(f3=f2+f1,i=2;in;++i){f1=f2;f2=f3;f3=f1+f2;}returnf3;}做一次递推fnfn-1fn-2北京交通大学计算机与信息技术学院教师:林友芳程序分析for(f3=f2+f1,i=2;in;++i){f1=f2;f2=f3;f3=f1+f2;}循环结束时i等于n,这时c的值是fn。要得到此结论,可设法证明:每次判断i的值时f3正是fi。北京交通大学计算机与信息技术学院教师:林友芳归纳证明第一次判断时i的值是2,f3的值2,正是fi(且f1的值是fi-1,f2的值是fi-2)若某次判断时i值是k(小于n),循环体中的语句使f1变成fk-1,f2变成fk,f3变成fk+1。i值增1使我们又有f1为fi-2,f2变成fi-1,f3变成fi根据归纳法,每次判断i的值时f3正是fi。北京交通大学计算机与信息技术学院教师:林友芳如何保证循环的正确执行循环实现重复性计算,循环体可能执行多次。如何保证对各种数据都能正确完成计算?循环中变量不断变化。写循环要考虑变量间的关系,保证某些关系在循环中不变:循环的不变关系。写循环时最重要的就是想清循环中应维持变量间的什么关系才能保证循环结束时变量能处在所需状态。写完循环后应仔细检查是否满足要求。循环不变关系(循环不变量)是理解循环、写好循环的关键。北京交通大学计算机与信息技术学院教师:林友芳问题本例中用循环的函数比用递归定义的好吗?新函数在计算时间上有极大优越性。计算时间由循环次数确定。循环体执行次数大致为n。fib(100)只需约100次循环,几乎察觉不到所花费时间。新函数定义较复杂,有复杂的循环。要理解程序意义,确认函数对任何参数都算出Fibonacci值,需要借助“循环不变关系”的概念和细致分析。注意:这个例子并不是说明递归比循环的效率低。完全可以写出计算fib的同样高效的递归定义的函数北京交通大学计算机与信息技术学院教师:林友芳最大公约数求两个整数的最大公约数(greatestcommondivisor,GCD),写函数longgcd(long,long)解法1从某个数开始,逐个判断当前数是否能同时整除m和n,在这个过程中记录下能同时整除m和n的最大整数。需要用一个辅助变量k记录当前需要判断的数。用一个循环实现k顺序取值初值设为1每次判断完后增1直到k大于m和n中其中的一个为止记下循环过程中出现的新的m和n的公约数,作为新的最大公约数用变量d表示当前的最大公约数初值1(是公约数),遇到新的公约数(一定更大)时记入d北京交通大学计算机与信息技术学院教师:林友芳程序有了d及其初值,k可以从2开始循环。函数定义longgcd(longm,longn){longd=1,k=2;for(;k=m&&k=n;k++)if(m%k==0&&n%k==0)d=k;returnd;}参数互素时初值1会留下来,能保证正确北京交通大学计算机与信息技术学院教师:林友芳计算过程示例mnkk=m&&k=nm%k==0&&n%k==0d2082是12是是23是否24是是45是否46是否47是否48是否4904北京交通大学计算机与信息技术学院教师:林友芳特殊情况处理一些特殊情况需要处理1)m和n都为0需特殊处理。令函数返回值0;2)若m和n中