題目

神炎皇烏利亞很喜歡數對，他想找到神奇的數對。
對于一個整數對(a,b)，若滿足a+b<=n且a+b是ab的因子，則成為神奇的數對。請問這樣的數對共有多少呢？

分析

設\(gcd(a,b)=d,a'd=a,b'd=b\)
那么\(a'+b'|a'b'd\)
因為\(gcd(a',b')=1\)
所以\(a'+b'|d\)。
又因為\((a'+b')d<=n\)
則\(a'+b'=\sqrt n\)
枚舉\(a'+b'=i\)
\(d就有\dfrac{n}{i^2}種情況\)
因為\(gcd(a',b')=gcd(a'+b',a')\)
所以\(a'和b'又有\varphi(i)種\)
線篩求\(\varphi()\)，時間復雜度\(O(\sqrt n)\)

#include <cmath> #include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <queue> const int maxlongint=2147483647; const int mo=1000000007; const int N=10000005; using namespace std; long long ans; long long n,qn,r,phi[N],p[N]; bool bz[N]; long long gcd(long long x,long long y) {for(;y;){r=x%y;x=y;y=r;}return x; } int main() {scanf("%lld",&n);qn=sqrt(n); phi[1]=1;memset(bz,true,sizeof(bz));for(long long i=2;i<=qn;i++){if(bz[i]) { bz[i]=false; p[++p[0]]=i; phi[i]=i-1; } for(long long j=1;j<=p[0] && i*p[j]<=qn;j++) {bz[i*p[j]]=false;if(i%p[j]) phi[i*p[j]]=phi[i]*(p[j]-1);else {phi[i*p[j]]=phi[i]*p[j]; break;} }ans+=(long long)phi[i]*(long long)(n/i/i); } printf("%lld",ans); }

轉載于:https://www.cnblogs.com/chen1352/p/9066604.html

總結

以上是生活随笔為你收集整理的【NOIP2017提高组模拟12.10】神炎皇的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。