[BZOJ4833] [Lydsy1704月赛]最小公倍佩尔数（Min-Max 容斥）

题意

• 令 $(1+\sqrt{2}{)}^n=e(n)+\sqrt{2}f(n),g(n)=\text{lcm}(f(1),f(2),...,f(n))$，求 $\sum _{i=1}^{n}g(i)\times i$。（$n\le 10^6$）

挺巧妙的数学题，全程对着题解抄的，要推出来脑洞得足够大。

首先可以知道： $f(n)=2f(n-1)+f(n-2)$

引理一：对于 $f(n)=af(n-1)+bf(n-2)$ 形式的递推式，如果有 $\gcd (a,b)=1$，那么 $\gcd (f(a),f(b))=f(\gcd (a,b))$，具体可以联系辗转相除的过程用归纳法证明。

引理二：
$$\text{lcm}(S)=\prod _{T\subset S}(-1{)}^{|T|+1}\gcd (T)(T̸=\varphi )$$
具体证明把 $\text{lcm}$ 看做指数取最大值，$\mathrm{g}\mathrm{c}\mathrm{d}$ 看做指数取最小值，套用 Min-Max 容斥即可。

定义 U = { 1 , 2 , . . . , n } U=\{1,2,...,n\} U={1,2,...,n}，$h$ 满足$f(n)=\prod _{d|n}h(d)$，那么 $h(n)=\frac{f(n)}{\prod _{d|n,d̸=n}h(d)}$，我们可以 $O(n\ln n)$ 枚举倍数计算。

由引理一和引理二可知：
$$\begin{gathered} g(n)=\prod _{T\subset U}(-1{)}^{|T|+1}f(\underset{i\in T}{\gcd }i) \\ g(n)=\prod _{T\subset U}(-1{)}^{|T|+1}\prod _{d|\underset{i\in T}{\gcd }i}{h}_d \end{gathered}$$
设 S d = { d , 2 d , . . . , ⌊ n d ⌋ } S_d=\{d,2d,...,\lfloor\frac n d \rfloor\} Sd​={d,2d,...,⌊dn​⌋}，那么：
$$g(n)=\prod _{d=1}^n{h}_d^{\sum _{T\subset S_d}(-1{)}^{|T|+1}}$$

又因为$\sum _{T\subset S_d}(-1{)}^{|T|+1}=\sum _{i=1}^{|S_d|}(-1{)}^{i+1}\left(\genfrac{}{}{0ex}{}{|S_d|}{i}\right)$，由二项式定理可知后面这个东西等于一，所以 $g(n)=\prod _{d=1}^n{h}_d$。

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 5;

int T, n, Mod, h[N], f[N];

inline int qpow(int a, int x)
{
	int ret = 1; 
	for (; x; x >>= 1, a = 1ll * a * a % Mod)
		if (x & 1) ret = 1ll * ret * a % Mod;
	return ret;
}

int main()
{
#ifdef ylsakioi
	freopen("4833.in", "r", stdin);
	freopen("4833.out", "w", stdout);
#endif
	f[0] = 0, f[1] = h[1] = 1;
	for (scanf("%d", &T); T -- ; )
	{
		scanf("%d%d", &n, &Mod);
		for (int i = 2; i <= n; ++ i)
			h[i] = f[i] = (2ll * f[i - 1] + f[i - 2]) % Mod;
		for (int i = 2; i <= n; ++ i)
		{
			int d = qpow(h[i], Mod - 2);
			for (int j = i << 1; j <= n; j += i) 
				h[j] = 1ll * h[j] * d % Mod;
		}
		int ans = 0, ret = 1; 
		for (int i = 1; i <= n; ++ i) 
		{
			ret = 1ll * ret * h[i] % Mod;
			(ans += 1ll * ret * i % Mod) %= Mod;
		}
		printf("%d\n", ans);
	}
	return 0;
}