GDUT寒假训练场4-Happy 2006

题目描述

Two positive integers are said to be relatively prime to each other if the Great Common Divisor (GCD) is 1. For instance, 1, 3, 5, 7, 9...are all relatively prime to 2006.

Now your job is easy: for the given integer m, find the K-th element which is relatively prime to m when these elements are sorted in ascending order.

输入格式

The input contains multiple test cases. For each test case, it contains two integers m (1 <= m <= 1000000), K (1 <= K <= 100000000).

输出格式

Output the K-th element in a single line.

输入样例

2006 1
2006 2
2006 3

输出样例

1
3
5

参考思路

题目意思是给定一个数让你求出和它互素的第k个数字。互素就是两者最大公约数是1。

我们发现对于公约数来说有gcd(a,b)=gcd(a+t*b,b)。所以在[1，m-1]中与m互素的个数与在[k*m+1，(k+1)*m-1]的互素的个数是一样的，即都等于pre(m)。

所以只要找出[1，m-1]中与m互素的数，就能找到对用的接下来的区间中互素的数：

那么答案就等于第k%pre(m)个与m互素的值p+m*(k/pre(m))。

我们可以直接用gcd一个个求放进数组里，然后根据上面的规律输出。也可以利用欧拉函数，这样快很多，毕竟欧拉函数就是小于等于n和n互素的数的个数，我们在函数模板里加入素数判断数组即可。

参考代码1

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int pri[1000000];
int gcd(int a, int b)
{
    return b?gcd(b,a%b):a;
}
int main()
{
    int m,k;
    while (~scanf("%d%d",&m,&k))
    {
        int i,t;
        t = 0;
        for (i = 1; i <= m; i++)
        {
            if (gcd(m, i) == 1)
                pri[t++] = i;
        }
        if (k % t)
            printf("%d\n",k/t*m+pri[k%t-1]);
        else
            printf("%d\n",(k/t-1)*m+pri[t-1]);
    }
    return 0;
}

 参考代码2

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N = 1000010;
bool prime[N];

int eluar(int n)
{
	int m = n, ans = n;
	for (int i = 2; i * i <= n; i++) {
		if (n % i == 0) {
			ans = ans / i * (i - 1);
			while (n % i == 0) n /= i;
			for (int j = 1; j * i <= m; j++) prime[j * i] = true;
		}
	}
	if (n > 1) {
		ans = ans / n * (n - 1);
		for (int j = 1; j * n <= m; j++) prime[j * n] = true;
	}
	return ans;
}
int main() {
	int m, k;
	while (cin >> m >> k)
	{
		memset(prime, 0, sizeof prime);
		int cnt = eluar(m);
		int t = k / cnt;
		if (k % cnt == 0) t--;//判断是否--
		k -= cnt * t;//判断区间剩余数字
		int sum = 0, i;
		for (i = 1; i <= m; i++)
		{
			if (!prime[i]) sum++;
			if (sum == k) break;
		}
		cout << t * m + i << endl;//i
	}
	return 0;
}

小楼一夜听春雨