bzoj #2813.奇妙的Fibonacci 题解

题目传送门

题目大意： 令 $f_i$ 为斐波那契数列第 $i$ 项，问有多少个 $j$ 满足 $f_j|f_i$，并求出 $j$ 的平方和。

题解

$f_j|f_i$ 等价于 $\gcd (f_i,f_j)=f_j$，而斐波那契有个神奇的结论：$\gcd (f_i,f_j)=f_{\gcd (i,j)}$，代入进去就是 $f_{\gcd (i,j)}=f_j$，即 $\gcd (i,j)=j$，也就是说，满足要求的 $j$ 就是 $i$ 的因子。

剩下就是线性筛了，约数个数 和 约数平方和 筛法本质是一样的，因为这两个东西就是本质一样的两个函数，即 ${\sigma }_0$ 和 ${\sigma }_2$，想必大家都很熟练了。

下面是那个神奇结论的证明：

引理1： $f_{n+m}=f_{m+1}{f}_n+f_m{f}_{n-1}$。

证明：
$$\begin{array}{rl}{f}_{n+m}& =f_{n+m-1}+f_{n+m-2}\\ & =2f_{n+m-2}+f_{n+m-3}\\ & =f_{k+1}{f}_{n+m-k}+f_k{f}_{n+m-k-1}\end{array}$$

令 $k=m$，代入得到：
$$=f_{m+1}{f}_n+f_m{f}_{n-1}$$

---

引理2： $\gcd (f_n,f_{n+1})=1$。

证明： 这个说实话应该挺众所周知的。

考虑归纳法，一开始有 $\gcd (f_1,f_2)=1$。

假设 $\gcd (f_{n-1},f_n)=1$，那么有 $\gcd (f_{n+1},f_n)=\gcd (f_{n+1}-f_n,f_n)=\gcd (f_{n-1},f_n)=1$。

---

引理3： $\gcd (f_{n+m},f_m)=\gcd (f_n,f_m)$。

证明： 根据引理1，有 $\gcd (f_{n+m},f_m)=\gcd (f_{m+1}{f}_n+f_m{f}_{n-1},f_m)=\gcd (f_{m+1}{f}_n,f_m)$。

根据引理1，有 $\gcd (f_{m+1},f_m)=1$，所以 $\gcd (f_{m+1}{f}_n,f_m)=\gcd (f_n,f_m)$。

---

注意到有了这个引理3，其实我们已经可以对下标做类似辗转相减的事情了，那么一开始的结论 $\gcd (f_n,f_m)=f_{\gcd (n,m)}$ 也就证明出来了。

证毕。

代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 10000010
#define mod 1000000007

int Q,Qx,A,B,C,ans1=0,ans2=0;
void reduce(int &x){x+=x>>31&mod;}
int prime[maxn],t=0;
bool v[maxn];
int d1[maxn],d2[maxn],pd[maxn];
void SieveInit(){
	d1[1]=d2[1]=1;
	for(int i=2;i<=maxn-10;i++){
		if(!v[i])prime[++t]=i,d1[i]=2,reduce(d2[i]=1+1ll*i*i%mod-mod),pd[i]=1;
		for(int j=1;j<=t&&i*prime[j]<=maxn-10;j++){
			v[i*prime[j]]=true;
			if(i%prime[j]==0){
				pd[i*prime[j]]=pd[i];
				reduce(d1[i*prime[j]]=d1[i]+d1[pd[i]]-mod);
				reduce(d2[i*prime[j]]=1ll*d2[i]*prime[j]%mod*prime[j]%mod+d2[pd[i]]-mod);
				break;
			}
			pd[i*prime[j]]=i;
			reduce(d1[i*prime[j]]=d1[i]+d1[i]-mod);
			reduce(d2[i*prime[j]]=1ll*d2[i]*prime[j]%mod*prime[j]%mod+d2[i]-mod);
		}
	}
}

int main()
{
	SieveInit();
	scanf("%d %d %d %d %d",&Q,&Qx,&A,&B,&C);
	while(Q--){
		reduce(ans1+=d1[Qx]-mod);
		reduce(ans2+=d2[Qx]-mod);
		if(Qx&1)reduce(ans1+=1-mod),reduce(ans2+=4-mod);
		Qx=(1ll*Qx*A+B)%C+1;
	}
	printf("%d\n%d",ans1,ans2);
}