阿申准备报名参加GT考试，准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字。他的不吉利数学A1A2...Am(0<=Ai<=9)有M位，不出现是指X1X2...Xn中没有恰好一段等于A1A2...Am. A1和X1可以为0
标签：
阿申准备报名参加GT考试，准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字。他的不吉利数学A1A2...Am(0<=Ai<=9)有M位，不出现是指X1X2...Xn中没有恰好一段等于A1A2...Am. A1和X1可以为0
第一行输入N,M,K.接下来一行输入M位的数。 100%数据N<=10^9,M<=20,K<=1000 40%数据N<=1000 10%数据N<=6
阿申想知道不出现不吉利数字的号码有多少种，输出模K取余的结果.
题解：做到这道题我很伤心。。。是ac(kmp)+dp
翻了这么多题解还是cjk讲的最好。。。请允许我转一下。。。
字符串上的动态规划：
按顺序处理准考证号每一位，
设f[i][j]表示：准考证号前i位中 后j位与不吉利数的前j位相同时，前i位的方案数
那么答案ans=f[n][0]+f[n][1]+…+f[n][m-1]
f[i][j]的准确含义：
1.f[i][j]表示的每种方案不仅与其后j位有关，还应保证不含不吉利数
2.为避免重复，f[i][j]表示的每种方案都不含长度大于j且与不吉利数的前缀相同 的后缀
否则就会出现：从1到m标号，不吉利数为123124时，f[i][2]计数的方案包含f[i][5]计数的方案 的情况
状态转移：
f[i][j]只能由f[i-1][k]得到，相当于填完第i-1位后，将其后缀k(长为k的后缀)后面新添一位num，之后这个i位数的 与不吉利数前缀相同的最长后缀是：后缀j
i>=1时：f[i][j]=f[i-1][0]*a[0][j]+f[i-1][1]*a[1][j]+…+f[i-1][m-1]*a[m-1][j]
比如：还是假设不吉利数为123124，那么 f[i][3]=f[i-1][2]+f[i-1][5]，因为 f[i-1][2]末尾的*****12不能是**12312，所以需要f[i-1][5]补充
但若不吉利数为123123，那么 f[i][3]=f[i-1][2]，因为 f[i][3]末尾的*****123不能是**123123
i==0时：f[0][0]=1,f[0][其他]=0
其中，a[k][j]就表示上面提到的num能取几个值，可以用kmp算法预处理出来，它是一个矩阵
这样就可以不重不漏地计数了
再来个矩阵加速：f[i][j]求法是个线性齐次递推式，可以构造成矩阵，然后加一下速。
这道题要留坑的。。。因为还不是很懂QAQ
1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<algorithm> 5 #include<queue> 6 #include<cstring> 7 #define PAU putchar(‘ ‘) 8 #define ENT putchar(‘\n‘) 9 using namespace std; 10 const int maxn=30; 11 int n,mod,len,nxt[maxn];char str[maxn]; 12 struct matrix{int f[maxn][maxn];matrix(){memset(f,0,sizeof(f));}}m,ans; 13 matrix operator*(const matrix&a,const matrix&b){ 14 matrix c; 15 for(int i=0;i<len;i++) 16 for(int j=0;j<len;j++) 17 for(int k=0;k<len;k++) 18 c.f[i][j]=(c.f[i][j]+a.f[i][k]*b.f[k][j])%mod; 19 return c; 20 } 21 void kmp(){ 22 int j=0;nxt[1]=0; 23 for(int i=2;i<=len;i++){ 24 while(j&&str[j+1]!=str[i]) j=nxt[j]; 25 if(str[j+1]==str[i])j++;nxt[i]=j; 26 }return; 27 } 28 void dp(){ 29 for(int i=0;i<len;i++) 30 for(int j=0;j<10;j++){ 31 int k; 32 for(k=i;k;k=nxt[k]) if(str[k+1]-‘0‘==j) break; 33 if(str[k+1]-‘0‘==j) k++; 34 m.f[i][k]=(m.f[i][k]+1)%mod; 35 } 36 return; 37 } 38 void pow(){ 39 for(int i=0;i<len;i++) ans.f[i][i]=1; 40 for(;n;n>>=1,m=m*m)if(n&1) ans=ans*m;return; 41 } 42 inline int read(){ 43 int x=0,sig=1;char ch=getchar(); 44 while(!isdigit(ch)){if(ch==‘-‘) sig=-1;ch=getchar();} 45 while(isdigit(ch)) x=10*x+ch-‘0‘,ch=getchar(); 46 return x*sig; 47 } 48 inline void write(int x){ 49 if(x==0){putchar(‘0‘);return;}if(x<0) putchar(‘-‘),x=-x; 50 int len=0,buf[15];while(x) buf[len++]=x%10,x/=10; 51 for(int i=len-1;i>=0;i--) putchar(buf[i]+‘0‘);return; 52 } 53 int main(){ 54 n=read();len=read();mod=read();scanf("%s",str+1);kmp();dp();pow(); 55 int res=0;for(int i=0;i<len;i++) res=(res+ans.f[0][i])%mod; 56 write(res); 57 return 0; 58 }
标签：
原文地址：http://www.cnblogs.com/chxer/p/4642446.html
