题意

xyq有\(n\)个骰子,第\(i\)个骰子有\(a_i\)面,每次xyq都会把\(n\)个骰子搞一遍,其中的最小值作为结果,问最终结果的期望\(\mod (10^9+7 )\)。

分析

lfx聚聚出的神仙题,瞎写了个\(dp\)解法

设\(dp[i][j]\)为前\(i\)个骰子的点数最小值为\(j\)的概率,有两种情况

  • \(dp[i][j]+=(\sum_{k=j+1}^{40}dp[i-1][k])\times \frac{1}{a[i]}\),(前\(i-1\)个骰子的最小值大于\(j\)时,第\(i\)个骰子的点数必须为\(j\))
  • \(dp[i][j]+=dp[i-1][j]\times \frac{a[i]-j+1}{a[i]}\),(前\(i-1\)个骰子的最小值等于\(j\)时,第\(i\)个骰子的点数的取值范围为\([~j,a[i]~]\) )

Code

#include<bits/stdc++.h>
#define fi first
#define se second
#define pb push_back
#define ll long long
using namespace std;
const ll inf=1e18;
const int mod=1e9+7;
const int maxn=1e5+10;
int n;
ll a[55],dp[55][55];
ll ksm(ll a,ll b){ll ret=1;while(b){if(b&1) ret=ret*a%mod;b>>=1;a=a*a%mod;}return ret;
}
int main(){//ios::sync_with_stdio(false);//freopen("in","r",stdin);cin>>n;for(int i=1;i<=n;i++){cin>>a[i];}sort(a+1,a+n+1);for(int i=1;i<=a[1];i++) dp[1][i]=ksm(a[1],mod-2);for(int i=2;i<=n;i++){for(int j=1;j<=a[i];j++){ll ret=0;for(int k=j+1;k<=40;k++) ret=(ret+dp[i-1][k])%mod;dp[i][j]=ret*ksm(a[i],mod-2)%mod+dp[i-1][j]*(a[i]-j+1)%mod*ksm(a[i],mod-2)%mod;dp[i][j]%=mod;}}ll ans=0;for(int i=1;i<=40;i++){ans=(ans+dp[n][i]*i%mod)%mod;}printf("%lld\n",ans);return 0;
}

转载于:https://www.cnblogs.com/xyq0220/p/11455537.html

HFUUOJ1023 闷声发大财 概率dp相关推荐

  1. ble 连接成功后找不到服务_闷声发大财的BLE芯片龙头

    来源:内容由半导体行业观察(ID:icbank)原创,作者:邱丽婷,谢谢. 你有多久没使用过蓝牙功能了?随着Wi-Fi和快传软件的普及,蓝牙(Bluethooth)这一每部手机都配备的数据传输工具,正 ...

  2. 有哪些闷声发大财的行业?

    闷声发大财的人,不会轻易告诉别人自己发财了,将自己发财的方法透露出来后,分一份羹给别人,这样的傻子确实有,但特别少. 有些人看似普普通通,没有高学历,人也不够聪明,但就是能够找到挣钱的方法.平时一块聊 ...

  3. 有什么毫不起眼却闷声发大财的行业或者生意?

    毫不起眼却闷声发大财的生意简直太多了,在2019年依然有很多暴利行业,可能投入的成本并不高,但是获得的利润却是非常可观的! 眼镜店 我舅妈就是开眼镜店的,所以对其盈利情况还是非常了解的.虽然近些年来眼 ...

  4. 大家都在积极发声的时候,这些互联网企业在闷声发大财!

    有人的地方就有江湖,就会有名利,只是有的企业乐于发声,例如BAT,而有的企业却隐于闹市,一直在默默的闷声发大财,本文就根据网友提供的种种线索,为大家将这些默默发财的公司展示出来. Viralnova: ...

  5. 闷声发大财的ARM凭什么这么牛

    http://tech.sina.com.cn/it/2017-03-30/doc-ifycwunr8116186.shtml?cre=techpagepc&mod=f&loc=1&a ...

  6. 思摩尔靠电子烟赚暴利:年利润达54亿 陈志平闷声发大财

    雷递网 雷建平 3月31日报道 闷声发大财的电子烟企业思摩尔国际控股有限公司(简称:"思摩尔国际",股票代码为:"6969")昨日发布财报.财报显示,思摩尔国际 ...

  7. 生意人告诉你闷声发大财的2个小技巧

    1.怎样找到闷声发大财的行业 虽然互联网信息足够发达,但很多重要信息并不是直接告诉你答案,而是通过各种现象表现出来.所以当你毫无头绪不知道做什么的时候,多去研究研究身边的小广告. 那些不知名的小广告大 ...

  8. 都说闷声发大财,这句话一点也没有错

    都说闷声发大财,这句话一点也没有错.不懂得低调.收敛,只会后悔不已. 堂弟以前是在菜市场卖蔬菜的,生意不咸不淡,比打工稍强点.这两年因为疫情,他改行蔬菜配送,反而生意好了不少. 去年应该是赚了一点钱吧 ...

  9. 年终奖8个月起,你肯定不知道蚂蚁金服有个闷声发大财的部门

    蚂蚁金服有个闷声发大财的神秘部门:月薪30K是超级白菜价,年终奖8个月起,上班时间普遍为1085,从去年年至今仍在扩招.同是进大厂享受福报,蚂蚁区块链部门的性价比真的非常高! 据知乎大V@半佛仙人透露 ...

最新文章

  1. bellman_ford寻找平均权值最小的回路
  2. C++ 虚函数个人理解
  3. Android Framework------之Property子系统
  4. 机电传动控制第三次作业
  5. MongoDB安装与副本集配置
  6. AIDL 发生异常的原因 Android java.lang.SecurityException: Binder invocation to an incorrect interface...
  7. 【bfs 反向建边】2016-2017 ACM-ICPC, Egyptian Collegiate Programming Contest (ECPC 16) Jumping
  8. OFF文件格式_拔剑-浆糊的传说_新浪博客
  9. Android O版本power按键锁屏亮屏流程
  10. QObject::moveToThread: Current thread(...) is not the object`s thread. Cannot move to target thread(
  11. 【C++005】结构体(值传递,地址传递,数组)
  12. linux下打开.mpp文件(微软project).
  13. Java进度条(excel文件解析)的实现
  14. Visual studio 无法下载安装
  15. org.jboss.deployment.DeploymentException: URL file:/X:/workspace/JBoss_AS/XXX.war/ deployment failed
  16. 前序遍历,中序遍历,后序遍历详解
  17. nmn吃第一天有什么感觉,吃完nmn的反应,一点点体会
  18. 吉林大学软件学院编译原理与实现习题(二) 期末复习用
  19. java咖啡_JAVA—咖啡馆
  20. 网络编程 TCP电子网络词库

热门文章

  1. 坚果O1S和当贝U1区别对比哪个好?3000元的档位超短焦投影怎么选?
  2. 计算机磁盘管理给优盘分区,给U盘分区的图文教程
  3. 评中级职称论文挂名可以快速通过职称评审吗?
  4. excel表格里的时间列杂乱无章,怎样快速统一格式?
  5. 得一微YS9082HP主控256G固态硬盘量产成功,简易教程,附YS9082HP开卡量产工具
  6. 好看的导航主页html,推荐|各种漂亮的网页导航,让浏览器靓起来
  7. Android Studio 安装、环境配置全家桶
  8. office2013零售版转vol版
  9. 苹果手机充电孔变松了_不久的将来,甚至下一代的苹果手机会取消充电接口吗?...
  10. 关于电流环带宽的一些理解