日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

P5319-[BJOI2019]奥术神杖【0/1分数规划,AC自动机,dp】

發布時間:2023/12/3 编程问答 32 豆豆
生活随笔 收集整理的這篇文章主要介紹了 P5319-[BJOI2019]奥术神杖【0/1分数规划,AC自动机,dp】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

正題

題目鏈接:https://www.luogu.com.cn/problem/P5319


題目大意

一個長度為nnn的串TTT,用0~90\sim 909填充所有的...

然后給出mmm個串和它們的價值。

一個填充方案的價值等于:若TTT中出現了ccc個給出的串,那價值等于它們的價值乘積開ccc次根。

1≤m≤1501,1≤Vi≤1091\leq m\leq 1501,1\leq V_i\leq 10^91m1501,1Vi?109


解題思路

ans=∏Vicans=\sqrt[c]{\prod V_i}ans=cVi??
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】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。