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;}