[51nod1668]非010串

题目描述

如果一个01字符串满足不存在010这样的子串，那么称它为非010串。
求长度为n的非010串的个数。(对1e9+7取模)

DP

写个DP式，然后矩阵乘法即可
因为比较懒现在懒得去推那个式子是啥了
看代码吧

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int mo=1000000007;
typedef long long ll;
struct dong{
    int a[3][3];
} ans,tmp,kqp;
int s[70];
int i,j,k,l,t,m,top;
ll n;
void mult(dong a,dong b,dong &c){
    dong d;
    int i,j,k;
    fo(i,0,2)
        fo(j,0,2)
            d.a[i][j]=0;
    fo(k,0,2)
        fo(i,0,2)
            fo(j,0,2)
                (d.a[i][j]+=(ll)a.a[i][k]*b.a[k][j]%mo)%=mo;
    fo(i,0,2)
        fo(j,0,2)
            c.a[i][j]=d.a[i][j];
}
int main(){
    scanf("%lld",&n);
    if (n==1) printf("2\n");
    else if (n==2) printf("4\n");
    else{
        kqp.a[0][0]=4;kqp.a[0][1]=1;kqp.a[0][2]=2;
        tmp.a[0][0]=2;
        tmp.a[0][2]=tmp.a[2][1]=1;
        tmp.a[1][0]=tmp.a[1][2]=-1;
        n-=2;
        while (n){
            s[++top]=n%2;
            n/=2;
        }
        ans.a[0][0]=ans.a[1][1]=ans.a[2][2]=1;
        while (top){
            mult(ans,ans,ans);
            if (s[top]) mult(ans,tmp,ans);
            top--;
        }
        mult(kqp,ans,ans);
        printf("%d\n",(ans.a[0][0]+=mo)%=mo);
    }
}