题意

给两个数n,m,问gcd(x,n)>=m (x<=n),问满足此式的x的数量(T<=100,m<=n<=1e9)

分析

枚举n的大于等于m的约数i,不难发现gcd(i,n/i * i)==i ,满足此式i的数量为(n/i)的欧拉函数,直接暴力枚举到√n即可

#include<bits/stdc++.h>
#define ll long long
using namespace std;const int maxn = 2e5+7;int t,n,m;int eular(int x)
{int ans=x;for(int i=2;i*i<=x;i++){if(x&&x%i==0){ans=ans-ans/i;while(x%i==0)x/=i;}}if(x>1)ans=ans-ans/x;return ans;
}int main()
{scanf("%d",&t);while(t--){scanf("%d%d",&n,&m);int nn=sqrt(n);int ans=0;for(int i=1;i<=nn;i++){if(n%i) continue;if(i>=m) ans+=eular(n/i);if(n/i>=m&&(n/i)!=nn) ans+=eular(i);}printf("%d\n",ans);}return 0;
}

View Code

转载于:https://www.cnblogs.com/Superwalker/p/8934296.html

HDU2588 欧拉函数相关推荐

  1. poj2154-color-polyan次二面体+欧拉函数优化

    N<=1e9,O(nlogn)的做法会超时.从枚举置换转变为枚举轮换长度,然后可以利用欧拉函数,把复杂度变为O(√n * logn) 1 /*-------------------------- ...

  2. hdu 1286( 欧拉函数 )

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1286 数学题真的是有点吃不消了... View Code 1 #include<iostream ...

  3. 费马定理中值定理_数论-欧拉函数、欧拉定理

    欧拉函数 积性函数 满足 ( 互质) 定义 对于正整数 ,欧拉函数是小于等于 的所有数中与 互质的数的 个数. 欧拉函数是积性函数(这个证明不是很显然,这个链接里面有很多种证明方法) 记作: 公式 , ...

  4. poj2154(Polya+欧拉函数优化模版)

    #include <cstdio> #include <cstring> #include<iostream> using namespace std; const ...

  5. 【数学专题】约数个数与欧拉函数

    整理的算法模板合集: ACM模板 目录 一.约数个数 1. AcWing 1291. 轻拍牛头 2. AcWing 1294. 樱花 2.1 AcWing 197. 阶乘分解 3. AcWing 19 ...

  6. 【数学知识】三种方法求 [1,n] 中所有数欧拉函数(线性筛欧拉函数优化至 O(n) )

    整理的算法模板合集: ACM模板 ①直接求小于或等于n,且与n互质的数个数(求[1,n]中所有数的欧拉函数时间复杂度:O(nn)O(n\sqrt{n})O(nn​)) ②求[1,n]之间每个数的质因数 ...

  7. bzoj 1409 Password 矩阵快速幂+欧拉函数

    可以发现,该数组的mi就是斐波那契数列 所以要矩阵快速幂搞出第n位 但是斐波那契数列上涨的很快,这就需要欧拉定理了 p^phi(q)%q=1(gcd(p,q)==1) p是素数,所以可以用 然后需要5 ...

  8. hdu1695(莫比乌斯)或欧拉函数+容斥

    题意:求1-b和1-d之内各选一个数组成数对.问最大公约数为k的数对有多少个,数对是有序的.(b,d,k<=100000) 解法1: 这个能够简化成1-b/k 和1-d/k 的互质有序数对的个数 ...

  9. POJ 2480 (约数+欧拉函数)

    题目链接: http://poj.org/problem?id=2480 题目大意:求Σgcd(i,n). 解题思路: 如果i与n互质,gcd(i,n)=1,且总和=欧拉函数phi(n). 如果i与n ...

最新文章

  1. 我花了三个小时写了一道题的六千字题解....(POJ 2888 Magic Bracelet)
  2. [Winform]检测exe是否已经运行,并将其置顶
  3. fzu 2139 久违的月赛之二
  4. hadoop综述_Hadoop书籍赠品–综述
  5. 将PostgreSQL PL / Java安装为PostgreSQL扩展
  6. Node.js抓取网页信息并展示(cheerio网络爬虫)
  7. SpringBoot : 注解@Resource
  8. 为提升在线语音识别效率,他创造了两种升级版算法模型
  9. 充一次电使用一年的手机_超级手机电池:充一次电用一年多?
  10. 44 年前的今天,改变世界的 TA 诞生了!
  11. [转载] python函数isdisjoint方法_Python中的isdisjoint()函数
  12. 数字图像处理-绪论(下)
  13. Swift5.2 Extension(工具类)
  14. GNU、GCC与G++的区别
  15. 可以用云服务器挂机传奇吗,蜂窝云挂机划算吗 游戏蜂窝云挂机怎么用
  16. 3D打印断电、断料续打,打印一半中断续打改代码,也适用于单喷头多色打印(用cura切片)
  17. 主板显卡服务器维修,显卡无显示怎么维修 显卡无显示解决方法【详细介绍】...
  18. 打开计算机 访问不到d盘,教大家电脑找不到d盘的原因及解决办法
  19. latex参考文献居中_LaTeX中参考文献管理以及注意事项
  20. iphone html阅读,iPhone如何使用Safari浏览器阅读列表功能

热门文章

  1. erlang注意(经典)
  2. 【SpringBoot 】 组件管理 + 属性注入
  3. 大厂HR:“不会数据分析,你还想干运营?”
  4. 中国3亿中年男人都在花钱买什么?数据分析后告诉你答案,很现实
  5. Python采集3000条北京二手房数据,看我都分析出了啥?
  6. python通过get方法获取key对应的值
  7. mysql 函数 数字转字串_mysql中字符串和数字的互转函数
  8. table合并单元格_element ui el-table 合并单元格
  9. 【论文笔记】Revisiting Graph based Collaborative Filtering: A Linear Residual Graph Convolutional Network
  10. 洛谷 P1983 [NOIP2013 普及组] 车站分级