以下内容为《算法:C语言实现》笔记。
1. 求对数:
/****************************************************************
求对数函数
这个函数更重要的意义在于:若用a进制表示n, 则函数的返回值为表示n所需要的位数
复杂度: O(lgn)
****************************************************************/
int loga(int n, int a)
{
assert(n > 0);
assert(a > 1);
int i;
for (i = 0; n > 0; i++, n /= a)
;
return i;
}
2. 求素数:
a. 基本方法:
/****************************************************************
基本法求素数
复杂度: Nsqrt(N)
****************************************************************/
int get_prime2(int n, int* primearray)
{
assert(n>1);
assert(NULL != primearray);
int cnt = 0;
int i, j;
for (i = 2; i <= n; i++)
{
for (j = 2; j <= sqrt(i); j++)
{
if (0 == (i%j)) break;
}
if (j > sqrt(i))
{
primearray[cnt++] = i;
}
}
return cnt;
}
b. 筛选法:
/****************************************************************
筛选法求小于等于n的所有素数:
算法:
首先,将所有数组元素设置为1,表示没有已知的非素数;
然后将已知为非素数(即为已知素数的倍数)的索引对应的数组元素设置为0;
如果所有较小素数的倍数都设置为0之后,a[i]仍然保持为1,则可判断它是素数
复杂度:
显然,程序的运行时间与下式成正比:
N + N/2 + N/3 + N/7 + N/11 + ...
< N + N/2 + N/3 + N/4 + ... = NH(N) (约为NlnN)
(H(n): 调和函数)
****************************************************************/
int get_prime(int n, int* primearray)
{
assert(n>1);
assert(NULL != primearray);
char *a = (char *)malloc(n+1);
assert(a != NULL);
int i, j;
int cnt = 0;
for (i = 2; i <= n; i++)
{
a[i] = 1;
}
for (i = 2; i <= n; i++)
{
if (1 == a[i])
{
for (j = i; j <= n / i; j++)
{
a[i*j] = 0;
}
}
}
for (i = 2; i <= n; i++)
{
if (1 == a[i])
{
primearray[cnt++] = i;
}
}
free(a);
return cnt;
}
注意: a的空间复杂度为O(1),b的空间复杂度为O(N)。