对"一维最大子数组和"问题的思考

时间:2023-03-09 05:13:04
对"一维最大子数组和"问题的思考

对"一维最大子数组和"问题的思考(homework-01)

 一维最大子数组和问题,即给定一个数组,在它所有的连续子数组的和中,求最大的那个和。“最大子数组和”是一个很好的IT面试考题,在《编程之美》一书中同时阐述了一维数组和二维数组的讨论。本篇博客将会更加细致的讨论一维部分。

  一. 最直观的 O(n3) 解法

  在课上看到这个问题,当然最直观的解决办法即穷举。我们通过对子数组的起点和重点进行二层循环,计算每一个子数组的和,取其最大值,这样当然能够解决。但是作为一种O(n3)的解法,显然是低效率的。

  二. 进一步思考——O(n2)解法

  为增加程序效率,在接触过的算法中,我们想到或许动态规划能够被运用。动态规划的思想是将前一步的结果存储起来并在后续的步骤中直接使用,减少重复性劳动。因为a[1,10]这样的子数组可以拆解为a[1,9] + a[10], 所以我们可以思考一下如何运用动态规划来增加效率。

  核心问题是,如何运用动态规划,在哪里用。因为是求和,动态规划就被用来存储已经求和的子数组和,供包含此子数组的子数组的求和来使用。那么,我们仅使用上述O(n3)算法的两层循环,求和过程就简化为线性的复杂度,因此此算法的时间复杂度能被降低为O(n2)。

  三. 惊艳!——O(n) 解法

  O(n)算法是我在博客上看到的。看到竟然有O(n)算法,不禁佩服起人类的智慧... 不过,尝试不看解答,我们先想想,是在哪一部分把数量级降低了1呢?要达到O(n)的时间复杂度,必须将循环减至1层。这是容易想到的。问题是:循环减到1层,又如何处理每一层、并且得到最大子数组和呢?下面我们来看看方法:

  对于一个数组a[1]...a[L-1],假设我们已经知道其最大子数组和为max_former。当数组增加一个元素a[L]在结尾,那么新数组和最大的子数组要么包含a[L+1]在末尾,要么没有a[L+1],仍保持原先的子数组不变。用max_updated表示以a[L+1]结尾的最大子数组的和,则增加1个元素后新数组的max_former新=max(max_updated, max_former旧).

  在max_former旧 已知情况下(动态规划存储的中间步骤)下面只需求max_updated。max_updated旧 表示以a[L-1]结束的最大子数组和。以显然增加一个元素后,max_updated新=max(max_updated旧+a[L+1], a[L+1])。

  上面两段的思考,我们得到:

  

max_updated = max(max_updated + a[L], a[L])
max_former = max(max_former, max_updated)

  四. O(n)代码示例及对比思考

   

对"一维最大子数组和"问题的思考
 1 /*
2 * Author: Shone JIN, 11061128
3 * Sept. 19, 2013. Version 1.1
4 */
5
6 #include <fstream>
7 #include <iostream>
8 using namespace std;
9
10 int max(int a, int b)
11 {
12 return a > b ? a : b;
13 }
14
15 int min(int *a, int l)
16 {
17 int tmp = a[0];
18 for (; l > 0; l--){
19 if (a[l] < tmp){
20 tmp = a[l];
21 }
22 }
    return tmp;
23 }
24 int sum_max(int * a, int l)
25 {
26 int max_former, max_updated;
27 max_former = min(a, l);
28 max_updated = max_former;
29
30 for (int i = 0; i <= l; i++){
31 max_updated = max(max_updated + a[i], a[i]);
32 max_former = max(max_former, max_updated);
33 }
34 return max_former;
35 }
36
37 int main()
38 {
39 int array[100];
40 int tmp;
41 int i = 0;
42
43 ifstream infile("input.txt", ios::in);
44 if(! infile){
45 cerr << "File Error" << endl;
46 exit(1);
47 }
48
49 while(! infile.eof()){
50 infile >> array[i++];
51 }
52
53 cout << sum_max(array, i - 1) << endl;
54 return 0;
55 }
对"一维最大子数组和"问题的思考

  前面讲到的三个算法,可以理解为循序渐进的优化过程。O(n3)最为直观,当我们尝试使用动态规划存储中间步骤求得的和时,我们将O(n3) 简化为(n2)。至于O(n)算法,就是在动态规划存储上一步结果的前提下,只对子数组末尾元素进行一层循环。这里面运用到了一个巧妙的思维避开了O(n2)算法的两层循环。

  四. a.程序架构和思路

    max()返回两个数的最大值,min()返回一个数组的最小值。整体的程序在从文件中读取数组后,将数组及长度传给函数sum_max()进行求最大子数组和。sum_max()函数的思想就是前面提到的解法,即对子数组末尾元素循环,找每一个a[0]..a[i]这样的子数组的最大子数组和并求最大值。

    b.心得与体会

    我在课上第一次接触这个问题时陷入了瓶颈,总觉得有很好的解法,但是迟迟无法落笔。经过上面的思考,我领悟到好的算法不总是一下子灵光一闪得到的,而更多的是在不断的改进中得到的。我不需要一次写出最高效的算法,而更应该先把最直观的解法写下,然后一层一层的运用传统的算法进行优化,最终能优化多少,就是多少。

    c.作业中的时间消耗和开发效率分析

    作业的任务量不大,50行,难点在于寻找高效率的解法。在本作业的过程中我使用约1个小时的时间在思考最优秀算法,使用1小时的时间查阅别人的资料并编程,20分钟的测试和修改,30分钟进行Github学习,并签入程序。最后,使用2小时的时间撰写本篇博文。我按时完成了任务,自认为代码质量和很好(代码很短也谈不上...)。总之虽然编程已经有点生疏、慢,但是态度自己还是满意的。

    d.测试结果及程序运行的截图

    首先是数组输入文件,我们举下面的例子:

对"一维最大子数组和"问题的思考

     其正确解应为60:

对"一维最大子数组和"问题的思考

    其他的测试数据,可以参考我的Github内容。测试中发现一个问题,无法对全负数数组正确求解,原因是前一个版本的max_former, max_updated初始化为0了,在这里已经修改为是数组的最小元素值。其他并无问题。

    By the way, 受助教老师启发,上面说的max_former, max_updated不用寻找数组最小值也可以得出正确解,而且效率还比较高。我们直接付一个比较小的值就可以,不妨在第一次赋值的时候都赋成a[0]的值,便不对结果产生影响。

  五. 其他

    我的Github账号: shoneJIN

    我选择的参考书是第二版的《代码大全》

面试:输出循环小数的循环节

Posted on 2013-09-20 15:46 DM张朋飞 阅读(451) 评论(2) 编辑 收藏

一.引入

偶然在网上看见的,手痒了,就试了试,哈哈。。。。。。。做出来了,实际就是模拟。

下面的来自*:循环小数是从小数部分的某一位起,一个数字或几个数字,依次不断地重复出现的小数。

循环小数即为有理数的小数表示形式,例:

对"一维最大子数组和"问题的思考

定理:一个分母为N的循环小数的循环节位数最多不超过N-1位。

(1) 化为分数的方法

我先想到用高中的等比数列的和求极限方式求取(a1/(1-q))。

对"一维最大子数组和"问题的思考 
         (2)利用长除法可以将分数(有理数)转化为循环小数。 
                         对"一维最大子数组和"问题的思考

对"一维最大子数组和"问题的思考

二.Java实现

import java.util.Arrays;

public class Xiaoshuo {

	//使用char类型可以省内存也可以直接存储小数点
	static char[] ch = new char[20];

	public static void main(String[] args) {
		// TODO Auto-generated method stub
		Arrays.fill(ch,'#');
		int a = 5;
		int b = 4;
		int index = 0;
		//在函数外部做处理方便多了
		if(a>=b) {
			if(a%b==0) {
				System.out.println(a/b+".0");
				return ;
			}
			System.out.print(a/b);
			ch[0] = '.';
			index = 1;
			a = a%b;
		}else {
			ch[0] = '0';
			ch[1] = '.';
			index = 2;
		}
		solve(index,a,b);
	}

	private static void solve(int index, int a, int b) {
		// TODO Auto-generated method stub
		//判断余数(你也可以叫商,好好看除法过程)是否出现过
		int from = -1;
		int to = -1;
		while(from==-1) {
			a *= 10;
			//必须加上48
			ch[index++] = (char)(a/b + 48);
			//不能再找到自身,所以需要第二个参数
			from = tell((char)(a/b+48),index-2);
			a = a%b;
		}
		to = index -2;
		for(int i=0; i<from; i++) {
			System.out.print(ch[i]);
		}
		System.out.print("{");
		for(int i=from; i<=to; i++) {
			System.out.print(ch[i]);
		}
		System.out.println("}");

	}

	private static int tell(char mod, int xiabiao) {
		// TODO Auto-generated method stub
		//先找到小数点后一位
		char c = ch[0];
		//index必须从0开始,因为有可能第一个就是小数点
		int index = 0;
		while(c!='.') {
			index++;
			c = ch[index];
		}
		//index是小数点的位置,所以从下一位开始
		index++;
		//采用while(mod1=‘#’)循环不行,这样是最好的,既有下标,又到非法字符处结束(前提是数组比较大,填充了非法字符)
		for(int i=index; i<=xiabiao; i++) {
			if(mod==ch[i])
				return i;
		}
		return -1;
	}

}

三.感悟 
        在函数外部先做数据处理化成统一格式方便多了。。。。。。。