Re: [问题] Interview street: zombie march

楼主: Leon (Achilles)   2012-10-31 23:58:22
※ 引述《DJWS (...)》之铭言:
: Leon提供了一个证明,
: 证明题目给定的矩阵存在stationary probability distribution,而且很容易证明。
: 但是Leon完全没有提到,
: 这矩阵究竟是几步之后才会收敛,只说它难解。
因为我知道这个一讲下去没人知道我再说什么啊,
我想这个版学过 Random process 的人应该不多.
: 至于原问题是问:k步之后,究竟state会变成如何?
: Leon提到的这两个命题,都无法完全解决原问题。
我想这是重点.
因为时间关系, 我写文章的时候并没有把所有条件列出来.
我所提供的方法, 是在考虑困难的情况下 (简单的我根本不想看)
也就是在 M, N, K 很大的时候,
我怎么去提供一个近似解.
: 所以Leon并没有解决原本的问题,只是提供了原问题的其中一些性质。
: 各位就算跟Leon继续辩论下去,也得不到原问题的答案的。
: 至于Leon回文的第一句:“其实是.. 有的, 而且快到你难以想像.”
: 想必这句话是说错了,其实他并没有解决原问题。
在研究的路程上, 通常你只能得到 heurisitic solution,
并且你只能解一步状况下的问题..
我想, 公平的说, 我不用解 eignvalue, eign-vector
就能得到 stationary probability.
在这点的进展, 并且引起讨论这是我乐于见到的.
我倒是很期待有人能放上更好的解法.
你要不要试试看?
: 照Leon的文章内容来看,
: 我想原问题应该是满难解的,
: 原因是牵扯到 Markov Chain mixing time,
: 如果要继续讨论的话,从这关键字下手,大家应该会比较有交集。
: 另外我想请教的是,
: Leon推文说到这是图论问题,
: 那么目前市面上有没有哪一本图论书籍,有涵盖到这个主题的?
这个问题夯的很呢,
去 google Markov Chain mixing time, random walk, gossip algorithm
连 S. Boyd 都在这个领域发过文章.
这个人写的论文比较浅显,
也有例子, 你也可以在 page 15 以后找到他推论 mixing time 的公式
http://keithbriggs.info/documents/Min_Chen_MSc.pdf
教科书, 通常都是要晚个好几年....
作者: DJWS (...)   2011-01-01 00:11:00
我想推文的人都应该有看懂你在说什么 只是你可能中文不好一直没办法理解大家想讨论一个确切解 而非近似解至于你提供近似解也是一个方向 不过也应该一开始就讲明了免得大家搞不清楚你究竟想不想解决原问题
楼主: Leon (Achilles)   2011-01-01 00:58:00
我会写那些东西, 是因为有人说可以观察 eignvalue..Based on what you said on ``eignvalue'', I suppose youalready accept the assumption on Markov Chainanyway, I don't think it's good to discuss this topicsince not many people know about the concept
作者: DJWS (...)   2011-01-01 07:53:00
是我提出来的 我当然知道呀...求eigenvalue是NP-hard所以我才说这个方法参考看看就好还有你eigenvalue打错字了 XD
楼主: Leon (Achilles)   2011-01-01 08:22:00
你.. 是112的吗? eigenvalue 是 NP-hard ????我.. 竟然浪费了时间和这群人讨论....我除了这篇留给你看之外, 其他的我全砍了..
作者: singlovesong (~"~)   2011-01-01 09:25:00
DJ大口误吧XDD
作者: neutrino (十年一梦)   2011-01-01 10:57:00
"提到eigenvalue"是指我吗? 我听过(在你给的这份文件Ch4也有提到, 不过我还没时间细看这份)那个second largeeigenvalue与mixing rate的方法 才提起来的eigen decompose我只知道O(n^3), 所以我想请教mixingtime有什么好进展. 毕竟自己不是做随机这方面领域的,自然希望有更多关键字指引, 不然MC,mixing time的文章对不是钻这领域的人来说, 可太多了...
作者: DJWS (...)   2011-01-01 15:44:00
就因式分解是NP-hard级别的 然后求eigenvalue可以化作多项式因式分解 所以求出所有eigenvalue是NP-hard因为很难算 所以可以用内插法、牛顿法等等的近似算法求根所以我们可以用P-time求得一个根 目前所有的eigen-decomposition的算法,都是基于近似算法的这段是之前从别人blog看到的 如果我的理解有误请告诉我另外就是 我并不觉得跟你讨论是浪费时间 还是能学到新观念只是你讲的内容实在偏离原问题太多 使得大家没共识而已...抱歉...第一句话开头请改成“多项式因式分解是NP-hard级别”
作者: suhorng ( )   2011-01-01 19:24:00
可不可以题外话问一下XD 像这种求解出来可能是什么根号什么奇怪的实数的问题 我们说计算时间一般是指什么的的时间? 比如 (-b+-sqrt(b^2-4ac))/(2a), 是忽略那sqrt时间吗?
作者: DJWS (...)   2011-01-01 22:03:00
如果是说big-O notation的话 只留下成长最快的那一项就可以
作者: suhorng ( )   2011-01-01 22:05:00
噢, 我说得不太清楚, 这样问好了我们说计算这个的时间复杂度 是指计算近似值到某特定精度的算法的复杂度吗?ZY
作者: DJWS (...)   2011-01-01 22:06:00
然后学理上的计算时间 就是算法的步骤数目喔喔 我开始理解你想问什么了 XD
作者: suhorng ( )   2011-01-01 22:07:00
耶XD
作者: DJWS (...)   2011-01-01 22:09:00
时间复杂度和复杂度分类,"应该"是指计算确切值吧!至于近似算法也有自己的复杂度分类例如PTAS的就是指可以P-time得到近似解 误差小于一定范围
作者: suhorng ( )   2011-01-01 22:14:00
喔喔~感谢!
作者: DJWS (...)   2011-01-01 22:14:00
我不是很懂复杂度分析 所以有错还请板友指正 谢谢!
作者: suhorng ( )   2011-01-01 22:15:00
有的则是可以依要求 花不同时间作到近似到任意精度大概是这样?
作者: DJWS (...)   2011-01-01 22:15:00
所以把它分类成P问题 好像也是可以 (?)对阿 我的理解是这样!
作者: suhorng ( )   2011-01-01 22:17:00
谢谢:D
楼主: Leon (Achilles)   2011-01-02 00:45:00
你如果不懂的话还是少说点话.. 免得误导别人另外, 我认为你连complexity基础都没有, 讲什么也无用我建议你把你上面解释 ``eigenvalue is NP-hard'' 那段话写下来, 拿去给 CSIE 教算法的老师看看,试着说服他你上面那段话
作者: DJWS (...)   2011-01-02 08:55:00
我的确没有complexity基础 有机会我会请教专业人士如果你能讲解的话 那就更好了!另外你应该是学者吧 想请教你是从事哪个领域的研究知道领域范畴的话 讨论起来比较有交集 / 我自己并不是学者
楼主: Leon (Achilles)   2011-01-03 03:42:00
去修课比较快
作者: DJWS (...)   2011-01-03 09:25:00
那么楼上修过complexity的课吗?我想知道complexity有哪本教科书有提到eigenvalue的时间复杂度...我找满久的都找不到 orz
楼主: Leon (Achilles)   2011-01-03 13:00:00
Take class first, because you even don't know the def
作者: suhorng ( )   2011-01-03 15:32:00
定义不是查一查就知到了吗... DJWS大显然一定知道呀
作者: stimim (qqaa)   2011-01-03 15:35:00
我觉得 DJWS 说的应该是求 exact value ,而不是数值解也就是答案如果是 sqrt(2) ,就不能输出 1.414...
作者: suhorng ( )   2011-01-03 15:51:00
那可能是个满吊诡的地方. sqrt(2)可以看成一个算法, 输出 x^2 - 2 = 0 的正根到任意精度位而一般任意五次或以上方程式, 其解没有通用的方法用系数的加减乘除及某些次方根来表示所以是必要用其他的表示法 就看接不接受
楼主: Leon (Achilles)   2011-01-03 15:58:00
我认为DJWS写的东西, 与Matrix computation谈的complexity定义完全不同, 所以我建议他去修课, 看看别人怎么定义问题不然就把自己的定义写成文章投稿 解释`eigenvalue is NP-hard顺便给楼上两位: finite state machine 怎么表示 irrationalnumber? 有办法 exactly 表示 [0,1] 间的 irrational number?
作者: stimim (qqaa)   2011-01-03 16:09:00
我对电脑中的exact value的理解是要像mathematica那样,1/3 他就写 1/3 ,sqrt(2) 他就写 sqrt(2),pi 就是 pi
作者: DJWS (...)   2011-01-03 18:46:00
所以你也没修过complexity囉? 那么等我问到,我再教你。如果你有管道可以问的话 也请你帮忙问一下 谢谢!
作者: suhorng ( )   2011-01-03 19:41:00
话说, Mathematica 能对于任意正整数 n, 都给出x^5 + x - n = 0 的解的 exact value 吗?我手边只有 Maple, 可是不知道怎么下指令..|||
作者: FRAXIS (喔喔)   2011-01-03 20:30:00
我觉得现在好像是在讨论eigenvalue是不是computable numberhttp://en.wikipedia.org/wiki/Computable_real
作者: suhorng ( )   2011-01-03 20:34:00
是好像有这种味道XDDD
楼主: Leon (Achilles)   2011-01-04 01:12:00
我认为 DJWS 写的东西还是不要看比较好, 对知识长进有害你先去念 CLRS 的 algorithm, 再去看 Golomb 的 matrixComputation, 另外, 敢上来呛人就得先学会 google, 不是嘛?
作者: DJWS (...)   2011-01-04 10:35:00
那我大概有个底了 我想你应该是做讯号方面的!CLRS我整本都读过了 至于矩阵运算我真的不熟还有就是关于eigenvalue的问题 我昨天努力google了很久多项式因式分解是P 包括有理数/代数扩展 http://ppt.cc/vnOA特征多项式是GapL http://ppt.cc/xBLC特征多项式事实上也有P算法 http://ppt.cc/T_E3求eigenvalue其实就是对特征多项式作因式分解 故推论为P至于求eigenvector、求eigenbasis的部分我还没找到资料所以我的推文确实推错了 求eigenvalue是P才对! 不好意思...另外想问 Golomb 的 matrix 是哪一本书? 因为找不到书名里有 martix 的...
作者: FRAXIS (喔喔)   2011-01-04 20:46:00
Eigenvalue的复杂度可以参考这篇论文The complexity of the matrix eigenproblem, STOC 1999至于书的话 我猜他是在说Matrix Computations, by Golub
作者: DJWS (...)   2011-01-04 22:39:00
谢谢楼上! 有空来去图书馆翻一翻

Links booklink

Contact Us: admin [ a t ] ucptt.com