• ! BJOI2019奥术神杖


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

    大菜鸡又去翻题解

    两边取对数!!!

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

    [lnans=frac1csum_{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);
    } 
    
  • 相关阅读:
    1+X云计算(中级) 单节点部署应用商城系统(gpmall)
    1+X云计算 应用商城系统(gpmall)-遇到的问题以及解决办法
    1+X云计算平台运维与开发(中级)eNSP A~E卷 试题+答案
    vi&vim 基本使用方法
    yum针对软件包操作的常用命令
    本地yum源配置
    SpringCloud微服务初体验
    SpringBoot自定义注解拦截器,实现登录token验证
    MySQL建立SSL连接问题,设置useSSL=false显式禁用SSL,或者设置useSSL=true
    TPL事务
  • 原文地址:https://www.cnblogs.com/aurora2004/p/12750601.html
Copyright © 2020-2023  润新知