标签:blog io os for sp 2014 log amp bs
题目:找到整数区间[1,n]中所有的互质数对。
分析:数论,筛法,欧拉函数。在筛素数的的同时,直接更新每个数字的欧拉函数。
每个数字一定会被他前面的每个素数筛到,而欧拉函数的计算是n*π(1-1/pi);
其中,pi是n的素数因子,所以可以利用筛法来计算欧拉函数,然后求和;
注意,这时求出的欧拉函数为所有小于n的数m与n的互质对,所以要乘以2;
(1,1)只有一个,所以还要减掉1。
说明:打表计算,查询输出。
#include <iostream>
#include <cstdlib>
using namespace std;
int euler[50001];
int sums[50001];
int used[50001];
int main()
{
for (int i = 0 ; i < 50001 ; ++ i) {
used[i] = 0;
euler[i] = i;
}
for (int i = 2 ; i < 50001 ; ++ i)
if (!used[i]) {
for (int j = i ; j < 50001 ; j += i) {
used[j] = 1;
euler[j] = euler[j]/i*(i-1);
}
}
sums[0] = 0;
for (int i = 1 ; i < 50001 ; ++ i)
sums[i] = sums[i-1]+euler[i];
int n;
while (cin >> n && n)
cout << 2*sums[n]-1 << endl;
return 0;
}
标签:blog io os for sp 2014 log amp bs
原文地址:http://blog.csdn.net/mobius_strip/article/details/40506089