p2257 yy的gcd 莫比乌斯反演入门题

a n s = ∑ d = 1 m i n ( n , m ) ⌊ n / d ⌋ ⌊ m / d ⌋ ( ∑ x ∣ d μ ( d / x ) ) ans=\sum_{d=1}^{min(n,m)}\lfloor n/d \rfloor \lfloor m/d \rfloor (\sum_{x|d}\mu(d/x)) ans=d=1min(n,m)n/dm/d(xdμ(d/x))
就是前面就是整除分块,后面就是前缀和预处理解决多组询问
式子是由莫比乌斯反演推出来的,不难,但LaTex太费事了,以后拍照上传把

看着大佬推的式子 敲代码a了,感觉还得自己练一下推式子的能力

#include 
#include 
#include 

using namespace std;
typedef long long ll;
const int maxn=1e7+7;
int prime[maxn],mu[maxn],g[maxn],tot;
ll sum[maxn];
bool p[maxn];
int T,n,m;

void get_mu(){
    mu[1]=1;
    for(int i=2;im)swap(n,m);ll ans=0;
        for(int i=1,gi;i<=n;i=gi+1){
            gi=min(n/(n/i),m/(m/i));
            ans+=(ll)(n/i)*(m/i)*(sum[gi]-sum[i-1]);
        }
        printf("%lld\n",ans);
    }
    return 0;
}

你可能感兴趣的:(数论)