传送门
题解:
对代价取一个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;
}