Description
给出n,求小于等于n的数中互素的整数对数
Input
多组输入,每组用例输入一整数n,以n=0结束输入
Output
对于每组输入,输出小于等于n的数中互素的整数对数
Sample Input
2
3
4
5
0
Sample Output
1
3
5
9
Solution
问题转化为前n项欧拉函数求和,注意此时euler[1]=0
Code
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<string>
using namespace std;
#define maxn 1000001
#define INF 0x3f3f3f3f
typedef long long ll;
int euler[maxn],prime[maxn],cnt;
void Get_Euler(int n)//求n以内所有数的欧拉函数值
{
memset(euler,0,sizeof(euler));
euler[1]=1;
cnt=0;
for(int i=2;i<=n;i++)
{
if(!euler[i])
{
euler[i]=i-1;
prime[cnt++]=i;
}
for(int j=0;j<cnt&&prime[j]*i<=n;j++)
{
if(i%prime[j])
euler[prime[j]*i]=euler[i]*(prime[j]-1);
else
{
euler[prime[j]*i]=euler[i]*prime[j];
break;
}
}
}
}
ll dp[maxn]={0ll};
int main()
{
Get_Euler(maxn);
euler[1]=0;//注意此处第一项欧拉函数值为0
for(int i=1;i<maxn;i++)
dp[i]=dp[i-1]+euler[i];
int n;
while(scanf("%d",&n),n)
printf("%lld\n",dp[n]);
return 0;
}