斐波那契数列的和

斐波那契数列的和

题目链接：luogu U145243（自己创的题目）

前言

在你谷上找了一会，发现好像没有这道题目。
然后就决定创一下这个题目。

题目大意

就是求斐波那契数列前 $n$ 项的和 $\mathrm{m}\mathrm{o}\mathrm{d}10^9+7$。

思路

这道题用矩阵乘法来做，然后又两种做法。
（设 $Q_i$ 为斐波那契数列前 $i$ 项的和）

第一种

我们直接构造一个 $1\times 3$ 的矩阵：

1(F[n-2])	1(F[n-1])	1(Q[n-2])

然后可以推出转移矩阵是这个：

0	1	0
1	1	1
0	0	1

然后就乘 $n-1$ 次转移矩阵，最后输出最右边的即可。

第二种

我们可以推出 $Q_i=F_{i+2}-1$。
原因看这里：
$F_{i+2}-1=F_{i+1}+F_i-1$
$=F_i+F_{i-1}+F_i-1$
$=F_{i-1}+F_{i-2}+F_{i-1}+F_i-1$
$=F_2+F_1+F_2+...+F_{i-2}+F_{i-1}+F_i-1$
又因为 $F_2=1$，所以
原来的式子$=1+F_1+F_2+...+F_{i-2}+F_{i-1}+F_i-1$
$=F_1+F_2+...+F_{i-2}+F_{i-1}+F_i=Q_i$

然后就是直接求斐波那契数列的模板题了。

代码

方法1

#include<cstdio>
#define mo 1000000007
#define ll unsigned long long

using namespace std;

struct matrix {
	ll n, m, a[5][5];
}a, b, ans, re;
ll n;

matrix operator *(matrix x, matrix y) {
	re.n = x.n;
	re.m = y.m;
	for (ll i = 0; i < re.n; i++)
		for (ll j = 0; j < re.m; j++)
			re.a[i][j] = 0;
	
	for (ll k = 0; k < x.m; k++)
		for (ll i = 0; i < re.n; i++)
			for (ll j = 0; j < re.m; j++)
				re.a[i][j] = (re.a[i][j] + (x.a[i][k] * y.a[k][j]) % mo) % mo;
	
	return re;
}

void jzksm(ll now) {
	b = a;
	now--;
	while (now) {
		if (now & 1) b = b * a;
		a = a * a;
		now /= 2;
	}
}

int main() {
//	freopen("sjs.txt", "r", stdin);
//	freopen("1.txt", "w", stdout);
	
	scanf("%lld", &n);
	
	if (n == 1) {
		printf("1");
		return 0;
	}
	
	ans.n = 1;
	ans.m = 3;
	ans.a[0][0] = 1;
	ans.a[0][1] = 1;
	ans.a[0][2] = 1;
	
	a.n = 3;
	a.m = 3;
	a.a[1][0] = 1;
	a.a[0][1] = 1;a.a[1][1] = 1;
	a.a[1][2] = 1;a.a[2][2] = 1;
	
	jzksm(n - 1);
	
	ans = ans * b;
	
	printf("%lld", ans.a[0][2]);
	
	fclose(stdin);
	fclose(stdout);
	
	return 0;
}

方法2

#include<cstdio>
#include<cstring>
#define mo 1000000007
#define ll unsigned long long

using namespace std;

ll n, ans[2][2] = {{1, 1}, {0, 0}}, a[2][2] = {{0, 1}, {1, 1}}, b[2][2], re[2][2];

void times(ll x[2][2], ll y[2][2]) {
	memset(re, 0, sizeof(re));
	for (ll k = 0; k < 2; k++)
		for (ll i = 0; i < 2; i++)
			for (ll j = 0; j < 2; j++)
				re[i][j] = (re[i][j] + (x[i][k] * y[k][j]) % mo) % mo;
	for (ll i = 0; i < 2; i++)
		for (ll j = 0; j < 2; j++)
			b[i][j] = re[i][j];
	return ;
}

void jzksm(ll now) {
	if (!now) return ;
	if (now == 1) {
		for (ll i = 0; i < 2; i++)
			for (ll j = 0; j < 2; j++)
				b[i][j] = a[i][j];
		return ;
	}
	
	jzksm(now / 2);
	times(b, b);
	if (now & 1) times(b, a);
	
	return ;
}

int main() {
//	freopen("sjs.txt", "r", stdin);
//	freopen("2.txt", "w", stdout);
	
	scanf("%lld", &n);
	
	if (n == 1) {
		printf("1");
		return 0;
	}
	
	jzksm((n + 2) - 1);
	
	times(ans, b);
	
	printf("%lld", (b[0][0] - 1 + mo) % mo);
	
	fclose(stdin);
	fclose(stdout);
	
	return 0;
}