1. 程式人生 > 實用技巧 >P3311 [SDOI2014] 數數

P3311 [SDOI2014] 數數

P3311 [SDOI2014] 數數

感覺我太不自信了,其實思路都是對的,就是不敢寫,想到了也不寫。

這題綜合了一些東西。

首先非常套路對於所有模式串建AC自動機,維護出 \(Trie\) 圖上哪些點不能到達

如何限制 \(\le n\) ?這就很像數位dp了,dp狀態多一維 \(0/1\) 表示是否頂到上界

我開心地寫完了這個dp卻發現WA40

到討論區發現被下面那組hack了

2000
1
001

很明顯是前導零的問題,但是我懷疑是別的問題懶得改。最後就是這個問題。

如果還有前導零那麼當前點設為 \(0\) ,否則走 \(Trie\) 圖上的出邊。

乾脆換成記搜寫數位dp了

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef double db;
#define x first
#define y second
#define sz(v) (int)v.size()
#define pb(x) push_back(x)
#define mkp(x,y) make_pair(x,y)
inline int read(){
	int x=0,f=1;char c=getchar();
	while(!isdigit(c)){if(c=='-')f=0;c=getchar();}
	while(isdigit(c))x=x*10+c-'0',c=getchar();
	return f?x:-x;
}
#define mod 1000000007
#define N 1505
char n[N];
int m,Len,ch[N][10],tot,dp[N][N][2][2],fail[N];
bool en[N];
void fmod(int&x){x-=mod,x+=x>>31&mod;}
void insert(char*str){
	int u=0,len=strlen(str);
	for(int i=0;i<len;++i){
		int c=str[i]-'0';
		if(!ch[u][c])ch[u][c]=++tot;
		u=ch[u][c];
	}
	en[u]=1;
}
void build_fail(){
	queue<int>q;
	for(int i=0;i<10;++i)if(ch[0][i])q.push(ch[0][i]);
	while(!q.empty()){
		int u=q.front();q.pop();
		en[u]|=en[fail[u]];
		for(int i=0;i<10;++i)
			if(ch[u][i])fail[ch[u][i]]=ch[fail[u]][i],q.push(ch[u][i]);
			else ch[u][i]=ch[fail[u]][i];
	}
}
int dfs(int u,int pos,int limit,int lead){
	if(!pos)return !en[u];
	if(en[u])return 0;
	if(~dp[u][pos][limit][lead])return dp[u][pos][limit][lead];
	int up=limit?n[pos-1]-'0':9,res=0;
	for(int i=0;i<=up;++i)fmod(res+=dfs(lead&&!i?0:ch[u][i],pos-1,limit&&i==up,lead&&!i));
	return dp[u][pos][limit][lead]=res;
}
signed main(){
	scanf("%s%d",n,&m),Len=strlen(n);
	for(int i=1;i<=m;++i){
		static char str[N];
		scanf("%s",str),insert(str);
	}
	build_fail();
	reverse(n,n+Len);
	memset(dp,-1,sizeof(dp));
	printf("%d\n",(dfs(0,Len,1,1)+mod-1)%mod);
	return 0;
}