整理的算法模板合集: ACM模板


题目传动门

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <cmath>
#include <cstring>
#include <vector>using namespace std;
typedef long long ll;
typedef pair<int , int> PII;
const int N = 5000007, M = 500007, T = 1000, Mod = 998244353, INF = 0x3f3f3f3f;int phi[N];
bool vis[N];
ll f[N], g[N];
int primes[N], cnt;void get_phi(int n)
{phi[1] = 1;f[1] = 1;for(int i = 2; i <= n; ++ i) {if(vis[i] == 0) {primes[ ++ cnt] = i;phi[i] = i - 1;}f[i] = 1ll * phi[i] * i / 2;for(int j = 1; j <= cnt && i * primes[j] <= n; ++ j) {vis[i * primes[j]] = true;if(i % primes[j] == 0) {phi[i * primes[j]] = phi[i] * primes[j];break;}else phi[i * primes[j]] = phi[i] * (primes[j] - 1);}}//Dirichelet前缀和优化for(int j = 1; j <= cnt; ++ j) {for(int i = 1; i * primes[j] <= n; ++ i) {f[i * primes[j]] += f[i];}}
}int n, t;void slove()
{scanf("%d", &n);printf("%lld\n", f[n] * n);
}int main()
{scanf("%d", &t);get_phi(N - 1);while(t -- ) {slove();}
}

SP5971 LCMSUM - LCM Sum(莫比乌斯反演 ,推柿子,经典)相关推荐

  1. SP5971 LCMSUM - LCM Sum

    SP5971 LCMSUM - LCM Sum 思路 ∑i=1nlcm(i,n)\sum_{i = 1}^{n}lcm(i, n)i=1∑n​lcm(i,n) =>∑i=1ningcd(i,n) ...

  2. 【BZOJ2226】LCM SUM,数论之一维LCM(莫比乌斯反演)

    Time:2016.06.18 Author:xiaoyimi 转载注明出处谢谢 传送门 思路: 一开始的我只能推出一个O(n√⋅n√=n)O(\sqrt n·\sqrt n=n)的式子(两次分块), ...

  3. 【BZOJ2154】Crash的数字表格,数论练习之二维LCM(莫比乌斯反演)

    Time:2016.06.18 Author:xiaoyimi 转载注明出处谢谢 看不见的分割线 传送门 思路: 设n≤m ∑i=1n∑j=1mlcm(i,j)\displaystyle \sum^n ...

  4. 【香蕉OI】GCD 和 LCM (莫比乌斯反演)

    文章目录 题意 思路 代码 题意 给出 TTT 组询问,每组询问求 ∑i≤n∑j≤m[gcd(i,j)≤a]lcm(i,j)\sum_{i\le n}\sum_{j\le m}[gcd(i,j)\le ...

  5. SPOJ LCMSUM - LCM Sum

    题意是求: $\sum_{i = 1}^{n}lcm(i, n)$ $= \sum_{i = 1}^{n}\frac{ni}{gcd(i, n)}$ $= n\sum_{i = 1}^{n}\frac ...

  6. 解题报告(十五)莫比乌斯反演与积性函数(ACM / OI)

    整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 好久没写反演了,写些水题慢慢捡起来 目录 A. (BZOJ 2440: 中山市选2011 )完全平方数 ...

  7. 我也不知道什么是莫比乌斯反演和杜教筛

    upd:发现这篇文章里面似乎有很多细节上的错误,如果还有错误的话在cnblogs下评论吧,我会改的QwQ,cnblogs戳这里. upd:正在写一篇复习向的文章,之后贴链接,可以作为这篇文章的一个补充 ...

  8. 莫比乌斯反演 做题记录

    来自Peterwuyihong 的题单. 前置知识 前置芝士1 数论分块 UVA11526 H(n) P2261 [CQOI2007]余数求和 P2260 [清华集训2012]模积和 其中有一个式子需 ...

  9. 数论 —— 莫比乌斯反演

    [反演] 假设我们手头有个数列 F,通过某种变换 H,可以得到函数 G.,即: 但现在只有函数 G,需要求 F,那么我们就需要寻找一种变换 ,使得 G 在经过这种变换后能够获得 F,这个过程即为反演, ...

最新文章

  1. Java入门时的一些基本概念的理解(j2ee,j2se,j2me,jdk,sdk,jre,jvm,跨平台)
  2. [Java] SSH框架笔记_Struts2配置问题
  3. 函数节流与防抖的实现
  4. python基础教程:装饰器
  5. 额外篇 | basemap(上)
  6. mysql 事务 myisam,InnoDB和MyISAM的差别(mysql事务处理)
  7. oracle数据库sqlloader,sql loader ---ORACLE SQLLDR
  8. 在子线程中更改主线程中的控件的信息,在子线程中用toast
  9. 【实践驱动开发3-003】TI WL1835MODCOM8 在android的移植 - 软件获取2
  10. Android 系统(268)---native保活5.0以下方案推演过程以及代码详述
  11. 错误未找到引用源_你好,C++(77)12.1 用右值引用榨干C++的性能
  12. 下载centos(阿里云)
  13. vmware7序列号
  14. dota2中断与服务器连接,dota2断开服务器连接
  15. 【盘点】2014年跨境进口型电商十大模式
  16. win7旗舰版升级win10
  17. java实现京东云第三方登录
  18. mac笔记本电脑外接显示器没有声音
  19. html立方体旋转展开,html5—旋转立方体
  20. JDBC中的statement、executeQuery()、SQL injection attacks注入式攻击、preparedStatement

热门文章

  1. 一文读懂图像局部特征点检测算法
  2. 【C++简明教程】Python和C++指定元素排序比较
  3. 用于道路目标检测的少镜头学习
  4. 【从零学习OpenCV 4】Image Watch插件的使用
  5. Factory-pattern 三种工厂模式
  6. 给你的网站穿上外衣 - HTTPS 免费部署指南
  7. Zabbix-03-3 告警变量值
  8. Altium Designer同一个工程里不同原理图导入到不同的PCB
  9. 推荐65个以自然风光为背景的UI设计
  10. php checkbox美化,纯CSS3美化radio和checkbox_html/css_WEB-ITnose