欧拉心算(bzoj4804)

问题 B: 欧拉心算(bzoj4804)
时间限制: 2 Sec 内存限制: 256 MB
提交: 74 解决: 28
[提交][状态]
题目描述
Description
给出一个数字N

Input
第一行为一个正整数T，表示数据组数。
接下来T行为询问，每行包含一个正整数N。
T<=5000,N<=10^7
Output
按读入顺序输出答案。
Sample Input
1

10
Sample Output
136
输入
输出
提示
林肯是大头：
http://www.accoders.com/problem.php?cid=1966&pid=1
bzoj
https://www.lydsy.com/JudgeOnline/problem.php?id=4804
思路：
隔了一天就忘了啊，，，
看看代码8
看了眼题解缓过来了
就是枚举谁为最大公约数->phi[i]
再计算有多少对i,j除以最大公约数后互质(利用phi前缀和计算，上限为n/i)
分块见代码8

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=10000000+10;
int phi[maxn];
int p[700000];
int num[maxn];
long long pre[maxn];
long long ans[maxn];
int main()
{
    for(int i=1;i<=maxn;i++)
    phi[i]=i;
    pre[1]=1;
    for(int i=2;i<=maxn;i++)
    {
        if(phi[i]==i)
        {
            p[++p[0]]=i;
            for(int j=i;j<=maxn;j+=i)
            phi[j]=phi[j]/i*(i-1);
        }
        num[i]=p[0];
        pre[i]=pre[i-1]+1ll*phi[i];
    }
    int t;
    scanf("%d",&t);
    for(int l=1;l<=t;l++)
    {
        int n;
        scanf("%d",&n);
        long long ans=0;
        int nxt;//nxt为无比玄学的分块操作 
        /*
        for(int i=1;i<=n;i++)//枚举最大公约数为i 
        {
            ans+=(pre[n/i]*2-1)*phi[i];
        }*///超时做法 
        for(int i=1;i<=n;i=nxt+1)
        {
            nxt=n/(n/i);//i到nxt之内，n/i余数相等 
            ans+=(pre[nxt]-pre[i-1])*(pre[n/i]*2-1);
            //此步操作把 phi[i]分块优化为前缀和之差 
        }
        printf("%lld\n",ans);
    }
    return 0;
}