クーポンコレクター時間の厳密な下限とは何ですか?


20

古典的なクーポンコレクターの問題では、ランダムに選択されたクーポンのセットを完了するのに必要な時間が、および満たすことがよく知られています。。TnE[T]∼nln⁡nVar(T)∼n2Pr(T>nln⁡n+cn)<e−c

この上限は、おおよそ1 / c ^ 2になるチェビシェフ不等式によって与えられる上限よりも優れてい 1/c2ます。

私の質問は次のとおりです。よりも良い-チェビシェフが対応している下限のためのT?(たとえば、Pr(T<nln⁡n−cn)<e−c)?


明白な下限があるPr(T<n)=0 ...が、私はそれをあなたがしている認識して推測する
— ワンストップ

回答:


14

分析は完全に初歩的であり、正確に目的の結果を提供するため、これを2番目の回答として提供します。

命題についてはおよび、 c>0n≥1

P(T<nlog⁡n−cn)<e−c.

証明の背後にある考え方は簡単です:

  1. すべてのクーポンが収集されるまでの時間をとして。ここで、は番目(以前)の一意のクーポンが収集される時間です。、平均時間を有する幾何ランダム変数である。T=∑i=1nTiTiiTinn−i+1
  2. チェルノフ限界のバージョンを適用し、単純化します。

証明

任意のおよび任意場合、 S > 0 P(T < T )= P(E - S T > E - S T)≤ E S T E E - S Tt s>0

P(T<t)=P(e−sT>e−st)≤estEe−sT.

以来、と独立している、我々は書くことができ T i E e − s T = n ∏ i = 1 E e − s T iT=∑iTiTi

Ee−sT=∏i=1nEe−sTi

以来、今、幾何学的である、のは、成功の確率で言わせて、その後、簡単な計算ショー p i E e − s T i = p iTipi

Ee−sTi=pies−1+pi.

この問題のは、、などです。したがって、 p 1 = 1 p 2 = 1 − 1 / n p 3 = 1 − 2 / n n ∏ i = 1 E e − s T i = n ∏ i = 1 i / npip1=1p2=1−1/np3=1−2/n

∏i=1nEe−sTi=∏i=1ni/nes−1+i/n.

いくつかのに対しておよびを選択し ましょう。次に、 およびで、 t = n logs=1/n、C > 0 、E S Tは = N E - C E S = E 1 / N ≥ 1 + 1 / N N Π iが= 1のI / Nをt=nlog⁡n−cnc>0

est=ne−c
es=e1/n≥1+1/n
∏i=1ni/nes−1+i/n≤∏i=1nii+1=1n+1.

これをまとめると、

P(T<nlog⁡n−cn)≤nn+1e−c<e−c

望んだ通りに。


それはとてもいいことで、医者が注文したものです。ありがとうございました。
— デビッド

好奇心が強い@David:意図したアプリケーションは何ですか?
— 枢機

長い話。私は、興味のあるアルゴリズムの実行時間を分析するために調理したマルコフ連鎖の混合時間の下限を証明しようとしています-これは、c .collectorの問題。ところで、私はまさにこの種のチェルノフスタイルの限界を見つけようとしていましたが、でその製品を取り除く方法を理解していませんでした。選択する良いコール:-)。s = 1 / nis=1/n
— デビッド

@David、はほぼ確実に準最適ですが、を与えたため、試してみるのは明らかなことのようでした。上限。s=1/nest=ne−c
— 枢機

1
リクエスト:上記の証明は私自身のものです。問題が私を興味をそそったので、私は楽しみからそれで働きました。しかし、私は目新しさを主張しません。実際、同様の手法を使用した同様の証明がまだ文献に存在していないとは想像できません。誰かが参照を知っているなら、ここにコメントとして投稿してください。私は次のようになり非常に 1を知って興味を持って。
— 枢機

9

@cardinalは、私が探していた境界を正確に与える答えをすでに与えていますが、より強い境界を与えることができる同様のChernoffスタイルの引数を見つけました。

命題: (これはより強力です)C > π 2

Pr(T≤nlog⁡n−cn)≤exp⁡(−3c2π2).
c>π23

証明:

の答えのように、は成功確率の独立した幾何学的ランダム変数合計であるという事実を使用できます。その次及び。T i p i = 1 − i / n ETTipi=1−i/n E [ T ] = ∑ n i = 1 E [ T i ] = n ∑ n i = 1 1E[Ti]=1/piE[T]=∑i=1nE[Ti]=n∑i=1n1i≥nlog⁡n

新しい変数、およびます。その後、 S := ∑Si:=Ti−E[Ti]のPr (T ≤ N ログN - C N )≤ のPr (T ≤ E [ T ] - のC 、N )= Prの(Sを≤ - Cをn )= Pr (exp (− s SS:=∑iSi

Pr(T≤nlog⁡n−cn)≤Pr(T≤E[T]−cn)=Pr(S≤−cn)
=Pr(exp⁡(−sS)≥exp⁡(scn))≤e−scnE[e−sS]

平均を計算すると、

ES-1≥SEのZ

E[e−sS]=∏iE[e−sSi]=∏ies/pi1+1pi(es−1)≤e12s2∑ipi−2
ここで、不等式はおよびのための。es−1≥sZ≥0ez1+z≤e12z2z≥0

したがって、、書くことができます のPr(T≤NログN-CN)≤ E 1∑ipi−2=n2∑i=1n−11i2≤n2π2/6

Pr(T≤nlog⁡n−cn)≤e112(nπs)2−scn.

で最小化すると、最終的に のPr (T ≤ N ログN - C N )≤ E - 3 、C 2s>0

Pr(T≤nlog⁡n−cn)≤e−3c2π2

1
(+1)いくつかの小さなタイプミスを法として、これは素晴らしいことです。あなたがやったように、平均に近いものを中心に拡大することは、しばしばよりうまくいきます。漸近的な結果に照らして、高次の収束を見て驚くことではありません。ここで、同様の上限を示すと、はVershyninの用語で準指数関数的であることが証明されます。これは、測定濃度に関して多くの意味を持ちます。(T−nlog⁡n)/n
— 枢機

1
引数は上限に直接一般化しないようです。を(およびを)と交換すると、。ただし、この時点でできることは、、それはまだと私はしませんこれをどうするかわからない- C S - S E [ E S S ] ≤ Π I E - S / PのIc−cs−sE[esS]≤∏ie−s/pi1−spie−z1−z≤exp⁡(z22(1−z))
E[esS]≤e12s2∑ipi2(1−s/pi)
— デビッド

2
ただし、興味深いことに、(下限の)引数全体は、クーポンコレクターの問題だけでなく、有界分散のある同一ではない独立した幾何学的変数の合計でも機能するようです。具体的に:与えられた場合、各は成功確率独立GVであり、ここで、その後 T i p i ∑ i pT=∑iTiTipi∑ipi−2≤A<∞
Pr(T≤E[T]−a)≤e−a22A
— デビッド

4

重要な注意:この回答で元々提供した証拠を削除することにしました。それはより長く、より計算的で、より大きなハンマーを使用し、私が与えた他の証拠と比較してより弱い結果を証明しました。あらゆる点で、劣ったアプローチ(私の意見では)。あなたが本当に興味があるなら、私はあなたが編集を見ることができると思います。

私が最初に引用した漸近的な結果は、この回答でまだ見つかっていますが、として、すべての当てはまる他の回答で証明された限界よりも少し良くできることを示しています。n→∞ n


次の漸近的な結果が成り立ちます

P(T>nlog⁡n+cn)→1−e−e−c

そして

P(T≤nlog⁡n−cn)→e−ec.

定数および制限はとして取得されます。2つの結果に分かれていますが、どちらの場合もが非負に制約されないため、結果はほとんど同じであることに注意してください。 N → ∞ Cc∈Rn→∞c

例えば、Motwani and Raghavan、Randomized Algorithms、pp。60--63を参照してください。


また、 Davidは、この回答へのコメントで述べられている上限の証拠を親切に提供します。


n

1
Aiinlog⁡n+cn1 - N - 1 I Σ IP(τ>nlog⁡n+cn)=P(∪iAi)≤∑iP(Ai)1−n−1i∑i(1−1/n)nlog⁡n+cn≤nexp⁡(nlog⁡n+cnn)=e−c

@David、それは素晴らしく簡単です。私はすぐに別の用語で包含/除外式を拡張してみましたが、すぐにどこにも到達せず、さらに詳しく調べる時間もありませんでした。イベントは、後にクーポンが残らないイベントと同等です{T<tn}tn

@David、上記の証明にはサインの反転がありますが、他の読者にも明らかであると確信しています。
— 枢機卿

@David、あなたの質問に対する他の投稿された答えを見てください。この方法は、指定した上限とは異なりますが、ここで示した答えとは対照的に、使用するツールはほとんど基本的です。
— 枢機

2

Benjamin Doerr は、本「ランダム化検索ヒューリスティックの理論」の「ランダム化検索ヒューリスティックの分析:確率理論からのツール」の章で、オンラインPDFのリンクを参照してください。

TPr[T≤(1−ϵ)(n−1)ln⁡n]≤e−nϵ

nϵ

これが証明スケッチです。

XiitPr[Xi=1]=(1−1/n)tXiI⊆[n]Pr[∀私∈私、バツ私=1]≤∏私∈私Pr[バツ私=1]私tjt

私Pr[∀私∈私、バツ私=1|バツj=1]≤Pr[∀私∈私、バツ私=1]j∉私Pr[∀私∈私、バツ私=1|バツj=0]≥Pr[∀私∈私、バツ私=1]j私私k|私|−kn−1|私|−knj私

Pr[T≤(1−ϵ)(n−1)ln⁡n]≤(1−(1−1/n)t)nt=(1−ϵ)(n−1)ln⁡n

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.