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

$\text{lcm}$不好处理，考虑转为$\gcd$，主要是因为$(f_n,f_m)=f_{(n,m)}$

$\text{lcm}$的本质是指数取$\max$，$\gcd$的本质是指数取$\min$，所以我们可以min-max容斥

设给定的集合为$S$

$\begin{align*}\mathop{\text{lcm}}\limits_{x\in S}f_x&=\prod\limits_{T\subset S}\left(\gcd\limits_{x\in T}f_x\right)^{(-1)^{|T|-1}}\\&=\prod\limits_{T\subset S}f_{\gcd(T)}^{(-1)^{|T|-1}}\end{align*}$

构造$g$使$f_n=\prod\limits_{d|n}g(d)$

$\begin{align*}\prod\limits_{T\subset S}f_{\gcd(T)}^{(-1)^{|T|-1}}&=\prod\limits_{T\subset S}\left(\prod\limits_{d|\gcd(T)}g(d)\right)^{(-1)^{|T|-1}}\\&=\sum\limits_dg(d)^{\sum\limits_{T\subset S,d|\gcd(T)}(-1)^{|T|-1}}\\&=\sum\limits_{\exists kd\in S}g(d)\end{align*}$

这里的$\sum\limits_{T\subset S,d|\gcd(T)}(-1)^{|T|-1}$，如果$S$中没有$d$的倍数，那么它是$0$，如果有，那么我们把$S$中非$d$倍数的数全部删掉形成$S'$，即是$\sum\limits_{T\subset S',T\ne\varnothing}(-1)^{|T|-1}=1$

#include<stdio.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int mod=1000000007,T=1000000;
int mul(int a,int b){return a*(ll)b%mod;}
int ad(int a,int b){return(a+b)%mod;}
int pow(int a,int b){
	int s=1;
	while(b){
		if(b&1)s=mul(s,a);
		a=mul(a,a);
		b>>=1;
	}
	return s;
}
int pr[T+10],mu[T+10];
bool np[T+10];
void sieve(){
	int i,j,M=0;
	mu[1]=1;
	for(i=2;i<=T;i++){
		if(!np[i]){
			pr[++M]=i;
			mu[i]=-1;
		}
		for(j=1;j<=M&&i*pr[j]<=T;j++){
			np[i*pr[j]]=1;
			if(i%pr[j]==0)break;
			mu[i*pr[j]]=-mu[i];
		}
	}
}
int f[T+10],g[T+10];
bool ex[T+10];
int main(){
	sieve();
	int n,m,i,j,s;
	scanf("%d",&n);
	m=0;
	for(i=1;i<=n;i++){
		scanf("%d",&j);
		m=max(m,j);
		ex[j]=1;
	}
	f[1]=1;
	for(i=2;i<=m;i++)f[i]=ad(f[i-1],f[i-2]);
	for(i=1;i<=m;i++)g[i]=1;
	for(i=1;i<=m;i++){
		s=pow(f[i],mod-2);
		for(j=i;j<=m;j+=i){
			if(mu[j/i]==1)
				g[j]=mul(g[j],f[i]);
			else if(mu[j/i]==-1)
				g[j]=mul(g[j],s);
		}
	}
	s=1;
	for(i=1;i<=m;i++){
		for(j=i;j<=m;j+=i){
			if(ex[j]){
				j=-1;
				break;
			}
		}
		if(j==-1)s=mul(s,g[i]);
	}
	printf("%d",ad(s,mod));
}

转载于:https://www.cnblogs.com/jefflyy/p/9408198.html