【題目來源】：https://vjudge.net/problem/HDU-2604
【題意】
f，m分別是female與male的縮寫，假設有一個佇列裡面是這些字母縮寫，長度為L，那麼共有2^L種，如果含有fmf或者fff這種子佇列的佇列被稱為0佇列，其餘的被稱為E佇列，問，長度為k的有多少種E佇列。
【思路】
先推出前幾項，得到：
f[0]=0,f[1]=2,f[2]=4,f[3]=6,f[4]=9,f[5]=15,f[6]=25;
故得到關係式：f[x]=f[x-1]+f[x-3]+f[x-4]。
其中x>=5.
所以（矩陣1）^(k-4)*矩陣2=矩陣3
我們要求矩陣3，首先要構造矩陣1和矩陣2。
那麼怎麼構造呢？
首先，構造矩陣1：
在f[x]=f[x-1]+f[x-3]+f[x-4]式子裡：（一般遞推的矩陣2都是由關係式得來的）
f[x]=f[x-1]++0*f[x-2]+f[x-3]+f[x-4]。
f[x-1]=1*f[x-1]；
f[x-2]=1*f[x-2]；
f[x-3]=1*f[x-3]；
f[x-4]=1*f[x-4]；
得到矩陣：
1 0 1 1
1 0 0 0
0 1 0 0
0 0 1 0
接著就是構造矩陣2：
f[x-1] 0 0 0
f[x-2] 0 0 0
f[x-3] 0 0 0
f[x-4] 0 0 0
……
沒啦。。。質疑可以發評論討論一下0.0 交流0.0
【程式碼】

#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<limits.h>
#include<algorithm>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace
 
 std;
const int mod=9973;
typedef unsigned long long ll;
typedef long long LL;
int k,m;
struct mat
{
    int a[5][5];
    mat()
    {
        mem(a,0);
        for(int i=1; i<=4; i++)
            a[i][i]=1;
    }
};
mat operator*(mat s,mat t)
{
    mat r;
    mem(r.a,0);
    for(int i=1; i<=4; i++)
    {
        for
 
(int j=1; j<=4; j++)
        {
            for(int p=1; p<=4; p++)
            {
                r.a[i][j]=r.a[i][j]+s.a[i][p]*t.a[p][j];
                if(r.a[i][j]>=m)
                    r.a[i][j]%=m;
            }
        }
    }
    return r;
}
void print(mat t)
{
    for(int i=1; i<=4; i++)
    {
        for(int j=1; j<=4; j++)
        {
            printf("%d ",t.a[i][j]);
        }
        printf("\n");
    }
}
mat pow_mat(mat &ans)
{
    mat base,temp;
    mem(base.a,0);
    base.a[1][1]=1;
    base.a[1][3]=1;
    base.a[1][4]=1;
    for(int i=2; i<=4; i++)
        base.a[i][i-1]=1;
    k-=4;
    while(k)
    {
        if(k&1)
            temp=temp*base;
        base=base*base;
        k>>=1;
    }
//    print(temp);
    ans=temp*ans;
}
int main()
{
    while(~scanf("%d%d",&k,&m))
    {
        if(k<5)
        {
            int p=0;
            switch(k)
            {
            case 0:
                p=0;
                break;
            case 1:
                p=2;
                break;
            case 2:
                p=4;
                break;
            case 3:
                p=6;
                break;
            case 4:
                p=9;
                break;
            }
            printf("%d\n",p%m);
        }
        else
        {
            mat ans;
            mem(ans.a,0);
            ans.a[1][1]=9;
            ans.a[2][1]=6;
            ans.a[3][1]=4;
            ans.a[4][1]=2;
            pow_mat(ans);
            printf("%d\n",(ans.a[1][1])%m);
        }
    }

}