简述 真没啥好说的，随便无脑化最后就出来了
ans=12∑d=1n?nd?dφ(d)+n2【#|51nod1227 平均最小公倍数】
代码
//杜教筛
#include
#include
#define maxn 4700000
#define mod 1000000007ll
using namespace std;
typedef long long ll;
ll N, phi[maxn+10], f[maxn+10], _2, _6;
int prime[maxn+10];
bool mark[maxn+10];
void m(ll &x){if(x>mod or x<-mod)x%=mod;
}
ll sqr(ll x){m(x);
return x*x%mod;
}
ll s1(ll x){m(x);
return x*(1+x)%mod*_2%mod;
}
ll s2(ll x){m(x);
return x*(1+x)%mod*(2*x+1)%mod*_6%mod;
}
void shai()
{
ll i, j;
phi[1]=1;
for(i=2;
i 推荐阅读
- 数据结构和算法|LeetCode 的正确使用方式
- #|7.分布式事务管理
- #|算法设计与分析（Java实现）——贪心算法（集合覆盖案例）
- #|算法设计与分析（Java实现）—— 动态规划 （0-1 背包问题）
- #|阿尔法点亮LED灯(一）汇编语言
- #|Multimedia
- #|ARM裸机开发（汇编LED灯实验（I.MX6UL芯片））
- 基础课|使用深度优先搜索（DFS）、广度优先搜索（BFS）、A* 搜索算法求解 (n^2 -1) 数码难题，耗时与内存占用（时空复杂度）对比（附（(n^2 - 1) 数码问题控
- #|学习笔记 | Ch05 Pandas数据清洗 —— 缺失值、重复值、异常值
- win10|搏一搏 单车变摩托，是时候捣鼓一下家中的小米电视机啦。