hdu 3658 How many words

思路：

构造矩阵，矩阵快速幂！！！

代码如下：

#include<cstdio>
#include<vector>
#include<cmath>
#include<iostream>
#include<cstring>
#define ll __int64
#define mod 1000000007
using namespace std;
int let[100];
struct ma
{
    ll a[60][60];
}A,B;
ll sum(ma a)
{
    ll ans=0;
    for(int i=0;i<52;i++)
        for(int j=0;j<52;j++)
            ans=(ans+a.a[i][j])%mod;
    return ans;
}
void init()
{
    int i,j;
    for(i='a';i<='z';i++){
        let[i-'a']=i;
        let[i-'a'+26]=i-'a'+'A';
    }
    for(i=0;i<52;i++)
    for(j=0;j<52;j++){
        if(abs(let[i]-let[j])<=32) A.a[i][j]=1;
        if(abs(let[i]-let[j])<32) B.a[i][j]=1;
    }
}
ma mul(ma a,ma b)
{
    ma ans;
    for(int i=0;i<52;i++)
    for(int j=0;j<52;j++){
        ans.a[i][j]=0;
        for(int k=0;k<52;k++)
            ans.a[i][j]=(ans.a[i][j]+a.a[i][k]*b.a[k][j])%mod;
    }
    return ans;
}
ma pows(ma a,int b)
{
    ma ans;
    for(int i=0;i<52;i++)
        for(int j=0;j<52;j++)
            ans.a[i][j]=(i==j);
    while(b){
        if(b&1) ans=mul(ans,a);
        b>>=1;
        a=mul(a,a);
    }
    return ans;
}
int main()
{
    init();
    int m,t;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&m);
        ll res=sum(pows(A,m-1))-sum(pows(B,m-1));
        if(res<0) res=(res%mod+mod)%mod;
        printf("%I64d
",res);
    }
    return 0;
}