比赛 |
noip普及组模拟赛 |
评测结果 |
AAAAAAAAAA |
题目名称 |
AABB(无题面) |
最终得分 |
100 |
用户昵称 |
胡嘉兴 |
运行时间 |
0.003 s |
代码语言 |
C++ |
内存使用 |
0.31 MiB |
提交时间 |
2018-10-27 17:59:43 |
显示代码纯文本
#include<bits/stdc++.h>
#define ll long long
#define lowbit(x) (x&-x)
using namespace std;
const int N = 1e6+7, p = 1e9+7;
int mul(int a, int b)
{
return (ll)a*b%p;
}
void add(int & a, int b)
{
a += b;
if(a>=p)
{
a -= p;
}
return;
}
struct matrix
{
int a[3][3];
matrix()
{
memset(a, 0, sizeof(a));
a[2][1] = 2;
a[1][1] = a[1][2] = a[2][2] = 1;
return;
}
void init()
{
memset(a, 0, sizeof(a));
return;
}
void Init()
{
memset(a, 0, sizeof(a));
a[1][1] = a[2][2] = 1;
}
matrix operator * (const matrix & an)const
{
matrix ret;
ret.init();
for(int i = 1; i <= 2; i++)
{
for(int j = 1; j <= 2; j++)
{
for(int k = 1; k <= 2; k++)
{
add(ret.a[i][j], mul(a[i][k], an.a[k][j]));
}
}
}
return ret;
}
};
matrix qpow(int k)
{
matrix ret, x;
ret.Init();
while(k)
{
if(k&1)
{
ret = ret*x;
}
x = x*x;
k >>= 1;
}
return ret;
}
int main()
{
int n;
freopen("AABB.in", "r", stdin);
freopen("AABB.out", "w", stdout);
scanf("%d", &n);
matrix res = qpow(n-1);
printf("%lld\n", (res.a[1][1]+res.a[2][1]+(res.a[1][2]+res.a[2][2])*2ll)%p);
return 0;
}