【BJOI2019】奥术神杖（分数规划）（AC自动机）

传送门

---

题解：

对代价取一个Ln之后就是分数规划的裸题，匹配按照套路使用AC自动机即可。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

inline bool ckmax(double &a,cs double &b){return a<b?(a=b,true):false;}

cs int N=1.5e3+7;

int n,m;

int son[N][10],fail[N],ps[N],tot;
double x[N],v[N],dp[N][N];
char s[N],t[N];int len;

inline int ins(){
	int u=0;for(char *c=t+1;*c;++c)
	u=son[u][48^*c]?son[u][48^*c]:(son[u][48^*c]=++tot);
	return u;
}

int q[N],qn;
inline void build_AC(){
	for(int re i=0;i<10;++i)if(son[0][i])q[++qn]=son[0][i];
	for(int re i=1;i<=qn;++i){
		int u=q[i];for(int re c=0;c<10;++c)
		son[u][c]?(q[++qn]=son[u][c],fail[son[u][c]]=son[fail[u]][c]):(son[u][c]=son[fail[u]][c]);
	}
}
inline double calc(double x){
	std::fill(v,v+tot+1,0);
	for(int re i=1;i<=m;++i)v[ps[i]]+=::x[i]-x;
	for(int re i=1;i<=qn;++i)v[q[i]]+=v[fail[q[i]]];
	for(int re i=0;i<=n;++i)std::fill(dp[i],dp[i]+tot+1,-1e9);dp[0][0]=0;
	for(int re i=1;i<=n;++i){
		for(int re j=0;j<=tot;++j)if(dp[i-1][j]>-1e9)
		if(s[i]=='.'){
			for(int re c=0;c<10;++c)
			ckmax(dp[i][son[j][c]],dp[i-1][j]);
		}else {
			int c=s[i]^48;
			ckmax(dp[i][son[j][c]],dp[i-1][j]);
		}
		for(int re j=0;j<=tot;++j)if(dp[i][j]>-1e9)dp[i][j]+=v[j];
	}
	double res=*std::max_element(dp[n],dp[n]+tot+1);
	return res;
}

int pre[N][N];
char ch[N][N];

signed main(){
#ifdef zxyoi
	freopen("arcana.in","r",stdin);
#endif
	scanf("%d%d%s",&n,&m,s+1);
	for(int re i=1;i<=m;++i){
		int v;scanf("%s%d",t+1,&v);
		x[i]=log2(v);ps[i]=ins();
	}build_AC();double l=0,r=22,mid;
	while(l+1e-6<r)calc(mid=(l+r)*0.5)>0?l=mid:r=mid-1e-6;
	std::fill(v,v+tot+1,0);
	for(int re i=1;i<=m;++i)v[ps[i]]+=::x[i]-l;
	for(int re i=1;i<=qn;++i)v[q[i]]+=v[fail[q[i]]];
	for(int re i=0;i<=n;++i)std::fill(dp[i],dp[i]+tot+1,-1e9);dp[0][0]=0;
	for(int re i=1;i<=n;++i){
		for(int re j=0;j<=tot;++j)if(dp[i-1][j]>-1e9)
		if(s[i]=='.'){
			for(int re c=0;c<10;++c)
			if(ckmax(dp[i][son[j][c]],dp[i-1][j])){
				pre[i][son[j][c]]=j;
				ch[i][son[j][c]]=c;
			}
		}else {
			int c=s[i]^48;
			if(ckmax(dp[i][son[j][c]],dp[i-1][j])){
				pre[i][son[j][c]]=j;
				ch[i][son[j][c]]=c;
			}
		}
		for(int re j=0;j<=tot;++j)if(dp[i][j]>-1e9)dp[i][j]+=v[j];
	}
	int cur=std::max_element(dp[n],dp[n]+tot+1)-dp[n];
	for(int re i=n;i;--i){
		t[i]=ch[i][cur]^48;cur=pre[i][cur];
	}t[n+1]='\0';puts(t+1);
	return 0;
}