欧拉函数:(转自百度百科)
欧拉函数F(x)的作用:求从1到x与x互质的数的个数。
直接求欧拉函数
int euler(int n)//计算与n互质的数的个数
{
int res = n,maxn = n;
for(int i=2;i*i < maxn;i++)
{
if(n % i == 0)
{
res = res / i * (i-1);
while(n % i == 0)//去除n中所有的i的倍数
n /= i;
}
}
if(n > 1)//最后还剩一个因子,也要算上
res = res / n * (n-1);
return res;
}
对一段区间内的数进行求欧拉函数,并打表
int euler[maxn];
void init(int maxn)
{
euler[1] = 1;
for(int i=2;i < maxn;i++)
euler[i] = i;
for(int i=2;i < maxn;i++)
{
if(euler[i] == i)
{
for(int j=i;j < maxn;j++)//对所有包含质因数i的数进行求欧拉函数值
{
euler[j] = euler[j] / i * (i-1);//先进行除法,是为了防止中间有数据溢出
}
}
}
return ;
}
题意:给出一个数n,问从1到n有多少对互质的数
HINT:(1,1)是互质的;
连续卡了几天的欧拉函数,在今天终于顿悟,虽然然不是很透彻,赶紧记下来,哈哈哈哈。。
这个题找到一对互质的数后,记得将他们两个互换也是符合条件的。这样求出总的互质的对数后乘以2就OK了。
等等!!!有个小坑就是(1,1)互换后是一样的。。。。。额,有点小尴尬。没事最后的数减去1就OK了。
上代码:
#include <cstdio>
#include <iostream>
#include <iomanip>
#include <vector>
#include <cstring>
using namespace std;
const int maxn = 50010;
int buf[maxn];
void init()
{
memset(buf,0,sizeof(buf));
buf[1] = 1;//边筛选素数,边计算合数的欧拉函数
for(int i=2; i <= maxn; i++) //当前的数是被筛到的数的质因子,当buf[i]为0是,代表buf[i]是素数
{
if(buf[i] == 0)//如果该数是素数
{
for(int j=i; j <= maxn; j+=i)//找到质因子中包含该数的合数
{
if(buf[j] == 0)
buf[j] = j;
//计算合数的欧拉函数
buf[j] = buf[j]/i*(i-1);
}
}
}
for(int i=2;i <= maxn;i++)
{
//printf("buf[i-1]: %d ",buf[i-1]);
buf[i] += buf[i-1];
}
//printf("buf[2]:%d\n",buf[2]);
}
int main()
{
int n;
init();
while(scanf("%d",&n) && n)
{
printf("%d\n",2*buf[n]-1);
}
return 0;
}