《GMOJ-Senior-1240 Fibonacci sequence》题解

题目大意

求${\sum }_{i=x}^{y}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000$（$f(i)=\{{}_{f(i-1)+f(i-2),(i\ge 3)}^{1,(i=1,2)}$）。

分析

这是一道数论题。
显然，$f(i)$为斐波那契数列。因为
$$\sum _{i=x}^{y}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000=[(\sum _{i=1}^{y}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000)-(\sum _{i=1}^{x-1}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000)+10000]\mathrm{m}\mathrm{o}\mathrm{d}10000\text{，}$$
所以问题就被转化为求${\sum }_{i=1}^{n}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000$。如此一来，问题就有很多种解法了。

找规律

因为只是$\mathrm{m}\mathrm{o}\mathrm{d}10000$，所以${\sum }_{i=1}^{n}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000$可能存在循环。我们可以先写一个程序，找一下${\sum }_{i=1}^{n}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000$的规律：

#include<cstdio>
int mod=10000;//模数
int main()
{
	printf("1");
	int a=1/*递推斐波那契数列用*/,b=1/*递推斐波那契数列用*/,s=2/*递推sum用*/;
	for(int i=3;i<=300000;i++)
	{
		int c=(a+b)%mod/*求出模意义下下一个斐波那契数列的数*/,ns=(s+c)%mod/*递推原式*/;
		if(s==1&&ns==2)//找到循环
		{
			printf(" %d",i-1);//输出
		}
		a=b;//更新
		b=c;
		s=ns;
	}
	return 0;
}

运行程序得（竖着看）：

$1$	$75001$	$150001$	$225001$
$15001$	$90001$	$165001$	$240001$
$30001$	$105001$	$180001$	$255001$
$45001$	$120001$	$195001$	$270001$
$60001$	$135001$	$210001$	$285001$

我们可以发现，这些数都可以表示为$15000n+1$（$n$为非负整数）。由此，我们可以知道，${\sum }_{i=1}^{n}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000={\sum }_{i=1}^{n\mathrm{m}\mathrm{o}\mathrm{d}15000}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000$。于是，我们只要求出${\sum }_{i=1}^{k}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000$（$0\le k<15000$），输出
$$[(\sum _{i=1}^{y\mathrm{m}\mathrm{o}\mathrm{d}15000}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000)-(\sum _{i=1}^{(x-1)\mathrm{m}\mathrm{o}\mathrm{d}15000}f(i)\mathrm{m}\mathrm{o}\mathrm{d}10000)+10000]\mathrm{m}\mathrm{o}\mathrm{d}10000$$
的值就可以了。
代码如下：

#include<cstdio>
int mod=10000/*模数*/,f[15051]/*sum*/;
int main()
{
	int T;
	scanf("%d",&T);//读入数据组数
	f[0]=0;//初始化f
	f[1]=1;
	f[2]=2;
	int a=1,b=1;//递推斐波那契数列
	for(int i=3;i<15000;i++)
	{
		int c=(a+b)%mod;//求出模意义下下一个斐波那契数列的数
		f[i]=(f[i-1]+c)%mod;//递推f
		a=b;//更新
		b=c;
	}
	while(T--)//注意多组数据
	{
		int x,y;
		scanf("%d%d",&x,&y);//读入x,y
		printf("%d\n",(f[y%15000]-f[(x-1)%15000]+mod)%mod);//输出
	}
	return 0;
}

矩阵快速幂

我们想到可以用矩阵快速幂求出斐波那契数列的任意一项，那么我们能否用矩阵快速幂求出斐波那契数列前$n$项的和呢？（注：这里不讨论取模的情况）

方法1

想要求出斐波那契数列前$n$项的和，我们先要求出斐波那契数列的第$n$项，于是我们要记录$f(n-1),f(n)$。又因为要求和，所以我们还要记录一项$su{m}_n$表示${\sum }_{i=1}^{n}f(i)$。由此我们得到一个$3\times 1$的矩阵$\left[\begin{array}{c}f(n-1)\\ f(n)\\ su{m}_n\end{array}\right]$。那我们要怎么从$\left[\begin{array}{c}f(n-1)\\ f(n)\\ su{m}_n\end{array}\right]$推出$\left[\begin{array}{c}f(n)\\ f(n+1)\\ su{m}_{n+1}\end{array}\right]$呢？显然要被一个$3\times 3$的矩阵乘。那么怎么推出这个矩阵呢？
因为
$$\left[\begin{array}{ccc}{a}_{1,1}& a_{1,2}& a_{1,3}\\ a_{2,1}& a_{2,2}& a_{2,3}\\ a_{3,1}& a_{3,2}& a_{3,3}\end{array}\right]\times \left[\begin{array}{c}{b}_1\\ b_2\\ b_3\end{array}\right]=\left[\begin{array}{c}{a}_{1,1}\times b_1+a_{1,2}\times b_2+a_{1,3}\times b_3\\ a_{2,1}\times b_1+a_{2,2}\times b_2+a_{2,3}\times b_3\\ a_{3,1}\times b_1+a_{3,2}\times b_2+a_{3,3}\times b_3\end{array}\right]\text{，}$$
所以我们可以通过$\left[\begin{array}{c}f(n-1)\\ f(n)\\ su{m}_n\end{array}\right]$和$\left[\begin{array}{c}f(n)\\ f(n+1)\\ su{m}_{n+1}\end{array}\right]$的关系求出这个$3\times 3$的矩阵：

1. 由$f(n-1)$到$f(n)$：因为$f(n)$在$\left[\begin{array}{c}f(n-1)\\ f(n)\\ su{m}_n\end{array}\right]$中的第二行已经出现了，我们可以直接用，所以$a_{1,1}=0,a_{1,2}=1,a_{1,3}=0$；
2. 由$f(n)$到$f(n+1)$：因为$f(n+1)=f(n-1)+f(n)$，而$f(n-1)$出现在$\left[\begin{array}{c}f(n-1)\\ f(n)\\ su{m}_n\end{array}\right]$的第一行，$f(n)$出现在$\left[\begin{array}{c}f(n-1)\\ f(n)\\ su{m}_n\end{array}\right]$的第二行，所以$a_{2,1}=1,a_{2,2}=1,a_{2,3}=0$；
3. 由$su{m}_n$到$su{m}_{n+1}$：因为$su{m}_{n+1}=su{m}_n+f(n+1)$，而$f(n+1)$的表示已经在上面得出，所以$a_{3,1}=1,a_{3,2}=1,a_{3,3}=1$。

于是，这个$3\times 3$的矩阵就被求出来了：$\left[\begin{array}{ccc}0& 1& 0\\ 1& 1& 0\\ 1& 1& 1\end{array}\right]$。有了这个矩阵，我们就可以用矩阵快速幂求出结果了。我们令$f(0)=0$，那么可得
$$\sum _{i=0}^{n}f(i)={\left[\begin{array}{ccc}0& 1& 0\\ 1& 1& 0\\ 1& 1& 1\end{array}\right]}^{n-1}\times \left[\begin{array}{c}f(0)\\ f(1)\\ su{m}_1\end{array}\right]\text{所得最后一项}={\left[\begin{array}{ccc}0& 1& 0\\ 1& 1& 0\\ 1& 1& 1\end{array}\right]}^{n-1}\times \left[\begin{array}{c}0\\ 1\\ 1\end{array}\right]\text{所得最后一项。}$$
因为$f(0)=0$，所以
$$\sum _{i=1}^{n}f(i)=\sum _{i=0}^{n}f(i)={\left[\begin{array}{ccc}0& 1& 0\\ 1& 1& 0\\ 1& 1& 1\end{array}\right]}^{n-1}\times \left[\begin{array}{c}0\\ 1\\ 1\end{array}\right]\text{所得最后一项。}$$
同时我们只要特判一下$n=0$的情况（结果为$0$）就好了。
如何进行矩阵快速幂呢？和数的快速幂类似，我们只要把数的乘法换成矩阵的乘法即可。如何进行矩阵的乘法呢？对于一个$n\times m$的矩阵$A$和一个$m\times k$的矩阵$B$，$A\times B$所得的矩阵$C$是$n\times k$的，且$C_{i,j}={\sum }_{t=1}^m(A_{i,t}\times B_{t,j})$。（此处$1\le i\le n,1\le j\le k$）
解题的代码如下：

#include<cstdio>
int mod=10000;//模数（注意取模）
struct matrix{int data[4][4];};//矩阵
matrix mul(matrix a,matrix b)//矩阵乘法
{
	matrix ans;
	for(int i=1;i<=3;i++)
	{
		for(int j=1;j<=3;j++)
		{
			ans.data[i][j]=0;
			for(int k=1;k<=3;k++)
			{
				ans.data[i][j]+=a.data[i][k]*b.data[k][j]%mod;
				ans.data[i][j]%=mod;
			}
		}
	}
	return ans;
}
int fsum(int n)
{
	if(n==0)//特判n=0
	{
		return 0;
	}
	matrix a/*推出的矩阵*/,ans/*最终矩阵*/;
	a.data[1][1]=0;//初始化a（推出的3*3矩阵）
	a.data[1][2]=1;
	a.data[1][3]=0;
	a.data[2][1]=1;
	a.data[2][2]=1;
	a.data[2][3]=0;
	a.data[3][1]=1;
	a.data[3][2]=1;
	a.data[3][3]=1;
	ans.data[1][1]=1;//初始化ans（单位矩阵）
	ans.data[1][2]=0;
	ans.data[1][3]=0;
	ans.data[2][1]=0;
	ans.data[2][2]=1;
	ans.data[2][3]=0;
	ans.data[3][1]=0;
	ans.data[3][2]=0;
	ans.data[3][3]=1;
	--n;//注意指数是n-1
	while(n>=1)//矩阵快速幂
	{
		if(n%2==1)
		{
			ans=mul(ans,a);
		}
		a=mul(a,a);
		n/=2;
	}
	return ans.data[3][2]+ans.data[3][3];//0*ans.data[3][1]+1*ans.data[3][2]+1*ans.data[3][3]=ans.data[3][2]+ans.data[3][3]
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)//注意多组数据
	{
		int x,y;
		scanf("%d%d",&x,&y);//读入x,y
		printf("%d\n",(fsum(y)-fsum(x-1)+mod)%mod);//输出
	}
	return 0;
}

方法2

方法$1$中的矩阵是$3\times 3$和$3\times 1$的，运算起来常数较大，能不能化简呢？
首先，我们需要知道斐波那契数列的一个性质：${\sum }_{i=1}^{n}f(i)=f(n+2)-f(2)$。证明如下：

$\because f(3)=f(1)+f(2),f(4)=f(2)+f(3),f(5)=f(3)+f(4),\cdots ,f(n+2)=f(n)+f(n+1)$

$\therefore f(1)=f(3)-f(2),f(2)=f(4)-f(3),f(3)=f(5)-f(4),\cdots ,f(n)=f(n+2)-f(n+1)$

$\because {\sum }_{i=1}^{n}f(i)=f(1)+f(2)+f(3)+\cdots +f(n)$

$\therefore {\sum }_{i=1}^{n}f(i)=f(3)-f(2)+f(4)-f(3)+f(5)-f(4)+\cdots +f(n+2)-f(n+1)$

得${\sum }_{i=1}^{n}f(i)=f(n+2)-f(2)$

有了这个性质后，我们求${\sum }_{i=1}^{n}f(i)$就变成了求$f(n+2)-f(2)=f(n+2)-1$,。于是我们就可以用可以用矩阵快速幂快速求出斐波那契数列的第$n+2$项，迅速得出答案了。
与方法$1$中的推到方法类似，我们可以确定出我们要记录的数值为$f(n-1),f(n)$，且$\left[\begin{array}{cc}0& 1\\ 1& 1\end{array}\right]\times \left[\begin{array}{c}f(n-1)\\ f(n)\end{array}\right]=\left[\begin{array}{c}f(n)\\ f(n+1)\end{array}\right]$，令$f(0)=0$，那么
$$\sum _{i=1}^{n}f(i)={\left[\begin{array}{cc}0& 1\\ 1& 1\end{array}\right]}^{n+2-1}\times \left[\begin{array}{c}f(0)\\ f(1)\end{array}\right]\text{所得最后一项}-1={\left[\begin{array}{cc}0& 1\\ 1& 1\end{array}\right]}^{n+1}\times \left[\begin{array}{c}0\\ 1\end{array}\right]\text{所得最后一项}-1\text{。}$$
解题代码如下：

#include<cstdio>
int mod=10000;//模数（记得取模）
struct matrix{int data[3][3];};//矩阵
matrix mul(matrix a,matrix b)//矩阵乘法
{
	matrix ans;
	for(int i=1;i<=2;i++)
	{
		for(int j=1;j<=2;j++)
		{
			ans.data[i][j]=0;
			for(int k=1;k<=2;k++)
			{
				ans.data[i][j]+=a.data[i][k]*b.data[k][j]%mod;
				ans.data[i][j]%=mod;
			}
		}
	}
	return ans;
}
int f(int n)
{
	matrix a/*推出的矩阵*/,ans/*最终矩阵*/;
	a.data[1][1]=0;//初始化a（推出的矩阵）
	a.data[1][2]=1;
	a.data[2][1]=1;
	a.data[2][2]=1;
	ans.data[1][1]=1;//初始化ans（单位矩阵）
	ans.data[1][2]=0;
	ans.data[2][1]=0;
	ans.data[2][2]=1;
	--n;//注意指数
	while(n>=1)//矩阵快速幂
	{
		if(n%2==1)
		{
			ans=mul(ans,a);
		}
		a=mul(a,a);
		n/=2;
	}
	return ans.data[2][2];//0*ans.data[2][1]+1*ans.data[2][2]=ans.data[2][2]
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)//注意多组数据
	{
		int x,y;
		scanf("%d%d",&x,&y);//读入x,y
		printf("%d\n",(f(y+2)-f(x+1/*x-1+2=x+1*/)/*(f(y+2)-1)-(f(x+1)-1)=f(y+2)-1-f(x+1)+1=f(y+2)-f(x+1)*/+mod)%mod);//输出
	}
	return 0;
}

总结

有时候一道题不止一种做法，我们要从多方面思考问题。