【矩阵乘法】【SSL 2514】 幼儿园数学题Ⅰ

【矩阵乘法】【SSL 2513】 幼儿园数学题Ⅰ

题目

某天，幼儿园学生LZH周测数学时吓哭了，一道题都做不出来。这下可麻烦了他马上就会成为垫底的0分啊。他的期望也不高，做出最简单的第一题就够了
题目是这样的，定义F(n)=((根号5+1)/2)^(n-1) ,当然为了凸显题目的简单当然不能是小数分数或无理数,F(x)因此需要向上取整,当然求F(n)是非常难的!因此幼儿园园长头皮决定简单一点,求下F(x)的前n项和就行了。

---

输入

输入 一个正整数n（保证0<=n<=2^31-1）

---

输出

输出 一个正整数S(n) 对1000000007 取余就好了

---

样例

input1
1

output1
1

input2
2

output 2
2

---

解题思路

额。。。
题目给少了一部分
其实就是比内数列
就是小学生都知道的斐波拉契数列的通项式啦
所以简单的说就是求斐波拉契数列前n项的和

---

代码

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const long long mo=1000000007;
struct lzf{
	long long n,m,h[5][5];
}a,b,x;
long long n;
lzf operator * (lzf l,lzf y)
{
	x.n=l.n,x.m=y.m;
	memset(x.h,0,sizeof(x.h));
	for (int k=1;k<=l.m;k++)
	    for (int i=1;i<=x.n;i++)
	        for (int j=1;j<=x.m;j++)
	            x.h[i][j]=(x.h[i][j]+l.h[i][k]*y.h[k][j]%mo)%mo;
	return x;
}
void power(long long n)
{
	 while (n)
	 {
	 	   if (n & 1) a=a*b; 
		   b=b*b;
		   n>>=1; 
	 }
}
int main()
{
	scanf("%lld",&n);
	a.n=1,a.m=3;
	a.h[1][1]=a.h[1][2]=a.h[1][3]=1;
	b.n=b.m=3;
	b.h[1][2]=b.h[2][1]=b.h[2][2]=1;
	b.h[2][3]=b.h[3][3]=1;
	power(n-1);
	printf("%lld",a.h[1][3]);
	return 0;
}