神炎皇——解题报告

题目链接：

http://132.232.5.128:666/problem/110

题目大意：

对于一个整数对 $(a,b)$，若满足 $a+b\le n$且 $a+b$是$a\times b$的因子，则称为神奇的数对。问这样的数对共有多少。

题目分析：

1.整除一类的题，我们很快能想到的就是$GCD$，我们可以利用$gcd$将式子化简，$a+b|a\times b$，设$d=gcd(a,b)$，$a^{\prime }=\frac{a}{d}$，$b^{\prime }=\frac{b}{d}$。那么$(a^{\prime }+b^{\prime })\times d|a^{\prime }\times b^{\prime }\times d\times d$，所以 $(a^{\prime }+b^{\prime })|a^{\prime }\times b^{\prime }\times d$，而根据辗转相除法，$gcd=(a^{\prime }+b^{\prime },a^{\prime })=gcd(a^{\prime }+b^{\prime },b^{\prime })=gcd(a^{\prime }+b^{\prime },a^{\prime }\times b^{\prime })=1$。所以我们得到了$(a^{\prime }+b^{\prime })|d$。
2.因为$(a^{\prime }+b^{\prime })\times d\le n$，所以$(a^{\prime }+b^{\prime })\le \sqrt{n}$。若$(a^{\prime }+b^{\prime })=k$，那么合法的$d$就有$n/k^2$个。而合法的$(a^{\prime },b^{\prime })$就有$\varphi (k)$个，所以答案就要加上$n/k^2\times \varphi (k)$。
3.于是我们枚举一下所有的$k$就行了。

正解程序：

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath> 

using namespace std;
typedef long long ll;
const ll maxn=10000010;
ll n,ans;
ll count1=0,prime[maxn],phi[maxn];
bool vis[maxn];
int main()
{
	scanf("%lld",&n);
	ll Max=sqrt(n); 
	for(ll i=2;i<=Max;i++)
	{
		if(!vis[i])
		{
			prime[++count1]=i;
			phi[i]=i-1;
		}
		for(ll j=1;j<=count1 && i*prime[j]<=Max;j++)
		{
			vis[i*prime[j]]=true;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			phi[i*prime[j]]=phi[i]*(prime[j]-1);
		}
		ans+=n/i/i*phi[i];
	}
	printf("%lld\n",ans);
	
	return 0;
}