矩阵加速

要点
若状态矩阵中的第x个数对下一单位时间状态矩阵中的第y个数产生影响，则把转移矩阵的第x行第y列赋值为适当的系数。
题目描述

小 C 数学成绩优异，于是老师给小 C 留了一道非常难的数学作业题：

给定正整数 NN 和 MM ，要求计算Concatenate (1 … N)Concatenate(1…N) ModMod MM
的值，其中 Concatenate (1 … N)Concatenate(1…N) 是将所有正整数 1, 2, …, N1,2,…,N
顺序连接起来得到的数。例如，N = 13N=13 , Concatenate (1 …
N)=12345678910111213Concatenate(1…N)=12345678910111213 .小C
想了大半天终于意识到这是一道不可能手算出来的题目，于是他只好向你求助，希望你能编写一个程序帮他解决这个问题。

输入格式

从文件input.txt中读入数据，输入文件只有一行且为用空格隔开的两个正整数N和M，其中30%的数据满足1≤N≤10000001≤N≤1000000
；100%的数据满足1≤N≤10^{18}1≤N≤10 18 且1≤M≤10^91≤M≤10 9

.

输出格式

输出文件 output.txt 仅包含一个非负整数，表示 Concatenate (1 … N)Concatenate(1…N)
ModMod MM 的值。

思路：f[n]=f[n-1]*10^k+n
根据数据范围所以矩阵加速
构造矩阵 {f[n-1] ,n-1,1}->{f[n],n,1}
{10^k,0,0 *f[n-1],n-1,1
1 ,1,0
1,1,1}
对10^k要分段讨论

#include <stdio.h>
#include <string.h>
#define min(a,b) a>b?b:a
#define max(a,b) a>b?a:b
#define ll long long
ll read()
{
    ll s=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){s=s*10+ch-'0';ch=getchar();}
    return s*w;
}
ll n,p;
struct mat
{
    ll m[4][4];
}a,ans,e;
ll t[20];
mat quick_mul(mat x,mat y)
{
    mat c;
    memset(c.m,0,sizeof(c.m));        //一定要赋值为0，不然c不知道为多少
    for(int i=1;i<=3;i++)
        for(int j=1;j<=3;j++)
            for(int k=1;k<=3;k++)
                c.m[i][j]=(c.m[i][j]+x.m[i][k]%p*y.m[k][j]%p)%p;
    return c;
}
mat quick_pow(mat x,ll y)
{
    mat ans=e;
    while(y)
    {
        if(y&1)ans=quick_mul(ans,x);
        x=quick_mul(x,x);
        y>>=1;
    }
    return ans;
}
int main()
{
	for(int i=1;i<=3;i++)e.m[i][i]=1;
    n=read(),p=read();
    ans.m[1][1]=ans.m[1][2]=ans.m[1][3]=1;
    a.m[2][1]=a.m[2][2]=a.m[3][1]=a.m[3][2]=a.m[3][3]=1;
    t[0]=1;
    for(int i=1;i<=18;i++)t[i]=t[i-1]*10;
    for(int i=1;;i++)
    {
        a.m[1][1]=t[i]%p;
        ll tem=min(n,t[i]-1);
        tem=tem-t[i-1]+1;
        if(i==1)tem--;
        ans=quick_mul(ans,quick_pow(a,tem));
        if(t[i]-1>=n)break;
    }
    printf("%lld\n",ans.m[1][1]);
}

描述

石头游戏在一个 n 行 m 列 (1≤n,m≤8)
的网格上进行，每个格子对应一种操作序列，操作序列至多有10种，分别用0~9这10个数字指明。
操作序列是一个长度不超过6且循环执行、每秒执行一个字符的字符串。每秒钟，所有格子同时执行各自操作序列里的下一个字符。序列中的每个字符是以下格式之一：
数字0_{9：表示拿0}9个石头到该格子。
NWSE：表示把这个格子内所有的石头推到相邻的格子，N表示上方，W表示左方，S表示下方，E表示右方。 D：表示拿走这个格子的所有石头。
给定每种操作序列对应的字符串，以及网格中每个格子对应的操作序列，求石头游戏进行了 t
秒之后，石头最多的格子里有多少个石头。在游戏开始时，网格是空

的。
输入格式

第一行4个整数n, m, t, act。 接下来n行，每行m个字符，表示每个格子对应的操作序列。
最后act行，每行一个字符串，表示从0开始的每个操作序列。

输出格式

一个整数：游戏进行了t秒之后，所有方格中最多的格子有多少个石头。

样例输入

1 6 10 3
011112
1E
E
0

样例输出

3

#include <stdio.h>
#include <string.h>
#define min(a,b) a>b?b:a
#define max(a,b) a>b?a:b
#define ll long long
ll read()
{
    ll s=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){s=s*10+ch-'0';ch=getchar();}
    return s*w;
}
int n,m;
struct mat
{
    ll m[65][65];
}a[61],e,b,f;
int num(int i,int j)
{
	
    return (i-1)*m+j;
}
mat mul(mat x,mat y)
{
    mat c;
    memset(c.m,0,sizeof(c.m));
    for(int i=0;i<=n*m;i++)
        for(int j=0;j<=n*m;j++)
            for(int k=0;k<=n*m;k++)
                c.m[i][j]+=x.m[i][k]*y.m[k][j];
    return c;
}
mat quick_pow(mat x,int y)
{
    mat ans=e;
    while(y)
    {
        if(y&1)ans=mul(ans,x);
        x=mul(x,x);
        y>>=1;
    }
    return ans;
}
char str[65],work[11][7],len[11];
int main()
{
   n=read(),m=read();
   int t=read(),act=read();
   for(int i=1;i<=n;i++)
        scanf("%s",str+(i-1)*m+1);
    for(int i=0;i<act;i++)
        scanf("%s",work[i]),len[i]=strlen(work[i]);
    for(int i=1;i<=60;i++)
    {
        for(int j=1;j<=n;j++)
            for(int k=1;k<=m;k++)
            {
                a[i].m[0][0]=1;
                int now=str[num(j,k)]-'0';
                char tem=work[now][(i-1)%len[now]];
                int pos=num(j,k);
                if(tem>='0'&&tem<='9')
                {
                    a[i].m[0][pos]=tem-'0';
                    a[i].m[pos][pos]=1;
                }
                else if(tem=='N'&&j>1)
                    a[i].m[pos][num(j-1,k)]=1;
                else if(tem=='W'&&k>1)
                    a[i].m[pos][num(j,k-1)]=1;
                else if(tem=='S'&&j+1<=n)
                    a[i].m[pos][num(j+1,k)]=1;
                else if(tem=='E'&&k+1<=m)
                    a[i].m[pos][num(j,k+1)]=1;
            }
            
    }
    for(int i=0;i<=n*m;i++)
        e.m[i][i]=1;
    b=e;
    f.m[0][0]=1;
    for(int i=1;i<=60;i++)
        b=mul(b,a[i]);
    mat ans=quick_pow(b,t/60);
    for(int i=1;i<=t%60;i++)
        ans=mul(ans,a[i]);
    ans=mul(f,ans);
    ll ma=0;
    for(int i=0;i<=n*m;i++)
        for(int j=0;j<=m*m;j++)
            ma=max(ans.m[i][j],ma);
    printf("%lld",ma);
}