最新要闻
- 三星Galaxy S24系列或取消Plus版本:销量太惨淡
- 【天天时快讯】大碗更尽兴!海福盛香辣牛肉面大促:每桶到手3块钱
- 今日视点:特斯拉2022年中国销量44万辆 还不敌比亚迪一个宋
- 《咬文嚼字》公布年度十大语文差错:连花清瘟?莲花清瘟?
- 世界观焦点:荣耀二代骁龙8新机来了!Magic 5系列入网:春节后发布
- 热门:别克GL8危险!腾势D9累计订单超5万:50%用户来自BBA
- 3岁男童反复呕吐被确诊癌症晚期:被称为儿童癌症之王
- 当前速读:《咬文嚼字》公布2022年度十大语文差错:天和核心舱、莘莘学子上榜
- 努比亚Z50限定版明天首销:搭载最纯净的定制系统 无广告
- 今日热议:2022年我国人均存款近1.3万元 网友:又拖后腿了
- 【报资讯】一家四口在三亚溺水全部遇难:官方科普“离暗流”危险性
- 天天微动态丨2022收官!合资车时代被终结 “迪王”养成 大票车企失去肥年
- 【环球速看料】《小美人鱼》真人电影周边童书曝光 黑美人鱼好可爱
- 天天即时:双形态不入耳!讯飞开放式办公耳机iFLYBUDS Air图赏
- 小米MIX Fold 2厚度与戴壳iPhone 14 Pro Max相当 雷军:惊艳
- 当前通讯!3秒复制任何人的嗓音!微软音频版DALL·E细思极恐 连环境背景音也能模仿
手机
iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
- 警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- 男子被关545天申国赔:获赔18万多 驳回精神抚慰金
- 3天内26名本土感染者,辽宁确诊人数已超安徽
- 广西柳州一男子因纠纷杀害三人后自首
- 洱海坠机4名机组人员被批准为烈士 数千干部群众悼念
家电
Codeforces 1278 F Cards 增强版 题解 (斯特林数,推式子)
原题链接
增强版链接
(相关资料图)
增强版中k=1e7
为啥网上题解的式子都那么长啊.jpg
首先令\(p=\frac 1m\)。求某个数的次幂的题通常都是无脑转下降幂:\(x^k=\sum_{i=0}^k S_2(k,i)x^{\underline i}\),其中\(S_2\)表示第二类斯特林数,\(x^{\underline i}\)表示下降幂,也就是\(\binom xi i!\)(i>x时值为0)。对于一种实际赢了\(x\)场的情况,\(S_2(k,j)\)对它的答案贡献应为\(\binom xjj!\)。因此我们可以把这个组合数中的每一种选取情况拿出来分别计算贡献,所以最终答案\(=\sum_{j=0}^k S_2(k,j)\binom njj!p^j\),这一步转化也可以列式子推,但是能用组合意义还是用吧。第二类斯特林数是可以\(O(klogk)\)求出一行的,如果NTT板子牛逼的话是可以过这题的,但有\(O(k)\)的做法。
有一个关于第二类斯特林数的公式:\(S_2(k,j)=\frac 1{j!}\sum_{i=0}^j (-1)^{j-i}\binom ji i^k\)。直接带入上面的式子得到:\(ans=\sum_{j=0}^k\sum_{i=0}^j (-1)^{j-i} \frac{n!}{i!(j-i)!(n-j)!}p^ji^k\)。这相当于是我们要把n个元素组成的集合分割成3份,大小分别为\(i,j-i,n-j\)(都可以为0,且前两部分大小之和\(\le k\)),其中第一部分的"权值"为\(i^kp^i\),第二部分的权值为\((-p)^{j-i}\),第三部分的权值为1,求所有分割方式的权值之积的和。我们枚举i,尝试把剩下的权值\(O(1)\)求出。改写一下答案:\(ans=\sum_{i=0}^k i^kp^i f(i)\),其中\(f(i)=\sum_{j=0}^{k-i}\binom{n-i}j(-p)^j\ \ (注意这里的j不是上面的j了)\),我们想要\(O(k)\)求出\(f_0\cdots f_k\)。
其实f是可以差分之后递推求的。注意到\(f_k=1\),所以我们反向差分:
\[\begin{align}f_i-f_{i+1}&=(\sum_{j=0}^{k-i}\binom{n-i}j(-p)^j)-(\sum_{j=0}^{k-i-1}\binom{n-i-1}j(-p)^j)\\&=\binom{n-i}{k-i}(-p)^{k-i}+\sum_{j=0}^{k-i-1}[\binom{n-i}j-\binom{n-i-1}j](-p)^j\\&=\binom{n-i}{k-i}(-p)^{k-i}+\sum_{j=1}^{k-i-1}\binom{n-i-1}{j-1}(-p)^j\ \ 考虑组合数的递推公式\\\end{align}\]再看看这个式子的后半部分等于什么:
\[\begin{align}&\sum_{j=1}^{k-i-1}\binom{n-i-1}{j-1}(-p)^j\\=&(-p)\sum_{j=1}^{k-i-1}\binom{n-i-1}{j-1}(-p)^{j-1}\\=&(-p)\sum_{j=0}^{k-i-2}\binom{n-i-1}j(-p)^j\\=&(-p)(f_{i+1}-\binom{n-i-1}{k-i-1}(-p)^{k-i-1})\end{align}\]因此只要预处理一下组合数就能\(O(k)\)求f了,也就是预处理\(g_i=\binom{n-i}{k-i}\),这个很容易\(O(k)\)求。
在\(ans=\sum_{i=0}^k i^kp^i f(i)\)中,求出了f还需要\(O(k)\)对所有i求\(i^kp^i\),这个东西是积性的,所以可以线性筛。但是快速幂常数不大所以应该也是可以过的。
时间复杂度\(O(k)\)。下面代码里用了快速幂,是\(O(klogk)\),没有特意去卡洛谷上的时间和空间。
点击查看代码
#include #define rep(i,n) for(int i=0;i#define fi first#define se second#define mpr make_pair#define pb push_backvoid fileio(){ #ifdef LGS freopen("in.txt","r",stdin); freopen("out.txt","w",stdout); #endif}void termin(){ #ifdef LGS std::cout<<"\n\nEXECUTION TERMINATED"; #endif exit(0);}using namespace std;const LL MOD=998244353;LL qpow(LL x,LL a){LL res=x,ret=1;while(a>0){if(a&1) (ret*=res)%=MOD;a>>=1;(res*=res)%=MOD;}return ret;}LL n,m,k,f[10000010],p,fac[10000010],inv[10000010],cval[10000010],mpp[10000010];int main(){ fileio(); freopen("sanhok.in","r",stdin); freopen("sanhok.out","w",stdout); fac[0]=1;repn(i,10000005) fac[i]=fac[i-1]*i%MOD; inv[10000003]=qpow(fac[10000003],MOD-2); for(int i=10000002;i>=0;--i) inv[i]=inv[i+1]*(i+1)%MOD; cin>>n>>m>>k; p=qpow(m,MOD-2); LL bas=(MOD-p)%MOD; mpp[0]=1;repn(i,10000005) mpp[i]=mpp[i-1]*bas%MOD; cval[k]=1;//cval[i]=C(n-i,k-i) LL cmul=1; for(int i=k-1;i>=0;--i) { (cmul*=(n-i))%=MOD; cval[i]=cmul*inv[k-i]%MOD; } f[k]=1; for(int i=k-1;i>=0;--i) { f[i]=(f[i+1]+cval[i]*mpp[k-i])%MOD; LL add=(f[i+1]-mpp[k-i-1]*cval[i+1]%MOD+MOD)%MOD; (add*=(MOD-p))%=MOD; (f[i]+=add)%=MOD; } LL pi=1,mul=1,ans=0; rep(i,k+1) { LL val=pi*qpow(i,k)%MOD*mul%MOD*inv[i]%MOD; (pi*=p)%=MOD;(mul*=(n-i))%=MOD; (ans+=val*f[i])%=MOD; } cout<
-
Codeforces 1278 F Cards 增强版 题解 (斯特林数,推式子)
原题链接增强版链接增强版中k=1e7为啥网上题解的式子都那么长啊 jpg首先令$p= frac1m$。求某个数的...
来源: -
热点在线丨sortablejs 列表拖拽排序,js vue2,解决拖拽排序乱序问题
功能:在列表中,需要给列表进行拖拽排序,并实时保存拖拽后的列表书序实现;运用js中的sortablejs库环境...
来源: Codeforces 1278 F Cards 增强版 题解 (斯特林数,推式子)
热点在线丨sortablejs 列表拖拽排序,js vue2,解决拖拽排序乱序问题
世界快资讯:用低代码这把“剑”之前,要先看定位,各取所需
当前短讯!开源动物行为分析实验箱(斯金纳箱)研发总结
三星Galaxy S24系列或取消Plus版本:销量太惨淡
【天天时快讯】大碗更尽兴!海福盛香辣牛肉面大促:每桶到手3块钱
今日视点:特斯拉2022年中国销量44万辆 还不敌比亚迪一个宋
《咬文嚼字》公布年度十大语文差错:连花清瘟?莲花清瘟?
世界观焦点:荣耀二代骁龙8新机来了!Magic 5系列入网:春节后发布
环球观速讯丨python之路 58 linux文件配置相关
学习笔记——MyBatis自动映射与自定义映射;Mybatis延迟加载
当前关注:C#、TS和Dart对比1:概述
即时焦点:SpringBoot Xss漏洞修复
软件开发入门教程网之Git 基本操作
热门:别克GL8危险!腾势D9累计订单超5万:50%用户来自BBA
3岁男童反复呕吐被确诊癌症晚期:被称为儿童癌症之王
当前速读:《咬文嚼字》公布2022年度十大语文差错:天和核心舱、莘莘学子上榜
努比亚Z50限定版明天首销:搭载最纯净的定制系统 无广告
环球观点:操作系统
环球播报:C++构造函数【cherno课程学习】
手机端H5 实现自定义拍照界面
今日热议:2022年我国人均存款近1.3万元 网友:又拖后腿了
【报资讯】一家四口在三亚溺水全部遇难:官方科普“离暗流”危险性
天天微动态丨2022收官!合资车时代被终结 “迪王”养成 大票车企失去肥年
【环球速看料】《小美人鱼》真人电影周边童书曝光 黑美人鱼好可爱
天天即时:双形态不入耳!讯飞开放式办公耳机iFLYBUDS Air图赏
小米MIX Fold 2厚度与戴壳iPhone 14 Pro Max相当 雷军:惊艳
当前通讯!3秒复制任何人的嗓音!微软音频版DALL·E细思极恐 连环境背景音也能模仿
每日热议!官宣:Intel发烧U回来了!350W 56核能打过AMD 280W 64核吗?
每日播报!人气爆棚!上美回应《中国奇谭》周边断货:已开足马力生产
Shell 命令奇淫技巧,就是有点短
当前头条:你买过哪些?苹果已售出23.2亿部手机 国人最爱iPhone 6
每日速递:美国加州风暴天气已致17死 有大树直接被连根拔起
环球速看:取消灵动岛!苹果iPhone 16 Pro将配备屏下Face ID
春运咋办?博主跑1千公里高速实测充电桩:有服务区一半都是坏的
【时快讯】《中国奇谭》口碑封神!仅上线三集 播放量突破5000万
【世界聚看点】CPU、显卡持续涨价!全球PC出货量暴跌 联想继续第一
简讯:(五)elasticsearch 源码之查询流程分析
环球热文:消息服务 + Serverless 函数计算如何助力企业降本提效?
世界滚动:el-table更新数据页面闪烁问题
全球今日报丨DJI这三个字母 是怎么占领你的背包的
环球热头条丨告别毛巾“一条恒久远”!金号纯棉抑菌毛巾大促:一条5块钱
今日快讯:哪吒汽车联手宁德时代共研“滑板底盘”:电池、底盘合体
世界信息:苹果加大降低中国工厂依赖程度:都要搬走?印度成香饽饽 出口激增
独占4K AMR 120帧高规格!《流浪地球2》发布CINITY海报
世界热点!指针知识点总结
每日热讯!TiDB 底层存储结构 LSM 树原理介绍
linux基础:2、前期必备知识、系统运行命令、快捷方式命令、目录结构相关命令、文件与文件夹相关命令、目录结构
环球今日报丨C# 循环给多个连续编号的控件赋值
网上银行怎么转账?网上银行转账限额是多少?
诺基亚5800xm当年多少钱?诺基亚5800XM手机参数
投影仪吊架怎么安装?吊式投影仪安装方法
华为gt2怎么设置相册表盘?华为gt2有血氧功能吗?
唐门鸟翔碧空在哪里学?唐门鸟翔碧空可以放什么技能?
雷龙鱼水温多少合适?雷龙吃什么饲料?
最新消息:三星Galaxy S23系列定档:2月2日登场 首发新版骁龙8 Gen2
《满江红》公布秦桧版预告:饰演者雷佳音狠辣狡诈
【全球新要闻】特斯拉大降价 其它车企跟不跟?乘联会秘书长发声
焦点速读:万物有灵 被收养流浪狗跳车拦住怀孕主人 下一秒山路塌方
每日聚焦:1208元!中国探月航天推出限量火箭碎片:运送嫦娥四号的长三乙
加减乘除是谁发明的?加减乘除混合运算100道
米亲韩语是什么意思?韩语shake it是什么意思?
全高清和超高清有什么区别?全高清和超高清4K哪个更护眼?
异丙醇的作用与用途有哪些?异丙醇和酒精的区别是什么?
Serverless 奇点已来,下一个十年将驶向何方?
每日热点:没电、没网也能支付 数字人民币全新功能上线:安卓先行
环球微资讯!用上比亚迪发动机 斯威大虎ED-i增程版亮相:油耗低至2.06升
关注:公司就给员工加薪50元致歉 已尽力对不起大家引热议:为何不知足?
Win7彻底停服 国产OS统信站出来了:打印机、软件轻松迁移
读编程与类型系统笔记04_类型安全
微动态丨【QtJson】用Qt自带的QJson,直接一步到位封装和解析一个类的实例对象!
还买iPhone 14/15?新iPhone准备中:苹果弃灵动岛 更完美
人民日报评电视收费乱象:从用户身上“薅羊毛” 广告还见缝插针
微头条丨今年首场寒潮来袭 最强雨雪下在哪?这些地区将迎暴雨、暴雪
天天滚动:2022年 我把比亚迪DM-i插混当纯电动车 开了1万公里
《王者荣耀》星会员正式上线:全新充值体系 打破贵10等级限制
当前速讯:神舟新款游戏本上架:13代i5+满血RTX 3050 4999元
最强高端!小米13 Ultra有望下月发:终于告别USB 3.0
动态:家长炮轰《中国奇谭》烂 吓哭孩子引热议:专家回应 给成年人看
环球观焦点:贾跃亭又摊上事了:法拉第未来接到纳斯达克摘牌警告
全球报道:JavaScript 扁平与树形数组数据的转换
热点在线丨层叠样式表(CSS)1
当前最新:ceph-3
世界速读:区块链特辑——solidity语言基础(四)
热文:华硕全球首秀四频段Wi-Fi 7路由器:峰值下载2.5万兆
滚动:惠普新款EliteBook 1040笔记本发布:13代酷睿、2K 120Hz屏
天天视点!为什么一个病毒株传着传着就没了?
Redmi K40S 12+256G顶配版不到1900元:骁龙870和OIS都有
各大新能源车企年度KPI出炉 特斯拉未达标 比亚迪称王
焦点速读:区块链特辑——solidity语言基础(二)
每日时讯!酷安最火骁龙8系手机诞生!一加11酷安热度第一
时隔两年 国美真快乐App重新更名国美
快看点丨大前端html学习06-宽高自适应
天天观察:SpringDataJPA 程序未配置乐观锁的情况下,报了乐观锁异常
Redis的客户端
每日时讯!关于19c RU补丁报错问题的分析处理
环球简讯:记录使用adb连接rn项目进行开发
焦点要闻:加油被惊喜到!一加11不杀后台:前一天打开的APP第二天还在
最新快讯!Python中高阶函数与装饰器教程
【新要闻】Python中的异常处理总结