[BZOJ 1009] HNOI 2008 GT考试 · 动态规划+矩阵乘法+KMP

KMP预处理，矩阵乘法加速DP。

刷了这么多矩阵乘法怎么感觉还是没什么进步。。。

这题就算是一条比较好的用矩阵加速DP的题目了，感觉还是有点难度的可是为什么AC率那么高。。。难道是我太弱了么。。。

还是就不要班门弄斧了，两个神犇的友链，写的很好：

/article/2352140.html /article/2649713.html

笔者写这题最大的收获就是顺便终于搞懂了KMP，如果不会的同学可以看下面这个博客，写的很详细：

/article/1362918.html

像笔者这种以前背过程序（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