最长不下降子序列.docVIP

求最长不下降序列 一．问题描述 设有由n个不相同的整数组成的数列，记为: a(1)、a(2)、……、a(n)且a(i)a(j) (ij) 例如3，18，7，14，10，12，23，41，16，24。 若存在i1i2i3 … ie 且有a(i1)a(i2) … a(ie)则称为长度为e的不下降序列。如上例中3，18，23，24就是一个长度为4的不下降序列，同时也有3，7，10，12，16，24长度为6的不下降序列。程序要求，当原数列给出之后，求出最长的不下降序列。 二．算法分析 （一）。根据动态规划的原理，由后往前进行搜索。 1·对a(n)来说，由于它是最后一个数，所以当从a(n)开始查找时，只存在长度为1的不下降序列； 2·若从a(n-1)开始查找，则存在下面的两种可能性： ①若a(n-1)a(n)则存在长度为2的不下降序列a(n-1)，a(n)。 ②若a(n-1)a(n)则存在长度为1的不下降序列a(n-1)或a(n)。 3·一般若从a(i)开始，此时最长不下降序列应该按下列方法求出: 在a(i+1),a(i+2),…,a(n)中，找出一个比a(i)大的且最长的不下降序列，作为它的后继。倒推公式为 F（I）=MAX{F（I+K）}+1 F[I]：表示以第I个位置为起点的最长不下降序列的长度。 K的选择范围：a(I+k)a(i) I+k≦n 最后从F[1]到F[N]中选取最大的即为最优解 当然也可以采用顺推的方法，顺推公式为 F（I）=MAX{F（I-K）}+1 F[I]：表示以第I个位置为终点的最长不下降序列的长度。 K的选择范围：a(I-k)a(i) I-k≧1 最后从F[1]到F[N]中选取最大的即为最优解 4·为算法上的需要，定义一个数组（倒推法） 整数类型二维数组d（N，3） d（I，1）表示点a（i） d（I，2）表示从I位置到达N的最长不下降序列长度 d（I，3）表示从I位置开始最长不下降序列的下一个数字的位置，以便打印。 初始化: FOR I = 1 TO N INPUT #1, D(I, 1) D(I, 2) = 1 D(I, 3) = 0 NEXT I 下面给出求最长不下降序列的算法: FOR I = N - 1 TO 1 STEP -1 L = 0: P = 0 FOR J = I + 1 TO N IF D(I, 1) D(J, 1) AND D(J, 2) L THEN L = D(J, 2) P = J END IF NEXT J IF L 0 THEN D(I, 2) = L + 1 D(I, 3) = P END IF NEXT I 下面找出最长不下降序列： DMAX = D(1, 2) L = 1 FOR I = 2 TO N - 1 IF D(I, 2) DMAX THEN DMAX = D(I, 2) L = I {记录最长不下降序列的起点位置} END IF NEXT I 最长不下降序列长度为D(L, 2) 序列 WHILE L 0 PRINT D(L, 1); L = D(L, 3) WEND （3）程序清单 PROGRAM MAX_RISE(INPUT,OUTPUT); CONST MAXN=100;FNAME=Q21.TXT; TYPE TLIST=ARRAY[1..MAXN,1..3] OF INTEGER; VAR D:TLIST;N:BYTE; PROCEDURE INIT; VAR I:INTEGER; F:TEXT; BEGIN ASSIGN(F,FNAME); RESET(F); READLN(F,N); FOR I:=1 TO N DO BEGIN READ(F,D[I,1]); D[I,2]:=1; D[I,3]:=0 END; CLOSE(F); END; PROCEDURE MAKE; VAR I,J,P:BYTE; L:INTEGER; BEGIN FOR I:=N-1 DOWNTO 1 DO BEGIN L:=0;P:=0; FOR J:=I+1 TO N DO IF (D[