标签：BJOI2019 sz 奥术 int 神杖 re while inline id

\(f_{i,u,v},i\)位，AC自动机上\(u\)位置，\(v\)次匹配的最大值

大菜鸡又去翻题解

两边取对数！！！

\[ams=\sqrt[c]{\prod_{i=1}^cw_i} \]

\[lnans=\frac1c\sum_{i=1}^clnw_i \]

分数规划，二分DP即可

#include<bits/stdc++.h>
using namespace std;
inline int read(){
	int x=0,f=1;char c=getchar();
	while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
	while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
	return f==1?x:-x;
}
const int N=1504;
const double eps=1e-8,inf=1e9;
vector<int>e[N];
int n,m,tot=1,sz[N],fail[N],ch[N][12],fi[N][N],fu[N][N];
double f[N][N],w[N];
void dfs(int x){
	for(auto v:e[x]){
		sz[v]+=sz[x];
		w[v]+=w[x];
		dfs(v);
	}
}
#define re register
inline void buildAC(){
	for(re int i=0;i<=9;++i)ch[0][i]=1;
	queue<int>q;
	q.push(1);
	while(!q.empty()){
		int x=q.front();q.pop();
		e[fail[x]].push_back(x); 
		for(re int i=0,u;i<=9;++i){
			u=ch[x][i];
			if(!u){
				ch[x][i]=ch[fail[x]][i];
				continue;
			} 
			fail[u]=ch[fail[x]][i];
			q.push(u);
		}
	}
	dfs(1);
} 
inline int cmp(double x,double y){
	if(fabs(x-y)<eps)return 0;
	return (x<y)?-1:1;
}
inline void update(int p,int x,int f1,int f2,double mid){
	double v=f[p-1][f1]+w[x]-mid*sz[x];
	if(v>f[p][x]){
		f[p][x]=v;
		fi[p][x]=f1;
		fu[p][x]=f2;
	}
} 
char s[N],c[N];
inline bool check(double mid){
	for(re int i=0;i<=n;++i)
		for(re int j=1;j<=tot;++j)f[i][j]=-inf;
	f[0][1]=0;
	for(re int i=0;i<n;++i)
		for(re int j=1;j<=tot;++j){
			if(!cmp(f[i][j],-inf))continue;
			if(s[i+1]=='.')for(int u=0;u<=9;u++)
				update(i+1,ch[j][u],j,u,mid);
			else update(i+1,ch[j][s[i+1]^48],j,s[i+1]^48,mid);
		}
	for(re int j=1;j<=tot;++j)
		if(cmp(f[n][j],0)==1)return 1;
	return 0;
}
void out(int p,int u){
	if(!p)return;
	out(p-1,fi[p][u]);
	putchar(fu[p][u]|48);
}
int main(){
	cerr<<log(inf);
	n=read();m=read();
	scanf("%s",s+1);
	for(re int i=1,p,u;i<=m;++i){
		scanf("%s",c+1);
		p=1;
		for(re int j=1,len=strlen(c+1);j<=len;++j){
			u=(c[j]^48);
			if(!ch[p][u])ch[p][u]=++tot;
			p=ch[p][u];
		}
		w[p]+=log(read());
		sz[p]++;
	}
	buildAC();
	double l=0,r=24,mid,ans=0;
	for(re int T=1;T<=40;++T){
		mid=(l+r)/2;
		if(check(mid))ans=l=mid;
		else r=mid;
	}
	check(ans);
	int id=1;
	for(re int i=2;i<=tot;++i)
		if(f[n][i]>f[n][id])id=i;
	out(n,id);
	return (0-0);
} 

标签：BJOI2019,sz,奥术,int,神杖,re,while,inline,id
来源： https://www.cnblogs.com/aurora2004/p/12750601.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。