コルモゴロフ複雑性のソースを表示

'''コルモゴロフ複雑性'''（コルモゴロフふくざつせい、{{lang-en|Kolmogorov complexity}}）とは、[[計算機科学]]において有限長のデータ列の[[複雑さ]]を表す[[指標]]のひとつで、出力結果がそのデータに一致する[[プログラム (コンピュータ)|プログラム]]の長さの最小値として定義される。'''コルモゴロフ複雑度'''、'''コルモゴロフ＝チャイティン複雑性''' ({{lang|en|Kolmogorov-Chaitin complexity}}) とも呼ばれる。

[[File:Mandelpart2.jpg|300px|right|thumb|この画像は[[フラクタル]]図形である[[マンデルブロ集合]]の一部である。この[[JPEG]]ファイルのサイズは17KB以上（約140,000ビット）ある。ところが、これと同じファイルは140,000ビットよりも遥かに小さいコンピュータ・プログラムによって作成することが出来る。従って、このJPEGファイルのコルモゴロフ複雑性は140,000よりも遥かに小さい。]]

コルモゴロフ複雑性の概念は一見すると単純なものであるが、[[停止性問題|チューリングの停止問題]]や[[ゲーデルの不完全性定理]]と関連する深遠な内容をもつ。コルモゴロフ複雑性やその他の文字列やデータ構造の[[複雑性]]の計量を研究する計算機科学の分野は[[アルゴリズム情報理論]]と呼ばれており、1960 年代末に[[アンドレイ・コルモゴロフ]]、[[レイ・ソロモノフ]]<!--Ray Solomonoff-->、[[グレゴリー・チャイティン]]によって創始された。

==定義==
厳密には、ある有限長の文字列 ''x'' として表されるデータ列のある万能計算機 ''u'' におけるコルモゴロフ複雑性は以下の式で定義される:
:<math>K_u (x) = \min_{ p : u(p)=x }\, l(p)</math>
ここで、''p'' は計算機 ''u'' のための[[プログラム (コンピュータ)|プログラム]]であり、''u''(''p'') はそれを実行したときに出力される文字列である。''l''(''p'') はプログラムの長さを表す。ただし、プログラムは入力をもたず、つねに決まった出力を返すようなものとする。 ここで万能計算機とは[[チューリングマシン|万能チューリングマシン]]と同等の能力をもつ計算機を意味し、例えば [[C言語|C]] など通常のプログラム言語を解釈し実行するような処理系だと考えてよい。

== 直観的な例 ==

ある文字列が別の文字列よりも複雑であると言うためにはどうすればよいかという問題を考えよう。 例えば、同じ 60 文字の長さの 0, 1 からなる以下の 2 種類の文字列が与えられたとする。 

 010101010101010101010101010101010101010101010101010101010101

 110010000110000111011110111011001111101001000010010101111001

これらを見比べると、直観的に前者の文字列はより簡単であって、後者の方はより複雑であると感じる。 この直観を明確にするひとつの方法として、文字列を言語で説明することを考えよう。 このとき前者は「01 の 30 回の繰り返し」と 11 文字で説明できる。 それに対して、後者はそのような簡単な説明が与えられそうもなく、説明するには文字列そのものを含めて 60 文字以上で記す他はないように思える。 

この例のように言語による記述によって短くできる文字列がある一方で、圧縮できないような文字列も存在しており、文字列の説明の長さは文字列そのものの複雑さと関係していると考えられる。 このことをより形式的に取り扱うために[[アルゴリズム]]という概念を用いて定式化されたものがコルモゴロフ複雑性である。 

コルモゴロフ複雑性を定めるためには、文字列に対する形式的な記述言語をまず指定しなければならない。 このような記述言語としては何かの具体的なプログラミング言語を選べばよい。 あるプログラム言語 ''L'' を固定したとき、 ''L'' のプログラム ''p'' が有限長の文字列 ''x'' を出力するなら、''p'' は ''x'' の「記述」であるということにする。 

このとき特に ''p'' として ''x'' 自身を明示的に含むような自明な記述がある。 例えば、上記の前者の文字列を[[標準出力]]に出力する [[Perl]] プログラムは、

 print "010101010101010101010101010101010101010101010101010101010101"

のようになるだろう。 このようなプログラムの長さは明らかに ''x'' の長さよりもある定数分だけ長くなる。 しかし文字列に何かの構造が発見できれば、適切なアルゴリズムを用いることによってより短いプログラムを作れるかもしれない。 例えば次のPerlプログラムは上と同じ文字列を出力する。

 print "01" x 30

このように記述言語 ''L'' における文字列 ''x'' のあらゆる記述のうちで最小の長さをもつ記述が見出せたとするなら、その長さが ''L'' における ''x'' のコルモゴロフ複雑性となる。

== 基本的な性質 ==

=== 自明な上限 ===

前述の例で示されたように明示的に文字列 ''x'' をプログラムに含めることができるので、すべての ''x'' に対してそのコルモゴロフ複雑性 ''K''<sub>''u''</sub>(''x'') は ''x'' 自体の長さ |''x''| を定数分以上越えることはない。 すなわち、

'''定理：''' 任意の ''u'' に対して、ある定数 ''c'' が存在して、任意の ''x'' に対して、
:<math>K_u(x) \leq |x| + c</math>
が成り立つ。

=== 記述言語による相対性 ===

定義から明らかなように、コルモゴロフ複雑性は記述言語あるいは万能計算機に依存する。 しかし、ある万能計算機 ''u''<sub>1</sub> から別の万能計算機 ''u''<sub>2</sub> にプログラムを移植しても、コルモゴロフ複雑性はたかだか（''u''<sub>1</sub> と ''u''<sub>2</sub> によって決まる） 定数分しか増えない。なぜなら2つの万能計算機は、必ずもう一方の計算機を[[エミュレーション|エミュレート]]できるからである。 ''u''<sub>2</sub> 上で ''u''<sub>1</sub>を模倣するエミュレーションプログラム &epsilon;<sub>1,2</sub> を作り、その上で ''u''<sub>1</sub> のためのプログラム ''p'' を動かせば、結果として ''u''<sub>2</sub> の上でプログラム ''p'' を動かせたことになる。 そしてこのエミュレーションプログラムはエミュレートするプログラムの大きさにかかわらずつねに一定である。 従って、 ''u''<sub>2</sub> 上でのコルモゴロフ複雑性はたかだか ''l''(''p'') + ''l''(&epsilon;<sub>1,2</sub>) である。 逆の場合も同様にエミュレートができるので、すなわち、

'''定理：''' 任意の万能計算機 ''u''<sub>1</sub>, ''u''<sub>2</sub> に対し、ある定数 ''c''<sub>1,2</sub> が存在して、任意の ''x'' に対し、
:<math>|K_1(x) - K_2(x)| \leq c_{1,2}.</math>
が成り立つ。

なおコルモゴロフ複雑性の議論では、記述言語の違いによりこのようなある定数分を除いて成立するという関係が頻出する。

=== 圧縮不能な文字列 ===

文字列 ''x'' が ''K''<sub>''u''</sub>(''x'') ≧ |''x''|  となるなら ''x'' は「圧縮不能」<!--incompressible-->であるという。 このような圧縮不能な文字列が存在するならどの程度存在するかを考えたい。 まず、簡単な[[濃度 (数学)|基数 (濃度)]] についての考察から、ある長さの文字列 ''x'' すべての中には「圧縮不能」な文字列が少なくとも 1 つは存在していることがわかる。 なぜなら、長さ ''n'' のバイナリ列 (2 種類のアルファベットからなる文字列) は 2<sup>''n''</sup> 個存在するが、 それより短いバイナリ列は長さ 0 も含めて 2<sup>''n''</sup> - 1 しか存在しないからである。さらに同じ理由で、長さ ''n'' 以下の 2<sup>''n''+1</sup> - 1 個のバイナリ列のうち、その半数より多い 2<sup>''n''</sup> が長さ ''n'' をもつのであるから、ある長さ以下のバイナリ列のうち圧縮不能な文字列は半数より多い。

一般に、圧縮不能の定義に自然数 ''c'' だけの余裕を入れ、''K''<sub>''u''</sub>(''x'') ≧ |''x''| - ''c'' を満たすとした ''x'' を「''c'' 圧縮不能」<!--''c''-incompressible--> という。 長さ ''n'' - ''c'' より短いバイナリ列の個数は、長さ ''n'' のバイナリ列と比べ、 (2<sup>''n''-''c''</sup> - 1) / 2<sup>''n''</sup> しかないので、 ''c'' 圧縮不能な文字列は全体の 1 - 1/2<sup>''c''</sup> より多い。 例えば、 ''c'' = 100 [[ビット]]だけ圧縮できるような文字列は、たかだか全体の 2<sup>-100</sup> &#8764; 10<sup>-30</sup> しかない。

よって、ほとんどの文字列は文字列の長さと大きく違わないコルモゴロフ複雑性しかもたないだろうと結論できる。

==コルモゴロフ複雑性の計算不能性==
コルモゴロフ複雑性に関する[[計算理論]]上の興味深い帰結のひとつは、コルモゴロフ複雑性が実効的に計算できないということである。

'''定理：''' ''K'' は[[計算可能関数]]ではない。

言い換えれば、任意の文字列 ''s'' を入力として整数 ''K''(''s'') を出力するようなプログラムは存在しない。証明は背理法による。以下、あるプログラムを構成し、そのプログラムから出力される文字列が、そのプログラムよりも長いプログラムからでなければ出力され得ないという矛盾を導く。次のようなプログラムがあるとしよう：

   '''function''' KolmogorovComplexity('''string''' ''s'')

これは文字列 ''s'' を入力とし ''K''(''s'') を返却する。ここで次のようなプログラムを考える。

   '''function''' GenerateComplexString('''int''' ''n'')
      '''for''' i = 1 '''to''' 無限大:
         '''for each''' 長さが丁度 i である string s の全集合に対して
            '''if''' KolmogorovComplexity(''s'') >= ''n''
               '''return''' ''s''
               '''quit'''

このプログラムは  KolmogorovComplexity() をサブルーチンとして呼び出す。このプログラムは長さが1から無限大に至るまでのありとあらゆる文字列を調べ、コルモゴロフ複雑性が少なくとも ''n'' 以上であるような文字列を見つけたらそれを返す。従って、任意の正の整数 ''n'' について、これはコルモゴロフ複雑性が少なくとも ''n'' 以上であるような文字列を生成する。このプログラム自身は固定的な長さを持つが、その長さを ''U'' と書こう。プログラム GenerateComplexString() に対する入力は整数 ''n'' だが、この大きさは ''n'' を表すのに必要なビット数で測ることとすると、これは log<sub>2</sub>(''n'') である。さて、ここで更に次のようなプログラムを考える：

   '''function''' GenerateParadoxicalString()
       '''return''' GenerateComplexString(''n''<sub>0</sub>)

このプログラムは GenerateComplexString() をサブルーチンとして呼び出すが、その際 ''n''<sub>0</sub> という任意の変数を引き渡す。このプログラムは コルモゴロフ複雑性が少なくとも ''n''<sub>0</sub> 以上であるような文字列 ''s'' を出力する。このとき、''n''<sub>0</sub> を適切に選ぶと矛盾を導くことができる。この値を選ぶに当っては、''s'' が GenerateParadoxicalString() というプログラムによって生成される点に着目すると良い。このプログラムの長さは、たかだか

:<math> U + \log_2(n_0) + C  \quad </math>

である（''C'' は GenerateParadoxicalString() が  GenerateComplexString() を呼び出す前後の固定部分の長さを表す）。''n'' と log<sub>2</sub>(''n'') では前者の方が急速に増大することが明らかなので、ある ''n''<sub>0</sub> が存在して次の関係を満たす。
:<math> U + \log_2(n_0) + C < n_0. \quad </math>
ところが、これは「少なくとも ''n''<sub>0</sub> 以上の複雑性を持つ」ということの定義に反してしまう。何故なら ''K''(''s'') の定義により、文字列 ''s'' を生成したプログラムの長さは少なくとも ''n''<sub>0</sub> 以上はある筈なのに、''s'' を生成した GenerateParadoxicalString() の長さは ''n''<sub>0</sub> よりも短いからである。以上より、"KolmogorovComplexity" というプログラムは、任意の文字列について複雑性を計算することは出来ないことが示された。

この証明は背理法によるが、ここで現れる矛盾は[[ベリーのパラドックス]]に似ている：「''n'' を30字未満では表現できない最小の正の整数としよう」。他の証明方法として、K が計算不能であることを[[停止性問題|停止問題]] H が計算不能であることに帰着して示すこともできる。K と H は[[チューリング次数|チューリング同値]]である<ref>{{cite web | url=http://www.daimi.au.dk/~bromille/DC05/Kolmogorov.pdf | title=Course notes for Data Compression - 2 Kolmogorov complexity Fall 2005 | last=Miltersen | first=Peter| accessdate=2009-07-19 | format=PDF | publisher=University of Aarhus | date=2005-09-29}}</ref>。

== 関連する問題 ==

=== データ圧縮の限界 ===

(可逆な) 圧縮プログラムは与えられた文字列 ''x'' に関して、文字列 ''y'' を返す関数とみることができる。 ただし ''y'' は対応する展開プログラムで ''x'' に復元できなければならない。 よって、コルモゴロフ複雑性の定義から、展開プログラムの長さと ''y'' の長さの和は ''x'' のコルモゴロフ複雑性以上でなければならない。 この意味でコルモゴロフ複雑性はその文字列データに対する究極的な圧縮を限界づけている。

通常、圧縮プログラムでは ''y'' の長さは ''x'' よりも小さくなることが期待されるが、上述の圧縮不能な文字列の議論は圧縮プログラムに関しても成立するので、実際には圧縮できない ''x'' が少なくとも半数存在し、ほとんどは大きな圧縮ができない。 圧縮プログラムが我々が扱う多くの文字列を圧縮できるようにみえるのは、我々が実際に扱う文字列が、可能なすべての文字列と比べて極めて偏ったものにすぎないからである。

また、コルモゴロフ複雑性が計算不能であるという事実は、すべての文字列 ''x'' に対してコルモゴロフ複雑性がしめす究極の圧縮を実現するプログラムを作ることが原理的に不可能であるということを意味している。もしすべての文字列をコルモゴロフ複雑性まで圧縮するプログラムが存在すれば、その出力結果の長さを数えることでコルモゴロフ複雑性が計算できてしまい、上述の事実と矛盾するからである。

=== ランダム性 ===

十分長い文字列において「圧縮不能」な文字列は、アルゴリズム化できるような何の規則性ももたないと考えられるので「ランダム」な文字列だとみなせるだろう。 コルモゴロフ複雑性を無限長の文字列に拡張したとき、圧縮できないような文字列は「[[アルゴリズム的にランダムな列]]」 (algorithmically random sequence) と呼ばれている。 正確には、与えられた無限列 ''x'' のすべての接頭部分列が ''c'' 圧縮不能であるような ''c'' が存在するとき、''x'' はアルゴリズム的にランダムである。 

有限列の場合と同様に濃度の議論から、プログラムは[[可算]]個しかないのに対して、無限長のバイナリ列は[[非可算]]個であるので、この意味で「ほとんどすべて」の無限列はアルゴリズム的にランダムである。 しかしこのようなランダム列をプログラムによって生成することは決してできない。 統計的にランダム性をもつように見える列、例えば[[疑似乱数]]列や、[[円周率]]のような計算可能な[[超越数]]、あるいは有限長で記述される初期状態をもつ{{仮リンク|記号力学系|en|Symbolic dynamics}}が示す[[カオス理論|カオス]]などは、すべてそれらを生成するプログラムが存在するので、この意味でのランダムではない。 ランダムな文字列はその特定の文字列の説明の複雑さという直観を元にしたコルモゴロフの意味ではもっとも複雑となるが、ランダムな列のアンサンブルは逆に統計的に特徴づけることがもっとも簡単な文字列でもある。

==チャイティンの不完全性定理==
殆どの文字列は、より「圧縮された」形では表現できないという意味で複雑である。しかしながら、文字列の長さがある閾値を超える場合、その文字列が複雑であることを形式的に証明することは出来ない、ということが言えてしまう。正確な定式化は以下の通り。まず、自然数を扱う特定の公理系 S を固定する。この公理系は次のことが出来る程度には強いものとする。即ち、'''S''' に含まれる式 '''F'''<sub>'''A'''</sub> を、文字列の複雑さに関する或る主張 '''A''' に関連付けることが出来る。この際、'''F'''<sub>'''A'''</sub>が '''S''' の公理から証明可能なら、これに対応する主張 '''A''' も真になるとする。この「定式化」に際しては、[[ゲーデル数化]]のような人工的な符号化を用いてもよいし、適用しようとする '''S''' をもっと明快に表現するような定式化を用いても良い。

'''定理：'''次の式（を'''S''' の中で定式化したもの）を公理系 '''S''' の中で証明できるように文字列 ''s'' を取れないような、定数 ''L''（具体的な値は公理系と記述言語にのみ依存する） が存在する。
: <math>K(s) \geq  L.</math>
圧縮不能に近い文字列は大量にあることから、殆ど全ての場合にこの式は真である筈なことに注意されたい。

この結果の証明は[[ベリーのパラドックス]]に似た自己言及的な構成を用いる。以下、背理法による。この定理が偽だと仮定すると、次のことが言える。

:'''主張(X)'''：任意の整数 ''n'' について、ある文字列 ''s'' が存在し、体系'''S'''の中で式 「''K''(''s'')≧''n''」（が'''S'''の中で定式化できると仮定して）を証明できる。

'''S'''内の形式的証明全てを実効的に枚挙する何らかの手段を見つけることができる。

   '''function''' NthProof('''int''' ''n'')

は入力として ''n'' を取り、適当な証明を出力する。この関数は全ての証明を枚挙する。その中には我々にとって差し当たり興味のない証明も混ざるだろう（NthProof()が枚挙する証明の中には例えば[[平方剰余の相互法則]]の証明、[[フェルマーの小定理]]の証明、[[フェルマーの最終定理]]の証明など、様々な既知の証明を'''S'''内の形式的言語に翻訳したものが現れるだろう）。この中の幾つかは ''K''(''s'')≧''n'' という形をした複雑性に関する式の証明である（''s'' と ''n'' は'''S'''内の言語における定数）。さて、次のプログラムが存在する。

   '''function''' NthProofProvesComplexityFormula('''int''' ''n'')

このプログラムは ''n'' 番目の証明が式 ''K''(''s'')≧''L'' を証明しているどうかを判定する。文字列 ''s'' と整数 ''L'' はそれぞれ以下のプログラムで計算可能である：

   '''function''' StringNthProof('''int''' ''n'')

   '''function''' ComplexityLowerBoundNthProof('''int''' ''n'')

ここで次のようなプログラムを考えよう。

   '''function''' GenerateProvablyComplexString('''int''' ''n'')
      '''for''' i = 1 '''to''' 無限大:
         '''if'''  NthProofProvesComplexityFormula(i) '''and''' ComplexityLowerBoundNthProof(i) >= ''n'' 
            '''return''' StringNthProof(''i'')
            '''quit'''

任意の ''n'' について、このプログラムは[[形式体系]]'''S'''内のありとあらゆる証明を調べ上げて、''K''(''s'')≧''L''（但し''L'' ≧ ''n'' ）を満たす文字列と証明を探し出そうとする。'''主張(X)'''によりこのプログラムは必ず停止する。このプログラムの長さを ''U'' としよう。
このとき、ある整数 ''n''<sub>0</sub> があって、''U'' + log<sub>2</sub>(''n''<sub>0</sub>) + ''C''  < ''n''<sub>0</sub> を満たす。ここで ''C'' は、次のプログラムがGenerateProvablyComplexString()を呼び出す前後の固定的な長さである。

    '''function''' GenerateProvablyParadoxicalString()
       '''return''' GenerateProvablyComplexString(''n''<sub>0</sub>)
       '''quit'''

プログラム GenerateProvablyParadoxicalString() は文字列 ''s'' を出力するが、このとき ''L'' が存在して ''K''(''s'')≧''L''（但し ''L'' ≧ ''n''<sub>0</sub>）を満たし、'''S'''内でこれを形式的に証明できる。特に ''K''(''s'')≧''n''<sub>0</sub> は真である。ところが、''s'' は長さが ''U''+log<sub>2</sub>(''n''<sub>0</sub>)+''C'' であるプログラムでも生成できるので、その複雑性は ''n''<sub>0</sub> よりも小さい。これは矛盾である。よって'''主張(X)'''は成り立たないことが証明された。

同様のアイディアは[[チャイティンの定数]]の性質を証明する際にも使われている。

==脚注==
<references/>

== 関連項目 ==
*[[チューリングマシン]]
*[[アンドレイ・コルモゴロフ]]
*[[可逆圧縮]]
*[[乱数列]]
*[[チャイティンの定数]]
*[[記述計算量]]

== 参考文献 ==
*Ming Li, P.M.B. Vitanyi, ''An Introduction to Kolmogorov Complexity and Its Applications'', Springer-Verlag, 1993:1997, ISBN 0387948686

{{データ圧縮}}
{{DEFAULTSORT:こるもころふふくさつせい}}
[[Category:計算理論]]
[[Category:アルゴリズム的情報理論]]
[[Category:数学に関する記事]]
[[Category:数学のエポニム]]