Codeforces900 D. Unusual Sequences（莫比乌斯反演）

题意：

给定x,y，问有多少种序列a,
满足a[i]>=1，gcd(a[])=x，且sum(a[])=y.

数据范围：x,y<=1e9

解法：

$$\begin{gathered} y\%x!=0时肯定无解 \\ 问题可以转换为满足和为k=y/x,gcd=1的序列数量 \\ 设g(k)为满足和为k的方案数 \\ 因为序列的长度可以不定,可以看作是k个1,在k-1个空位插入0到k-1个隔板 \\ 因此g(k)=\sum _{i=0}^{k-1}C(k-1,i)=2^{k-1} \\ g(k)没有考虑gcd,里面包含gcd=1到gcd=k的所有方案. \\ 设f(x)为gcd=1,和为y/x的方案数, \\ 因为对于g(k),所有可能的序列的gcd一定是k的约数 \\ 因此g(k)=\sum _{d|k}f(\frac{k}{d}). \\ 为什么是f(\frac{k}{d})呢?f(\frac{k}{d})的和是\frac{k}{d}不是k呀？ \\ 这是因为和为\frac{k}{d},gcd=1的序列,每个数都乘上d之后 \\ 就变成和为k，gcd=d的序列了 \\ 反演得:f(k)=\sum _{d|k}\mu (d)g(\frac{k}{d}) \\ O(sq(k))枚举因子d即可. \\ 预处理质数，然后\mu (d)就可以O(sq(d)/log(d))计算. \end{gathered}$$

code：

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxm=1e6+5;
const int mod=1e9+7;
int notprime[maxm];
int prime[maxm],cnt;
int x,y;
int ppow(int a,int b,int mod){
    int ans=1%mod;a%=mod;
    while(b){
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return ans;
}
int mu(int x){
    if(x==1)return 1;
    int ans=1;
    for(int j=0;j<cnt&&prime[j]*prime[j]<=x;j++){
        if(x%prime[j]==0){
            int cnt=0;
            while(x%prime[j]==0){
                x/=prime[j];
                cnt++;
            }
            if(cnt!=1)return 0;
            ans=-ans;
        }
    }
    if(x!=1)ans=-ans;
    return ans;
}
int g(int x){
    return ppow(2,x-1,mod);
}
void init(){
    for(int i=2;i<maxm;i++){
        if(!notprime[i]){
            prime[cnt++]=i;
        }
        for(int j=0;j<cnt;j++){
        	if(prime[j]*i>=maxm)break;
            notprime[prime[j]*i]=1;//素数的倍数标记为非素数
            if(i%prime[j]==0)break;//看外面解释
        }
    }
}
signed main(){
    init();
    cin>>x>>y;
    if(y%x!=0){
        cout<<0<<endl;
        return 0;
    }
    int k=y/x;
    int ans=0;
    for(int i=1;i*i<=k;i++){
        if(k%i==0){
            ans+=mu(i)*g(k/i);
            if(k/i!=i){
                ans+=mu(k/i)*g(i);
            }
            ans%=mod;
        }
    }
    ans=(ans%mod+mod)%mod;
    cout<<ans<<endl;
    return 0;
}

---