题解 \(by\;zj\varphi\)

一道 \(\varphi()\) 的题。

对于一个合法的数对,设它为 \((a*m,b*m)\) 则 \(((a+b)*m)|a*b*m^2\),所以 \((a+b)|a*b\),因为 \(\gcd(a,b)=1\),所以 \(a+b|m\)

那么设 \(a+b=k\),且 \(k*m\le n\),那么 \(k\) 最大为 \(\sqrt n\),所以枚举 \(k\) 即可,对于每个 \(k\),有 \(\rm \frac{n}{k^2}\) 个 \(m\) ,同时又有 \(\rm \varphi(k)\) 对 \(a+b\)

Code:
#include<bits/stdc++.h>
#define ri register signed
#define p(i) ++i
using namespace std;
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf,OPUT[100];
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?(-1):*p1++
template<typename T>inline void read(T &x) {
ri f=1;x=0;register char ch=gc();
while(!isdigit(ch)) {if (ch=='-') f=0;ch=gc();}
while(isdigit(ch)) {x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
x=f?x:-x;
}
template<typename T>inline void print(T x) {
if (x<0) putchar('-'),x=-x;
if (!x) return putchar('0'),(void)putchar('\n');
ri cnt(0);
while(x) OPUT[p(cnt)]=x%10,x/=10;
for (ri i(cnt);i;--i) putchar(OPUT[i]^48);
return (void)putchar('\n');
}
}
using IO::read;using IO::print;
namespace nanfeng{
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
typedef long long ll;
static const int N=1e7+7;
int vis[N],phi[N],prim[N],nm[N],cnt,sn;
ll n,ans;
inline void Getphi(int n) {
for (ri i(2);i<=n;p(i)) {
if (!vis[i]) phi[i]=i-1,vis[prim[p(cnt)]=i]=i;
for (ri j(1);j<=cnt&&prim[j]*i<=n;p(j)) {
int nw=i*prim[j];
vis[nw]=prim[j];
if (vis[i]==prim[j]) {
phi[nw]=phi[i]*prim[j];
break;
}
else phi[nw]=phi[i]*(prim[j]-1);
}
}
}
inline int main() {
//FI=freopen("nanfeng.in","r",stdin);
//FO=freopen("nf.out","w",stdout);
read(n);
sn=sqrt(n);
Getphi(sn);
for (ri i(2);i<=sn;p(i)) ans+=(ll)phi[i]*(n/i/i);
print(ans);
return 0;
}
}
int main() {return nanfeng::main();}

最新文章

  1. 也说php从mysql数据库通过服务器端json返回数据出现乱码问题
  2. java 24 - 5 GUI之 鼠标移动改变窗体背景颜色
  3. Java 集合系列11之 Hashtable详细介绍(源码解析)和使用示例
  4. jQuery对input select操作小结
  5. mars android视频学习笔记一:Activity生命周期
  6. 仿今日头条最强顶部导航指示器,支持6种模式-b
  7. ios 阻止GDB依附
  8. 关于org.openqa.selenium.ElementNotVisibleException
  9. POJ2299 Ultra-QuickSort(归并排序求逆序数)
  10. PYTHON线程知识再研习E---条件变量同步Condition
  11. ubuntu 16.04 Ubuntu 安装GDebi,从而安装deb文件
  12. K60 启动过程分析
  13. 一个小团队TDD游戏及实践
  14. 自己整理的openresty安装步骤
  15. 发短信utils
  16. Angular开发实践(八): 使用ng-content进行组件内容投射
  17. [Ahoi2005]LANE 航线规划
  18. js 实现论坛评论模块原理
  19. 【MIT-6.824】Lab 1: MapReduce
  20. P2659 美丽的序列

热门文章

  1. 利用漏洞破解win7密码
  2. ArchLinux安装步骤(一)
  3. [刘阳Java]_程序员Java编程进阶的5个注意点,别编程两三年还是增删改查
  4. nodejs 文本逐行读写功能的实现
  5. UIAutomator2 之 计算机积极拒绝
  6. [WinError 10013]以一种访问权限不允许的方式做了一个访问套接字的尝试
  7. Idea快捷键 累积大全
  8. Centos 7 安装mysql5.7 nginx tomcat
  9. 某学院m3u8视频解密获取分析实战分享
  10. 终于彻底搞清楚了spin-lock 之一次CPU问题定位过程总结