国际惯例不放题干

扯淡

其实题目翻译过来是杀戮尖塔,某steam上的卡牌游戏,我也曾热衷刷榜

题解

首先题目中要求的期望是假期望,结合题目中所给的阶乘就可以看出这其实是从$2*n$张牌中选择$m$张牌使用,并且所有情况都取最大值时的和

首先排序贪心最大

再说一个非常显然的结论,有强化牌就用,强化牌上数大于2,所以打出强化牌一定收益>=比打出一张攻击牌带来(保证最后一定要出至少一张攻击牌就行)

例如  你有强化牌  2 2 攻击牌 4 4 限制出2张

那么你先打一张2再打一张4和 先打4 再打4 收益一样

那么如果强化牌比k张要小,那么要把所有强化牌都用了,最后打输出.如果强化牌比k还要大,那么最后留一张攻击牌打出来就行了

我们可以假设$F(i,j)$表示你有i张强化牌时打出j张时收益,$G(i,j)$表示你有i张攻击牌时打出j张时收益

那么最后答案就是

$\sum\limits_{i=0}^{k-1} F(i,i)*G(m-i,k-i)+\sum\limits_{i=k}^{k=min(n,m)} F(i,k-1)*G(m-i,1)$

那么我们现在的任务就是快速求出来F和G

考虑直接求非常难求,我们这时一般要开几个别的数组

首先开$f[i][j]$ $g[i][j]$表示当你抽到i张牌必打第i张一共打了j张时收益

那么可以推出来$f[i][j]=w_i\times \sum f[l][j-1]$和$g[i][j]=C_{i-1}^{j-1}w_i + \sum g[l][j-1]$ (这里的组合数表示一共有$C_{i-1}^{j-1}$种方式转移过来,每种转移都要加一个w)

然后得到$F(i,j)=\sum_{l=j}^{n} C_{n-l}^{i-j} f[l][j]$

解释一下 F可以从l中使用j张(这j张都是大牌)再从剩余n-l中选择i-j张不打共同构成F的i张打了j张

G同理得到$G(i,j)=\sum{l=j}^{n} C_{n-l}^{i-j} g[i][j]$

然后转移就完了

本题稍卡常

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define A 1510
#define py printf("f**k\n")
#define mod 998244353
ll n,m,t,k,ans=0,f[1510][1510],g[1510][1510],C[3010][3010],sum1[A],sum2[A],wg[A],wq[A];
inline ll read()
{ll x=0; char c=getchar();while(isdigit(c)){x=x*10+c-'0';c=getchar();}return x;
}
ll gcd(ll x,ll y){if(y==0) return x;return gcd(y,x%y);
}
inline ll F(ll i,ll j){ll r=0;for(ll l=j;l<=n;l++){r=(r+C[n-l][i-j]*f[l][j]%mod)%mod;//表示含义为首先从前l张选出j张再从n-l中选出l-j张
    }return r;
}
inline bool cmp(const ll & a,const ll & b){return a>b;
}
inline ll G(ll i,ll j){ll r=0;for(ll l=j;l<=n;l++){r=(r+C[n-l][i-j]*g[l][j]%mod)%mod;//含义同上
    }return r;
}
inline void re(){memset(sum1,0,sizeof(sum1));memset(sum2,0,sizeof(sum2));
}
inline void init(){n=read(),m=read(),k=read();f[0][0]=1;for(ll i=1;i<=n;i++)wq[i]=read();sort(wq+1,wq+n+1,cmp);for(ll i=1;i<=n;i++)wg[i]=read();sort(wg+1,wg+n+1,cmp);for(ll i=1;i<=n;i++)f[i][1]=wq[i],g[i][1]=wg[i];for(ll i=1;i<=n;i++){for(ll j=2;j<=i;j++){sum1[j-1]=(sum1[j-1]+f[i-1][j-1])%mod;f[i][j]=wq[i]*sum1[j-1]%mod;sum2[j-1]=(sum2[j-1]+g[i-1][j-1])%mod;g[i][j]=(C[i-1][j-1]*wg[i]+(sum2[j-1]))%mod;}}
/*    for(ll i=1;i<=n;i++)for(ll j=1;j<=i;j++)printf("f[%lld][%lld]=%lld g[%lld][%lld]=%lld\n",i,j,f[i][j],i,j,g[i][j]);
*/}
inline void work(){ll ans=0;for(ll i=min(n,m);i>=0;i--){if(i<k)ans=(ans+F(i,i)*G(m-i,k-i)%mod)%mod/*,printf("F=%lld G=%lld\n",F(i,i),G(m-i,k-i))*/;else ans=(ans+F(i,k-1)*G(m-i,1)%mod)%mod/*,printf("F=%lld G=%lld\n",F(i,k-1),G(m-i,1))*/;
//        printf("ans=%lld\n",ans);
    }printf("%lld\n",ans);
}
int main()
{for(ll i=0;i<=3000;i++)C[i][0]=1;for(ll i=1;i<=3000;i++)for(ll j=1;j<=i;j++)C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;t=read();while(t--){re();init();work();}
}

转载于:https://www.cnblogs.com/znsbc-13/p/11238425.html

「PKUWC2018」Slay the Spire相关推荐

  1. loj2538 「PKUWC2018」Slay the Spire 【dp】

    题目链接 loj2538 题解 比较明显的是,由于强化牌倍数大于\(1\),肯定是能用强化牌尽量用强化牌 如果强化牌大于等于\(k\),就留一个位给攻击牌 所以我们将两种牌分别排序,企图计算\(F(i ...

  2. LOJ#2542. 「PKUWC2018」随机游走

    LOJ#2542. 「PKUWC2018」随机游走 题目描述 Solution 去过一个点集中所有节点的期望时间不好求,考虑min−maxmin-maxmin−max容斥,转化为求第一次到达某一个点集 ...

  3. LOJ 2537 「PKUWC2018」Minimax

    BZOJ 5461. 线段树合并优化$dp$. 假设所有离散之后的权值$\in [1, m]$,对于一个点$x$它的权值是$i$的概率是$f(x, i)$,那么 1.假如这个点只有一个儿子$y$,那么 ...

  4. 「PKUWC2018」Minimax

    传送门 Solution 发现叶子节点的值都不样,所以可以线段树合并. 然后因为我们要维护一个后缀,所以我们先合并右儿子,在合并左儿子 Code //2019.1.14 8:59~10:15 Pape ...

  5. 「PKUWC2018」随机游走

    Min-Max容斥真好用 LOJ #2542 题意:给一棵不超过1818个节点的树,50005000次询问,每次问从根随机游走走遍一个集合的期望步数 Solution: 考虑Min-Max容斥 有Ma ...

  6. 【LOJ】 #2540. 「PKUWC2018」随机算法

    题解 感觉极其神奇的状压dp \(dp[i][S]\)表示答案为i,然后不可选的点集为S 我们每次往答案里加一个点,然后方案数是,设原来可以选的点数是y,新加入一个点后导致了除了新加的点之外x个点不能 ...

  7. loj2537 「PKUWC2018」Minimax 【概率 + 线段树合并】

    题目链接 loj2537 题解 观察题目的式子似乎没有什么意义,我们考虑计算出每一种权值的概率 先离散化一下权值 显然可以设一个\(dp\),设\(f[i][j]\)表示\(i\)节点权值为\(j\) ...

  8. loj#2542. 「PKUWC2018」随机游走(MinMax容斥 期望dp)

    题意 题目链接 Sol 考虑直接对询问的集合做MinMax容斥 设\(f[i][sta]\)表示从\(i\)到集合\(sta\)中任意一点的最小期望步数 按照树上高斯消元的套路,我们可以把转移写成\( ...

  9. loj#2541. 「PKUWC2018」猎人杀

    传送门 思路太清奇了-- 考虑容斥,即枚举至少有哪几个是在\(1\)号之后被杀的.设\(A=\sum_{i=1}^nw_i\),\(S\)为那几个在\(1\)号之后被杀的人的\(w\)之和.关于杀了人 ...

最新文章

  1. PHP支持模板,让ThinkPHP5支持模板主题(参照ThinkPHP3.2)
  2. Sending e-mail with Spring MVC--转载
  3. 符号引用(typeglob,别名)与全局变量的修改
  4. 数据库元数据数据字典查询_6_列出给定表的外键引用
  5. javaweb学习总结(十五)——JSP基础语法
  6. 面试题:聊一聊设计模式的基本原则
  7. Java基础学习总结(76)——Java异常深入学习研究
  8. 浙江万里学院计算机专业宿舍,2020年浙江万里学院宿舍条件环境照片 宿舍空调相关配置介绍...
  9. winform前后端框架_ABP开发框架前后端开发系列(1)框架的总体介绍
  10. Home Assistant系列 -- 设置界面语言与地理位置
  11. Git以及Git小乌龟的安装使用
  12. scp构造端口_scp端口号(scp默认端口号)
  13. Hadoop Failover
  14. c# wia摄像_使用WIA捕获网络摄像机的视频-不可能吗?
  15. 战争军事历史纪录片收集
  16. new joiner
  17. Java实现快递管理系统一(views)
  18. GPS是如何定位你在哪的?
  19. CREO工程图学习笔记
  20. 使用WTSGetActiveConsoleSessionId()的VISTA服务与桌面交互

热门文章

  1. Maven的资源坐标
  2. Linux 命令之 nslookup 命令-查询域名 DNS 信息的工具
  3. 舱机器人尾巴毛茸茸_毛茸茸的(Fluffy)——欧美有史以来最伟大的机器人战争机器人...
  4. 归并排序 java_归并排序(Java实现)
  5. C语言 | 一维数组
  6. 什么是C语言中的隐式函数声明?
  7. 里氏替换原则_代码需要有单一职责,还要开闭,里氏替换又是什么鬼?
  8. ospf工作原理_OSPF动态路由配置经典案例
  9. C/C++头文件与变量的声明和定义
  10. 什么是DVI光端机?dvi光端机的优势有哪些?