首页 > 其他分享 >单位根反演小记

单位根反演小记

时间:2024-04-27 14:22:32浏览次数:25  
标签:le frac sum 单位根 反演 binom omega 小记

反演公式

\[[n | v] = \frac{1}{n}\sum_{0 \le j < n}(\omega_n^v)^j \]

证明很简单,等比数列求和即可。

应用

牛客Wannafly挑战赛11E 白兔的***难

题意:

给定 \(k \le 2^20, n \le 10^{16}, p = 998244353\),求 \(t \in [0, k)\),\(a_t =\sum_{k | i,0 \le i + t \le n} \binom{n}{i + t}\)。

思路:

直接单位根反演:

\[\begin{aligned} \sum_{k | i,0 \le i + t \le n} \binom{n}{i + t} &= \sum_{0 \le i\le n} \binom{n}{i}[k | i - t]\\ &= \sum_{0 \le i\le n} \binom{n}{i}\frac{1}{k}\sum_{j = 0}^{k-1}(\omega_{k}^{i -t})^j\\ &= \frac{1}{k}\sum_{j = 0}^{k-1}\omega_{k}^{-tj}\sum_{0 \le i\le n} \binom{n}{i}\omega_{k}^{ij}\\ &= \frac{1}{k}\sum_{j = 0}^{k-1}\omega_{k}^{-tj}(1 + \omega_k^j)^n\\ \end{aligned} \]

不妨令 \(y_j = (1 + \omega_k^j)^n\),则上式就是一个 IDFT 形式的东西,利用 NTT 解决即可。

P5293 [HNOI2019] 白兔之舞

设 \(W\) 是 \(n \times n\) 的邻接矩阵,则我们有走 \(m\) 步的方法数等于:

\[\binom{L}{m}W^m \]

所以我们要求的就是对于所有 \(t\):

\[\sum_{m=0}^L\binom{L}{m}W^m[k | m - t] \]

单位根反演可以化为:

\[\frac{1}{k}\sum_{j=0}^{k-1}\omega_k^{-tj}(1 + \omega_k^jW)^m \]

不妨设 \(A_j = (I + \omega_{k}^jW)^L\),\(a_j\) 为 \(W\) 中编号第 \(j\) 个的系数,则我们要求的就是 \(IDFT(a)\)。使用 MTT 即可解决。

【牛客Wannafly挑战赛23 F】计数

直接大力推式子:

\[\begin{aligned} ANS &= \sum_{T}[k | \sum_{e \in T}val_e]\\ &= \frac{1}{k}\sum_{T}\sum_{j=0}^{k-1}(\omega_{k}^{\sum_{e \in T}val_e})^j\\ &= \frac{1}{k}\sum_{j=0}^{k-1}\sum_{T}\prod_{e \in T}(\omega_k^{j})^{val_e} \end{aligned} \]

转化成标准的矩阵树定理形式求解,时间复杂度 \(O(kn^3)\)。

标签:le,frac,sum,单位根,反演,binom,omega,小记
From: https://www.cnblogs.com/rlc202204/p/18161998

相关文章

  • 点分治小记
    点分治是一类高效统计树上路径问题的算法,通过优化递归深度的方法来有效保证时间复杂度。具体操作一般是以下几步:找到当前子树的重心以重心为根计算经过根节点的路径对答案的贡献将根删去并递归处理它的所有子树因为我们每次都以树的重心来作为根节点,递归深度不会超过......
  • 网络流小记
    基本定义:网络:一张有向图。流量:经过一条边的流的大小,一条边\((u,v)\)的流量记为\(flow(u,v)\),一个网络的流量定义为\(∑f(s,x)\)。容量:一条边的流量上限,一条边\((u,v)\)的容量记为\(cap(u,v)\)。费用:经过一条边单位流量的所需费用,一条边\((u,v)\)的费用记为......
  • 小记 Demo
    定义领域模型:AggregateRoot:定义SalesOrder作为聚合根,其中包括订单明细、客户信息、订单总额等。Entity:定义如OrderItem(订单项)、Inventory(库存)等实体。ValueObject:定义如Address(地址)、Money(金额)等值对象。建立仓储接口:使用ABPvNext框架的仓储模式来实现数据的......
  • 博弈论小记
    以下我们都考虑这样一种游戏:两个人,轮流进行;游戏总是在有限步内结束;同一个状态不可能多次抵达,且没有平局;每个时刻的合法决策集合仅与当前局面有关,而与游戏者无关;不能操作者输。我们定义:必败态:无论如何先手必败的状态(局面)。必胜态:先手存在必胜策略的状态(局面)。......
  • hexo 折腾小记
    hexo是一套静态网页生成框架类似的还有JekyllGitHub默认推荐的框架/Hugo(......
  • windows平台vs2019编译Luabind小记
    前言写这篇文章的目的是Luabind这个库比较老旧,对于新编译器需要做一些代码上的兼容,参考资料又都有点过时,所以特写此篇,记录踩坑过程;参考资料用VS2010编译luabind如何编译luabind支持vs2010之后所有版本Lua官网LuabindRepo编译前准备准备相关前置组件基本编译依赖Des......
  • 1.0 多层感知机&BP传播 小记
    1.感知机与线性模型单层感知机的表达式和线性分类表达式等同,可以将一个单层感知机看作是一个线性分类器。单层感知机可以解决与、或、非的分类问题,但是不能解决异或分类(非线性)问题。howtosolvetheproblem:多个线性分类器解决线性不可分问题,即:多个单层感知机组合叠加解......
  • 13年 史密斯热水器全拆清洗 小记
    ~~~~~~~~~~~~~~~排水~~~~~~~~~~~~~~~拆机过程不复杂,提前断开电源,切断进水阀门.用过导水管从进水口排干净热水器中的水即可,注意排水过程需要打开热水器出水口,用来进气.排干净水之后就可以拆下来进出水管,这样热水器就从水路上分离开来了.把热水器从挂墙上拿下来,50L的......
  • PWN出题小记
    记录一下PWN出题的源码、环境部署。PWN题环境部署图方便的话可以使用pwn_deploy_chroot这个项目。如何安全快速地部署多道ctfpwn比赛题目也可以自己写dockerfile拉取镜像。将题目名命名为pwn,与Dockerfile、ctf.xinetd、start.sh三个文件放在同一目录下,下面提供相......
  • STM32外部中断小记
    一、EXTI配置步骤//1.配置RCC时钟RCC_APB2PeriphClockCmd(RCC_APB2Periph_GPIOx,ENABLE);RCC_APB2PeriphClockCmd(RCC_APB2Periph_AFIO,ENABLE);//开启AFIO时钟,AFIO:GPIO复用/重映射功能//2.配置EXTIGPIO端口及工作模式(输入模式)//3.配置EXTI中断线、模式(上升沿、下降沿......