bzoj 4804 欧拉心算 欧拉函数,莫比乌斯

Posted Kaiser

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj 4804 欧拉心算 欧拉函数,莫比乌斯相关的知识,希望对你有一定的参考价值。

欧拉心算

Time Limit: 15 Sec  Memory Limit: 256 MB
Submit: 408  Solved: 244
[Submit][Status][Discuss]

Description

 给出一个数字N

Input

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

 

Output

按读入顺序输出答案。

 

Sample Input

1
10

Sample Output

136

HINT

 

Source

 

 

 1 #include<cstring>
 2 #include<cmath>
 3 #include<iostream>
 4 #include<algorithm>
 5 #include<cstdio>
 6 
 7 #define ll long long
 8 #define N 10000007
 9 using namespace std;
10 inline int read()
11 {
12     int x=0,f=1;char ch=getchar();
13     while(ch<\'0\'||ch>\'9\'){if (ch==\'-\')f=-1;ch=getchar();}
14     while(ch>=\'0\'&&ch<=\'9\'){x=(x<<1)+(x<<3)+ch-\'0\';ch=getchar();}
15     return x*f;
16 }
17 
18 int n;
19 int cnt,pri[N],phi[N];
20 ll sum[N];
21 bool flag[N];
22 
23 void init_phi()
24 {
25     phi[1]=1;
26     for (int i=2;i<=10000000;i++)
27     {
28         if (!flag[i])pri[++cnt]=i,phi[i]=i-1;
29         for (int j=1;j<=cnt&&pri[j]*i<=10000000;j++)
30         {
31             flag[pri[j]*i]=1;
32             if (i%pri[j]==0)
33             {
34                 phi[pri[j]*i]=pri[j]*phi[i];
35                 break;
36             }
37             else phi[pri[j]*i]=phi[i]*phi[pri[j]];
38         }
39     }
40     for (int i=1;i<=10000000;i++)sum[i]=sum[i-1]+phi[i];
41 }
42 void solve(ll n)
43 {
44     ll last,ans=0;
45     for (int i=1;i<=n;i=last+1)
46     {
47         last=n/(n/i);//琛ㄧずn涓埌杈惧摢閲屼负姝㈤兘鏄痭/i鐨?
48         ans+=(sum[last]-sum[i-1])*((sum[n/i]<<1)-1);
49     }
50     printf("%lld\\n",ans);
51 }
52 int main()
53 {
54     init_phi();
55     int T=read();
56     while(T--)
57     {
58         n=read();
59         solve(n);
60     }
61 }

 

以上是关于bzoj 4804 欧拉心算 欧拉函数,莫比乌斯的主要内容,如果未能解决你的问题,请参考以下文章

bzoj4804欧拉心算 欧拉函数

[BZOJ4804]欧拉心算

bzoj4804: 欧拉心算 欧拉筛

BZOJ2226[Spoj 5971] LCMSum 莫比乌斯反演(欧拉函数?)

bzoj 4916: 神犇和蒟蒻欧拉函数+莫比乌斯函数+杜教筛

BZOJ.2705.[SDOI2012]Longge的问题(莫比乌斯反演 欧拉函数)