hdu 6053 TrickGCD （莫比乌斯反演）

题目链接：哆啦A梦传送门

题意：给出长度为n的序列a，找出长度也为n的序列b，
满足$b_i\le a_i,gcd(b_1,b_2,...,b_n)\ge 2$，问：有多少种不同的序列b？

题解：

我们先设：
$f(d)$为满足$gcd(b_1,b_2,...,b_n)=d$的方案数。

那么我们根据莫比乌斯第二描述：
$$F(n)=\sum _{n|d}f(d)$$
$$f(n)=\sum _{n|d}u(\frac{d}{n})\ast F(d)$$

易得：$F(n)$为满足$d|gcd(b_1,b_2,...,b_n)$的方案数。

参考博客：(https://blog.csdn.net/V5ZSQ/article/details/76381584)

即：$\begin{array}{r}F(d)=\prod _{i=1}^{mi}\lfloor \frac{a_i}{d}\rfloor \end{array}$ ，显然这样求复杂度有点大，我们再设$sum[i]表示a_i\le i的总数$，对于每个d，
它在每一段$a\in (i\ast d-1,(i+1)\ast d-1)$ ($i\ast d\le 1e5$)的值除d的值都是一样的，都为i。那么我们只需算下它们的幂$i^{sum[(i+1)\ast d-1]-sum[i\ast d-1]}$就好了。

那么此题的结果就为：

$\begin{array}{rl}ans& =F(1)-f(1)\\ & =F(1)-\sum _{d=1}^{mi}F(d)\\ & =-\sum _{d=2}^{mi}F(d)\end{array}$

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N=1e5+10;
const LL mod=1e9+7;

int prime[N],tot,mu[N];

bool vis[N];

void init()
{
    mu[1]=1;

    for(int i=2;i<N;i++)
    {
        if(!vis[i]){
            prime[++tot]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=tot&&i*prime[j]<N;j++){
            int x=i*prime[j];

            vis[x]=1;
            if(i%prime[j]==0){
                mu[x]=0;
                break;
            }
            mu[x]=-mu[i];

        }

    }

}


int sum[N];

LL F[N];

LL fast(LL x,LL y)
{
    LL re=1;
    while(y)
    {
        if(y&1) re=re*x%mod;
        x=x*x%mod;
        y>>=1;
    }
    return re;
}


int main()
{

    init();
    int n,T=0;
    int ncase;
    scanf("%d",&ncase);
    while(ncase--){

        int n;
        scanf("%d",&n);
        int mi=1e5,item;
        memset(sum,0,sizeof(sum));
        for(int i=1;i<=n;i++){
            scanf("%d",&item);
            mi=min(item,mi);
            sum[item]++;
        }

        for(int i=1;i<=1e5;i++){
            sum[i]+=sum[i-1];
        }

        for(int d=2;d<=mi;d++)
        {
            F[d]=1LL;
            for(int i=1;i*d<=100000;i++){
                F[d]=F[d]*fast(i,sum[min((i+1)*d-1,100000)]-sum[i*d-1])%mod;
//            printf("%d %d %d %d\n",i,d,mi,F[d]);
            }

        }


        LL ans=0;
        for(int d=2;d<=mi;d++)
        {
            ans=ans-mu[d]*F[d]+mod;
            ans%=mod;
        }
        printf("Case #%d: %lld\n",++T,ans);




    }
    return 0;
}