[杜教筛模板] 51Nod 1239 欧拉函数之和
来源:互联网 发布:淘宝买家订单打印 编辑:程序博客网 时间:2024/06/10 00:18
模板题
#include<cstdio>#include<cstdlib>#include<algorithm>#include <tr1/unordered_map>typedef long long ll;using namespace std;using namespace std::tr1;const int maxn=10000000;int prime[1000000],num;int vst[maxn+5],phi[maxn+5];const int P=1e9+7;const int inv=500000004;inline void Pre(){ phi[1]=1; for (int i=2;i<=maxn;i++){ if (!vst[i]) prime[++num]=i,phi[i]=i-1; for (int j=1;j<=num && (ll)i*prime[j]<=maxn;j++){ vst[i*prime[j]]=1; if (i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break; }elsephi[i*prime[j]]=phi[i]*phi[prime[j]]; } } for (int i=1;i<=maxn;i++) (phi[i]+=phi[i-1])%=P;}unordered_map<ll,int> S;inline int Sum(ll n){ if (n<=maxn) return phi[n]; if (S.find(n)!=S.end()) return S[n]; int tem=(ll)(n%P)*((n+1)%P)%P*inv%P; ll l,r; for (l=2;l*l<=n;l++) (tem+=P-Sum(n/l))%=P; for (ll t=n/l;l<=n;l=r+1,t--) r=n/t,(tem+=P-(ll)(r-l+1)*Sum(t)%P)%=P; return S[n]=tem;}int main(){ ll n; freopen("t.in","r",stdin); freopen("t.out","w",stdout); Pre(); scanf("%lld",&n); printf("%d\n",Sum(n)); return 0;}
0 0
- [杜教筛模板] 51Nod 1239 欧拉函数之和
- 【51Nod 1239】欧拉函数之和
- 51nod 1239欧拉函数之和
- 【51NOD 1239】欧拉函数之和
- 51nod 1239 欧拉函数之和
- [51NOD]1239 欧拉函数之和
- 51nod 1239 欧拉函数之和
- 51Nod-1239-欧拉函数之和
- 51Nod 1239 欧拉函数之和
- 51nod 1239 欧拉函数之和
- 51nod 1239 欧拉函数之和
- 51nod 1239 欧拉函数之和
- 51 NOD 1239 欧拉函数之和(杜教筛)
- [杜教筛] 51Nod 1239: 欧拉函数之和
- 51NOD 1136 欧拉函数模板
- 51nod 1136 欧拉函数 【模板 】
- 51nod-1040-最大公约数之和(欧拉函数)
- 51NOD 1040 最大公约数之和(分析 + 欧拉函数)
- J.U.C--locks--AQS分析
- Oracle命令---第一天
- MPI 数据结合体
- wordpress配置邮箱发送功能
- 代数问题3*3矩阵的逆矩阵
- [杜教筛模板] 51Nod 1239 欧拉函数之和
- 切换到24的SDK出现了ViewPager消失
- 用cygwin移植c写的基于cocos2dx引擎的游戏为android编译时android不支持libiconv
- [杜教筛模板] BZOJ 3944 Sum
- 2016-12-28 迪米特法则+依赖倒转原则+里氏替换原则
- JPA中的manytomany
- 2016-12-28 外观模式 + 适配器模式+迪米特法则+依赖倒转原则+里氏替换原则
- 【记录】第一次搭SSM框架
- angular实现商品筛选功能(过滤器)