P2568 GCD （欧拉函数+gcd）

传送门

1<=n<=1e7

首先考虑gcd(x,y)==1 (x,y<=n)的对数，答案就是2*Σphi[i]-1。(从1到n枚举x，欧拉函数加进答案，然后乘2，最终结果要-1因为phi[1]算了两次)

但是题目要求是gcd(x,y)==p (p为素数),其等价为gcd(a*p,b*p)==p,(a b互质)
于是我们可以令x=a*p,y=b*p，

那么对于每个p 都只需要求出满足1<=a*p,b*p<=n,gcd(a,b)==1的a,b对数即可，这个就转化成了我们起初考虑的问题。

具体做法就枚举素数j，求出2*Σphi[n/j]-1,再相加即可

可以推出来:gcd(x,y)==k (x<=n,y<=n)的x,y的对数为:
所以这题就枚举素数，然后求和就行。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pii pair<int,int>
const int maxn = 1e7+10;
const int mx = 40;
const int mod = 1e9+7;
const ll inf = 0x3f3f3f3f3f3f3f3f;
const int INF = 1e9+7;
//给定整数n  求满足1<=x,y<=n gcd(x,y)是素数的x y的对数 ,x y可以相同
int vis[maxn];//0表示素数
int book[maxn];//第i个素数
ll phi[maxn];//欧拉函数
ll pre[maxn];
int n;
void phi_primer()
{
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            book[++book[0]]=i;
            phi[i]=i-1;
        }
        for(int j=1;j<=book[0] && i*book[j]<=n;j++)
        {
            vis[i*book[j]]=1;
            if(i%book[j]==0)
            {
                phi[i*book[j]]=phi[i]*book[j];
                break;
            }
            else //i book[j]互质
            {
                phi[i*book[j]]=phi[i]*(book[j]-1);
            }
        }
    }
}

int main()
{
  	scanf("%d",&n);
    phi_primer();
    for(int i=1;i<=n;i++) pre[i]=pre[i-1]+phi[i];
    ll ans=0;
    for(int i=1;i<=book[0] && book[i]<=n;i++)
    {
        ans = ans + pre[n/book[i]]*2-1;
    }
    printf("%lld",ans);
    return 0;
}

补充一些欧拉函数的性质:

• 定义函数φ(k)表示数字1…(k-1)中，和k互质的数字的个数。 φ(k)就是欧拉函数

• φ(1)=1 如果i是质数 φ(i)=i-1

• 质数x的k次幂 φ(x^k) = (x-1)*(x^(k-1)) = φ(x^(k-1)) * x

• 如果m n互质 φ(mn)=φ(m)*φ(n)

• 如果m == 2 n为奇数, 由于φ(2)=1,φ(2n)=φ(n)

• 如果i%p==0 其p是质数 则 φ(i*p)=φ(i)*p

• 当n>2 φ(n)为偶数

• 小于n的数中 与n互质的数的和为 n*φ(n)/2
  理由:若gcd(n,i)=1 则gcd(n,n-i)=1 , 也就是n ,n-i成对出现