Problem D: qwb与神奇的序列

Description

qwb又遇到了一道题目：

有一个序列，初始时只有两个数x和y，之后每次操作时，在原序列的任意两个相邻数之间插入这两个数的和，得到新序列。举例说明：

初始：1 2
操作1次：1 3 2
操作2次：1 4 3 5 2
……
请问在操作n次之后，得到的序列的所有数之和是多少？

Input

多组测试数据，处理到文件结束（测试例数量<=50000）。

输入为一行三个整数x，y，n，相邻两个数之间用单个空格隔开。（0 <= x <= 1e10, 0 <= y <= 1e10, 1 < n <= 1e10）。

Output

对于每个测试例，输出一个整数，占一行，即最终序列中所有数之和。
如果和超过1e8，则输出低8位。（前导0不输出，直接理解成%1e8）

Sample Input

1 2 2

Sample Output

15

#include<bits/stdc++.h>
using namespace std;
const int inf=(1<<30)-1;
const int maxn=100010;
#define REP(i,n) for(int i=(0);i<(n);i++)
#define FOR(i,j,n) for(int i=(j);i<=(n);i++)
typedef long long ll;
ll x,y,n;
const int Mo=1e8;
struct matrix {
	ll a[2][2];
	matrix operator *(const matrix &b) const {
		matrix c;memset(c.a,0,sizeof(c.a));
		REP(i,2) REP(j,2) REP(k,2) 
		c.a[i][j]=(c.a[i][j]%Mo+(a[i][k]%Mo)*(b.a[k][j]%Mo)%Mo+Mo)%Mo;
		return c;
	}
};
matrix A,B;
void pow(ll x)
{
	while(x)
	{
		if(x&1) A=A*B;
		x>>=1;
		B=B*B;
	}
}
int main()
{
	//A.a[0][0]=A.a[1][1]=1;
	//B.a[0][0]=3;B.a[0][1]=-1;B.a[1][0]=0;B.a[1][1]=1;
	while(~scanf("%lld%lld%lld",&x,&y,&n))
	{
		memset(A.a,0,sizeof(A.a));
		memset(B.a,0,sizeof(B.a));
		A.a[0][0]=A.a[1][1]=1;
		B.a[0][0]=3;B.a[0][1]=-1;B.a[1][0]=0;B.a[1][1]=1;
		ll tt=x+y;
		pow(n);
		ll Ans=(A.a[0][0]+A.a[0][1])%Mo;
		Ans=(Ans+Mo)%Mo;
		printf("%lld\n",(Ans*tt%Mo)%Mo);
	}
	return 0;
}