【51nod1355】斐波那契的最小公倍数（Min-Max容斥）（莫比乌斯反演）

传送门

---

题解：

首先所有相同的数只保留一个。

然后由于$gcd$和$lcm$可以表示成各个质因子的次数取$\max$和$\min$的情况，我们由$Min-Max$容斥可以得到：

$$lcm(S)=\prod _{T\subset S}gcd(T{)}^{-1^{|T|+1}}$$

注意上面的式子那里有双重指数。

斐波那契数有一个性质：$gcd(f_n,f_m)=f_{gcd(n,m)}$

这里给出一个证明，由于$f_n=f_{n-1}+f_{n-2}$，归纳可知：$gcd(f_n,f_{n-1})=1$，同时有一个递推的性质：$f_{n+m}=f_{m+1}\ast f_n+f_m\ast f_{n-1}$，由以上两个结论可以方便地推出上面那个性质。

$$gcd(f_{n+m},f_n)=gcd(f_{n-1}\cdot f_m+f_n\cdot f_{m+1},f_n)=gcd(f_m,f_n)$$

下标上的迭代会在到达$gcd(n,m)$时候停止。归纳可以知道$gcd(f_n,f_m)=f_{gcd(n,m)}$

对我们要求的东西进行转化：

$$\begin{array}{rlrl}Ans& =& & lcm(f_S)\\ & =& & \prod _{T\subset S}gcd(f_T{)}^{-1^{|T|+1}}\\ & =& & \prod _{T\subset S}{f}_{gcd(T)}^{-1^{|T|+1}}\end{array}$$

感觉这个玩意不好算，观察到数据范围只有$1e6$，所以我们感觉可以上$O(A\log A)$的反演来一下。

考虑构造函数$g_n$满足：$f_n={\prod }_{d\mid n}{g}_d$

于是可以继续化简：$$\begin{array}{rlrl}Ans& =& & \prod _{T\subset S}\prod _{d\mid gcd(T)}{g}_d^{-1^{|T|+1}}\\ & =& & \prod _d{g}_d^{{\sum }_{T\subset S}[d\mid gcd(T)](-1^{|T|+1})}\end{array}$$

我们观察一下上面这个式子中的指数：$$\sum _{T\subset S}[d\mid gcd(T)](-1^{|T|+1})$$

看上去不好处理，但是由这种集合的性质（或许应该称为$gcd$的性质，但是总得来说还是来自容斥）我们知道：

$$\sum _{T\subset S}[d\mid gcd(T)](-1^{|T|+1})=\{\begin{array}{rlrl}1& & & \exists a\in S,[d\mid a]\\ 0& & & otherwise\end{array}$$

于是我们可以$O(n\log n)$处理$g$数组，然后随便做了。

但是实际上我们没必要算出$g$数组，我们直接算出在所有$g$的乘积中，$f_i$出现了多少次就行了。

这个也可以在调和级数时间内算出来。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

using std::cerr;
using std::cout;

cs int mod=1e9+7;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){ll r=(ll)a*b;return r>=mod?r%mod:r;}
inline void Inc(int &a,int b){(a+=b)>=mod&&(a-=mod);}
inline void Dec(int &a,int b){(a-=b)<0&&(a+=mod);}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
} 

cs int N=1e6+6;
int lim;
int f[N],cnt[N];

int n,ans=1;
bool a[N];

signed main(){
//	freopen("fiblcm.in","r",stdin);
	scanf("%d",&n);
	for(int re i=1;i<=n;++i){
		int x;scanf("%d",&x);
		a[x]=true;
		lim=std::max(lim,x);
	}
	f[1]=1;for(int re i=2;i<=lim;++i)f[i]=add(f[i-1],f[i-2]);
	for(int re i=lim;i>2;--i){
		for(int re j=i;j<=lim&&!a[i];j+=i)a[i]|=a[j];
		if(!a[i])continue;
		cnt[i]=1;
		for(int re j=i+i;j<=lim;j+=i)cnt[i]-=cnt[j];
		int tmp=cnt[i]%(mod-1);
		if(tmp<0)tmp+=mod-1;
		if(tmp)ans=mul(ans,power(f[i],tmp));
	}
	cout<<ans<<"\n";
	return 0;
}