51nod 1355 斐波那契的最小公倍数 莫比乌斯反演+数学

题意

F(0) = 0 F(1) = 1
F(n) = F(n-1) + F(n-2)
给出n个正整数a1, a2,…… an，求对应的斐波那契数的最小公倍数，由于数字很大，输出Mod 1000000007的结果即可。
2 <= N <= 50000，1 <= ai <= 1000000。

分析

跪zyz大爷写的题解。

首先有个结论就是

$$gcd(f_n,f_m)=f_{gcd(n,m)}$$
证明的话可以根据 $f_{n+m}=f_{n}*f_{m-1}+f_{n+1}*f_m$来推导。

还有一个结论就是

$$lcm\{S\}=\prod\limits_{T\subseteq S,T\not=\phi}gcd(T)^{(-1)^{|T|+1}}$$
证明的话可以考虑知道每个子集的最小值通过容斥来求全集的最大值。
$Max(S)=\sum\limits_{T\subseteq S,T\not=\phi}(-1)^{|T|+1}*Min(T)$
为什么这样是对的呢？如果真的无法把这个东西看做容斥的话可以这样想：
我们把全集从大到小排序后，第n个数出现的次数就是 $\sum\limits_{i=0}^{n-1}(-1)^i*C_{n-1}^i$
注意到这条式子只有当n=0时为1否则为0，因此得出来的是最大值。

从而可以得到

$$lcm(f_S)=\prod\limits_{T\subseteq S,T\not=\phi}gcd(f_T)^{(-1)^{|T|+1}}$$
$$=\prod\limits_{T\subseteq S,T\not=\phi}(f_{gcd(T)})^{(-1)^{|T|+1}}$$
然后我们构造数列 $\{g_d\}$满足 $f_n=\prod\limits_{d|n}g_d$
那么有
$$lcm(f_S)=\prod\limits_{T\subseteq S,T\not=\phi}(\prod\limits_{d|gcd(T)}g_d)^{(-1)^{|T|+1}}$$
$$=\prod\limits_d(g_d)^{\sum\limits_{T\subseteq S,T\not=\phi,d|gcd(T)}(-1)^{|T|+1}}$$
注意到其指数上的式子 $$\sum\limits_{T\subseteq S,T\not=\phi,d|gcd(T)}(-1)^{|T|+1}$$
当存在 $x\in S$满足 $d|x$时该式子为1，否则为0。
那么有
$$lcm(f_S)=\prod\limits_{\exists x\in S,d|x}g_d$$
在求出 $g_d$后这条式子显然可以在 $O(nlogn)$内计算完毕。
求 $g_d$的话可以通过对定义式进行变形得到 $$g_n=f_n*(\prod_{d|n,d\not= n}g_d)^{-1}$$
那么 $g_d$显然也可以在 $O(nlogn)$内预处理完毕，问题解决。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int N=1000005;
const int MOD=1000000007;

int n,f[N];
bool vis[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

int ksm(int x,int y)
{
    int ans=1;
    while (y)
    {
        if (y&1) ans=(LL)ans*x%MOD;
        x=(LL)x*x%MOD;y>>=1;
    }
    return ans;
}

int main()
{
    n=read();int mx=0;
    for (int i=1;i<=n;i++)
    {
        int x=read();
        mx=max(mx,x);vis[x]=1;
    }
    for (int i=1;i<=mx;i++)
        for (int j=i;j<=mx;j+=i)
            vis[i]|=vis[j];
    f[1]=f[2]=1;
    for (int i=3;i<=mx;i++) f[i]=f[i-1]+f[i-2],f[i]-=f[i]>=MOD?MOD:0;
    int ans=1;
    for (int i=1;i<=mx;i++)
    {
        int ny=ksm(f[i],MOD-2);
        for (int j=i*2;j<=mx;j+=i) f[j]=(LL)f[j]*ny%MOD;
        if (vis[i]) ans=(LL)ans*f[i]%MOD;
    }
    printf("%d",ans);
    return 0;
}