C++实例分析组合数的计算与排列组合的产生
组合数的计算
使用加法递推—O(n^2)
边界条件
int C[1001][1001]; // 根据实际需要开数组,必要时采用高精度类型
//……
memset(C,0,sizeof(C));
for (int i=0;i<=n;i++)
{
C[i][0]=1;
for (int j=0;j<=i;j++)
C[i][j]=C[i-1][j-1] + C[i-1][j];
}
使用乘法递推—O(n)
边界条件
必须先乘后除,否则除不开
一个小优化:
(m>n/2时用)
int C[1001]; // C[m]其实表示C[n][m],必要时采用高精度类型
……
C[0]=1;
if (m>n-m) m=n-m;
for (int i=1;i<=m;i++)
C[i] = (n-i+1) * C[i-1] / i; // 如果怕溢出,可以把中间结果转化成long long。
排列和组合的产生(无重集元素)
int item[N]; // 第i位要放置的数字
bool used[N];
int n,m;
全排列
将 n 个数字 1~n 进行排序,有多少种排序方法?
使用深度优先搜索,对 n 个位置逐个进行试探。时间复杂度为 O(n!)。
void full_ permutation(int depth)
{
if (depth==n)
{
// print(); // 输出结果
return;
}
for (int i=0; i<n; i++)
if (!used[i])
{
used[i]=true;
item[depth]=i+1;
try(depth+1);
used[i]=false; // 别忘记清除”使用”标记
}
}
一般组合
从 n 个元素 1~n 中任取 m 个元素,有多少种取法?
一个合法的组合有这样一个特点:排在右面的数字一定严格大于左面的数字。比如说某一位上取了 3,那么从 4 开始搜索下一位就可以了。
void combination(int depth, int p)
{
if (depth==m)
{
// print(); // 输出结果
return;
}
for (int i=p+1 ; i<n-(m-depth) ; i++)
{
// 由于后面的元素一定前面的大,所以不需要标记used了。
item[depth]=i;
try(depth+1);
}
}
combination(0,0);
全组合
输入 n 个数,求这 n 个数构成的集合的所有非空子集。
和一般组合不同,这次只要产生一个解,就马上输出。
void full_combination(int l, int p)
{
for (int i=0; i<l; i++) // 每次进入递归函数都输出
cout<<item[i]<<" ";
cout<<endl;
for (int i=p; i<n; i++)
{
item[l] = i; // 在l位置放上该数
full_combination(l+1, i+1); // 填下一个位置
}
}
full_combination(0, 0);
注意:对于一个整数,每一位不是 0 就是 1,所以可以用整数来表示一个集合。具体实例可参见 “2.8Healthy Holsteins”。
由上一排列产生下一排列
① 从右往左寻找第一个小于右边的数,位置为 j。
② 在 j 位置的右边寻找大于 aj的最小数字 ak(位置 k)
③ 将 aj 与 ak的值进行交换
④ 将数列的 j+1 位到 n 位倒转。
int a[N]; // 初始化:a[i]是字典序最小的排列, 0≤i<N
int j,k, p,q, temp;
j=(n-1) - 1;
while ((j>=0)&&(a[j]>a[j+1])) j--; // 从右往左寻找第一个小于右边的数,位置为j。
if (j>=0) // 如果j<0说明已经排完了。
{
k=n-1;
while (a[k]<a[j]) k--; // 在j位置的右边寻找大于aj的最小数字ak(位置k)
swap(a[j], a[k]); // 将aj与ak的值进行交换
for (p=j+1,q=n-1; p<q; p++,q--) // 将数列的j+1位到n位倒转
swap(a[p], a[q]);
}
STL 中有与此相同的算法。头文件为<alGorithm>。
next_permutation(序列第一项的地址, 序列最后一项的地址+1):产生下一排列。
prev_permutation(序列第一项的地址, 序列最后一项的地址+1):产生上一排列。
这两个函数能够用于可重集的排列。
由上一组合产生下一组合
① 从右向左寻找可以往下取一个元素的数,位置为 j。
(举个例子:从 7 个数中取 4 个数,有一个组合为 1367,那么 6、7 就不能再往下取了)
② 数列的 j 位到 n 位重新取元素。
注意:
① 从 N 个连续元素中取 M 个元素。如果元素序号不连续,就需要修改下面的“+1”。
② 右侧的数字一定严格大于左侧的数字。
int a[M]; // 初始化:a[i]是字典序最小的排序, 0≤i<M,1≤a[i]≤N
//……
int j=m-1;
while ((j>=0)&&(a[j]==n-(m-1 -j))) j--;
if (i>=0)
{
a[j]++;
for (int k=j+1; k<m; k++) a[k]=a[k-1]+1;
}
到此这篇关于c++实例分析组合数的计算与排列组合的产生的文章就介绍到这了,更多相关C++组合数的计算内容请搜索以前的文章或继续浏览下面的相关文章希望大家以后多多支持!
相关文章