HDU6030 Happy Necklace[矩阵快速幂]

Happy Necklace

  HDU - 6030

 

题意:

小Q要弄一条项链给女朋友，而弄这条项链的要求，就是要每素数长度里，红色的珠子的个数都不能少于蓝色珠子的个数(就是说红色个数小于等于蓝色个数)，然后给你n个珠子(红色蓝色都可以)，问弄成符合这样条件的一条项链，能有多少种方案。

题解:

要保证每个素数长度都符合这个条件，那么只要保证最小的两个素数满足就可以了，就是说，必须保证两个珠子的时候，红色跟蓝色个数相同，或者只有红色珠子，三个珠子的时候，红色的个数要比蓝色的多。

用o表示红色珠子，用x表示蓝色珠子。

那么满足条件的只有 ox xo oo这三种情况。

然后接下来要继续满足这三种情况的话。

① ox  的后面，只能接 o. 会构成 oxo 

因为如果接的是x，会构成 oxx 。本身红色跟蓝色的个数相同，还要加上一个蓝色，就会使得蓝色比红色多，不符合条件。

② xo 的后面，也只能接 o 会构成 xoo

因为如果接的是x，会构成 xox ，跟上面的情况一样，本身红色跟蓝色的个数相同，还要加上一个蓝色，就会使得蓝色比红色多，不符合条件。

③ oo 的后面，既可以接 o ，也可以接 x。

因为接 o 的话 构成的是 ooo ， 接 x 的话，构成 oox ，无论是长度为2 还是为3 都是符合要求的。

最后根据这三种情况，我们可以列出一个矩阵。

初始的时候 F2的情况

最后结果就是

(渣渣的人生第一次一血啊T_T[在女生专场里]  哇 刚看到  来补补博客纪念一下)

#pragma comment(linker, "/STACK:102400000,102400000")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<string>
#include<algorithm>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
const int N=5;
ll ans[N][N],a[N][N],tmp[N][N];
void initial()
{
    memset(ans,0,sizeof(ans));
    memset(a,0,sizeof(a));
    for (int i=0 ; i<3 ; ++i)
        ans[i][i]=1;
    a[0][2]=a[1][0]=a[2][1]=a[2][2]=1;
}
void Pow(ll a[][N],ll b[][N],int n)
{
    memset(tmp,0,sizeof(tmp));
    for (int i=0 ; i<n ; ++i)
        for (int j=0 ; j<n ; ++j)
            if (a[i][j])
                for (int k=0 ; k<n ; ++k)
                    tmp[i][k]=(tmp[i][k]+a[i][j]*b[j][k]%mod)%mod;
    for (int i=0 ; i<n ; ++i)
        for (int j=0 ; j<n ; ++j)
            a[i][j]=tmp[i][j];
}
void Pow(ll a[][N],ll m)
{
    while (m)
    {
        if (m&1)
            Pow(ans,a,3);
        Pow(a,a,3);
        m>>=1;
    }
}
int main()
{
	int T;
	scanf("%d",&T);
	while (T--)
    {
        ll n;
        scanf("%lld",&n);
        initial();
        ll res=0;
        Pow(a,n-2);
        for (int i=0 ; i<3 ; ++i)
            for (int j=0 ; j<3 ; ++j)
                res=(res+ans[i][j])%mod;
        printf("%lld\n",res);
    }
	return 0;
}