CF585F Digits of Number Pi

题目链接
如果子串必须是 S S S的前缀,那么是一个经典的单个字符串匹配问题,用KMP+数位DP记录目前匹配到的位置即可。
而任意子串就相当于 S S S的所有后缀的前缀,是多个字符串匹配问题,又注意到有用的只有每个后缀的前 d / 2 d/2 d/2个字符,对它们建立 A C AC AC自动机即可,效率 O ( 10 d 2 s ) O(10d^{2}s) O(10d2s)
代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll P=1e9+7;
const int N=55,M=50005,K=10;
char ch[M];
void read(int* a,int& n){
	scanf("%s",ch+1); n=strlen(ch+1);
	for(int i=1;i<=n;i++) a[i]=ch[i]-'0';
}
int n,m,a[M],to[M][K],fa[M],ct=1;
bool is[M];
void add(int l,int r){
	int u=1;
	for(int i=l;i<=r;i++){
		int x=a[i];
		if(!to[u][x]) to[u][x]=++ct;
		u=to[u][x];
	}
	is[u]=1;
}
void bfs(){
	queue<int>q; fa[1]=1;
	for(int i=0;i<K;i++){
		int& u=to[1][i];
		if(u) fa[u]=1,q.push(u);
		else u=1; 
	}
	while(!q.empty()){
		int u=q.front(); q.pop();
		for(int i=0;i<K;i++){
			int& v=to[u][i];
			if(v) fa[v]=to[fa[u]][i],q.push(v);
			else v=to[fa[u]][i];
		}
	}
}
ll f[N][M];
ll cnt(int* s){
	memset(f,0,sizeof(f));
	int u=1;
	for(int i=0;i<m;i++){
		for(int k=0;k<s[i+1];k++){
			int v=to[u][k]; if(is[v]) v=0;
			(f[i+1][v]+=1)%=P;
		}
		for(int j=0;j<=ct;j++){
			for(int k=0;k<K;k++) if(f[i][j]){
				int v=to[j][k]; if(is[v]) v=0;
				(f[i+1][v]+=f[i][j])%=P;
			}
		}
		u=to[u][s[i+1]]; if(is[u]) u=0;
	}
	//cout<<(f[m][0]+(u==0))%P<<endl;
	return (f[m][0]+(u==0))%P;
}
int x[N],y[N];
int main()
{
	//freopen(".in","r",stdin);
	//freopen(".out","w",stdout);
	read(a,n),read(x,m),read(y,m);
	for(int i=m;i;i--){
		if(!x[i]) x[i]=9;
		else{
			x[i]--; break;
		}
	}
	for(int i=1;i+m/2-1<=n;i++) add(i,i+m/2-1);
	bfs(); //return 0;
	/*for(int i=1;i<=ct;i++){
		cout<<"u="<<i<<" fa="<<fa[i]<<endl;
		for(int j=0;j<K;j++) cout<<to[i][j]<<" ";
		puts("");
	}*/
	cout<<(cnt(y)+P-cnt(x))%P<<endl;
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值