51nod 1355 斐波那契的最小公倍数

原题链接.

###题解：
这题肯定是我做过的最简单的八级算法题了。

设$Fi$为斐波拉契序列的第i项。

通过打表可得：
$gcd(Fn,Fm)=F_{gcd(n,m)}$

接下来尝试去证明这个结论。

首先斐波拉契序列有一个递推式：
$F_{n+m}=F_{n-1}\ast F_m+F_n\ast F_{m+1}$

还有一个结论就是$gcd(F_{n-1},F_n)=1$。

这两个结论都是显然的，就不证明了。

不妨设$n<m$，那么
$gcd(Fn,Fm)$
$=gcd(F_n,F_{n-1}\ast F_{m-n}+F_n\ast F_{m-n+1})$
因为$F_n\ast F_{m-n+1}$一定是$F_n$的倍数，因此可以直接去掉。
$=gcd(F_n,F_{n-1}\ast F_{m-n})$
因为$gcd(F_n,F_{n-1})=1$，所以可以去掉
$=gcd(F_n,F_{m-n})$
。
发现下标在更相减损，归纳可得$gcd(Fn,Fm)=F_{gcd(n,m)}$

$lcm$有一个容斥求法：

$lcm(a1,a2,\dots ,an)$
$={\prod }_{b\in a且b̸=\varnothing }gcd(b1,b2,\dots ,b_{|b|}{)}^{-1^{|b|+1}}$

证明可通过对每一个质因子单独讨论得到。

设$q_i$为$a_i$含质因子p的幂次。

原式等价于：
$={\prod }_{b\in q且b̸=\varnothing }min(b_i{)}^{-1^{|b|+1}}$
假设已经把q从小到大排序了
$={\sum }_{k=1}^n{q}_k\ast {\sum }_{i=0}^{n-k}{C}_{n-k}^i\ast (-1{)}^i$
$={\sum }_{k=1}^n{q}_k\ast [n-k=0]$
$=q_n$
$=排序前的max(q)$

把每个质因子合并起来即得证。

设$f_i$表示$gcd$恰好为$i$的每种方案×这种方案的容斥系数（就是那个-1的多少次方）的和。

直接算不好算，于是反演一下，$g_i$定义类似，为$gcd$是$i$的倍数的XXX

设读入的n个数中有k个数是i的倍数，则$g(k)=[k>0]$

由$g_i={\sum }_{i|j}{f}_j$得
$f_i={\sum }_{i|j}{g}_j\ast \mu (j/i)$

至此此题完成，写起来真的好写，就是证明有点多。

Code：

#include<cstdio> 
#define ll long long
#define fo(i, x, y) for(int i = x; i <= y; i ++)
using namespace std;

const int mo = 1e9 + 7;

const int N = 1e6;

int n, x;
bool bz[N + 1];
int p[N], mu[N + 1], s[N + 1], f[N + 1];

ll ksm(ll x, ll y) {
	if(y < 0) x = ksm(x, mo - 2), y = -y;
	ll s = 1;
	for(; y; y /= 2, x = x * x % mo)
		if(y & 1) s = s * x % mo;
	return s;
}

int main() {
	f[0] = 0; f[1] = 1; mu[1] = 1;
	fo(i, 2, N) {
		f[i] = (f[i - 1] + f[i - 2]) % mo;
		if(!bz[i]) p[++ p[0]] = i, mu[i] = -1;
		for(int j = 1; i * p[j] <= N; j ++) {
			int k = i * p[j]; bz[k] = 1;
			if(i % p[j] == 0) {mu[k] = 0; break;}
			mu[k] = -mu[i];
		}
	}
	scanf("%d", &n);
	fo(i, 1, n) scanf("%d", &x), s[x] ++;
	fo(i, 1, N) {
		fo(j, 2, N / i) s[i] += s[i * j];
		s[i] = s[i] > 0;
	}
	fo(i, 1, N) fo(j, 2, N / i) s[i] += s[i * j] * mu[j];
	ll ans = 1; fo(i, 1, N) ans = ans * ksm(f[i], s[i]) % mo;
	printf("%lld", ans);
}