问题定义:
输入一个数组a[1,2,...,n](假设均可通过下标访问,包括a[n]),输出一个随机顺序的新序列a[1',2',...,n']
比如输入{1,2,3},则该数组一共有6种排列{1,2,3}{1,3,2}{2,1,3}{2,3,1}{3,1,2}{3,2,1},则程序输出以上任何一种排列的概率均为1/6.
算法(假设数组下标从1到n):
for i=1:n-1
index=random(i:n)
swap(a[i],a[index])
n个数的排列一共有n!种组合,那么输出的序列的概率为1/n!
从第i步开始,某个特定的数被选中的概率为1/(n-i+1),则最后输出数组选定的概率为p=1/n * 1/(n-1) *...*1/(n-(n-1))
输入一个数组a[1,2,...,n](假设均可通过下标访问,包括a[n]),输出一个随机顺序的新序列a[1',2',...,n']
比如输入{1,2,3},则该数组一共有6种排列{1,2,3}{1,3,2}{2,1,3}{2,3,1}{3,1,2}{3,2,1},则程序输出以上任何一种排列的概率均为1/6.
算法(假设数组下标从1到n):
for i=1:n-1
index=random(i:n)
swap(a[i],a[index])
end
n个数的排列一共有n!种组合,那么输出的序列的概率为1/n!
从第i步开始,某个特定的数被选中的概率为1/(n-i+1),则最后输出数组选定的概率为p=1/n * 1/(n-1) *...*1/(n-(n-1))
即p=1/n!;
源码:
int _tmain(int argc, _TCHAR* argv[])
{
srand((unsigned int)time(0));
int a[4]={1,2,3,4};
RandSequence(a,4);
return 0;
}
void RandSequence(int a[], int length)
{
for(int i=0;i<length-1;i++)
{
int toSwap=rand()%(length-i)+i;
Swap(&a[i],&a[toSwap]);
}
}
void Swap(int *a, int *b)
{
int temp=*a;
*a=*b;
*b=temp;
}