版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/sugar_free_mint/article/details/87256433
题目
分析
原式=
重点就是
那么原式
那么就可以用整除分块求解,然而也快不了多少
代码
#include <cstdio>
#define rr register
#define min(a,b) ((a)<(b)?(a):(b))
using namespace std;
typedef long long ll;
ll n,m,t,ans,phi[100001],v[100001],prime[10001],cnt;
inline void init(int N){
phi[1]=1;
for (rr int i=2;i<=N;++i){
if (!v[i]) phi[i]=(prime[++cnt]=v[i]=i)-1;
for (rr int j=1;j<=cnt&&prime[j]*i<=N;++j){
v[i*prime[j]]=prime[j];
if (i%prime[j]) phi[i*prime[j]]=phi[i]*(prime[j]-1);
else {phi[i*prime[j]]=phi[i]*prime[j]; break;}
}
}
for (rr int i=2;i<=N;++i) phi[i]+=phi[i-1];
}
signed main(){
scanf("%lld%lld",&n,&m);
init(t=min(n,m));
for (rr ll l=1,r;l<=t;l=r+1){
r=min(n/(n/l),m/(m/l));
ans+=(phi[r]-phi[l-1])*(n/l)*(m/l);
}
printf("%lld",2*ans-n*m);
return 0;
}