POJ 2773 欧拉函数， 素数表

这题很很很不懂，谁能详细讲解下。

素数表第一次用，用了不熟练，欧拉函数求了好久出来，不太懂。

#include <stdio.h>
#include <string.h>
#include <math.h>
#define MAXN 1000010
int prime[MAXN], phi[MAXN];
void Init()
{
int i, j, k = ceil(sqrt(1000010));
//求素数表
memset(prime, 0, sizeof(prime));
prime[0] = prime[1] = 1;
for (i = 2; i < k; i++)
{
if (prime[i] == 0)
{
for (j = 2 * i; j < MAXN; j += i)
prime[j] = 1;
}
}
//求欧拉函数值
for (i = 1; i < MAXN; i++)
phi[i] = i;
for (i = 2; i < MAXN; i++)
{
if (prime[i] == 0)
{
for (j = i; j < MAXN; j += i)
//phi(n) = n * (1 - 1/p1) * (1 - 1/p2)....[pi是与n的质因素]
phi[j] = phi[j] / i * (i - 1);
}
}
}
int GCD(int a, int b)
{
int t;
while (b)
{
t = b;
b = a % b;
a = t;
}
return a;
}
int main()
{
int m, K, cnt, i, flag, mul;
Init();
while (scanf("%d%d", &m, &K) != EOF)
{
//判断m内第几个与m互素
cnt = K % phi[m];
//判断分为了几组
mul = K / phi[m];
//正好分为整数组
if (cnt == 0)
{
//m的最后一个与其互素的
cnt = phi[m];
//分组减一
mul--;
}
flag = 0;
for (i = 1; i < m; i++)
//判断是否互素
if (GCD(m, i) == 1)
{
flag++;
if (flag == cnt)
break;
}
//printf("phi[m] = %d i = %d mul = %d\n", phi[m], i, mul);
printf("%d\n", mul * m + i);
}
return 0;
}

这个是我第n次RE后写的，超时了，有什么特别好的优化方法吗？ 请教下大家

#include <stdio.h>
int Coprime (int a, int b)
{
int t;
while (b)
{
t = b;
b = a % b;
a = t;
}
if (a == 1)
return 1;
return 0;
}
int main()
{
int i, m, cnt, K, flag, phi_m, mul;
while (scanf("%d%d", &m, &K) != EOF)
{
phi_m = flag = 0;
for (i = 1; i <= m; i++)
if (Coprime(m, i))
phi_m++;
cnt = K % phi_m;
mul = K / phi_m;
if (cnt == 0)
{
cnt = phi_m;
mul--;
}
for (i = 1; i <= m;  i++)
if (Coprime(m, i))
{
flag++;
if (flag == cnt)
break;
}
printf("%d\n", mul * m + i);
}
return 0;
}

今天又练了到素数表的题, 不过时间好慢200+MS 不知道别人的0MS怎么弄出来的

有谁知道怎么优化比较好, 发个链接, 谢谢。

POJ 3006 筛素数

#include <stdio.h>
#include <string.h>
int prime[1000010];
void init()
{
int i, j;
memset(prime, 0, sizeof(prime));
prime[0] = prime[1] = 1;
for (i = 2; i < 1000; i++)
{
if (prime[i] == 0)
for (j = 2 * i; j < 1000000; j += i)
prime[j] = 1;
}
}
int main()
{
int a, d, n, num, i;
init();
while (scanf("%d %d %d", &a, &d, &n) != EOF)
{
if (!(a || d || n)) break;
num = i = 0;
while (num != n)
{
if (prime[a + i * d] == 0) num++;
i++;
}
printf("%d\n", a + (i - 1) * d);
}
return 0;
}