蓄水池抽樣Reservior Sampling

標籤：style   blog   http   color   io   檔案   for   ar   

編程珠璣第12章練習題10：

如何從n個對象（可以依次看到這n個對象，但事先不知道n的值）中隨機播放一個？具體說來，如何在事先不知道文字檔行數的情況下讀取檔案，從中隨機播放並輸出一行?

解答：我們總選擇 第1行，並以機率1/2選擇第2行，以機率1/3選擇第3行，依次類推，在這一過程結束時，每一行選中的機率是相等的（都是1/n，其中n是檔案的總行數）

i = 0;while more input lines    with probability 1.0/++i            choice = this input lineprint choice

推廣：

如何從未知或者很大樣本空間隨機地取k個數？

for i=k+1:n   m = random(1,i);  //inclusive   if(m<=k)     swap(array[m],array[i]);end print array[1:k]

 

                       --------------引用http://www.cnblogs.com/HappyAngel/archive/2011/02/07/1949762.html