HDU 6053 TrickGCD

HDU 6053 TrickGCD

原题连接：
http://acm.hdu.edu.cn/showproblem.php?pid=6053
（对于不了解，莫比乌斯反演的同学不需要纠结。直接用结论就好了）

莫比乌斯反演公式

若对于 定义域 为正整数，值域为复数子集 的函数$f$：

定义它的和函数

$$F(n)=\sum_{d|n}f(d)=\sum_{d|n}f(\frac{n}{d})$$

则对于任意正整数n：

$$f(n)=\sum_{d|n}\mu (d)F(\frac{n}{d})=\sum_{d|n}\mu (\frac{n}{d})F(d)$$

同时 。也有一种无穷形式的反演公式：

$$F(n)=\sum_{n|d}f(d)$$

这里，d取n的倍数。 所以这是一个无穷和式子。有：

$$f(n)=\sum_{n|d}\mu (\frac{d}{n})F(d)$$

反演的形式很多。都是很巧妙的（可以参阅组合数学，有更一般的反演与容斥）。

上面的的莫比乌斯反演的定义是很强的。- -！一切算术函数

————————————————————————————————————————————————

对于这一题。$gcd(B_l,B_{l+1},...B_{r})>=2,对于任意的l<=r,l,r\in[1,n]$

设：$gcd(B)=n$的B数列个数为$f(n)$

设

$$F(n)=\sum_{n|d}^{n}f(d)$$

由 第二个莫比乌斯反演得：

$$f(1)=\sum_{d>=1}\mu (d)F(d)$$

这里$F(n)$就是$gcd=n$的倍数的所有方案数量。

虽然 得到的式子是一个无穷和式子。

但。当$d>min(A_1,A_2,A_3,...A_n)$时，$F(d)=0$

我么可以通过预处理与快速幂快速计算出F函数

对于$F(k)$那么任意B都应该含有k这个因子。

那么$B_i$可选择的方案数就是 $\lfloor\frac{A_i}{k}\rfloor$

那么乘法原理得：

$$F(k)=\prod_{i=1}^n\big\lfloor\frac{A_i}{k}\big\rfloor$$

令$m=min(A_1,A_2,..A_n)$

令$h=max(A_1,A_2,...A_n)$

暴力算所有的$F$从$1$到$m$的时间就是$O(nm)$

设$cnt[u]$是$A$中$A[i]<=u$的元素个数。

对于$A$中满足$\big\lfloor\frac{A_i}{k}\big\rfloor=t$的所有的元素个数为：

$$cnt[(t+1)k-1]-cnt[tk-1]$$

那么 $A$中与$k$的商等于$t$ 的元素对$F(k)$的贡献为：

$$t^{cnt[(t+1)k-1]-cnt[tk-1]}$$

那么计算$F(k)$的复杂度为$O(\frac{h}{k}log(cnt[(t+1)k-1]-cnt[tk-1]))$

计算出$1$到$m$的$F(k)$所用时间为：

$O(h\sum_{k=1}^m\frac{log(cnt)}{k})=O(h*logm*logh)<O(10^5log^2(10^5))$

利用反演公式：

$$f(1)=\sum_{d>=1}\mu(d)F(d)$$

计算出$f(1)$

$$answer=F(1)-f(1)$$

#include <algorithm>
#include <string.h>
#include <stdio.h>
#include <set>
#include <math.h>
#define MAXN 100005
using namespace std;

typedef long long LL;

const LL mod=1e9+7;

int cnt[MAXN];
LL f[MAXN];
const int N=MAXN;
int mu[N];
void getMu()
{
    for(int i=1;i<=N;i++)
    {
        int target=i==1?1:0;
        int delta=target-mu[i];
        mu[i]=delta;
        for(int j=i+i;j<=N;j+=i)
        {
            mu[j]+=delta;
        }
    }
}

LL pow(LL a,int b)
{
    LL tmp=1;
    while(b)
    {
        if(b&1)
            tmp=tmp*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return tmp;
}

int main ()
{
    int T,ca=1;
    getMu();
    scanf("%d",&T);
    while(T--)
    {
        int A;
        LL ans=0;
        memset(cnt,0,sizeof cnt);
        int n,m=MAXN,ma=-1;
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&A);
            m=min(A,m);
            ma=max(A,ma);
            cnt[A]++;
        }
        for(int i=1;i<=ma;i++) cnt[i]+=cnt[i-1];
        for(int k=1;k<=m;k++)
        {
            f[k]=1;
            for(int l=-1,j=k-1,t=1; ;t++)
            {
                j+=k;
                l+=k;
                if(j>ma)j=ma;
                f[k]*=pow((LL)t,cnt[j]-cnt[l]);
                f[k]%=mod;
                if(j==ma)break;
            }
        }
        for(int i=1;i<=m;i++)
            ans=(ans+f[i]*mu[i]+mod)%mod;
        printf("Case #%d: %lld\n",ca++,(f[1]-ans+mod)%mod);
    }
    return 0;
}