2017广西邀请赛 Covering（矩阵快速幂）

题意：给你4*n的矩形，用1*2的方格铺满。

题解：经典问题，我们可以先用公式找出前几项，然后求递推式即可，但是n很大，所以要快速幂搞搞。

我们可以由前几项推得递推公式为：f[i]=f[i-1]+5*f[i-2]+f[i-3]-f[i-4];

因此我们可以构造矩阵为：

1 5 1 -1

1 0 0 0

0 1 0 0

0 0 1 0       然后矩阵快速幂求解就好啦。

附带公式代码：

#include<set>
#include<map>
#include<stack>
#include<queue>
#include<vector>
#include<string>
#include<math.h>
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<functional>
using namespace std;
typedef long long ll;
#define inf 1000000000
#define mod 1000000007
#define maxn 2600005
#define PI acos(-1.0)
#define lowbit(x) (x&-x)
#define eps 1e-9
struct node
{
ll a[5][5];
}a, b;
int main(void)
{
double n, m, i, j, sum;
while (scanf("%lf", &m) != EOF)
{
sum = 1;n = 4;
for (i = 1;i <= (m + 1) / 2;i++)
for (j = 1;j <= (n + 1) / 2;j++)
sum *= 4 * cos(PI*i / (m + 1))*cos(PI*i / (m + 1)) + 4 * cos(PI*j / (n + 1)) * cos(PI*j / (n + 1));
printf("%.0f\n", sum);
}
return 0;
}


AC代码：
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<vector>
#include<string>
#include<math.h>
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<functional>
using namespace std;
typedef long long ll;
#define inf 1000000000
#define mod 1000000007
#define maxn 2600005
#define PI acos(-1.0)
#define lowbit(x) (x&-x)
#define eps 1e-9
struct node
{
ll a[5][5];
}a, b;
ll c[10];
void init()
{
int i, j;
c[1] = 1;c[2] = 5;
c[3] = 11;c[4] = 36;
b.a[1][1] = 1;b.a[1][2] = 5;b.a[1][3] = 1;b.a[1][4] = -1;
for (i = 2;i <= 4;i++)
{
for (j = 1;j <= 4;j++)
{
if (i - j == 1)
b.a[i][j] = 1;
else
b.a[i][j] = 0;
}
}
}
node qq(node a, node b)
{
node c;
int i, j, k;
for (i = 1;i <= 4;i++)
for (j = 1;j <= 4;j++)
{
c.a[i][j] = 0;
for (k = 1;k <= 4;k++)
{
a.a[i][k] %= mod;
b.a[k][j] %= mod;
c.a[i][j] += a.a[i][k] * b.a[k][j];
c.a[i][j] %= mod;
}
}
return c;
}
node q(node a, ll y)
{
node tmp;
int i, j;
for (i = 1;i <= 4;i++)
for (j = 1;j <= 4;j++)
{
if (i == j)
tmp.a[i][j] = 1;
else
tmp.a[i][j] = 0;
}
while (y)
{
if (y % 2)
tmp = qq(tmp, a);
a = qq(a, a);
y /= 2;
}
return tmp;
}
int main(void)
{
init();
ll n, i, j, ans;
while (scanf("%lld", &n) != EOF)
{
ans = 0;
if (n <= 4)
{
printf("%lld\n", c
);
continue;
}
a = q(b, n - 4);
for (i = 1;i <= 4;i++)
{
a.a[1][4 - i + 1] %= mod;
ans += c[i] * a.a[1][4 - i + 1];
ans %= mod;
while (ans < 0)
ans += mod;
}
printf("%lld\n", ans);
}
return 0;
}