冠群's profile成长的见证PhotosBlogLists Tools Help

成长的见证

冠群 张

Occupation
Location
Photo 1 of 10
November 21

Verapamil

一点都不乖。。。让我想想。。。
November 16

11/15/09——无

本周无所事事,本篇无事可记,本人无可救药。
November 13

被迫消极怠工

Maximum number of users for Signal_Toolbox reached.
Try again later.

Maximum number of users for Statistics_Toolbox reached.
Try again later.
November 08

11/08/09——娘死人不偿命

不知不觉的,又到了该写日志的时候。这个星期写些什么呢?让我想一想啊。。。想来又想去,想到了CANOEE上周来访密歇根。密歇根人民可开心了,很久没有这么开心过了。后来,他回去了,一切又都恢复了平静,就跟以前一样一样的。在CANOEE临走前,给我看了一个人写的几篇文章,文风很是新颖,看得我很是受启发,觉得原来日志也可以这样写啊,要是能早点见识到那该多好啊。于是呢,这个周末,天气这样好,我就仿效此文风,写下了上述的一段话。

======================华丽的分割线===========================
皑皑,以后要一直这么写东西,从此也不用再世上混了。。。
所以说,发PAPER一定要认真严谨准确,一不小心有失水准就别想翻身了。

此外,人人都有把代码弄乱的潜质,人人也都有把文章写恶心的潜质。
关于如何弄乱代码,请移步:

6个变态的C语言Hello World程序: http://coolshell.cn/?p=914

如何加密/混乱C源代码: http://coolshell.cn/?p=933

大家领会精神,适可而止。

November 05

Talking about [转] P问题、NP问题、NPC问题、NP-Hard问题

 

Quote

[转] P问题、NP问题、NPC问题、NP-Hard问题
这或许是众多OIer最大的误区之一。

    你会经常看到网上出现“这怎么做,这不是NP问题吗”、“这个只有搜了,这已经被证明是NP问题了”之类的话。你要知道,大多数人此时所说的NP问题其实 都是指的NPC问题。他们没有搞清楚NP问题和NPC问题的概念。NP问题并不是那种“只有搜才行”的问题,NPC问题才是。好,行了,基本上这个误解已 经被澄清了。下面的内容都是在讲什么是P问题,什么是NP问题,什么是NPC问题,你如果不是很感兴趣就可以不看了。接下来你可以看到,把NP问题当成是 NPC问题是一个多大的错误。

    还是先用几句话简单说明一下时间复杂度。时间复杂度并不是表示一个程序解决问题需要花多少时间,而是当问题规模扩大后,程序需要的时间长度增长得有多快。 也就是说,对于高速处理数据的计算机来说,处理某一个特定数据的效率不能衡量一个程序的好坏,而应该看当这个数据的规模变大到数百倍后,程序运行时间是否 还是一样,或者也跟着慢了数百倍,或者变慢了数万倍。不管数据有多大,程序处理花的时间始终是那么多的,我们就说这个程序很好,具有O(1)的时间复杂 度,也称常数级复杂度;数据规模变得有多大,花的时间也跟着变得有多长,这个程序的时间复杂度就是O(n),比如找n个数中的最大值;而像冒泡排序、插入 排序等,数据扩大2倍,时间变慢4倍的,属于O(n^2)的复杂度。还有一些穷举类的算法,所需时间长度成几何阶数上涨,这就是O(a^n)的指数级复杂 度,甚至O(n!)的阶乘级复杂度。不会存在O(2*n^2)的复杂度,因为前面的那个“2”是系数,根本不会影响到整个程序的时间增长。同样地,O (n^3+n^2)的复杂度也就是O(n^3)的复杂度。因此,我们会说,一个O(0.01*n^3)的程序的效率比O(100*n^2)的效率低,尽管 在n很小的时候,前者优于后者,但后者时间随数据规模增长得慢,最终O(n^3)的复杂度将远远超过O(n^2)。我们也说,O(n^100)的复杂度小 于O(1.01^n)的复杂度。

    容易看出,前面的几类复杂度被分为两种级别,其中后者的复杂度无论如何都远远大于前者:一种是O(1),O(log(n)),O(n^a)等,我们把它叫 做多项式级的复杂度,因为它的规模n出现在底数的位置;另一种是O(a^n)和O(n!)型复杂度,它是非多项式级的,其复杂度计算机往往不能承受。当我 们在解决一个问题时,我们选择的算法通常都需要是多项式级的复杂度,非多项式级的复杂度需要的时间太多,往往会超时,除非是数据规模非常小。

    自然地,人们会想到一个问题:会不会所有的问题都可以找到复杂度为多项式级的算法呢?很遗憾,答案是否定的。有些问题甚至根本不可能找到一个正确的算法 来,这称之为“不可解问题”(Undecidable Decision Problem)。The Halting Problem就是一个著名的不可解问题,在我的MSN Space上有过专门的介绍和证明。再比如,输出从1到n这n个数的全排列。不管你用什么方法,你的复杂度都是阶乘级,因为你总得用阶乘级的时间打印出结 果来。有人说,这样的“问题”不是一个“正规”的问题,正规的问题是让程序解决一个问题,输出一个“YES”或“NO”(这被称为判定性问题),或者一个 什么什么的最优值(这被称为最优化问题)。那么,根据这个定义,我也能举出一个不大可能会有多项式级算法的问题来:Hamilton回路。问题是这样的: 给你一个图,问你能否找到一条经过每个顶点一次且恰好一次(不遗漏也不重复)最后又走回来的路(满足这个条件的路径叫做Hamilton回路)。这个问题 现在还没有找到多项式级的算法。事实上,这个问题就是我们后面要说的NPC问题。

    下面引入P类问题的概念:如果一个问题可以找到一个能在多项式的时间里解决它的算法,那么这个问题就属于P问题。P是英文单词多项式的第一个字母。哪些问 题是P类问题呢?通常NOI和NOIP不会出不属于P类问题的题目。我们常见到的一些信息奥赛的题目都是P问题。道理很简单,一个用穷举换来的非多项式级 时间的超时程序不会涵盖任何有价值的算法。

    接下来引入NP问题的概念。这个就有点难理解了,或者说容易理解错误。在这里强调(回到我竭力想澄清的误区上),NP问题不是非P类问题。NP问题是指可以在多项式的时间里验证一个解的问题(但未能找到在多项式的时间内求解的算法,即“猜得到算不出”)。NP问题的另一个定义是,可以在多项式的时间里猜出一个解的问题。比方说,我RP很好,在程序中需要枚举时,我可以 一猜一个准。现在某人拿到了一个求最短路径的问题,问从起点到终点是否有一条小于100个单位长度的路线。它根据数据画好了图,但怎么也算不出来,于是来 问我:你看怎么选条路走得最少?我说,我RP很好,肯定能随便给你指条很短的路出来。然后我就胡乱画了几条线,说就这条吧。那人按我指的这条把权值加起来 一看,嘿,神了,路径长度98,比100小。于是答案出来了,存在比100小的路径。别人会问他这题怎么做出来的,他就可以说,因为我找到了一个比100 小的解。在这个题中,找一个解很困难,但验证一个解很容易。验证一个解只需要O(n)的时间复杂度,也就是说我可以花O(n)的时间把我猜的路径的长度加 出来。那么,只要我RP好,猜得准,我一定能在多项式的时间里解决这个问题。我猜到的方案总是最优的,不满足题意的方案也不会来骗我去选它。这就是NP问 题。当然有不是NP问题的问题,即你猜到了解但是没用,因为你不能在多项式的时间里去验证它(即“无解”)。下面我要举的例子是一个经典的例子,它指出了一个目前还没有 办法在多项式的时间里验证一个解的问题。很显然,前面所说的Hamilton回路是NP问题,因为验证一条路是否恰好经过了每一个顶点非常容易。但我要把 问题换成这样:试问一个图中是否不存在Hamilton回路。这样问题就没法在多项式的时间里进行验证了,因为除非你试过所有的路,否则你不敢断定它“没 有Hamilton回路”。

    之所以要定义NP问题,是因为通常只有NP问题才可能找到多项式的算法。我们不会指望一个连多项式地验证一个解都不行的问题存在一个解决它的多项式级的算 法。相信读者很快明白,信息学中的号称最困难的问题——“NP问题”,实际上是在探讨NP问题与P类问题的关系。

    很显然,所有的P类问题都是NP问题。也就是说,能多项式地解决一个问题,必然能多项式地验证一个问题的解——既然正解都出来了,验证任意给定的解也只需 要比较一下就可以了。关键是,人们想知道,是否所有的NP问题都是P类问题。我们可以再用集合的观点来说明。如果把所有P类问题归为一个集合P中,把所有 NP问题划进另一个集合NP中,那么,显然有P属于NP。现在,所有对NP问题的研究都集中在一个问题上,即究竟是否有P=NP?通常所谓的“NP问 题”,其实就一句话:证明或推翻P=NP。

    NP问题一直都是信息学的巅峰。巅峰,意即很引人注目但难以解决。在信息学研究中,这是一个耗费了很多时间和精力也没有解决的终极问题,好比物理学中的大统一和数学中的歌德巴赫猜想等。

    目前为止这个问题还“啃不动”。但是,一个总的趋势、一个大方向是有的。人们普遍认为,P=NP不成立,也就是说,多数人相信,存在至少一个不可能有多项 式级复杂度的算法的NP问题。人们如此坚信P≠NP是有原因的,就是在研究NP问题的过程中找出了一类非常特殊的NP问题叫做NP-完全问题,也即所谓的 NPC问题。C是英文单词“完全”的第一个字母。正是NPC问题的存在,使人们相信P≠NP。下文将花大量篇幅介绍NPC问题,你从中可以体会到NPC问 题使P=NP变得多么不可思议。

    为了说明NPC问题,我们先引入一个概念——约化(Reducibility,有的资料上叫“归约”)。

    简单地说,一个问题A可以约化为问题B的含义即是,可以用问题B的解法解决问题A,或者说,问题A可以“变成”问题B。《算法导论》上举了这么一个例子。 比如说,现在有两个问题:求解一个一元一次方程和求解一个一元二次方程。那么我们说,前者可以约化为后者,意即知道如何解一个一元二次方程那么一定能解出 一元一次方程。我们可以写出两个程序分别对应两个问题,那么我们能找到一个“规则”,按照这个规则把解一元一次方程程序的输入数据变一下,用在解一元二次 方程的程序上,两个程序总能得到一样的结果。这个规则即是:两个方程的对应项系数不变,一元二次方程的二次项系数为0。按照这个规则把前一个问题转换成后 一个问题,两个问题就等价了。同样地,我们可以说,Hamilton回路可以约化为TSP问题(Travelling Salesman Problem,旅行商问题):在Hamilton回路问题中,两点相连即这两点距离为0,两点不直接相连则令其距离为1,于是问题转化为在TSP问题 中,是否存在一条长为0的路径。Hamilton回路存在当且仅当TSP问题中存在长为0的回路。

    “问题A可约化为问题B”有一个重要的直观意义:B的时间复杂度高于或者等于A的时间复杂度。也就是说,问题A不比问题B难。这很容易理解。既然问题A能 用问题B来解决,倘若B的时间复杂度比A的时间复杂度还低了,那A的算法就可以改进为B的算法,两者的时间复杂度还是相同。正如解一元二次方程比解一元一 次方程难,因为解决前者的方法可以用来解决后者。

    很显然,约化具有一项重要的性质:约化具有传递性。如果问题A可约化为问题B,问题B可约化为问题C,则问题A一定可约化为问题C。这个道理非常简单,就不必阐述了。

    现在再来说一下约化的标准概念就不难理解了:如果能找到这样一个变化法则,对任意一个程序A的输入,都能按这个法则变换成程序B的输入,使两程序的输出相同,那么我们说,问题A可约化为问题B。

    当然,我们所说的“可约化”是指的可“多项式地”约化(Polynomial-time Reducible),即变换输入的方法是能在多项式的时间里完成的。约化的过程只有用多项式的时间完成才有意义。

    好了,从约化的定义中我们看到,一个问题约化为另一个问题,时间复杂度增加了,问题的应用范围也增大了。通过对某些问题的不断约化,我们能够不断寻找复杂 度更高,但应用范围更广的算法来代替复杂度虽然低,但只能用于很小的一类问题的算法。再回想前面讲的P和NP问题,联想起约化的传递性,自然地,我们会想 问,如果不断地约化上去,不断找到能“通吃”若干小NP问题的一个稍复杂的大NP问题,那么最后是否有可能找到一个时间复杂度最高,并且能“通吃”所有的 NP问题的这样一个超级NP问题?答案居然是肯定的。也就是说,存在这样一个NP问题,所有的NP问题都可以约化成它。换句话说,只要解决了这个问题,那 么所有的NP问题都解决了。这种问题的存在难以置信,并且更加不可思议的是,这种问题不只一个,它有很多个,它是一类问题。这一类问题就是传说中的NPC 问题,也就是NP-完全问题。NPC问题的出现使整个NP问题的研究得到了飞跃式的发展。我们有理由相信,NPC问题是最复杂的问题。再次回到全文开头, 我们可以看到,人们想表达一个问题不存在多项式的高效算法时应该说它“属于NPC问题”。此时,我的目的终于达到了,我已经把NP问题和NPC问题区别开 了。到此为止,本文已经写了近5000字了,我佩服你还能看到这里来,同时也佩服一下自己能写到这里来。

    NPC问题的定义非常简单。同时满足下面两个条件的问题就是NPC问题。首先,它得是一个NP问题;然后,所有的NP问题都可以约化到它。证明一个问题是 NPC问题也很简单。先证明它至少是一个NP问题,再证明其中一个已知的NPC问题能约化到它(由约化的传递性,则NPC问题定义的第二条也得以满足;至 于第一个NPC问题是怎么来的,下文将介绍),这样就可以说它是NPC问题了。

    既然所有的NP问题都能约化成NPC问题,那么只要任意一个NPC问题找到了一个多项式的算法,那么所有的NP问题都能用这个算法解决了,NP也就等于P 了。因此,给NPC找一个多项式算法太不可思议了。因此,前文才说,“正是NPC问题的存在,使人们相信P≠NP”。我们可以就此直观地理解,NPC问题 目前没有多项式的有效算法,只能用指数级甚至阶乘级复杂度的搜索。

    顺便讲一下NP-Hard问题。NP-Hard问题是这样一种问题,它满足NPC问题定义的第二条但不一定要满足第一条(就是说,NP-Hard问题要比 NPC问题的范围广)。NP-Hard问题同样难以找到多项式的算法,但它不列入我们的研究范围,因为它不一定是NP问题。即使NPC问题发现了多项式级 的算法,NP-Hard问题有可能仍然无法得到多项式级的算法。事实上,由于NP-Hard放宽了限定条件,它将有可能比所有的NPC问题的时间复杂度更 高从而更难以解决。

    不要以为NPC问题是一纸空谈。NPC问题是存在的。确实有这么一个非常具体的问题属于NPC问题。下文即将介绍它。

    下文即将介绍逻辑电路问题。这是第一个NPC问题。其它的NPC问题都是由这个问题约化而来的。因此,逻辑电路问题是NPC类问题的“鼻祖”。

    逻辑电路问题是指的这样一个问题:给定一个逻辑电路,问是否存在一种输入使输出为True。

    什么叫做逻辑电路呢?一个逻辑电路由若干个输入,一个输出,若干“逻辑门”和密密麻麻的线组成。看下面一例,不需要解释你马上就明白了。
  ┌───┐
  │ 输入1├─→┐    ┌──┐
  └───┘    └─→┤    │
                     │ OR ├→─┐
  ┌───┐    ┌─→┤    │    │    ┌──┐
  │ 输入2├─→┤    └──┘    └─→┤    │
  └───┘    │               ┌─→┤AND ├──→输出
               └────────┘ ┌→┤    │
  ┌───┐    ┌──┐            │  └──┘
  │ 输入3├─→┤ NOT├─→────-┘
  └───┘    └──┘
    这是个较简单的逻辑电路,当输入1、输入2、输入3分别为True、True、False或False、True、False时,输出为True。

    有输出无论如何都不可能为True的逻辑电路吗?有。下面就是一个简单的例子。
  ┌───┐
  │输入1 ├→─┐    ┌──┐
  └───┘    └─→┤    │
                     │AND ├─→┐
               ┌─→┤    │    │
               │    └──┘    │  ┌──┐
               │               └→┤    │
  ┌───┐    │                   │AND ├─→输出
  │输入2 ├→─┤  ┌──┐      ┌→┤    │
  └───┘    └→┤NOT ├→──┘  └──┘
                   └──┘
    上面这个逻辑电路中,无论输入是什么,输出都是False。我们就说,这个逻辑电路不存在使输出为True的一组输入。

    回到上文,给定一个逻辑电路,问是否存在一种输入使输出为True,这即逻辑电路问题。

    逻辑电路问题属于NPC问题。这是有严格证明的。它显然属于NP问题,并且可以直接证明所有的NP问题都可以约化到它(不要以为NP问题有无穷多个将给证 明造成不可逾越的困难)。证明过程相当复杂,其大概意思是说任意一个NP问题的输入和输出都可以转换成逻辑电路的输入和输出(想想计算机内部也不过是一些 0和1的运算),因此对于一个NP问题来说,问题转化为了求出满足结果为True的一个输入(即一个可行解)。

    有了第一个NPC问题后,一大堆NPC问题就出现了,因为再证明一个新的NPC问题只需要将一个已知的NPC问题约化到它就行了。后来,Hamilton 回路成了NPC问题,TSP问题也成了NPC问题。现在被证明是NPC问题的有很多,任何一个找到了多项式算法的话所有的NP问题都可以完美解决了。因此 说,正是因为NPC问题的存在,P=NP变得难以置信。P=NP问题还有许多有趣的东西,有待大家自己进一步的挖掘。攀登这个信息学的巅峰是我们这一代的 终极目标。现在我们需要做的,至少是不要把概念弄混淆了。

November 02

主题周记11/01/09——CANOEE

不早了,简记:

1 CANOEE驱车七小时来到了美丽富饶的MICHIGAN,连声惊讶这里居然有MACY‘S。CANOEE, EXCOOL, SAVVY三巨头在WESTLAND召开了对世界格局影响深远的西大地会议。

2 冬令时,又迎来了伸手不见五指的日子。

3 上周四跟老板MEET,老板口述文章的修改方案。当时没有记下来,现在已经完全想不起来他当时讲过什么了。估计他也想不起来了

4 看了一场阿童木,差点两个人包场。

5 ITHOUCH很不错,要不要入一个。

6 周济这个____终于出事了,我很高兴。作为一个律己甚严的人,我一向提醒自己:别人出事了,我不应该幸灾乐祸,因为这是不对的。

但确实TMD很高兴。

7 WINDOWS 7 和 UBUNTU 9.10上周发布,善。
October 27

主题周记10/25/09——福特堡

稍微晚了一点,现在补上。。。

这周末实验室组织了一年一度的活动:全体成员在著名的SPARTAN STADIUM观看举世瞩目的FOOTBALL比赛。
参赛双方是以EXCOOL为代表的主场应战CANOEE为代表的IOWA。当天蒙蒙细雨,但是在晚上七点比赛的时候雨过天晴。

比赛的大部分时间双方进行友好和谐的比赛,双方未有一次TOUCHDOWN,都是三分三分交替提高着比分。直到终场前两分钟,比分还是6:9,其中以CANOEE为代表的IOWA暂时领先。

就在这最后两分钟里,MSU率先打破僵局,以一记漂亮的TD拉开比分,变成13:9,整个赛场成了绿色的海洋。
IOWA在极其不利的条件下,发扬拼搏第一,比分第二的高尚风格。他们的努力在最后两秒得到了回报:以迅雷不及掩耳盗铃之势对MSU进行报复,以其人之道还制其人之身。比分最后锁定在13:16

引用THE STATE NEWS的话:

Hawkeye heartbreak——Spartan Stadium goes from hysteria to horror in 2-minute turn


这是据说自1995年MSU第一次败给IOWA,也是IOWA历史上第一次八连胜。我们实验室继续保持每看必败的记录。

October 21

吼一吼(非常boring,慎入)

已知有这样一组数据,产生方式:y = x^2 + 4x + 25

现在要求用y = (ax+b)^2来做最小二乘拟合。

我知道可以通过回归的办法可以求的a b的数值解。

问题是:有没有可能求得a b的解析表达式?

x的范围是[0, 200]

————————————————————————————————
多谢大家的回复,这个想出了解决办法,跟光浪头的思路基本上是一致的

这个问题实际上就是求(ax+b)^2与(x^2+4x+25)在[0 200]上面误差的平方和最小。
也就是:误差的平方 = ((ax+b)^2-(x^2+4x+25))^2
然后在[0 200]上面对x进行定积分,这样会得到一个关于a b的四次二元函数。
然后找一个最小值,对应的a b就是满足要求的解。

October 19

主题周记10/18/09——车

好久没有成长的见证了。。。自从回了一趟国以后。。。

今天被PLATAS鄙视了,说我居然不在实验室。着实惭愧,我周末一般是不去实验室的,在家从事形式丰富多彩的主题周末活动。

为了证明时光没有虚度,为了留下时光的美好记忆,特此开辟主题周记活动。

上周概要:子非鱼及其家属冒着生命危险前来参观访问,参观了DOWNTOWN老区和OKEMOS新区。以及EAST LANSING新老交替区。

本周概要:盼望着,盼望着,西北风来了,冬天的脚步近了。决定从此周开始偃旗息鼓,韬光养晦,把征战一夏的小破车翻翻新,升升值。顺便简单科普了一下车子的基本结构以及修车的理论常识,这种XXX FOR DUMMIES的书一般都会让人心里面多少痒痒的,觉得自己无所不能,战无不胜。
October 10

读《论东北育才东楼的倒掉》有感

暑假回国的时候,还顺道回访了一下母校。这么多年过去了,母校的一些臭规矩依然没改,于是我只能扒在铁门上向里面望了望,又望了望。分明望见那个时候东楼还是在的。

在育才的东楼只呆过初三的下半年,可以说是六年里最悠闲的时光。不过对于东楼感情无多,只是觉得每天早上要从东门走到最西边存车然后再走回最东边颇为麻烦。

怎么说呢,总而言之:

这是一段不能常常有,却十分有意义的时光,
有着一个明确的目标,并且为之不计一切的努力着。
追求的目标和实现目标的过程都是快乐的,
我愿意这每一分钟都是不虚度的。

不过做题比看PAPER有意思多了。。。