CF183D-T-shirtx【dp,贪心】
生活随笔
收集整理的這篇文章主要介紹了
CF183D-T-shirtx【dp,贪心】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/CF183D
題目大意
nnn個人,mmm種衣服,給出每個人喜歡某件衣服的概率,你可以選擇nnn件衣服帶過去(可以重復款式)。求最大化能拿到喜歡衣服人的期望數量。
1≤n≤3000,1≤m≤3001\leq n\leq 3000,1\leq m\leq 3001≤n≤3000,1≤m≤300
解題思路
考慮暴力的dpdpdp,設fi,j,kf_{i,j,k}fi,j,k?表示對于前kkk個人種類為jjj的衣服選擇了iii件。
這樣顯然過不了。
但是考慮答案,假設我們第iii種衣服選擇了kkk件那么產生的貢獻就是
∑j=0ki×fi,j,n+k∑j=k+1nfi,j,n\sum_{j=0}^k i\times f_{i,j,n}+k\sum_{j=k+1}^nf_{i,j,n}j=0∑k?i×fi,j,n?+kj=k+1∑n?fi,j,n?
然后對于k?>k+1k->k+1k?>k+1會多產生的貢獻就是1?∑j=1kfi,j,n1-\sum_{j=1}^kf_{i,j,n}1?∑j=1k?fi,j,n?。考慮到這個值肯定是單調遞減的,所以貢獻函數是一個關于kkk的上凸函數。
然后就是很經典的方法了,每次暴力選擇一個能擴展的最大的擴展即可。
時間復雜度O(n(n+m))O(n(n+m))O(n(n+m))
code
#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int M=310,N=3100; int n,m,k[M];double s[M],f[2][M][N],a[M][N],ans; int main() {scanf("%d%d",&n,&m);for(int i=1;i<=m;i++)f[0][i][0]=1;for(int i=1;i<=n;i++)for(int j=1;j<=m;j++){scanf("%lf",&a[j][i]);a[j][i]/=1000.0;f[0][j][i]=f[0][j][i-1]*(1-a[j][i]);}for(int i=1;i<=m;i++){for(int j=1;j<=n;j++)f[1][i][j]=f[1][i][j-1]*(1-a[i][j])+f[0][i][j-1]*a[i][j];k[i]=1;s[i]=f[0][i][n];}for(int p=1;p<=n;p++){int pos=1;for(int i=2;i<=m;i++)if(s[i]<s[pos])pos=i;ans=ans+(1-s[pos]);s[pos]=s[pos]+f[k[pos]][pos][n];k[pos]^=1;int o=k[pos];for(int i=0;i<=n;i++)f[o][pos][i]=0;for(int i=1;i<=n;i++)f[o][pos][i]=f[o][pos][i-1]*(1-a[pos][i])+f[!o][pos][i-1]*a[pos][i];}printf("%.12lf\n",ans);return 0; }總結
以上是生活随笔為你收集整理的CF183D-T-shirtx【dp,贪心】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 立个flag是啥意思 立flag是个什么
- 下一篇: CF235D-Graph Game【LC