【bzoj4804】欧拉心算 狄利克雷卷积

Description

给出一个数字N

Input

第一行为一个正整数T，表示数据组数。
接下来T行为询问，每行包含一个正整数N。
T<=5000,N<=10^7
Output

按读入顺序输出答案。
Sample Input

1

10
Sample Output

136

题解
给YZH dalao 跪烂烂烂烂烂


F(D)为欧拉函数和莫比乌斯函数的卷积

代码

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#define ll long long
#define inf 1000000000
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
const int N=10000005;
ll F[N],ans;
int n,tot,p[N],pos;
bool flag[N];
void prepapre()
{
    F[1]=1;
    for (int i=2;i<N;i++)
    {
        if (!flag[i]){p[++tot]=i;F[i]=i-2;}
        for (int j=1;j<=tot&&p[j]*i<N;j++)
        {
            flag[i*p[j]]=1;
            if (i%p[j]==0)
            {
                if (i/p[j]%p[j]==0) F[i*p[j]]=F[i]*p[j];
                else F[i*p[j]]=F[i/p[j]]*(p[j]-1)*(p[j]-1);
                break;
            }
            F[i*p[j]]=F[i]*F[p[j]];
        }
    }
    for (int i=2;i<N;i++) F[i]+=F[i-1];
}
int main()
{
    prepapre();int Case=read();
    while (Case--)
    {
        n=read();ans=0;
        for (int i=1;i<=n;i=pos+1)
        {
            pos=n/(n/i);
            ans+=(F[pos]-F[i-1])*(n/i)*(n/i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}