结论:如果p是n的约数,那么满足gcd(i,n)==p的i的个数是Φ(n/p)

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1e5+11;
typedef long long ll;
ll euler(int n){
ll ans=n;
for(ll i = 2; i*i <= n; i++){
if(n%i==0){
ans=ans/i*(i-1);
while(n%i==0) n/=i;
}
}
if(n>1) ans=ans/n*(n-1);//note
return ans;
}
int main(){
int T; scanf("%d",&T);
while(T--){
int n,m; ll ans=0;
scanf("%d%d",&n,&m);
for(ll i = 1; i*i<=n; i++){
if(n%i==0){
if(i>=m) ans+=euler(n/i);
if(n/i>=m&&i*i!=n) ans+=euler(i);
}
}
printf("%lld\n",ans);
}
}

最新文章

  1. 练习题(登陆-进度条-微信接口判断qq-微信接口判断列车时刻表-)
  2. JAVA编程讲座-吴老
  3. 【DP】组合数字
  4. 正确的选择log级别
  5. 12. javacript高级程序设计-DOM2和DOM3
  6. 争夺 &amp; KM思想
  7. 小议安卓定位伪造-实战足不出户畅玩Pokemon Go
  8. shell脚本的入参
  9. javascrip笔记——图片加载
  10. android 控制手机的体积的大小 切换音效模式
  11. IDEA2017版本打可运行jar包并运行
  12. 30.Django CSRF 中间件
  13. HTTP传输编码增加了传输量,只为解决这一个问题 | 实用 HTTP
  14. 解决 VS2019 打开 edmx 文件时没有 Diagram 视图的 Bug
  15. ext整合highcharts实现饼图
  16. Mac 装机必备软件推荐
  17. 【转】stm32CubeMx上移植自己的printf()和scanf()函数
  18. [Codeforces721E]Road to Home
  19. flume-拦截器、channel选择器、sink组合sink处理器
  20. Flex的一些小实例

热门文章

  1. 2-javascript::笔记
  2. 面试题:java实例变量,局部变量,类变量 背1
  3. UOJ#46. 【清华集训2014】玄学
  4. Windows下用Nginx配置遇到的问题
  5. glib hash库GHashTable的使用实例
  6. 关于pycharm字体大小的调整
  7. http请求和返回的head字段
  8. javascript总结7:算术运算符
  9. PAT 1017 Queueing at Bank (25) (坑题)
  10. C# Excel 操作