P5319-[BJOI2019]奥术神杖【0/1分数规划,AC自动机,dp】
生活随笔
收集整理的這篇文章主要介紹了
P5319-[BJOI2019]奥术神杖【0/1分数规划,AC自动机,dp】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/P5319
題目大意
一個長度為nnn的串TTT,用0~90\sim 90~9填充所有的...。
然后給出mmm個串和它們的價值。
一個填充方案的價值等于:若TTT中出現了ccc個給出的串,那價值等于它們的價值乘積開ccc次根。
1≤m≤1501,1≤Vi≤1091\leq m\leq 1501,1\leq V_i\leq 10^91≤m≤1501,1≤Vi?≤109
解題思路
ans=∏Vicans=\sqrt[c]{\prod V_i}ans=c∏Vi??
ln?ans=1c∑(ln?Vi)\ln ans=\frac{1}{c}\sum (\ln V_i)lnans=c1?∑(lnVi?)
然后就是一個0/10/10/1分數規劃問題了,因為要匹配,所以要先跑一個ACACAC自動機上dpdpdp就好了。
時間復雜度O(n2log?109)O(n^2\log 10^9)O(n2log109)
code
#include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #include<queue> using namespace std; const int N=1510; const double eps=1e-6; int n,m,cnt,ch[N][10],fail[N],num[N]; double w[N],v[N],f[N][N]; pair<int,int>pre[N][N];queue<int> q; char T[N],S[N]; void Insert(char *s,double val){int l=strlen(s),x=0;for(int i=0;i<l;i++){int c=s[i]-'0';if(!ch[x][c])ch[x][c]=++cnt;x=ch[x][c];}v[x]+=val;num[x]++;return; } void Build(){for(int i=0;i<10;i++)if(ch[0][i])q.push(ch[0][i]);while(!q.empty()){int x=q.front();q.pop();v[x]+=v[fail[x]];num[x]+=num[fail[x]];for(int i=0;i<10;i++){if(!ch[x][i])ch[x][i]=ch[fail[x]][i];else{fail[ch[x][i]]=ch[fail[x]][i];q.push(ch[x][i]);}}}return; } int check(double mid){for(int i=0;i<=cnt;i++)w[i]=v[i]-mid*num[i];for(int i=0;i<=n;i++)for(int j=0;j<=cnt;j++)f[i][j]=-1e100;f[0][0]=0;for(int i=1;i<=n;i++){for(int j=0;j<=cnt;j++){if(f[i-1][j]<=-1e99)continue;if(T[i]=='.'){for(int k=0;k<10;k++){int y=ch[j][k];if(f[i][y]<f[i-1][j]+w[y])f[i][y]=f[i-1][j]+w[y],pre[i][y].first=j,pre[i][y].second=k;}}else{int k=T[i]-'0',y=ch[j][k];if(f[i][y]<f[i-1][j]+w[y])f[i][y]=f[i-1][j]+w[y],pre[i][y].first=j,pre[i][y].second=k;}}}int root=0;for(int j=0;j<=cnt;j++)if(f[n][j]>f[n][root])root=j;return root; } void print(int k,int x){if(k<1)return; print(k-1,pre[k][x].first);printf("%d",pre[k][x].second);return; } int main() {scanf("%d%d",&n,&m);scanf("%s",T+1);for(int i=1;i<=m;i++){int x;scanf("%s",S);scanf("%d",&x);Insert(S,log(x));}Build();double l=0,r=21;while(r-l>eps){double mid=(l+r)/2.0;int root=check(mid);if(f[n][root]>eps)l=mid;else r=mid;}int root=check(l);print(n,root);return 0; }總結
以上是生活随笔為你收集整理的P5319-[BJOI2019]奥术神杖【0/1分数规划,AC自动机,dp】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: P4585-[FJOI2015]火星商店
- 下一篇: YbtOJ#463-序列划分【二分答案,