最新要闻
- 每日动态!2023春运车票24日开售 除夕车票要等到1月7日
- 《妮姬》首月收入突破6.9亿!腾讯海外收入占比提升达12.5%
- 今日热门!超越电竞机!Redmi K60要榨干第二代骁龙8:画质、帧率、亮度三不降
- 每日精选:效果堪比镀铬 2.2万元的特斯拉Model Y新配色值不值?
- Intel拆分GPU部门 一把手重回技术岗 累计亏损超20亿美金
- 天天通讯!iPhone 14 Pro爆出“闪线门” :屏幕出现诡异的绿色和黄色细横线
- 全球热点评!当ChatGPT遇上弱智吧:全程爆笑
- 夫妻的世界翻拍哪部电视剧?夫妻的世界最后结局是什么意思?
- 滕王阁为什么叫阁不叫楼?滕王阁为什么是三大名楼之首?
- 情非情砸车是第几集?情非情盖总和保姆的结局是什么?
- 小昭去波斯是哪一集?小昭去波斯后她母亲去哪儿了?
- 男人是大猪蹄子是什么意思?男人是大猪蹄子女人是什么?
- 每日消息!性能超越电竞手机!Redmi K60 Pro综合跑分达135万
- 信息:千万别强忍 20岁小伙憋气压抑咳嗽导致昏厥
- 特斯拉今年股价累计暴跌超60%!马斯克透露大跌原因
- 收购动视暴雪遇阻 微软哭弱:根本打不过索尼、任天堂
手机
iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- iphone11大小尺寸是多少?苹果iPhone11和iPhone13的区别是什么?
- 警方通报辅警执法直播中被撞飞:犯罪嫌疑人已投案
- 男子被关545天申国赔:获赔18万多 驳回精神抚慰金
- 3天内26名本土感染者,辽宁确诊人数已超安徽
- 广西柳州一男子因纠纷杀害三人后自首
- 洱海坠机4名机组人员被批准为烈士 数千干部群众悼念
家电
Codeforces 1630 E Expected Components 题解 (组合数学)
(相关资料图)
题目链接
首先明确概念:
- 排列。指的就是一个把数组a重排得到的序列,两个排列相等当且仅当它们对应位全都相等
- 环形排列。指的是把数组a重排得到的序列首尾相接得到的环形数组,两个环形排列相等当且仅当它们通过位移能够变得对应位全部相等
这题其实是要对所有不同的环形排列,求其中连续相同的数构成的连通块个数的期望值。特判a中数全相同的情况,其余情况其实只要求出相邻两个数不同的位置个数的期望值就行了,把这个期望值称为一个环形排列的"权值"。一个排列的权值指的是把它视为环形排列之后的权值(1,n也算相邻)。对于一个排列p,有多少个与它不同的排列在作为循环排列时与它相同?其实这与p的最短循环节长度有关,令这个长度为\(len\)。这里的循环节必须是完整的,也就是\(|p|\ mod \ len=0\),且把\(p\)平均切成\(\frac{|p|}{len}\)段之后每段都相同。观察发现对于p,有\(len-1\)个与它不同的排列在作为循环排列时与它相同,也就是等价类大小为\(len\)。如果能对每个\(len\)求出所有最小循环节大小为\(len\)的不同排列的权值和\(f_{len}\)与总个数\(cnt_{len}\),那么题目要求的答案就是\(\frac xy\),其中\(x=\sum_{len|n}\frac{f_{len}}{len},y=\sum_{len|n}\frac{cnt_{len}}{len}\)。
令\(a\)中数\(i\)出现的次数为\(b_i\)。现在枚举\(len\),尝试计算出\(f_{len}与cnt_{len}\)。令\(B=\frac n{len}\),显然对于任意i,必须满足\(B|b_i\),否则就不存在最短循环节为len的排列。然后想到用容斥,先计算出实际最小循环节长度为\(len\)的约数的排列的权值和与总个数。令\(b"_i=\frac{b_i}B\)。显然总个数为\(\frac{len!}{\prod{b"_i!}}\)。对于权值和,注意到整个排列被分成了\(B\)段,每一段都相等,所以可以把每一段都视为首尾相接的一个环形排列计算权值和,然后乘\(B\)。在每一段构成的小环形排列中,任意相邻两个位置贡献1的概率都相同,因为这个排列是随机的。这个概率的计算方法是先枚举靠前一个位置的值为\(i\),其发生的概率为\(\frac{b"_i}{len}\);后面一个位置的值必须与其不同,概率为\(\frac{len-b"_i}{len-1}\)。枚举所有的i,把这个概率相加即可。先暂时让\(f_{len}=\)此时算出的权值和,\(cnt_{len}=\)此时算出的总个数,然后进行这样一步容斥就能算出真正的\(f和cnt\):
for(int i=1;i<=n;++i) for(int j=i+i;j<=n;j+=i) (f[j]+=MOD-f[i])%=MOD,(cnt[j]+=MOD-cnt[i])%=MOD;
时间复杂度\(O(能过)\)。不知道怎么证明,但是看上去感觉不会被卡。
点击查看代码
#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 t,n,a[1000010],b[1000010],f[1000010],cnt[1000010],fac[1000010],inv[1000010],rinv[1000010];vector v;void solve(LL len){ LL dv=n/len; rep(i,v.size()) v[i]/=dv; cnt[len]=fac[len];rep(i,v.size()) (cnt[len]*=inv[v[i]])%=MOD; LL possi=0,mul=rinv[len]*rinv[len-1]%MOD; rep(i,v.size()) (possi+=v[i]*(len-v[i])%MOD*mul)%=MOD; f[len]=possi*n%MOD*cnt[len]%MOD; rep(i,v.size()) v[i]*=dv;}int main(){ fileio(); fac[0]=1;repn(i,1000005) fac[i]=fac[i-1]*i%MOD; rep(i,1000003) inv[i]=qpow(fac[i],MOD-2),rinv[i]=qpow(i,MOD-2); cin>>t; rep(tn,t) { scanf("%lld",&n); rep(i,n+3) b[i]=f[i]=cnt[i]=0; rep(i,n) scanf("%d",&a[i]),++b[a[i]]; v.clear(); repn(i,n) if(b[i]) v.pb(b[i]); if(v.size()==1) { puts("1"); continue; } repn(len,n) if(n%len==0) { bool bad=false; LL dv=n/len; rep(i,v.size()) if(v[i]%dv!=0){bad=true;break;} if(bad) continue; solve(len); } //cout<
-
Codeforces 1630 E Expected Components 题解 (组合数学)
题目链接首先明确概念:排列。指的就是一个把数组a重排得到的序列,两个排列相等当且仅当它们对应位全都...
来源: -
头条:Java基础项目:超市管理项目
JavaSE项目:超市管理系统项目:上面项目所用的表结构:实体类:@Setter@Getter@NoArgsConstructor@AllArgsConstructorpubliccl
来源: Codeforces 1630 E Expected Components 题解 (组合数学)
头条:Java基础项目:超市管理项目
每日动态!2023春运车票24日开售 除夕车票要等到1月7日
《妮姬》首月收入突破6.9亿!腾讯海外收入占比提升达12.5%
今日热门!超越电竞机!Redmi K60要榨干第二代骁龙8:画质、帧率、亮度三不降
每日精选:效果堪比镀铬 2.2万元的特斯拉Model Y新配色值不值?
环球观察:高可用 Canal集群 实操( 秒懂 + 史上最全)
微头条丨认证管理(锐捷业软篇)
Intel拆分GPU部门 一把手重回技术岗 累计亏损超20亿美金
天天通讯!iPhone 14 Pro爆出“闪线门” :屏幕出现诡异的绿色和黄色细横线
全球热点评!当ChatGPT遇上弱智吧:全程爆笑
夫妻的世界翻拍哪部电视剧?夫妻的世界最后结局是什么意思?
滕王阁为什么叫阁不叫楼?滕王阁为什么是三大名楼之首?
情非情砸车是第几集?情非情盖总和保姆的结局是什么?
小昭去波斯是哪一集?小昭去波斯后她母亲去哪儿了?
男人是大猪蹄子是什么意思?男人是大猪蹄子女人是什么?
排序算法模板(更新中)
当前速读:机器学习——果蔬分类
每日消息!性能超越电竞手机!Redmi K60 Pro综合跑分达135万
信息:千万别强忍 20岁小伙憋气压抑咳嗽导致昏厥
特斯拉今年股价累计暴跌超60%!马斯克透露大跌原因
收购动视暴雪遇阻 微软哭弱:根本打不过索尼、任天堂
到手9袋!良品铺子坚果礼盒1440 仅44元包邮
每日讯息!教你用JavaScript实现背景图像滑动
户外运动有哪些项目?户外运动品牌排行榜
什么鱼营养价值最高?什么鱼只会逆流而上?
金木水火土命怎么算出来的?金木水火土哪个腿长?
玉面小飞龙是什么意思?玉面小飞龙出自哪里?
Redmi K60系列上架:三颗口碑最好的芯片都拿到了 12月27日发
每日聚焦:最快闪充旗舰!真我GT Neo5充电头曝光:支持240W充电
环球热讯:紫米裁员80%并入小米?官方澄清:ZMI品牌将继续存在
全球新资讯:9.99万元遭疯抢 五菱宏光MINI EV敞篷版下线:能跑280km
苯胺皮是什么皮?苯胺皮和纳帕皮有什么区别?
世界新动态:CloudCanal实战-五分钟搞定Oracle到StarRocks数据迁移与同步
(一)elasticsearch 编译和启动
【速看料】马斯克辞任CEO,产品经理如何用项目协作软件武装自己?
焦点速讯:字节鏖战美团的关键一役
重点聚焦!糗事百科宣布将关闭服务 自侃“享年17岁”
全球观点:神似苹果AirPower!特斯拉推出无线充电板:最高功率15W
手慢无 民族品牌两面针牙膏大促:四支到手20元还送牙刷
又一新能源品牌官宣涨价:最少涨5千 今年买车还剩最后一周“窗口期”
全球速看:盘点适合《战神》奎爷的演员:道恩·强森、杰森·莫玛等
新型复兴号CR200J首次亮相:Wi-Fi全覆盖 充电插口增加
环球微动态丨比亚迪DM-i再外放 东风小康风光蓝电E5官图发布:综合续航1150km
霍乱疫情卷土重来:已致马拉维国410人死亡
环球今热点:随身咖啡馆 精神X小时:Nevercoffee咖啡1.99元(京东5元)
天天微头条丨什么是 HTML5?
每日消息!Ubuntu:Docker 容器操作
天天关注:苹果降低中国工厂依赖:真要搬走?iPhone 14制造难度降低
全球聚焦:不装了!日本万亿重新发展核能:新一代核反应堆准备中
【热闻】冬至湖南浏阳全城燃放烟花 满城烟花一河诗画:网友羡慕哭
焦点简讯:顺丰又上热搜!买Chanel耳钉顺丰运掉五颗珍珠
焦点热门:修复RX 7900显卡功耗异常 AMD新驱动实测:有用 但没什么大用
天天简讯:比iPhone 14 Pro Max更轻更便宜 OPPO Find N2首销:7999元
4插槽怪兽 华硕、猫头鹰合作打造最安静、最冷静的RTX 4090/4080显卡
动态:5.2万亿财富没了 特斯拉股东喊话马斯克:别只顾着推特了
世界微速讯:小岛秀夫:只有Xbox懂我
天天通讯!本田思域Type R各国/地区售价曝光 在日本才卖20多万?
每日短讯:负债585.68亿:国美获黄光裕公司三笔贷款累计5亿港元
全球快看点丨新能源车国补退场倒计时!车企打响价格战:现金立减、保险补贴
时隔半年 终于不寂寞!讯景发布全球第二款RX 6700
中国哪里的羊肉最好吃?这5个地方 你最爱谁?
后壳质感堪比玉石!vivo S16 Pro图赏
微软重构资源管理器进程:Windows 11运行速度大提升
支付宝接入技术
Python requests库指定IP请求,并使用HTTPS证书验证
世界今热点:MAUI新生4.5-字体图像集成Font&Image
精彩看点:Codeforces 1654 G Snowy Mountain 题解 (重心分治)
美国遭史上最严重禽流感疫情:鸡蛋价格创纪录 真吃不起节奏
环球速看:FreeSWITCH学习笔记:Lua脚本
每日短讯:剪映上线团队剪辑“神技”:异地多端一起剪视频成为可能
3299元起 vivo S16 Pro手机发布:首发双面柔光人像拍摄
环球信息:童年的味道 大白兔奶糖促销:1斤20元到手
环球聚焦:自拍绝了!vivo发布新机S16e:2099元起、行业首创“玉质玻璃”工艺
软链接和硬链接
世界热消息:渗透实录-02
雷军宣布小米人事调整:总裁王翔退休 卢伟冰晋升
环球热消息:特斯拉北美大降价5.2万 超高折扣只为保住销量?
vivo S16系列亮相:标准版搭载骁龙870 Pro版搭载天玑8200
vivo S16系列出厂预装OriginOS 3 虚拟内存提升8G
vivo S16系列7.36mm机身塞进4600mAh:苹果iPhone都没做到
快资讯:教你用JavaScript实现鼠标特效
【天天聚看点】男子开宝马专挑外地牌照车碰瓷 套路防不胜防:扔石头制造声响
迪士尼神话剧《美生中国人》新剧照:杨紫琼饰演观音 吴彦祖变身孙悟空
环球滚动:FIFA年终国家队排名:国足降至第80 美加墨世界杯出线希望增加
【报资讯】RTX 40系列移动显卡参数曝光:价格可不低
世界通讯!一加首款键盘曝光:全铝机身、自带USB-C/A接口
当前动态:基于Netty的IM聊天加密技术学习:一文理清常见的加密概念、术语等
JDK源码分析实战系列-PriorityBlockingQueue
资讯:Altium Designer v23.0.1.38图文详解
焦点热议:历史总是惊人地相似:复古主机Atari VCS宣告停产 一个时代终结
网友称考研民宿房费暴涨近20倍 店家:每年都一样
中国企业站稳全球LCD市场!李东生:TCL部分技术领先三星
消息!超可爱!《王者荣耀》梦奇赛年皮肤来了 特效贼棒
环球热讯:《王者荣耀》《合金弹头》联动:联名首发新英雄莱西奥
【世界播资讯】高能吸水 洁丽雅纯棉毛巾:15.9元/3条
今日报丨Intel显卡事业部突然解散!掌门人Raja回归首席架构师
全球新消息丨国内油价要止步“三连跌”!元旦后或迎新一轮价格上调
世界微头条丨比iPhone 14 Pro Max还轻 OPPO Find N2明天首销:7999元
Go 快速入门指南 - 环境安装