鱼C论坛

 找回密码
 立即注册
查看: 10176|回复: 58

[技术交流] 最优?步步贪婪,误入歧途!

[复制链接]
发表于 2014-12-8 15:05:47 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
看到一个算法问题,说的是:


有一个m*n的矩阵格子,里面每个格子有一个正数,一个人从左上要走到右下,每次只能往右或者下走,每经过一个格子他就拿走其分数,拿了格子里面的分数就没有了(为0),求如何走让分数最大。以下是难度:这个人这么走了2次,求怎么走让2次的总和最大。


思路:
先来走一次的,这个比较好写。
三种方法因为只能走右或下,所以选择两个中最大的来走。
第一种也就是贪婪的写法。

乔巴.jpg
def tl(seq):
        count1,count2 = [],[]
        results = []
        def move(seq,a = [0,0],result = 0,remove = []):
                if a == [len(seq)-1,len(seq[0])-1]:
                        count1.append(result)
                        count2.append(remove)
                else:
                        if a[0] == len(seq)-1:
                                move(seq,a = [a[0],a[1]+1],result = result+seq[a[0]][a[1]] ,remove = remove + [a])
                        elif a[1] == len(seq[0])-1:
                                move(seq,a = [a[0]+1,a[1]],result = result+seq[a[0]][a[1]],remove = remove + [a])
                        else:
                                if seq[a[0]][a[1]+1] < seq[a[0]+1][a[1]]:
                                        move(seq,a = [a[0]+1,a[1]],result = result+seq[a[0]][a[1]],remove = remove + [a])
                                else:
                                        move(seq,a = [a[0],a[1]+1],result = result+seq[a[0]][a[1]],remove = remove + [a])
        move(seq)
这样做确保的是每一步都是最优,但有时候不会整体最优。
[
[0,3,4,5,6],
[1,2,3,4,5],
[5,6,7,8,1],
[3,5,6,7,0]
]
根据贪婪,走法是0-3-4-5-6-5-1-0,这样虽然确保了每一步的最优,但是无法保证了整体的最优,所以贪婪往往得到的是次优或次次,次次次- -。

这也就是只贪图眼前,而不够长远的做法,所以想解决问题还是不要贪婪!
贪婪.jpg

第二种方法,暴力随机数。
想拼人品就把次数弄少点,想全点就把次数弄多点。

此方法太过暴力,慎用!
想试试的看一下这个问题。当时我用的随机数,那效率,没的说~,低到爆- -。

http://bbs.fishc.com/forum.php?m ... d=54578&ctid=12

第三种方法,动态规划。
动态规划是比较好的解决最优问题的方法,
也就是每走一步将结果暂存再将结果的可能性都列出,之后求最大。(不知道理解的对不对。)

将贪婪的修改一下就可以用啦。
def dp(seq):
        count1,count2 = [],[]
        results = []
        def move(seq,a = [0,0],result = 0,remove = []):
                if a == [len(seq)-1,len(seq[0])-1]:
                        count1.append(result)
                        count2.append(remove)
                else:
                        if a[0] == len(seq)-1:
                                move(seq,a = [a[0],a[1]+1],result = result+seq[a[0]][a[1]] ,remove = remove + [a])
                        elif a[1] == len(seq[0])-1:
                                move(seq,a = [a[0]+1,a[1]],result = result+seq[a[0]][a[1]],remove = remove + [a])
                        else:
                                move(seq,a = [a[0],a[1]+1],result = result+seq[a[0]][a[1]],remove = remove + [a])
                                move(seq,a = [a[0]+1,a[1]],result = result+seq[a[0]][a[1]],remove = remove + [a])
        move(seq)
利用递归,最后返回最大值就好了。

当然也可以用迭代来完成,用迭代可能要嵌套多个循环,能理出头绪就试试吧~。

---------------------------------------凶残的化身----------------------暴力分割线--------------------------------------

回到问题,上面只是解决了一步,

P.S(我是问题:  有一个m*n的矩阵格子,里面每个格子有一个正数,一个人从左上要走到右下,每次只能往右或者下走,每经过一个格子他就拿走其分数,拿了格子里面的分数就没有了(为0),求如何走让分数最大。以下是难度:这个人这么走了2次,求怎么走让2次的总和最大。)

让第一次的最优了,那剩下的直接替换掉路径,再来一次最优那不就完成了!OK,说干就干。

加上
        a = max(count1)
        b = count2[count1.index(a)]
        for i in range(len(b)):
                seq[b[i][0]][b[i][1]] = 0
        results.append(a)
        count1,count2 = [],[]
        move(seq)
        results.append(max(count1))
栗子:
[
[0,3,0,2],
[0,3,0,0],
[0,3,0,0],
[0,0,0,4],
[0,0,0,4],
[0,3,0,0]]
结果:
[17, 3]

也就是20,方法是动态规划本身是求的最优,应该是最优解了吧~。那他是怎么得来的呢?
第一次 360截图20141208143147136.jpg ,那第二次要么取3,要么取2,所以取了3。这样应该就最大了吧。

但是还有一种可以全取到的走法。(不唯一。)
360截图20141208143712209.jpg


这样的结果是 22。显然最优。

现在问题来了:
20140819091757514.png

我要说的显然不是如何编写最优解让挖掘鸡获取最大的效率T~T。(中国山东找蓝翔!)

动态规划还不能最优?
[url=]代码[/url]
其实是犯了一个错误,动态规划是整体上解决,但是我们是分步来看,局部动态规划,整体上就成了贪婪。谨慎!

不过,我没写出来直接求解的答案,只会后面进行处理。

        for i in range(len(count1)):                    
                for c in range(len(count1[i:])):
                        temp = count1[i]+count1[c]
                        #将每一步的结果两两相加。
                        for j in range(len(count2[0])):
                                if count2[i][j] == count2[c][j]:
                                        temp -= seq[count2[i][j][0]][count2[i][j][1]]
                        #重复的地方减少一个
                        results.append(temp)
这样问题就解决了,整体上最优解,效率也变得有些慢。。

少鱼币的童鞋,
练思路的童鞋,
来捧场的童鞋,
只有鱼币奖励- -。



360截图20141208150256962.jpg

是时候展示你犀利的思路,风骚的代码了!。
(more and more fish money!)




游客,如果您要查看本帖隐藏内容请回复

玉:


游客,如果您要查看本帖隐藏内容请回复



评分

参与人数 1荣誉 +5 鱼币 +5 收起 理由
bevin + 5 + 5

查看全部评分

本帖被以下淘专辑推荐:

想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2014-12-8 15:20:19 | 显示全部楼层

回帖奖励 +5 鱼币

不明觉厉,顶一个!

评分

参与人数 1荣誉 +2 鱼币 +2 收起 理由
wei_Y + 2 + 2 - -,我写的有这么差。。

查看全部评分

想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 17:18:01 | 显示全部楼层

回帖奖励 +5 鱼币

很久没来鱼C了,看到了小甲鱼的微博,进来瞧一瞧,顶一顶!
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 17:53:01 | 显示全部楼层

回帖奖励 +5 鱼币

刚刚开始学还不太明白学习学习!
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 20:08:21 | 显示全部楼层

回帖奖励 +5 鱼币

卧槽 不当版主,隐藏的东西都看不了了
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2014-12-8 20:12:01 | 显示全部楼层
bevin 发表于 2014-12-8 20:08
卧槽 不当版主,隐藏的东西都看不了了

叫你不当,吃亏了吧。
话说我隐藏的东西一般没啥作用,就是听说隐藏回复的人多哈哈。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 20:34:26 | 显示全部楼层
wei_Y 发表于 2014-12-8 20:12
叫你不当,吃亏了吧。
话说我隐藏的东西一般没啥作用,就是听说隐藏回复的人多哈哈。

算法太久没碰基本忘得差不多了,寒假有空准备补补,到时可以交流交流 哈哈哈
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 20:35:35 | 显示全部楼层
wei_Y 发表于 2014-12-8 20:12
叫你不当,吃亏了吧。
话说我隐藏的东西一般没啥作用,就是听说隐藏回复的人多哈哈。

感觉主要还是深度和广度搜索用的比较多,没事多做做欧拉计划上的题,会有提高的~
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2014-12-8 20:42:48 | 显示全部楼层
bevin 发表于 2014-12-8 20:35
感觉主要还是深度和广度搜索用的比较多,没事多做做欧拉计划上的题,会有提高的~

我还没学到嘛,欧拉计划我感觉要不能直接写出来要不就看半天没思路,没思路的还是大部分!
我没事就看兜兜给我说的那个外国站- -,全是算法的,还是英文。。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 20:44:59 | 显示全部楼层
wei_Y 发表于 2014-12-8 20:42
我还没学到嘛,欧拉计划我感觉要不能直接写出来要不就看半天没思路,没思路的还是大部分!
我没事就看兜 ...

哈哈哈 哪个网站?  我最近太忙了,这几周要汇报三四篇论文
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2014-12-8 20:51:37 | 显示全部楼层
bevin 发表于 2014-12-8 20:44
哈哈哈 哪个网站?  我最近太忙了,这几周要汇报三四篇论文

http://www.checkio.org

我这里晚上打不开- -。

论文辛苦啊,等你不忙了得教教我,我这一个问题要搞好久好久。。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 21:00:48 | 显示全部楼层
wei_Y 发表于 2014-12-8 20:51
http://www.checkio.org

我这里晚上打不开- -。

嗯,我也在群里的,到时在qq上说
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-8 22:13:50 | 显示全部楼层

回帖奖励 +5 鱼币

看看砖和玉
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2014-12-9 10:53:47 | 显示全部楼层
bevin 发表于 2014-12-8 21:00
嗯,我也在群里的,到时在qq上说

嗯。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-10 16:32:27 | 显示全部楼层

回帖奖励 +5 鱼币

阔以哦{:1_1:}
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-10 19:13:40 | 显示全部楼层

回帖奖励 +5 鱼币

有点吃力。。。
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-10 19:50:34 | 显示全部楼层

回帖奖励 +5 鱼币

正在努力构思……
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-10 20:19:55 | 显示全部楼层

回帖奖励 +5 鱼币

非常支持!!!!!!!
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-10 21:49:33 | 显示全部楼层

回帖奖励 +5 鱼币

好厉害
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2014-12-10 22:18:10 | 显示全部楼层

回帖奖励 +5 鱼币

喵 我也恢复一个看看
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2024-11-23 21:29

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表