CodeForces 593E

按时间分段，
在每个时间段内，进行矩阵快速幂

#include
#include
#include
using namespace std;
const int mod=1000000007;
int M[44][44],Me[44][44],V[44],Ve[44],n;
void pow_M(int e){
    while(1)
    {
        if(e&1)
        {
            memset(Ve,0,sizeof(Ve));
            for(int i=0;ifor(int j=0;jlong long)V[j])%mod;
            for(int i=0;i>=1;
        if(e==0)break;

        memset(Me,0,sizeof(Me));
        for(int i=0;ifor(int j=0;jfor(int k=0;klong long)M[k][j])%mod;
        for(int i=0;ifor(int j=0;jint C[44][44],x,y,q;
void build()
{
    for(int i=0;ifor(int j=0;jint a=i/y,b=i%y;
        int A=j/y,B=j%y;
        int k=abs(A-a)+abs(B-b);
        if(k<=1&&C[a][b]==0&&C[A][B]==0)M[i][j]=1;
        else M[i][j]=0;
    }
}

int main()
{
    #ifdef DouBi
    freopen("input.cpp","r",stdin);
    #endif // DouBi
    while(scanf("%d%d%d",&x,&y,&q)!=EOF){
        n=x*y;
        int t=1;
        memset(V,0,sizeof(V));V[0]=1;
        memset(C,0,sizeof(C));
        while(q--)
        {
            build();
            int l,a,b,k;scanf("%d%d%d%d",&l,&a,&b,&k);--a;--b;
            pow_M(k-t);
            t=k;
            if(l==3){V[a*y+b]=0;C[a][b]=0;}
            if(l==2){V[a*y+b]=0;C[a][b]=1;}
            if(l==1)printf("%d\n",V[a*y+b]);
        }
        //printf("line\n");
    }
    return 0;
}