【集训队作业2018】复读机【指数型生成函数】【单位根反演】【二项式定理】

传送门

单位根反演听着高级,其实没啥技术含量……

本文是篇几乎没有证明的佛系讲解


单位根反演的式子长这样:

1n∑i=0n−1ωnik=[k∣n]\frac{1}{n}\sum_{i=0}^{n-1}\omega_n^{ik}=[k|n]n1i=0n1ωnik=[kn]

其实本质是IFFT

感觉懵的?

那直接记结论:对于一个多项式f(x)f(x)f(x),有

∑k∣i[xi]f(x)=1k∑i=0k−1f(ωix)\sum_{k|i}[x^i]f(x)=\frac{1}{k}\sum_{i=0}^{k-1}f(\omega^ix)ki[xi]f(x)=k1i=0k1f(ωix)

请参考下面的题解以更好地理解这个式子


d=1d=1d=1的时候直接快速幂走人

d=2d=2d=2的时候

考虑生成函数 带标号所以用EGF

得到一位群友的EGF为

f(x)=∑2∣i1i!xif(x)=\sum_{2|i}\frac{1}{i!}x^if(x)=2ii!1xi

也就是只有偶数项有值

可以构造出

f(x)=ex+e−x2f(x)=\frac{e^x+e^{-x}}{2}f(x)=2ex+ex

kkk位群友就是fk(x)f^k(x)fk(x)

强行二项式定理拆开,可以得到

12k∑i=0k(ki)[xn]ex(2i−k)\frac{1}{2^k}\sum_{i=0}^k\binom{k}{i}[x^n]e^{x(2i-k)}2k1i=0k(ik)[xn]ex(2ik)

对于[xn]ex(2i−k)[x^n]e^{x(2i-k)}[xn]ex(2ik),可以考虑它的组合意义,相当于用2i−k2i-k2ik种颜色涂nnn个格子,即(2i−k)n(2i-k)^n(2ik)n,然后得到的是EGF,也就是除以了一个n!n!n!;最终算答案会乘上一个n!n!n!,刚好抵消。

当然很不严谨,因为有负数的情况。严谨证明可以用泰勒展开。

然后暴力枚举即可

对于d=3d=3d=3

f(x)=∑3∣i1i!xif(x)=\sum_{3|i}\frac{1}{i!}x^if(x)=3ii!1xi

仿照上面的栗子,考虑到1,−11,-11,1111的二次方根,那么现在应该用333的单位根

大胆猜想

f(x)=ex+eωx+eω2x3f(x)=\frac{e^x+e^{\omega x}+e^{\omega^2x}}{3}f(x)=3ex+eωx+eω2x

后面就没啥区别了

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#define MAXN 500005
using namespace std;
const int MOD=19491001;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
typedef long long ll;
inline int qpow(int a,int p)
{int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD;p>>=1;}return ans;
}
#define inv(x) qpow(x,MOD-2)
const int w1=18827933,w2=663067;
int fac[MAXN],finv[MAXN];
inline void init(const int& N)
{fac[0]=1;for (int i=1;i<=N;i++) fac[i]=(ll)fac[i-1]*i%MOD;finv[N]=inv(fac[N]);for (int i=N-1;i>=0;i--) finv[i]=(ll)finv[i+1]*(i+1)%MOD;
}
inline int C(const int& n,const int& m){return (ll)fac[n]*finv[m]%MOD*finv[n-m]%MOD;}
int main()
{int n,k,d;scanf("%d%d%d",&n,&k,&d);init(k);if (d==1) return printf("%d\n",qpow(k,n)),0;if (d==2){int ans=0;for (int i=0;i<=k;i++) ans=add(ans,(ll)C(k,i)*qpow(dec(2*i,k),n)%MOD);ans=(ll)ans*inv(qpow(2,k))%MOD;printf("%d\n",ans);return 0;}int ans=0;for (int i=0;i<=k;i++)for (int j=0;j<=k-i;j++)ans=add(ans,(ll)C(k,i)*C(k-i,j)%MOD*qpow(add(add(i,(ll)j*w1%MOD),(ll)(k-i-j)*w2%MOD),n)%MOD);ans=(ll)ans*inv(qpow(3,k))%MOD;printf("%d\n",ans);return 0;
}

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

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

相关文章

ASP.NET Core on K8S学习初探(2)

“ [LOG] ASP.NET Core on K8S Starting...”在上一篇《单节点环境搭建》中&#xff0c;通过Docker for Windows在Windows开发机中搭建了一个单节点的K8S环境&#xff0c;接下来就是动人心弦的部署ASP.NET Core API到K8S了。但是&#xff0c;在部署之前&#xff0c;我还是把基本…

Educational Codeforces Round 96 E. String Reversa 线段树模拟序列交换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 与上一篇题解大同小异&#xff0c;无非就是不需要枚举排列了。 // Problem: E. String Reversal // Contest: Codeforces - Educational Codeforces Round 96 (Rated for Div. 2) // URL:…

【LOJ6363】「地底蔷薇」【点双】【指数型生成函数】【扩展拉格朗日反演】【多项式幂函数】

传送门 题意&#xff1a;给定nnn和集合SSS,求出nnn个点的「所有极大点双连通分量的大小都在SSS 内」的不同简单无向连通图的个数 模 998244353998244353998244353。 n,∑i∈Si≤105n,\sum_{i\in S}i \leq10^5n,∑i∈S​i≤105 道理我都懂&#xff0c;可为啥我百度搜地灵殿ex终…

ASP.NET Core on K8S学习初探(1)

“ [LOG] ASP.NET Core on K8S Starting...”01—写在之前当近期的一个App上线后&#xff0c;发现目前的docker实例&#xff08;应用服务BFF中台服务工具服务&#xff09;已经越来越多了&#xff0c;而我司目前没有专业的运维人员&#xff0c;发现运维的成本逐渐开始上来&#…

AtCoder Regular Contest 120 C - Swaps 2 线段树模拟

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个序列a,ba,ba,b&#xff0c;每次可以执行一个操作&#xff1a;将a[i]a[i]a[i]与a[i1]a[i1]a[i1]交换&#xff0c;且让交换后的a[i]1,a[i1]−1a[i]1,a[i1]-1a[i]1,a[i1]−1&#xff0c;问将aaa变成bbb…

【BZOJ 4671】异或图 【斯特林反演】【线性基】【贝尔数复杂度】

传送门 题意&#xff1a;定义两个图的异或的边集为在两张图中恰出现一次的边。给sss张nnn个点的图的集合&#xff0c;求异或和为连通图的子集数。 s≤60,n≤10s \leq 60,n \leq 10s≤60,n≤10 设GiG_iGi​表示异或出iii个连通块的子集数&#xff0c;答案就是G1G_1G1​ GGG并不…

.NET Core IdentityServer4实战 第Ⅳ章-集成密码登陆模式

回顾下ClientCredentials模式&#xff0c;在ReSourceApi中定义了我们公开服务&#xff0c;第三方网站想要去访问ReSourceApi则需要在身份验证服务中获取toekn&#xff0c;根据token的内容&#xff0c;硬编码去访问公开服务&#xff08;ResApi&#xff09;,这个还是非常简单的&a…

Codeforces Round #656 (Div. 3) F. Removing Leaves 贪心 + 模拟

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先有一个贪心策略就是每次都找一个叶子节点最多的点&#xff0c;让后删掉他的kkk个叶子节点&#xff0c;现在我们就来考虑如何模拟这个过程。 我们整一个vector<set<int>>ve…

.NetCore中三种注入生命周期的思考

.NetCore彻底诠释了“万物皆可注入”这句话的含义&#xff0c;在.NetCore中到处可见注入的使用。因此core中也提供了三种注入方式的生命周期使用&#xff0c;分别是&#xff1a;AddTransient&#xff1a;每次请求&#xff0c;都获取一个新的实例。即使同一个请求获取多次也会是…

兰州大学第一届 飞马杯 体育课排队 二分 + 最大流 + 输出路径

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 读懂题就会辣&#xff0c;经典模型了&#xff0c;二分时间&#xff0c;让后将其转换成二分图&#xff0c;左边是人&#xff0c;右边是位置&#xff0c;能在规定时间到的连边&#xff0c;跑…

从零开始实现ASP.NET Core MVC的插件式开发(一) - 使用Application Part动态加载控制器和视图...

如果你使用过一些开源CMS的话&#xff0c;肯定会用过其中的的插件化功能&#xff0c;用户可以通过启用或者上传插件包的方式动态添加一些功能&#xff0c;那么在ASP.NET Core MVC中如何实现插件化开发呢&#xff0c;下面我们来探究一下。本系列只是笔者的一些尝试&#xff0c;并…

感性理解Berlekamp-Massey算法

引入 BM算法主要解决的是根据数列求最短线性齐次递推式的问题在OI中主要辅助打表使用 即&#xff1a;已知FFF&#xff0c;求序列AAA使得 Fn∑i1mAiFn−i(n>m)F_n\sum_{i1}^mA_iF_{n-i} \quad(n>m)Fn​i1∑m​Ai​Fn−i​(n>m) 其中mmm尽量小 算法流程 文中的所有…

兰州大学第一届 飞马杯 ★★飞马祝福语★★ 线段树维护dp(动态dp)

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个串&#xff0c;每次将区间都修改为某一个字母&#xff0c;问最终包含多少个FeiMaFeiMaFeiMa子序列。 思路&#xff1a; 首先暴力修改肯定是不行的&#xff0c;复杂度nqnqnq。 如果没有修改操作&am…

兰州大学第一届 飞马杯 ★★快乐苹果树★★ 树链剖分 + 懒标记 + 树状数组

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 第一次听说树链剖分能在fa[top[i]]fa[top[i]]fa[top[i]]的地方加懒标记&#xff0c;学到了学到了。 首先不能被题目吓住&#xff0c;这个题目仔细剖析一下不难发现一些性质&#xff1a; 以…

【BZOJ3328】PYXFIB【矩阵快速幂】【单位根反演】【二项式定理】

传送门 题意&#xff1a; ∑i0⌊nk⌋(nik)Fik\sum_{i0}^{\lfloor\frac nk\rfloor}\binom n{ik}F_{ik}i0∑⌊kn​⌋​(ikn​)Fik​ FFF为斐波拉契数列 n≤1e18,k≤2e4,p≤1e9n\leq 1e18,k\leq2e4,p\leq1e9n≤1e18,k≤2e4,p≤1e9且为质数且模kkk余111 显然就是求 ∑i0n[k∣i](ni…

.NET分布式框架 | Orleans 知多少

引言公司物联网项目集成Orleans以支持高并发的分布式业务&#xff0c;对于Orleans也是第一次接触&#xff0c;本文就分享下个人对Orleans的理解。这里先抛出自己的观点&#xff1a;Orleans 是一个支持有状态云生应用/服务水平伸缩的基于Virtual Actor 模型的.NET分布式框架。下…

nowcoder 牛牛的最大兴趣组 质因子 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先nnn很小的话可以暴力连边&#xff0c;让后染个色求一个颜色最多的即可。但是这个题显然不行&#xff0c;由于是三次方&#xff0c;所以考虑质因子入手。 首先很容易就能想到将所有的数…

ASP.NET Core on K8S学习初探(3)部署API到K8S

“ 终于可以部署ASP.NET Core到K8S中了...”在上一篇《基本概念快速一览》中&#xff0c;我们把基本的一些概念快速地简单地不求甚解地过了一下&#xff0c;本篇开始我们会将ASP.NET Core WebAPI部署到K8S&#xff0c;从而结束初探的旅程。01—准备一个WebAPI这里准备一个空的A…

K-D Tree学习笔记

引入 K-D Tree 是一种处理高维空间的数据结构。 支持O(nk−1k)O(n^{\frac {k-1}k})O(nkk−1​)查询给定超矩形内的点的信息&#xff0c; kkk 为维数。 可以用替罪羊树的思想实现动态插入。 但其实用的最多的是在错误的复杂度内查询奇奇怪怪的点信息。 构建 K-D Tree 是平…

Codeforces Round #633 B. Edge Weight Assignment 结论题 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先考虑最小值&#xff0c;如果从一个叶子结点出发到任意叶子的距离都为偶数&#xff0c;那么只需要一个值就可以满足条件。如果有奇数的&#xff0c;考虑111 ^ 222 ^ 303030&#xff0c;…