51nod1355 斐波那契的最小公倍数

最新推荐文章于 2021-01-25 12:49:25 发布

sdfzyhx

最新推荐文章于 2021-01-25 12:49:25 发布

阅读量443

点赞数

分类专栏：数学其他oj

本文链接：https://blog.csdn.net/sdfzyhx/article/details/74370399

版权

数学同时被 2 个专栏收录

229 篇文章 0 订阅

订阅专栏

其他oj

66 篇文章 0 订阅

订阅专栏

首先我们可以知道 $f_{n+m}=f_{n-1}f_m+f_nf_{m+1}$ ，这一点可以用数学归纳法证明。于是我们有 $\gcd(f_{n+m},f_m)=\gcd(f_n,f_m)$ 。进一步地， $\gcd(f_n,f_m)=f_{\gcd(n,m)}$ 。
令 $g=f\times 1$ ，那么根据最大最小值反演我们有

lcm (f S) = \prod T \in S gcd (f T) (- 1) | T | + 1 = \prod T \in S f gcd (T) (- 1) | T | + 1 = \prod d g d \sum T \in S, d | gcd (T) (- 1) | T | + 1 = \prod \exists s i, d | s i g d

$\begin{align} \text{lcm}(f_S) & =\prod_{T\in S}\gcd(f_T)^{(-1)^{|T|+1}} \\ & = \prod_{T\in S}{f_{\gcd(T)}}^{(-1)^{|T|+1}} \\ & =\prod_d {g_d}^{\sum_{T\in S,d|\gcd(T)}(-1)^{|T|+1}} \\ & = \prod_{\exists s_i,d|s_i} g_d \end{align}$
复杂度

O(nm−−√+mlogm) $O(n\sqrt m+m\log m)$ 。

#include<cstdio>
#include<algorithm>
using namespace std;
#define LL long long
const int maxm=1000010,maxn=50010,p=1000000007;
int f[maxm],g[maxm],a[maxn],vis[maxm],n,m;
int inc(int x,int y)
{
    x+=y;
    return x>=p?x-p:x;
}
int pow(int b,int k)
{
    int r=1;
    for (;k;k>>=1,b=(LL)b*b%p)
        if (k&1) r=(LL)r*b%p;
    return r;
}
int main()
{
    int x,ans=1;
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d",&a[i]),m=max(m,a[i]);
    f[1]=g[1]=1;
    for (int i=2;i<=m;i++) f[i]=g[i]=inc(f[i-1],f[i-2]);
    for (int i=1;i<=m;i++)
    {
        x=pow(g[i],p-2);
        for (int j=i*2;j<=m;j+=i)
            g[j]=(LL)g[j]*x%p;
    }
    for (int i=1;i<=n;i++)
        for (int j=1;(LL)j*j<=a[i];j++)
            if (a[i]%j==0) vis[j]=vis[a[i]/j]=1;
    for (int i=1;i<=m;i++)
        if (vis[i]) ans=(LL)ans*g[i]%p;
    printf("%d\n",ans);
}