【CF 185A】Plant

解析：

       矩阵快速幂。

       令表示年后的上三角形个数，为年后的三角形个数于是得到递推式：

       

       于是可转化成矩阵乘法用矩阵快速幂解决。

 

代码：

#include <bits/stdc++.h>
using namespace std;

const int mod=1e9+7;
const int Max=3;
long long n;
int ans[Max]={0,1,1},a[Max][Max];

inline void mul(int ans[],int a[Max][Max])
{
	int c[Max];
	memset(c,0,sizeof(c));
	for(int j=1;j<=2;j++)
	  for(int k=1;k<=2;k++)
	    c[j]=(c[j]+1ll*ans[k]*a[k][j]%mod)%mod;
	memcpy(ans,c,sizeof(c));
}

inline void mulself(int a[Max][Max])
{
	int c[Max][Max];
	memset(c,0,sizeof(c));
	for(int i=1;i<=2;i++)
	  for(int j=1;j<=2;j++)
	    for(int k=1;k<=2;k++)
	      c[i][j]=(c[i][j]+1ll*a[i][k]*a[k][j]%mod)%mod;
	memcpy(a,c,sizeof(c));
}

inline void solve()
{
	while(n)
	{
	  if(n&1) mul(ans,a);
	  n>>=1;
	  mulself(a);
	}
}

int main()
{
	scanf("%I64d",&n);
	a[1][1]=2,a[2][1]=1,a[2][2]=4;
	solve();
	printf("%d",ans[1]);
	return 0;
}