[莫比乌斯反演] 51Nod 1355 斐波那契的最小公倍数

现在我自己推不出来，看题解的。
http://blog.csdn.net/ZLH_HHHH/article/details/77587832
大概是考虑用经典的容斥写出答案的式子，考虑每个数对答案的贡献，然后需要反演一下算出贡献的式子。

#include<cstdio>#include<algorithm>using namespace std;typedef long long LL;const int N=1000000,maxv=1000005,MOD=1000000007;int n,fib[maxv],cnt[maxv],ctb[maxv],a[maxv],mu[maxv],p[maxv];LL ans;bool vis[maxv];void get_mu(){    mu[1]=1;    for(int i=2;i<=N;i++){        if(!vis[i]) p[++p[0]]=i, mu[i]=-1;        for(int j=1;j<=p[0]&&(LL)i*p[j]<=N;j++){            vis[i*p[j]]=true;            if(i%p[j]==0){ mu[i*p[j]]=0; break; }            mu[i*p[j]]=-mu[i];        }    }}int Pow(LL a,LL b){    LL res=1;    for(;b;b>>=1,a=a*a%MOD) if(b&1) res=(res*a)%MOD;    return res; }LL get(LL a,LL b){    if(b>=0) return Pow(a,b);    return Pow(Pow(a,-b),MOD-2);}int main(){    scanf("%d",&n);    for(int i=1;i<=n;i++){        int x; scanf("%d",&x);        if(x>=3) a[x]++;    }     for(int i=3;i<=N;i++)     for(int j=1;i*j<=N;j++) cnt[i]+=a[i*j];    get_mu();    for(int i=3;i<=N;i++)     for(int d=1;i*d<=N;d++) if(cnt[i*d]) ctb[i]+=mu[d];        fib[1]=1; for(int i=2;i<=N;i++) fib[i]=(fib[i-1]+fib[i-2])%MOD;    ans=1; for(int i=3;i<=N;i++)     ans=(ans*get(fib[i],ctb[i]))%MOD;    printf("%lld\n",ans);    return 0;}