cdqz2017-test10-柚的策略(期望DP & 组合数学)
根据期望的可加性,我们可以算出每一位客人的期望等待时间,将他们累加
即 每一位客人所有可能情况的时间之和 / n!
设S= 每一位客人所有可能情况的时间之和
如果有f(i,p)种方案使客人i是恰好第p个进入花亭的,那对S的贡献为(n-p+1)* t[i] * f(i,p)
所以问题转变为计算f(i,p),客人i是恰好第p个进入花亭的方案数
这个恰好很难算
所以转化为在前p个进入花亭的客人中有i,最后f(i,p)减f(i,p-1)就得到了恰好是第p个
若前k+p-1个客人中至少有k-1个客人的用时比第i个客人的用时多,那么客人i可以在前p个进入花亭
所以问题又转化为了计算dp(i,j,l),在i个客人里至少有j个客人的用时比第l个客人用时多,且第i个客人一定在这i个客人里的方案数
换个状态定义会更好算:
用时相同的客人谁先进谁后进对答案没有影响
将客人的用时映射到1——n
dp(i,j,l),在1——n里选i个数至少有j个数比l大 且 不能选l的方案数
可以得到方程:(把f换成dp)
有了dp(i,j,l),再来算f(i,p)
前面说了若前k+p-1个客人中至少有k-1个客人的用时比第i个客人的用时多,那么客人i可以在前p个进入花亭
假设t[i]映射到了c
所以f(i,p)= dp[k+p-1-1][k-1][c]*(k+p-1)*(n-(k+p-1))!
因为客人i要占据一个位置,这个位置有(k+p-1)种选择,
在客人i之后的客人可以随意组合,所以是全排列
#include<cstdio>
#include<iostream>
#include<algorithm> using namespace std; #define N 301
#define M 1000001 const int mod=1e9+; int t[N];
int v[M],rk[N]; int C[N][N];
int fac[N]; int dp[N][N][N];
int f[N][N]; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} int Pow(int a,int b)
{
int res=;
for(;b;a=1LL*a*a%mod,b>>=)
if(b&) res=1LL*a*res%mod;
return res;
} int main()
{
freopen("strategy.in","r",stdin);
freopen("strategy.out","w",stdout);
int n,k;
read(n); read(k);
for(int i=;i<=n;++i) read(t[i]),v[t[i]]++;
for(int i=;i<M;++i) v[i]+=v[i-];
for(int i=;i<=n;++i) rk[i]=v[t[i]]--;
C[][]=;
fac[]=;
for(int i=;i<=n;++i)
{
C[i][]=;
fac[i]=1LL*fac[i-]*i%mod;
for(int j=;j<=i;++j) C[i][j]=(C[i-][j-]+C[i-][j])%mod;
}
for(int i=;i<=n;++i)
{
for(int j=;j<=i;++j)
for(int l=;l<=n;++l)
dp[i][j][l]=1LL*C[n-l][j]*C[l-][i-j]%mod*fac[i]%mod;
//选i个数恰好有j个数比l大的方案数
for(int j=i;j>=;--j)
for(int l=;l<=n;++l)
{
dp[i][j][l]+=dp[i][j+][l];
dp[i][j][l]-=dp[i][j][l]>=mod ? mod : ;
}
//选i个数至少有j个数比l大的方案数
}
int c,m,h;
int ans=;
for(int i=;i<=n;++i)
{
c=rk[i];
for(int j=;j<=n;++j) // 第i个人是前j个进入花亭的
{
m=min(k+j-,n);
h=k--max(,k+j--n);
f[i][j]=1LL*dp[m-][h][c]*m%mod*fac[n-m]%mod;
}
for(int j=n;j;--j)
{
f[i][j]-=f[i][j-];
if(f[i][j]<) f[i][j]+=mod;
ans=(ans+1LL*f[i][j]*(n-j+)%mod*t[i]%mod)%mod;
// printf("%d\n",ans);
}
}
ans=1LL*ans*Pow(fac[n],mod-)%mod;
printf("%d",ans);
}
cdqz2017-test10-柚的策略(期望DP & 组合数学)的更多相关文章
- [CSP-S模拟测试]:chemistry(期望DP+组合数学)
题目传送门(内部题27) 输入格式 第一行有$4$个整数$n,k,p,q$.第二行有$n$个整数$a_i$.接下来有$n-1$行,每行有两个整数$u,v$,表示$u$与$v$之间通过化学单键连接. 输 ...
- 概率期望dp
对于概率dp,我一直都弄得不是特别明白,虽然以前也有为了考试去突击过,但是终究还是掌握得不是很好,所以决定再去学习一遍,把重要的东西记录下来. 1.hdu4405 Description 在一个 \( ...
- 2018.08.30 花园(期望dp)
题目背景 SCOI2017 DAY2 T1 题目描述 小 A 的花园的长和宽分别是 L,H .小 A 喜欢在花园里做游戏.每次做游戏的时候,他都先把花园均匀分割成 L×H 个小方块,每个方块的长和宽都 ...
- 【bzoj5197】[CERC2017]Gambling Guide 期望dp+堆优化Dijkstra
题目描述 给定一张n个点,m条双向边的无向图. 你要从1号点走到n号点.当你位于x点时,你需要花1元钱,等概率随机地买到与x相邻的一个点的票,只有通过票才能走到其它点. 每当完成一次交易时,你可以选择 ...
- [BZOJ4872][六省联考2017]分手是祝愿(期望DP)
4872: [Shoi2017]分手是祝愿 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 516 Solved: 342[Submit][Statu ...
- 【BZOJ4872】[Shoi2017]分手是祝愿 数学+期望DP
[BZOJ4872][Shoi2017]分手是祝愿 Description Zeit und Raum trennen dich und mich. 时空将你我分开.B 君在玩一个游戏,这个游戏由 n ...
- 【bzoj4872】[Shoi2017]分手是祝愿 数论+期望dp
题目描述 Zeit und Raum trennen dich und mich. 时空将你我分开. B 君在玩一个游戏,这个游戏由 n 个灯和 n 个开关组成,给定这 n 个灯的初始状态,下标为从 ...
- P3750 [六省联考2017]分手是祝愿 期望DP
\(\color{#0066ff}{ 题目描述 }\) Zeit und Raum trennen dich und mich. 时空将你我分开. B 君在玩一个游戏,这个游戏由 \(n\) 个灯和 ...
- BZOJ4872 [六省联考2017]分手是祝愿 【期望dp】
题目 Zeit und Raum trennen dich und mich. 时空将你我分开.B 君在玩一个游戏,这个游戏由 n 个灯和 n 个开关组成,给定这 n 个灯的初始状态,下标为 从 1 ...
随机推荐
- Spring Framework: @RestController vs @Controller
https://docs.spring.io/spring-framework/docs/current/javadoc-api/org/springframework/web/bind/annota ...
- Linux Network Commands
https://www.tecmint.com/linux-network-configuration-and-troubleshooting-commands/ http://www.tldp.or ...
- Which path should be used jdk or jre for JAVA_HOME environment variable?
https://stackoverflow.com/questions/17601827/which-one-should-java-home-to-point-jdk-or-jre 临时变更JAVA ...
- linux命令学习head和tail
linux命令head和tail是一对:more和less是一对. head和tail https://www.2cto.com/os/201507/414753.html 一个头,一个尾. tail ...
- FreeMarker boolean Issue
FreeMarker template error:Can't convert boolean to string automatically, because the "boolean_f ...
- [转帖] 数据库用优化方案 https://segmentfault.com/a/1190000006158186
Mysql大表优化方案 当MySQL单表记录数过大时,增删改查性能都会急剧下降,可以参考以下步骤来优化: 单表优化 除非单表数据未来会一直不断上涨,否则不要一开始就考虑拆分,拆分会带来逻辑.部 ...
- Test Scenarios for Filter Criteria
1 User should be able to filter results using all parameters on the page2 refine search functionalit ...
- luogu P2644 树上游戏
一道点分难题 首先很自然的想法就是每种颜色的贡献可以分开计算,然后如果你会虚树就可以直接做了 点分也差不多,考虑每个分治重心的子树对它的贡献以及它对它子树的贡献 首先,处理一个\(cnt\)数组,\( ...
- 【ZJOI 2019】麻将(dp of dp)
这是我第一次写$dp \; of \; dp$,大致思路参考了xyx的做法,可能有些地方不太一样,但也许会详细一点. 考虑给你一副牌,如何判断这副牌是否是胡的. 容易发现相同的顺子不会选三个以上,于是 ...
- ZJOI 2019 一试记
没有心思做太多的铺述. $Day1$考完已经过去了几天,逐渐变得冷静下来,细碎的思绪也已经飘远,剩下一些难以忘怀的东西.$WC$以后一段时间的训练,本以为把状态找了回来,豪情壮志地出发,结果竟又是落魄 ...