标签：des 约数 space 意义 inpu www 关系 sea 表示
对于100%的数据满足0 < A, B < 10^ 6
好久不写数学了...水一发...
$\sum_ {i=1}^{n} \sum_ {j=1}^{m} [gcd(i,j)==k]$
$n=n/k,m=m/k(n<=m)$
$\sum_ {i=1}^{n} \sum_ {j=1}^{m} [gcd(i,j)==1]$
$\sum_ {i=1}^{n} \sum_ {j=1}^{m} \sum_ {d\mid gcd(i,j)} μ(d)$
$\sum_ {d=1}^{n} μ(d) \left\lfloor \frac{n}{d} \right \rfloor \left\lfloor \frac{m}{d} \right \rfloor$
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
//by NeighThorn
#define int long long
using namespace std;
const int maxn=1000000+5;
int n,m,d,k,cnt,mu[maxn],pri[maxn],vis[maxn];
inline void prework(void){
mu[1]=1;
for(int i=2;i<=1000000;i++){
if(!vis[i])
pri[++cnt]=i,mu[i]=-1;
for(int j=1;j<=cnt&&pri[j]*i<=1000000;j++){
vis[pri[j]*i]=1;
if(i%pri[j]==0){
mu[i*pri[j]]=0;
break;
}
mu[i*pri[j]]=-mu[i];
}
mu[i]+=mu[i-1];
}
}
signed main(void){
scanf("%lld%lld%lld",&n,&m,&d);
n/=d,m/=d;if(n>m) swap(n,m);prework();
int ans=0,r;
for(int i=1;i<=n;i=r+1){
r=min(n/(n/i),m/(m/i));
ans+=(mu[r]-mu[i-1])*(n/i)*(m/i);
}
printf("%lld\n",ans);
return 0;
}
By NeighThorn
标签：des 约数 space 意义 inpu www 关系 sea 表示
原文地址：http://www.cnblogs.com/neighthorn/p/6414862.html
