数学与程序设计..docVIP

数学与程序设计 数学是研究现实世界的空间形式和数量关系的科学,是处理客观问题的强有力的工具,几乎在一切自然科学领域中都起着基础性的作用。数学的特点不仅在于概念的抽象性、逻辑的严密性、结论的明确性，还在于它应用的广泛性。 数学方法在程序设计中的运用包括两个方面：化简题目和直接解决问题。 应用数学方法化简题目是解决问题必不可少的重要步骤，也是分析题目的基本方法。应用数学方法化简题目，发掘题目中的隐含条件，寻求更多的“已知”条件，从而为建立数学模型提供依据。而用数学方法直接解题，其效率更是一般算法所不可及的。下面以具体的问题为例，介绍有关的数学知识和数学方法，体会利用数学方法解决问题时的乐趣。 一．数论基础 1．欧几里德a，b的最大公约数function gcd(a,b:longint):longint;?begin?if b=0 then gcdd:=a?else gcd:=gcd(b,a mod b);?end; 思考：如何把上述算法写成迭代形式？ 2．扩展的欧几里德算法gcd(a,b)=d，一定存在整数x和y满足gcd(a,b)=ax+by满足gcd(a,b)=ax+by的整数x，y）的算法如下： function exgcd(a,b:longint;var x,y:longint):longint;vart:longint;beginif b=0 then??begin? exgcd:=a;? x:=1;? y:=0;?endelse?begin?exgcd:=exgcd(b,a mod b,x,y);?t:=x;?x:=y;?y:=t-(a div b)*y;?end;end; 注释1 思考：1）如何把上述算法写成迭代形式？2）满足gcd(a,b)=ax+by的整数x，y）是否是唯一的？ 3. 求解二元一次不定方程ax+by=c整数解 我们的任务是解二元一次不定方程 ax+by=c ① 其中a,b,c都是整数，所求的解(x,y)也是整数 关于方程①的可解性，有下面的两个重要的结论： （1）设gcd(a,b)表示整数a,b的最大公约数。方程①有解的充分必要条件是gcd(a,b)|c。（记号“x|y”表示x能整除y，即存在整数k，使y=kx）。 （2）如果(x0,y0)是方程①的一组解，则对任何整数t,(x0+bt,y0-at)也都是方程①的解。 下面我们讨论具体求解的方法。 为了避免计算中对负数和0的讨论，我们假定a0，b0，并且a=b。 假定方程①有解，即系数满足：gcd(a,b)|c，这时，c’=c/gcd(a,b)一定是个整数。我们先讨论下面的方程： ax+by= gcd(a,b) ② 根据上述扩展的欧几里德算法整数x和y满足ax+by =gcd(a,b)(x0,y0)是方程②的一组解，则(c’x0,c’y0)也是方程①的一组解，即 a(c’x0)+b(c’y0)=(c’f)=c。 下面给出求二元一次不定方程ax+by=c一组整数解（x0,y0）的算法： procedure equation(a,b,c:longint;var x0,y0:longint);var d,x,y:longint;begin?d:=exgcd(a,b,x,y);{参见扩展的欧几里德算法?if c mod d0 then?begin? writeln(no answer);? halt;?end else?begin? x0:=x*(c div d);? y0:=y*(c div d);?end;end; 说明： （1）如果a,b中有一个小于0，例如a0，可以令x’=-x，解方程 ax’+by=c。 求解后，再令x=-x’就可以了。 （2）利用前面讲过的性质：“如果(x0,y0)是方程①的一组解，则对任何整数t,(x0+bt,y0-at)也都是方程①的解”。可以通过任何整数t得到方程①的其余解。 方程ax+by=c整数解的应用 例1：有三个分别装有a升水、b升水和c升水的量筒(cba0),现c筒装满水，问能否在c筒量出d升水(cd0)。若能，请列出一种方案。算法分析：?量水过程实际上就是倒来倒去，每次倒的时候总有如下几个持点：?1.总有一个筒中的水没有变动；?2．不是一个筒被倒满就是另一个筒被倒光；?3．c筒仅起中转作用，而本身容积除了必须足够装下a筒和b筒的全部水外，别无其它限制。 a筒b筒倒来倒去量出d升水升量出d升水能否在c筒量出d升水方程ax+by=整数解程序如下: program mw;typenode=array[0..1] of longint;vara,b,c:no