【BZOJ3530】数数(AC自动机,动态规划)

题面

BZOJ

题解

很套路的\(AC\)自动机+\(DP\)

首先,如果长度小于\(N\)

就不存在任何限制

直接大力\(DP\)

然后强制限制不能走到带有标记的点上面

如果长度恰好为\(N\)的长度

那么,要考虑是否恰好卡在范围里面

于是\(DP\)状态多记一维

表示是否卡在范围里面

最后求一下和就行啦

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define MAX 2000
#define MOD 1000000007
struct Node
{
int vis[10];
int lt,fail;
}t[MAX];
int tot,m;
char N[MAX],ch[MAX];
int f[MAX][MAX][2],g[MAX][MAX];
void insert(char *s)
{
scanf("%s",s+1);
int l=strlen(s+1),now=0;
for(int i=1;i<=l;++i)
{
if(!t[now].vis[s[i]-48])
t[now].vis[s[i]-48]=++tot;
now=t[now].vis[s[i]-48];
}
t[now].lt=1;
}
void GetFail()
{
queue<int> Q;
for(int i=0;i<=9;++i)
if(t[0].vis[i])Q.push(t[0].vis[i]);
while(!Q.empty())
{
int u=Q.front();Q.pop();
t[u].lt|=t[t[u].fail].lt;
for(int i=0;i<=9;++i)
if(t[u].vis[i])
t[t[u].vis[i]].fail=t[t[u].fail].vis[i],Q.push(t[u].vis[i]);
else t[u].vis[i]=t[t[u].fail].vis[i];
}
}
int main()
{
scanf("%s",N+1);
scanf("%d",&m);
while(m--)insert(ch);
GetFail();
int l=strlen(N+1);
long long ans=0;
g[0][0]=1;
for(int i=0;i<l;++i)
for(int u=0;u<=tot;++u)
if(!t[u].lt)
{
for(int k=0;k<=9;++k)
if(!t[t[u].vis[k]].lt)
{
if(!i&&!k)continue;
(g[i+1][t[u].vis[k]]+=g[i][u])%=MOD;
}
}
for(int i=1;i<l;++i)
for(int j=0;j<=tot;++j)
ans=(ans+g[i][j])%MOD;
f[0][0][1]=1;
for(int i=0;i<l;++i)
for(int u=0;u<=tot;++u)
if(!t[u].lt)
{
for(int k=0;k<=9;++k)
if(!t[t[u].vis[k]].lt)
{
if(!i&&!k)continue;
(f[i+1][t[u].vis[k]][0]+=f[i][u][0])%=MOD;
if(k<N[i+1]-48)(f[i+1][t[u].vis[k]][0]+=f[i][u][1])%=MOD;
if(k==N[i+1]-48)(f[i+1][t[u].vis[k]][1]+=f[i][u][1])%=MOD;
}
}
for(int i=0;i<=tot;++i)
ans=(ans+f[l][i][0])%MOD,ans=(ans+f[l][i][1])%MOD;
printf("%lld\n",ans);
return 0;
}

最新文章

  1. thread.join 从异步执行变成同步
  2. Unattend.xml应答文件制作(WISM)
  3. Mysql 创建用户 授权
  4. Sencha Touch 手机移动开发框架 HTML5 项目压缩方案;
  5. 5分钟弄懂Docker!
  6. e2e 自动化集成测试 架构 实例 WebStorm Node.js Mocha WebDriverIO Selenium Step by step (六) 自动化测试结构小节
  7. nyoj 17 单调递增最长子序列
  8. TCP 的那些事儿(下)
  9. Merge into的使用详解-你Merge了没有【转】
  10. python学习day3------列表、元组、字符串操作
  11. 两层fragment嵌套时出现空白,(收藏别人的)
  12. windows CMD常用命令
  13. Shiro入门 - 通过ini文件进行认证
  14. 1.HTTP协议|web框架
  15. C# WPF 进度条,根据读取数据显示进度条进度,根据Excel文件读取数据,进度条样式
  16. HDU 1166 敌兵布阵(线段树单点更新,区间查询)
  17. hdu 4169 二分匹配最大独立集 ***
  18. WPF中的3D变换PlaneProjection
  19. byte字节数组的压缩
  20. C#存取数据库图片

热门文章

  1. Burp_用户名密码爆破
  2. python学习:使用正则收集ip信息
  3. file_put_contents写入文字换行
  4. 《InsideUE4》UObject(六)类型系统代码生成重构-UE4CodeGen_Private
  5. 读取Pdm文件内容(含源码)
  6. python2 =&gt; python3 踩坑集合
  7. Sublime Text编辑远程Linux服务器上的文件
  8. nyoj720 项目安排 二分+dp
  9. Nginx 调优经验记录
  10. C++ 监测磁盘空间