coin
### Decsription
数据范围:\(n<=3000,m<=300\),保证\(\forall i,\sum\limits_{j}p_{ij}=1000\)
Solution
日常期望算不对。。
首先比较明显的一点是,每种类型是独立的,我们可以分开来考虑
先来想一个比较直接的dp
用\(a[i][j]\)表示第\(i\)个人最喜欢\(j\)的概率,记\(f[c][i][j]\)表示前\(i\)个人里面恰好有\(j\)个人最喜爱类型\(c\),转移显然:
\]
再考虑用\(g[i][j]\)表示第\(i\)种钞票携带了\(j\)张,期望拿走的数量,那么有转移:
g[i][j]&=\sum\limits_{k=0}^mmin(k,j)*f[i][n][k]\\
&=\sum\limits_{k=0}^j k*f[i][n][k]+\sum\limits_{k=j+1}^mj*f[i][n][k]\\
\end{aligned}
\]
然后我们就可以用一个背包求出答案了,但是这样无论是空间还是时间都很爆炸
考虑\(\nabla g[i][j]=g[i][j]-g[i][j-1]\),也就是\(\nabla g[i][j]=\sum\limits_{k=j}^m f[i][n][k]\),首先注意到\(f[i][n][k]\)在\(i\)确定的情况下大小一定是随着\(k\)的增大而增大的,然后再看回这个\(\nabla g[i][j]\),显然这个东西的值应该是非负的,并且当\(i\)一定,\(\nabla g[i][j]\)的值随着\(j\)的增大而单调不升
所以我们可以得到一个结论,当\(i\)确定的时候,\(g[i][j]\)是不降的,并且增幅逐渐下降
然后又因为每种钞票是独立的,我们可以考虑一个贪心
枚举\(n\)张钞票的每一张都选什么,对于每一种钞票\(c\),我们维护一个隐性的\(p_c\)(说是“隐性”是因为在实现的时候你并不需要真的去维护这么一个东西),表示钞票\(c\)当前已经取了\(p_c\)张,那么从贪心的角度来看,我们当前枚举到的这张钞票,肯定希望新加入这张钞票之后,对应的\(g[c]\)的增幅最大,所以就可以得到一个大致的流程:从\(1\)到\(n\)枚举每一张钞票选什么,对于第\(i\)张钞票,\(O(m)\)求得最大的\(g\)的增幅,加入答案,然后对于选择的这个种类\(c\),更新其\(g\)值
现在的问题就是怎么维护增幅,为了方便下面只针对确定的\(c\)类钞票进行表述
一开始的时候\(p_c=0\),\(g[c][0]\)到\(g[c][1]\)的增幅是很好计算的,\(\nabla g[c][1]=1-\prod\limits_{i=1}^n(1-a[i][c])\),但是\(\nabla g[c][2]\)看起来就不能那么直接地进行计算了,所以我们还是看回这个式子:\(\nabla g[i][j]=\sum\limits_{k=j}^mf[i][n][k]=1-\sum\limits_{k=0}^{j-1}f[i][n][k]\),那所以我们只要对于第\(c\)种钞票维护\(f[c][i][j]\)就好,然后维护一个\(sum[c]=\sum\limits_{k=0}^{p_c}f[c][n][k]\),每次更新完\(f[c]\)之后把\(f[c][n][p_c]\)加到\(sum[c]\)里面就好了
接下来就是空间怎么解决:对于\(f[c][i][j]\)来说,我们是按照\(j\)进行dp的,显然\(j\)那维可以滚动掉,\(f[c][][j]\)到\(f[c][][j+1]\)的更新是\(O(n)\)的,然后我们就获得了一个\(O(nm)\)的算法
mark:(弱智操作)谁告诉你听算期望的时候可以钦定顺序的???
Code
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=3010,M=310;
double sum[N],f[2][M][N];
double a[N][M];
int now[M],pre[M];
int n,m;
double ans;
void dp(int x){
swap(now[x],pre[x]);
int Now=now[x],Pre=pre[x];
f[Now][x][0]=0;
for (int i=1;i<=n;++i)
f[Now][x][i]=f[Now][x][i-1]*(1-a[i][x])+f[Pre][x][i-1]*a[i][x];
}
int main(){
#ifndef ONLINE_JUDGE
freopen("a.in","r",stdin);
#endif
int x;
scanf("%d%d",&n,&m);
for (int i=1;i<=n;++i)
for (int j=1;j<=m;++j)
scanf("%d",&x),a[i][j]=1.0*x/1000;
for (int i=1;i<=m;++i){
f[0][i][0]=1;
for (int j=1;j<=n;++j)
f[0][i][j]=f[0][i][j-1]*(1-a[j][i]);
sum[i]=f[0][i][n];
}
for (int i=1;i<=m;++i) now[i]=0,pre[i]=1;
int which;
double mx;
for (int j=1;j<=n;++j){
mx=0; which=-1;
for (int i=1;i<=m;++i)
if (1-sum[i]>mx)
which=i,mx=1-sum[i];
ans+=mx;
dp(which);
sum[which]+=f[now[which]][which][n];
}
printf("%.10lf\n",ans);
}
coin的更多相关文章
- [LeetCode] Coin Change 硬币找零
You are given coins of different denominations and a total amount of money amount. Write a function ...
- 洛谷P2964 [USACO09NOV]硬币的游戏A Coin Game
题目描述 Farmer John's cows like to play coin games so FJ has invented with a new two-player coin game c ...
- [luogu2964][USACO09NOV][硬币的游戏A Coin Game] (博弈+动态规划)
题目描述 Farmer John's cows like to play coin games so FJ has invented with a new two-player coin game c ...
- LeetCode Coin Change
原题链接在这里:https://leetcode.com/problems/coin-change/ 题目: You are given coins of different denomination ...
- ACM Coin Test
Coin Test 时间限制:3000 ms | 内存限制:65535 KB 难度:1 描述 As is known to all,if you throw a coin up and let ...
- HDOJ 2069 Coin Change(母函数)
Coin Change Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- leetcode:Coin Change
You are given coins of different denominations and a total amount of money amount. Write a function ...
- UVa 674 Coin Change【记忆化搜索】
题意:给出1,5,10,25,50五种硬币,再给出n,问有多少种不同的方案能够凑齐n 自己写的时候写出来方案数老是更少(用的一维的) 后来搜题解发现,要用二维的来写 http://blog.csdn. ...
- Epic - Coin Change
Something cost $10.25 and the customer pays with a $20 bill, the program will print out the most eff ...
- UVA 674 Coin Change (DP)
Suppose there are 5 types of coins: 50-cent, 25-cent, 10-cent, 5-cent, and 1-cent. We want to make c ...
随机推荐
- 20155216 Exp8 WEB基础实践
Exp8 WEB基础实践 实践内容 Apache环境配置 apache是kali下的web服务器,通过访问 ip地址+端口号+文件名称 打开对应的网页. 输入命令 vi /etc/apache2/po ...
- 11、可扩展MySQL+12、高可用
11.1.扩展MySQL 静态分片:根据key取hash,然后取模: 动态分片:用一个表来维护key与分片id的关系: 11.2.负载均衡 12. 12.2导致宕机得原因: 35%环境+35%性能+2 ...
- POJ 1068&&2632&&1573&&2993&&2996
这次的题目是著名的模拟(水逼)专题 题目难度都不大,思维也不深,就是一个字——烦 同时很多题目都有很多坑点 1608 题意是告诉你一个只有()的字符串(且匹配正确),每个)的左边有几个(,让你求每一对 ...
- Luogu P1198 [JSOI2008]最大数
我会用高级(???)的单调栈来打这道题吗? 线段树即可水过. 假设这个数列刚开始所有数都是0,然后我们每次只要进行一个点的修改和区间求和即可. 这不就是 线段树大法. 只要用一个len记录一下当前数列 ...
- LeetCode 3Sum (Two pointers)
题意 Given an array S of n integers, are there elements a, b, c in S such that a + b + c = 0? Find all ...
- POJ1094——拓扑排序和它的唯一性
比较模板的topological-sort题,关键在于每个元素都严格存在唯一的大小关系,而一般的拓扑排序只给出一个可能解,这就需要每趟排序的过程中监视它是不是总坚持一条唯一的路径. 算法导论里面的拓扑 ...
- 通过C#的HttpClient模拟form表单请求
post提交表单一般无非是一般text文本和文件类型,如下 <input type="file"/> <input type="text"/& ...
- Asp.Net_后台代码访问前台html标签
//单击按钮后批量改变li元素的内联文本值及样式 ; i <= ; i++) { HtmlGenericControl li = this.FindControl("li" ...
- SCRUM 12.14
由于最近的课业较多,大家平时有些力不从心,对于工作完成得有限. 最近课业压力小了一些,我们决定从今天起,投入精力. 以下为我们的任务分配情况: 人员 任务 高雅智 清除缓存 彭林江 统计活跃用户数量 ...
- M1阶段事后分析
M1阶段的开发结束了,在周四的课上我们组也进行了alpha阶段的汇报.我们的努力得到了应有的回报,下面我们将针对M1阶段产生的一些问题进行分析和反思. 一.设想和目标 1.我们的app更像是一款针对北 ...