hdu 6053 trick gcd 容斥

http://acm.hdu.edu.cn/showproblem.php?pid=6053

题意：给定一个数组，我们定义一个新的数组b满足bi<ai 求满足gcd(b1,b2....bn)>=2的数组b的个数

题解：利用容斥定理。我们先定义一个集合f(x)表示gcd(b1,b2...bn)为x倍数的个数（x为质数），我们在定义一个数mi为数组中的最小值，那么集合{f(2)Uf(3)....f(n)}就是我们想要的答案。f(x)=(a1/x)*(a2/x)*.....(ai/x),直接累加肯定是有重复的，我们得用容斥定理筛一下，如果x是奇数个不同素数因子的乘积最后的结果要加上f(x);如果x为偶数个不同素数因子的乘积，最后的结果要减去f(x)，其他情况贡献为0。是不是和莫比乌斯函数的情况正好相反？这里筛值的时候，用0(n)求到的莫比乌斯函数筛时间复杂度还是可以的。光这样还是不够，因为数组的长度为1e5,我们求f(x)的时候也得优化，怎么优化呢。我们把用一个权值数组把a[i]的值离散上去，为啥要用权值的形式存放a数组的值呢，我们把权值数组分成x段，每段的贡献由1开始递增到x（这段比较抽象，具体看下代码）

ac代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
typedef long long ll;
const int mod=1e9+;
ll mu[];
int prime[],vis[];
ll a[],sum[];
// 比较大的数组还是定义在外面比较好
ll qpow(ll a,ll b)
{
    ll f=;
    while(b)
    {
        if(b%==) f=(f*a)%mod;
        a=(a*a)%mod;
        b/=;
    }
    return f;
}

void init()
{
    mu[]=;
    memset(mu,,sizeof(mu));
    memset(prime,,sizeof(prime));
    memset(vis,,sizeof(vis));
    int ret=;
    for(int i=;i<;i++)
    {
        if(!vis[i])
        {
           prime[ret++]=i;
           mu[i]=-1LL;
        }
        for(int j=; j<ret && i*prime[j] < ;j++)
        {
            int temp=i*prime[j];
            vis[temp]=;
            if(i%prime[j]) mu[temp]=-mu[i];
            else
            {
                mu[temp]=;
                break;
            }
        }
    }
}

int main()
{
    int t;
    scanf("%d",&t);
    init();

    int Case=;
    while(t--)
    {
        ll n;
        cin>>n;
        ll mi=;
        memset(sum,,sizeof(sum));
        for(int i=;i<=n;i++)
        {
            scanf("%d",&a[i]);
            mi=min(a[i],mi);
            sum[a[i]]++;
        }
        for(int i=;i<=;i++) sum[i]+=sum[i-];
        ll zz=;
        for(int i=;i<=mi;i++)
        {
            ll ans=;
            if(mu[i]==) continue;
            for(int j=;j*i<=;j++)// 平铺分段的思想  枚举贡献的思想吧，对于f(x)来说，x/n相同的值比较多，这样就可以把问题的规模变小 这个思维比较常见
            {
                ans=(ans*qpow(j,sum[i*(j+)-]-sum[j*i-])%mod)%mod;
            }
            zz=(zz-mu[i]*ans%mod+mod)%mod;// ! 取模的时候 如果有减法 要注意
        }
        printf("Case #%d: ",++Case);
        cout<<zz<<endl;
    }
    return ;
}