依旧看人代码写,不过我觉得自己慢慢写一个也可以写成?

原题:

阿申准备报名参加GT考试,准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字。
他的不吉利数学A1A2...Am(0<=Ai<=9)有M位,不出现是指X1X2...Xn中没有恰好一段等于A1A2...Am. A1和X1可以为
0

N<=10^9,M<=20,K<=1000

一开始完全是懵逼的,后来找到两个解释结合一下就理解了

用a(i,j)表示原串走到i,不吉利的串走到j的方案数

设一个数组b(i,j),表示从(i,j)转移到(i+1,k)的方案数,这个可以用kmp处理,先kmp出next,然后枚举i和j,根据kmp求出k并将b(i+1,k)++

然后这个b表示的就是a中的元素下一步会贡献到哪里

这个东西比较玄,只能勉强意会……

NOIP吧里有一种解释,虽然写法似乎和我的不太一样,引导思路效果不错:

“构造转移矩阵A和列向量B,B=(f[0][0],f[0][1],...,f[0][m]),A可以由DP得到,那么A*B的结果就是(f[1][0],f[1][1],...,f[1][m]),所以A^n*B的结果就是(f[n][0],f[n][1],...,f[n][m])”

然后矩阵快速乘即可

小技巧:使用a&1判断奇偶

代码:

 #include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
int n,m,mo,s[];
int next[];
int b[][],a[][],c[][];
int ans=;
void kmp(){
int temp=; next[]=;
for(int i=;i<=m;i++){
while(temp && s[temp+]!=s[i]) temp=next[temp];
if(s[temp+]==s[i]) temp++;
next[i]=temp;
}
}
void get_b(){
int temp;
for(int i=;i<m;i++){//注意j枚举的是下一位,因为这里j是对下一位的转移,所以矩阵乘法从0开始写,比较方便
a[i][i]=;
for(int j=;j<=;j++){
temp=i;
while(temp && s[temp+]!=j) temp=next[temp];
if(j==s[temp+]) b[i][temp+]+=;
else b[i][]+=;
}
}
}
void fast_mi(){
while(n){
if(n&){//快速判断奇偶
for(int i=;i<m;i++)
for(int j=;j<m;j++){
c[i][j]=;//反正也要枚举,就不用memset了
for(int k=;k<m;k++)
c[i][j]=(c[i][j]+b[i][k]*a[k][j])%mo;
}
for(int i=;i<m;i++)
for(int j=;j<m;j++)
a[i][j]=c[i][j];
}
n>>=;
for(int i=;i<m;i++)
for(int j=;j<m;j++){
c[i][j]=;
for(int k=;k<m;k++)
c[i][j]=(c[i][j]+b[i][k]*b[k][j])%mo;
}
for(int i=;i<m;i++)
for(int j=;j<m;j++)
b[i][j]=c[i][j];
/*for(int i=0;i<m;i++){
for(int j=0;j<m;j++)
cout<<a[i][j]<<" ";
cout<<endl;
}*/
}
}
int main(){//freopen("ddd.in","r",stdin);
memset(b,,sizeof(b));
cin>>n>>m>>mo;
for(int i=;i<=m;i++){
scanf("%c",&s[i]); while(s[i]<''||s[i]>'') scanf("%c",&s[i]);
s[i]-='';
}
kmp(); get_b();
fast_mi();
for(int i=;i<m;i++) ans=(ans+a[][i])%mo;
cout<<ans<<endl;
return ;
}

【BZOJ1009】【HNOI2008】GT考试的更多相关文章

  1. [BZOJ1009] [HNOI2008] GT考试(KMP+dp+矩阵快速幂)

    [BZOJ1009] [HNOI2008] GT考试(KMP+dp+矩阵快速幂) 题面 阿申准备报名参加GT考试,准考证号为N位数X1X2-.Xn,他不希望准考证号上出现不吉利的数字.他的不吉利数学A ...

  2. BZOJ1009 [HNOI2008]GT考试 矩阵

    去博客园看该题解 题目 [bzoj1009][HNOI2008]GT考试 Description 阿申准备报名参加GT考试,准考证号为N位数X1X2….Xn(0<=Xi<=9),他不希望准 ...

  3. bzoj1009 [HNOI2008] GT考试 矩阵乘法+dp+kmp

    1009: [HNOI2008]GT考试 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 4542  Solved: 2815[Submit][Statu ...

  4. [Bzoj1009][HNOI2008]GT考试(KMP)(矩乘优化DP)

    1009: [HNOI2008]GT考试 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 4309  Solved: 2640[Submit][Statu ...

  5. bzoj1009: [HNOI2008]GT考试(kmp+矩阵乘法)

    1009: [HNOI2008]GT考试 题目:传送门 题解: 看这第一眼是不是瞬间想起组合数学??? 没错...这样想你就GG了! 其实这是一道稍有隐藏的矩阵乘法,好题! 首先我们可以简化一下题意: ...

  6. [bzoj1009](HNOI2008)GT考试 (kmp+矩阵快速幂加速递推)

    Description 阿 申准备报名参加GT考试,准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字.他的不吉利数学 A1A2...Am(0&l ...

  7. [bzoj1009][HNOI2008]GT考试

    Description 阿申准备报名参加考试,准考证号为位数,他不希望准考证号上出现不吉利的数字. 他的不吉利数学有位,不出现是指中没有恰好一段等于. 可以为. Input 第一行输入.接下来一行输入 ...

  8. [BZOJ1009] [HNOI2008] GT考试 (KMP & dp & 矩阵乘法)

    Description 阿申准备报名参加GT考试,准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字. 他的不吉利数学A1A2...Am(0< ...

  9. bzoj1009: [HNOI2008]GT考试 ac自动机+矩阵快速幂

    https://www.lydsy.com/JudgeOnline/problem.php?id=1009 阿申准备报名参加GT考试,准考证号为N位数X1X2....Xn(0<=Xi<=9 ...

  10. BZOJ1009:[HNOI2008]GT考试(AC自动机,矩乘DP)

    Description 阿申准备报名参加GT考试,准考证号为N位数X1X2....Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字. 他的不吉利数学A1A2...Am(0< ...

随机推荐

  1. android webview实战

    webSettings = wvShowProduce.getSettings();//设置WebView属性,能够执行Javascript脚本webSettings.setJavaScriptEna ...

  2. ubuntu 添加源

    edit file :  /etc/apt/sources.list add: deb http://mirrors.163.com/ubuntu/ intrepid main restricted ...

  3. 自定义圆的半径layout

    <RelativeLayout xmlns:android="http://schemas.android.com/apk/res/android"    xmlns:too ...

  4. 在shell脚本中进行条件控制以及使用循环

    转载请标明:http://www.cnblogs.com/winifred-tang94/ if条件语句语法: if [ 条件表达式 ] then 代码 else 代码 fi 注意:在上面的if条件语 ...

  5. web api post传一个参数时 值永远是null

    这个问题纠结了我一个早上,不管用什么样的传参方法,走到控制器中,那个参数永远不变的等于null 在网上找了很多解决方案 上面这个是从网上截图的,第一:要将参数标记为[FromBody],变为简单参数 ...

  6. java generic type

    java generic type: 类型安全.类型参数化.不需要强制转换

  7. C#操作txt问件,进行清空添加操作

    //把txt清空 FileStream stream = File.Open(Adr,FileMode.OpenOrCreate,FileAccess.Write); stream.Seek(, Se ...

  8. TVB-Gone 红外编码方法

    每次都重新推算一遍,年纪大了还是记录一下吧(硬件版的TV-B-Gone的压缩编码跳过此段往下看) N900上有个TVB-Gone的程序,但是程序的红外编码往往都是对应的国外的电视,好多国产电视都不支持 ...

  9. flex css 布局

    http://www.w3cplus.com/css3/flexbox-basics.html

  10. 使用httputils上传图片到服务器

    //创建httpUtils对象 HttpUtils mRegHttpUtils = new HttpUtils(); //图片路径 String path = "/sdcard/Downlo ...