bzoj 4805: 欧拉函数求和

Time Limit: 15 Sec Memory Limit: 256 MB
Submit: 539 Solved: 304
[Submit][Status][Discuss]

Description

给出一个数字N，求sigma(phi(i)),1<=i<=N

Input

正整数N。N<=2*10^9

Output

输出答案。

Sample Input

Sample Output

HINT

Source

By FancyCoder

杜教筛模板题，好久没打了练练手hhh

(一开始把i和j写反了WA了一墙hhhh)

/**************************************************************Problem: 4805User: JYYHHLanguage: C++Result: AcceptedTime:1524 msMemory:64772 kb
****************************************************************/#include<bits/stdc++.h>
#define ll unsigned long long
#define maxn 5000000
using namespace std;
map<ll,ll> mmp;
int zs[maxn+5],t=0;
ll phi[maxn+5],n;
bool v[maxn+5];inline void init(){phi[1]=1;for(int i=2;i<=maxn;i++){if(!v[i]) zs[++t]=i,phi[i]=i-1;for(int j=1,u;j<=t&&(u=zs[j]*i)<=maxn;j++){v[u]=1;if(!(i%zs[j])){phi[u]=phi[i]*zs[j];break;}phi[u]=phi[i]*(zs[j]-1);}}for(int i=1;i<=maxn;i++) phi[i]+=phi[i-1];
}inline ll sum(ll x){if(x<=maxn) return phi[x];if(mmp.count(x)) return mmp[x];ll ans=x*(x+1)>>1ll,j;for(ll i=2;i<=x;i=j+1){j=x/(x/i);ans-=sum(x/i)*(j-i+1);}mmp[x]=ans;return ans;
}int main(){init();scanf("%llu",&n);printf("%llu\n",sum(n));return 0;
}

转载于:https://www.cnblogs.com/JYYHH/p/8287379.html

bzoj 4805: 欧拉函数求和相关推荐

欧拉函数的相关应用 noj欧拉函数求和+noj 最大公约数求和
注意求欧拉函数之和是每个因子的欧拉函数之和不是质因子.而欧拉函数的值是它本身与它的因子件事互质的关系,这样的因子有多少个.点击打开链接 #include<stdio.h> #include ...
BZOJ - 2186 欧拉函数
题目的意思大概是求1~N!中和M!互质的数的个数因为对欧拉函数理解不够深刻所以我是分析得到结果的: 当N<=M的时候显然符合要求的数的个数为0: 当N>M的时候我们要求的是1~N!中不含 ...
bzoj 2818 欧拉函数
思路:就是对于某个数q,跟他互质的数p,kp和kq的最大公约数是k,那么这个数能组成的答案的数量就是phi[i]乘以某个质数,且乘积小于n 基于这种思路写下这个代码 #include <cstd ...
BZOJ 4802: 欧拉函数（大数因数分解算法 Pollard_rho 和素数测试算法 Miller_Rabin）
Description 已知N,求phi(N) Input 正整数N.N<=10^18 Output 输出phi(N) Sample Input 8 Sample Output 4 Soluti ...
BZOJ2818-莫比乌斯反演/欧拉函数
这道题之前没有看数论函数的时候搞懂了,想到直接用欧拉函数做,现在再来看第一个想法就是这不是莫比乌斯反演嘛. 但还是能用简单数论知识直接做出来的还是尽量做简单一点. 两种方法想到后都写的差不多对了,都爆 ...
AcWing 874. 筛法求欧拉函数（线性筛法求欧拉函数模板）
给定一个正整数 n,求 1∼n 中每个数的欧拉函数之和. 输入格式共一行,包含一个整数 n. 输出格式共一行,包含一个整数,表示 1∼n 中每个数的欧拉函数之和. 数据范围 1≤n≤10^6 输入 ...
bzoj 1409 Password 矩阵快速幂+欧拉函数
可以发现,该数组的mi就是斐波那契数列所以要矩阵快速幂搞出第n位但是斐波那契数列上涨的很快,这就需要欧拉定理了 p^phi(q)%q=1(gcd(p,q)==1) p是素数,所以可以用然后需要5 ...
bzoj 3884: 上帝与集合的正确用法（欧拉函数）
3884: 上帝与集合的正确用法 Time Limit: 5 Sec Memory Limit: 128 MB Submit: 2574 Solved: 1151 [Submit][Status] ...
bzoj 2705: [SDOI2012]Longge的问题（欧拉函数）
2705: [SDOI2012]Longge的问题 Time Limit: 3 Sec Memory Limit: 128 MB Submit: 3077 Solved: 1914 [Submit ...