1349. 最大公约数

Description

　　小菜的妹妹小诗就要读小学了！正所谓计算机要从娃娃抓起，小菜决定在幼儿园最后一段轻松的时间里教妹妹编程。
　　小菜刚教完gcd即最大公约数以后，一知半解的妹妹写了如下一段代码：
　　 sum:=0;
　　 for i:=1 to n-1 do
　　 for j:=i+1 to n do sum:=sum+gcd(i,j)

　　显然这个程序的效率是很低的，小明打算写一个更强的程序，在求出sum的同时比妹妹跑的更快。

Input

　　第一行一个整数t，即表示有t组数据
　　接下来t行，每行一个整数n

Output

　　t行，每行一个整数，表示n所对应的sum值

Sample Input

2
10
100

Sample Output

67
13015

Data Constraint

Hint

【数据规模】
　　20%数据t≤100,n≤100
　　40%数据t≤1000,n≤2000
　　100%数据t≤10000,n≤1000000

题解

设x为1~n-1的一个数. 若gcd(n,i)=k, 则gcd(n/k,i/k)=1. 设x=i/k,则求1~n-1中满足gcd(n,i)=k的k之和等价于满足gcd(n/k,x)=1的x个数*k之和,即∑ϕ(n/k)*k.

线性筛Euler函数
f[i]=k*f[i/k];
sum[i]=s[i-1]+f[i];

#include<cstdio>
#define N 1000005
#define ll long long
using namespace std;
int t;
ll p[N],f[N],sum[N];
void phi()
{
    p[1]=1;
    for(int i=2;i<N;i++)
        p[i]=i;
    for(int i=2;i<N;i++)
        if(p[i]==i)
            for(int j=i;j<N;j+=i)
                p[j]=p[j]/i*(i-1);
}
int main()
{
    phi();
    scanf("%d",&t);
    for(int i=1;i<N;i++)
        for(int j=i<<1;j<N;j+=i) f[j]+=i*p[j/i];
    sum[1]=0;
    for(int i=2;i<N;i++) sum[i]=sum[i-1]+f[i];
    while(t--)
    {
        int n;
        scanf("%d",&n);
        printf("%lld\n",sum[n]);
    }
    return 0;
}