P5591-小猪佩奇学数学【单位根反演】

正题

题目链接:https://www.luogu.com.cn/problem/P5591


题目大意

给出n,p,kn,p,kn,p,k
(∑i=0n(ni)pi⌊ik⌋)mod998244353\left(\sum_{i=0}^n\binom{n}{i}p^i\left\lfloor\frac{i}{k}\right\rfloor \right)\mod 998244353(i=0n(in)piki)mod998244353

1≤n,p<998244353,k=2w(w∈[0,20])1\leq n,p<998244353,k=2^w(w\in[0,20])1n,p<998244353,k=2w(w[0,20])


解题思路

开始以为推错了,结果是要特判
推出了看上去不是我能推的式子

∑i=1n(ni)pi∑j=1i[k∣j]\sum_{i=1}^n\binom{n}{i}p^i\sum_{j=1}^i[k|j]i=1n(in)pij=1i[kj]
然后单位根反演
∑i=1n(ni)pi∑j=1i1k∑l=0k−1ωkl×j\sum_{i=1}^n\binom{n}{i}p^i\sum_{j=1}^i\frac{1}{k}\sum_{l=0}^{k-1}\omega_k^{l\times j}i=1n(in)pij=1ik1l=0k1ωkl×j
系统整理一下
1k∑l=0k−1∑i=1n(ni)pi∑j=1iωkl×j\frac{1}{k}\sum_{l=0}^{k-1}\sum_{i=1}^n\binom{n}{i}p^i\sum_{j=1}^i\omega_k^{l\times j}k1l=0k1i=1n(in)pij=1iωkl×j
然后等比数列通项公式拆开
1k∑l=0k−1∑i=1n(ni)piωkl−ωkliωkl1−ωkl\frac{1}{k}\sum_{l=0}^{k-1}\sum_{i=1}^n\binom{n}{i}p^i\frac{\omega_k^l-\omega_{k}^{li}\omega^l_k}{1-\omega_k^l}k1l=0k1i=1n(in)pi1ωklωklωkliωkl
1k∑l=0kωkl1−ωkl(∑i=1n(ni)pi−∑i=1n(ni)piωkli)\frac{1}{k}\sum_{l=0}^k\frac{\omega_k^l}{1-\omega_k^l}\left(\sum_{i=1}^n\binom{n}{i}p^i-\sum_{i=1}^n\binom{n}{i}p^i\omega_k^{li}\right)k1l=0k1ωklωkl(i=1n(in)pii=1n(in)piωkli)
1k∑l=0kωkl1−ωkl((p+1)n−(pωkl+1)n)\frac{1}{k}\sum_{l=0}^k\frac{\omega_k^l}{1-\omega_k^l}\left((p+1)^n-(p\omega_k^l+1)^n\right)k1l=0k1ωklωkl((p+1)n(pωkl+1)n)

然后写出来会愉快的发现没有过样例,仔细看我们式子里面有一个ωkl1−ωkl\frac{\omega_k^l}{1-\omega_k^l}1ωklωkl

l=0l=0l=0 的时候1−ωkl=01-\omega_k^l=01ωkl=0,所以不能直接这么求,我们这个得分开考虑。

就是
∑i=1n(ni)pii⇒∑i=1nn!(i−1)!(n−i)!pi\sum_{i=1}^n\binom{n}{i}p^ii\Rightarrow \sum_{i=1}^n\frac{n!}{(i-1)!(n-i)!}p^ii=1n(in)piii=1n(i1)!(ni)!n!pi
np∑i=1n(n−1)!(i−1)!(n−i)!pi−1⇒n∑i=1n(n−1i)pi−1⇒np(p+1)n−1np\sum_{i=1}^n\frac{(n-1)!}{(i-1)!(n-i)!}p^{i-1}\Rightarrow n\sum_{i=1}^n\binom{n-1}{i}p^{i-1}\Rightarrow np(p+1)^{n-1}npi=1n(i1)!(ni)!(n1)!pi1ni=1n(in1)pi1np(p+1)n1
就好了

时间复杂度:O(klog⁡P):O(k\log P):O(klogP)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll P=998244353;
ll n,p,k,ans;
ll power(ll x,ll b){ll ans=1;x%=P;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
signed main()
{scanf("%lld%lld%lld",&n,&p,&k);ll d=power(3,(P-1)/k);ans=n*p%P*power(p+1,n-1)%P;for(ll i=1,w=d;i<k;i++,w=w*d%P){ll inv=power(P+1-w,P-2)*w%P;ans+=power(p+1,n)*inv%P;ans-=power(w*p+1,n)*inv%P;ans=(ans+P)%P;}printf("%lld\n",ans*power(k,P-2)%P);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318310.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

YBTOJ:字符串题(KMP)

文章目录题目描述解析1.p[i]>0:2.p[i]0&#xff1a;代码题目描述 解析 我个人做起来很费劲的一道题&#xff0c;用vector&#xff0c;并查集等等东西搞来搞去过掉了&#xff08;竟然只WA了一次 &#xff09; 看题解思路就一下子清晰了&#xff0c;还是对KMP的理解不到位。 …

Strange List CodeForces - 1471B

题意&#xff1a; 对于长度为n的数组a&#xff0c;从第一位开始如果可以整除x&#xff0c;就将x个a/x的结果加到数组最后&#xff0c;然后对下一位进行一样的操作&#xff0c;直到第x位不可以整除x&#xff0c;到此结束&#xff0c;然后计算此时数组的总和 题解&#xff1a; …

Asp.Net Core 轻松学-正确使用分布式缓存

前言本来昨天应该更新的&#xff0c;但是由于各种原因&#xff0c;抱歉&#xff0c;让追这个系列的朋友久等了。上一篇文章 在.Net Core 使用缓存和配置依赖策略 讲的是如何使用本地缓存&#xff0c;那么本篇文章就来了解一下如何使用分布式缓存&#xff0c;通过本章&#xff0…

计算几何初步

两点之间距离 欧氏距离 即欧几里得距离。 平面内两点的距离为 \[\sqrt{(x_1-x_2)^2(y_1-y_2)^2} \]立体空间内两点的距离为 \[\sqrt{(x_1-x_2)^2(y_1-y_2)^2(z_1-z_2)^2} \]\(\dots\) \(n\) 维空间内两点的距离为 \[\sqrt{\sum_{i1}^{n}{(x_1-x_2)^2}} \]曼哈顿距离 二维空间内…

树链剖分概念及模板 + 例题 [POJ3237 tree + 软件包管理器]

文章目录概念模板例题1&#xff1a;软件包管理器题目题解代码实现例题2&#xff1a;POJ3237 tree题目题解代码实现概念 树链剖分主要是用于解决以下这两个问题。 1、更改树上点x到点y的最短路径上的所有结点的值 2、查询树上点x到点y的最短路径上的所有结点的和。 在讲树链剖分…

P3706-[SDOI2017]硬币游戏【高斯消元,字符串hash】

正题 题目链接:https://www.luogu.com.cn/problem/P3706 题目大意 给出 nnn 个长度为 mmm 的 H/TH/TH/T 串。 开始一个空序列&#xff0c;每次随机在后面加一个 H/TH/TH/T &#xff0c;求每个串第一次出现的概率。 1≤n,m≤3001\leq n,m\leq 3001≤n,m≤300 解题思路 数据范…

YBTOJ:字符串匹配(KMP)

文章目录题目描述解析代码题目描述 解析 看了题解。。。 这题的关键在于可以变换匹配的一个充要条件&#xff1a; 每个字符与前一个相同字符的距离相同 这个搞出来之后就可以以它为关键字进行KMP了 注意&#xff01; 当与前一个字符的距离超过匹配长度时&#xff0c;是没有意义…

Strange Birthday Party CodeForces - 1471C

Strange Birthday Party CodeForces - 1471C 题意&#xff1a; 我有n个朋友&#xff0c;商店有m种商品&#xff0c;这m种商品按序号价格从小到大排列&#xff0c;对于每一个朋友我给出一个序号k&#xff0c;我可以直接给朋友序号k的商品价格的金钱或给朋友买一个序号小于k的商…

TransactionScope事务处理方法介绍及.NET Core中的注意事项

作者&#xff1a;依乐祝原文链接&#xff1a;https://www.cnblogs.com/yilezhu/p/10170712.html今天在写CzarCms的UnitOfWork的使用使用到了这个TransactionScope事务&#xff0c;因此对它进行了相关资料的查阅并记录如下&#xff0c;希望对大伙在.NET Core中使用有所帮助。写在…

P4548-[CTSC2006]歌唱王国【概率生成函数,KMP】

正题 题目链接:https://www.luogu.com.cn/problem/P4548 题目大意 ttt次询问&#xff0c;给出一个长度为mmm的串SSS和一个空串TTT&#xff0c;每次在TTT后面随机加入1∼n1\sim n1∼n的字符&#xff0c;直到TTT中出现SSS为止&#xff0c;求期望次数。 1≤n≤105,t≤50,1≤m≤10…

YBTOJ:字符匹配(KMP)

文章目录题目描述解析代码题目描述 解析 显然应该要尝试套kmp的板子 关键是如何套 也就是那个判断匹配的条件是什么的问题 本题的关键是当kmp匹配时&#xff0c;匹配位之前的所有位大小关系的顺序都是匹配的&#xff0c;所以我们只需要看当前位即可 考虑对b预处理出3个数组&a…

历史上的今天(history)+ 勇者斗恶龙(dragon)

朋友们我来了&#xff0c;好久都没有更新了&#xff0c;手实在痒痒&#xff0c;不撸两道&#xff0c;内心过意不去 文章目录A&#xff1a;历史上的今天&#xff08;history&#xff09;题目题解代码实现B&#xff1a; 勇者斗恶龙&#xff08;dragon&#xff09;题目描述题解代…

Codeforces Round #694 Div. 2

Codeforces Round #694 Div. 2 CodeForces 1471A Strange Partition CodeForces 1471B Strange List CodeForces 1471C Strange Birthday Party 签到题 稍有难度 CodeForces 1471D Strange Definition CodeForces 1471E Strange Shuffle CodeForces 1471F Strange Housing

领域驱动设计,让程序员心中有码(四)

#领域驱动设计&#xff0c;让程序员心中有码&#xff08;四&#xff09;----------------------追忆有关分层的古老往事我一直认为&#xff0c;程序员也是艺术家&#xff0c;他们撰写的每一行代码&#xff0c;是献给这大好世界的优美诗篇。不同的人&#xff0c;写的代码也许风格…

YBTOJ洛谷P4551:最长异或路径(trie树)

洛谷传送门 文章目录题目描述解析代码题目描述 解析 本题关键就在于一点&#xff1a; 若把每个点的深度dep[i]定义为从根到节点边权的异或和 那么i到j的路径异或和可以表示为&#xff1a; dep[i] ^ dep[j] 首先要是i、j在不同子树上显然成立 如果他们在同一子树上&#xff0c;…

P5074-Eat the Trees【插头dp】

正题 题目链接:https://www.luogu.com.cn/problem/P5074 题目大意 给出一个nmn\times mnm的网格&#xff0c;有的必须铺线有的不能&#xff0c;铺成若干条闭合回路&#xff0c;求方案数。 1≤n,m≤121\leq n,m\leq 121≤n,m≤12 解题思路 考虑插头dpdpdp&#xff0c;因为可以…

初赛—算法复杂度

对于 \(T(n) a\times T(\frac{n}{b})c\times n^k\) 这样的递归关系&#xff0c;有这样的结论&#xff1a; \(a>b^k\) 则有 &#xff1a; \(T(n) O(c\times n^{\log_{b} a})\)\(ab^k\) 则有 &#xff1a; \(T(n) O(c\times n^k\times log n)\)\(a<b^k\) 则有 &#xff…

1414. 牛异或

题意&#xff1a; 区间长度为n的数组a&#xff0c;取连续一段子序列&#xff0c;问哪段子序列的异或和最大&#xff0c;如果存在多个这样的序列&#xff0c;那么选择序列末端整数对应的编号更小的那个序列。 如果仍然存在多个可选的序列&#xff0c;那么选择长度最短的那个序列…

基于Asp.Net Core打造轻量级内部服务治理RPC(一)

继之前的《Asp.Net Core Docker 搭建》文章末尾说过的&#xff0c;将陆续编写基于asp.net core 打造一个内部服务治理的rpc框架。不过前端时间较忙&#xff0c;所以搁置了一段时间。闲话不多说&#xff0c;下面就来讲讲为什么需要去做一个该框架&#xff0c;以及想法的来源和设…

P7046-「MCOI-03」诗韵【SAM,倍增,树状数组】

正题 题目链接:https://www.luogu.com.cn/problem/P7046 题目大意 给出一个长度为 nnn 的字符串&#xff0c;然后 mmm 次把它的一个子串加入集合。如果一个字符串在这个集合中作为字符串的后缀出现次数大于 kkk 那么这个字符串就会被计入贡献。 每次求计入贡献的字符串数和最…