KMP预处理,矩阵乘法加速DP。
刷了这么多矩阵乘法怎么感觉还是没什么进步。。。
这题就算是一条比较好的用矩阵加速DP的题目了,感觉还是有点难度的可是为什么AC率那么高。。。难道是我太弱了么。。。
还是就不要班门弄斧了,两个神犇的友链,写的很好:
http://blog.csdn.net/popoqqq/article/details/40188173 http://blog.csdn.net/cjk_cjk/article/details/43038377
笔者写这题最大的收获就是顺便终于搞懂了KMP,如果不会的同学可以看下面这个博客,写的很详细:
http://blog.csdn.net/v_july_v/article/details/7041827
像笔者这种以前背过程序(orz蒟蒻)有一定基础但是不怎么理解的可以看这个比较简洁的:
http://kb.cnblogs.com/page/176818/
#include <stdio.h>
#include <algorithm>
#include <string.h>
#include <iostream>
#include <string>
using namespace std;
int n,m,k,p[100];
char s[100];
struct M{
int num[100][100];
void init(){
memset(num,0,sizeof num);
}
friend M operator *(M a,M b){
M c;c.init();
for (int i=0;i<m;i++)
for (int j=0;j<m;j++)
for (int t=0;t<m;t++)
c.num[i][j]=(c.num[i][j]+a.num[i][t]*b.num[t][j])%k;
return c;
}
}a,b;
void KMP(){
int j=0;p[1]=0;
for (int i=2;i<=m;i++){
while (j && s[j+1]!=s[i]) j=p[j];
if (s[j+1]==s[i]) j++;
p[i]=j;
}
j=0;
for (int i=0;i<m;i++)
for (int t=0;t<=9;t++){
j=i;
while (j && s[j+1]!=(t+'0')) j=p[j];
if (s[j+1]==(t+'0')) a.num[i][j+1]++;
else a.num[i][0]++;
}
}
int main(){
cin>>n>>m>>k;
scanf("%s",s+1);
KMP();
for (int i=0;i<m;i++) b.num[i][i]=1;
for (;n;n>>=1,a=a*a)
if (n&1) b=b*a;
int ans=0;
for (int i=0;i<m;i++)
ans=(ans+b.num[0][i])%k;
cout<<ans<<endl;
return 0;
}
突然发现这篇博客写得好敷衍。。。各位不要见怪。。。XD