bzoj2818(欧拉函数)

题目链接：https://www.lydsy.com/JudgeOnline/problem.php?id=2818

2818: Gcd

Time Limit: 10 Sec   Memory Limit: 256 MB
Submit: 7551   Solved: 3362
[ Submit][ Status][ Discuss]

Description

给定整数N，求1<=x,y<=N且Gcd(x,y)为素数的
数对(x,y)有多少对.

Input

一个整数N

Output

如题

Sample Input

4

Sample Output

4

思路：一道入门级别的欧拉函数题（对我是final级别）。。。

由于当gcd(i,j)=1时，必有gcd(i*p,j*p)=p（其中p为素数），所以我们先算出sumphi,即phi的前缀和。又由于这道题（x,y）和（y,x）是两种，所以在最后加的时候还需乘2，再减掉多加 的1个phi[1]，最后答案就是程序里写的那个式子，这里就不写了。

详见代码：

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 10000000;
typedef long long ll;
bool check[MAXN + 10];
int phi[MAXN + 10];
int prime[MAXN + 10];
int tot;    //  素数个数
ll sumphi[MAXN+10];
void phi_and_prime_table(int N)
{
    memset(check, false, sizeof(check));
    phi[1] = 1;
    tot = 0;
    for (int i = 2; i <= N; i++)
    {
        if (!check[i])
        {
            prime[tot++] = i;
            phi[i] = i - 1;
        }
        for (int j = 0; j < tot; j++)
        {
            if (i * prime[j] > N)
            {
                break;
            }
            check[i * prime[j]] = true;
            if (i % prime[j] == 0)
            {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
            else
            {
                phi[i * prime[j]] = phi[i] * (prime[j] - 1);
            }
        }
    }
    return ;
}
int main()
{
	int n;
	scanf("%d",&n);
	phi_and_prime_table(n);
	sumphi[1]=phi[1];
	for(int i=2;i<=n;i++)
	sumphi[i]=sumphi[i-1]+phi[i];
	ll ans=0;
	for(int i=0;i<tot;i++)
	{
		 if(prime[i]<n)
		 ans+=(2*sumphi[n/prime[i]]-1);
		 else
		 break;
	}
	printf("%lld\n",ans);
	return 0;
}