Happy 2006

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.

Input：

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

Output：

Output the K-th element in a single line.

Sample Input：

Smaple Output：

题目大意：

题目给出了m和k，找出第k个与m互素（最大公约数都为1）的数，这道题考察了欧几里得定理。

思路分析：

gcd（a，b），如果a与b互素，那么就有gcd（tb+a，b），tb+a与b也一定互素，这道题也可以参考求模运算(n=k*q+r），k=n/q,r=n%q。

代码块：

#include<stdio.h>
int gcd(int m,int i)//找出与m互素的数 
{
	return i==0 ? m : gcd(i,m%i);
}
int main()
{
	int m,k,i;
	int a[1000005];//这里数组要大一点，要不然就WA， 
	while(scanf("%d %d",&m,&k)!=EOF)
	{
		int j=0;
		for(i=1;i<=m;i++)
		{
			if(gcd(m,i)==1)
			{
				a[j++]=i;
			}
		}
		if(k%j!=0)//这里分开讨论，如果余数不0，说明他们不是倍数关系。 
		{		//这里一般都是小于j的数，所以k/j是为0，后面就直接在a数组里面直接寻找
		        //当k大于j时候，根据周期性来找，例：2017，可以当做是2016+a[0]. 
			printf("%d\n",(k/j)*m+a[k%j-1]);//因为下标从0开始，所以减去1 
		}
		else
		printf("%d\n",(k/j-1)*m+a[j-1]);//如果余数为0，说明是倍数，K/j得出是整数，
										//并且k是大于j，不是在a数组中找出
										//只能根据周期性来分析.
	}
 }