原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ450.html

题解

首先有一个东西叫做“单位根反演”,它在 FFT 的时候用到过:

$$\frac 1 n \sum_{i=0}^{n-1} \omega_n ^{d\cdot i} = [n|d]$$

其中 $\omega_n$ 表示 $n$ 次单位根。

接下来我们回到本题。

我们来搞一个指数生成函数,第 i 项表示总共复读 i 次,使得一个复读机开心的方案。

$$f(x) = \sum_{i\geq 0} [d|i] \frac{x^i} {i!}$$

那么我们要求的东西就是:

$$f^k(x)[n]$$

我们来给 $f(x)$ 推一推式子:

$$f(x) = \sum_{i\geq 0} [d|i] \frac{x^i} {i!} \\ = \sum_{i\geq 0} \left(  \frac 1 d \sum_{j=0}^{d-1} \omega_d^{ij} \right) \frac{x^i}{i!}\\ = \frac 1 d \sum_{j=0}^{d-1} \sum_{i\geq 0} \frac{x^i (\omega _d ^j)^i}{i!}\\ = \frac 1 d \sum_{i=0}^{d-1} e^{\omega_d ^i x}$$

在 NTT 的时候,我们用原根的幂来代替单位根。

我们发现 19491001 是个质数,它的最小原根是 7 ,而且 19491000 = 2*2*2*3*5*5*5*73*89,含有因子 2 和 3,这说明能找到原根来分别代替 $\omega _2$ 和 $\omega_3$ 。

接下来我们分情况讨论:

d = 1 : ans = $k^n$ 。

d = 2 :

$$f(x) = \frac 12 (e^x + e^{-x})$$

$$f^k(x)[n] = (\frac 1 2 )^k\sum_{i=0}^{k} \binom k i c^n $$

由于 $k\leq 500000$,直接爆算就好了。

d = 3 :

$$f(x) = \frac 13 (e^x + e^{\omega_3 x } + e^{\omega_3^2 x })$$

注意到由于 d = 3 时, $k\leq 1000$ ,所以和 $d = 2 $ 的情况差不多,暴力展开 2 层就好了。

具体怎么做直接看代码吧。懒得码式子了。

时间复杂度 $O(k^{d-1})$ 。

代码

#pragma GCC optimize("Ofast","inline")
#include <bits/stdc++.h>
#define clr(x) memset(x,0,sizeof (x))
#define For(i,a,b) for (int i=a;i<=b;i++)
#define Fod(i,b,a) for (int i=b;i>=a;i--)
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define _SEED_ ('C'+'L'+'Y'+'A'+'K'+'I'+'O'+'I')
#define outval(x) printf(#x" = %d\n",x)
#define outvec(x) printf("vec "#x" = ");for (auto _v : x)printf("%d ",_v);puts("")
#define outtag(x) puts("----------"#x"----------")
#define outarr(a,L,R) printf(#a"[%d...%d] = ",L,R);\
For(_v2,L,R)printf("%d ",a[_v2]);puts("");
using namespace std;
typedef long long LL;
LL read(){
LL x=0,f=0;
char ch=getchar();
while (!isdigit(ch))
f|=ch=='-',ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return f?-x:x;
}
const int N=500005,mod=19491001,G=7;
int Pow(int x,int y){
int ans=1;
for (;y;y>>=1,x=(LL)x*x%mod)
if (y&1)
ans=(LL)ans*x%mod;
return ans;
}
int n,k,d;
int Fac[N],Inv[N];
void prework(){
int n=N-1;
for (int i=Fac[0]=1;i<=n;i++)
Fac[i]=(LL)Fac[i-1]*i%mod;
Inv[n]=Pow(Fac[n],mod-2);
Fod(i,n,1)
Inv[i-1]=(LL)Inv[i]*i%mod;
}
int C(int n,int m){
if (m>n||m<0)
return 0;
return (LL)Fac[n]*Inv[m]%mod*Inv[n-m]%mod;
}
int main(){
prework();
n=read(),k=read(),d=read();
if (d==1){
cout<<Pow(k,n)<<endl;
}
else if (d==2){
int ans=0;
For(i,0,k){
int c=(i-(k-i)+mod)%mod;
ans=((LL)C(k,i)*Pow(c,n)+ans)%mod;
}
ans=(LL)ans*Pow(2,mod-1-k)%mod;
cout<<ans<<endl;
}
else {
int ans=0;
int w0=1,w1=Pow(G,(mod-1)/3),w2=(LL)w1*w1%mod;
For(i,0,k)
For(j,0,k-i){
int c=((LL)w0*i+(LL)w1*j+(LL)w2*(k-i-j))%mod;
ans=((LL)C(k,i)*C(k-i,j)%mod*Pow(c,n)+ans)%mod;
}
ans=(LL)ans*Pow(3,mod-1-k)%mod;
cout<<ans<<endl;
}
return 0;
}

  

UOJ#450. 【集训队作业2018】复读机 排列组合 生成函数 单位根反演的更多相关文章

  1. uoj #450[集训队作业2018]复读机

    传送门 \(d=1\),那么任何时刻都可以\(k\)个复读机的一种,答案为\(k^n\) \(d>1\),可以枚举某个复读机的复读次数(必须是\(d\)的倍数),然后第\(i\)个复读时间为\( ...

  2. 【UOJ#450】[集训队作业2018] 复读机

    题目链接 题目描述 群里有\(k\)个不同的复读机.为了庆祝平安夜的到来,在接下来的\(n\)秒内,它们每秒钟都会选出一位优秀的复读机进行复读.非常滑稽的是,一个复读机只有总共复读了\(d\)的倍数次 ...

  3. [2018集训队作业][UOJ450] 复读机 [DP+泰勒展开+单位根反演]

    题面 传送门 思路 本文中所有$m$是原题目中的$k$ 首先,这个一看就是$d=1,2,3$数据分治 d=1 不说了,很简单,$m^n$ d=2 先上个$dp$试试 设$dp[i][j]$表示前$i$ ...

  4. UOJ 449 【集训队作业2018】喂鸽子 【生成函数,min-max容斥】

    这是第100篇博客,所以肯定是要水过去的. 首先看到这种形式的东西首先min-max容斥一波,设\(f_{c,s}\)表示在\(c\)只咕咕中,经过\(s\)秒之后并没有喂饱任何一只的概率. \[ \ ...

  5. UOJ 422 [集训队作业2018] 小Z的礼物 min-max容斥 期望 轮廓线dp

    LINK:小Z的礼物 太精髓了 我重学了一遍min-max容斥 重写了一遍按位或才写这道题的. 还是期望多少时间可以全部集齐. 相当于求出 \(E(max(S))\)表示最后一个出现的期望时间. 根据 ...

  6. 【UOJ#450】【集训队作业2018】复读机(生成函数,单位根反演)

    [UOJ#450][集训队作业2018]复读机(生成函数,单位根反演) 题面 UOJ 题解 似乎是\(\mbox{Anson}\)爷的题. \(d=1\)的时候,随便怎么都行,答案就是\(k^n\). ...

  7. uoj450 【集训队作业2018】复读机(生成函数,单位根反演)

    uoj450 [集训队作业2018]复读机(生成函数,单位根反演) uoj 题解时间 首先直接搞出单个复读机的生成函数 $ \sum\limits_{ i = 0 }^{ k } [ d | i ] ...

  8. UOJ #449. 【集训队作业2018】喂鸽子

    UOJ #449. [集训队作业2018]喂鸽子 小Z是养鸽子的人.一天,小Z给鸽子们喂玉米吃.一共有n只鸽子,小Z每秒会等概率选择一只鸽子并给他一粒玉米.一只鸽子饱了当且仅当它吃了的玉米粒数量\(≥ ...

  9. 【UOJ#422】【集训队作业2018】小Z的礼物(min-max容斥,轮廓线dp)

    [UOJ#422][集训队作业2018]小Z的礼物(min-max容斥,轮廓线dp) 题面 UOJ 题解 毒瘤xzy,怎么能搬这种题当做WC模拟题QwQ 一开始开错题了,根本就不会做. 后来发现是每次 ...

随机推荐

  1. Django 后台定制自己的选择框删除函数

    from django.contrib import admin from .models import Article,Category from datetime import datetime ...

  2. 洛谷红名+AC150祭

    emmmm没什么想说的,随便放个图吧23333(逃~

  3. Asp.net Core导出Excel

    本篇文章是在MVC设计模式下,基于windows系统的Excel导出 1.前台的实现不用我多说了吧,加一个a标签链接地址跳到它所调用的方法里面,可以根据当前页面的查询条件去传值,从而查询出你想要的数据 ...

  4. 2019年最受欢迎IMX6系列开发板,资料全开源,助力产品研发不在话下

    迅为IMX6开发板: Android4.4系统  Linux + Qt5.7系统  Ubuntu12.04系统 部分真实案例:HMI:3D打印机:医疗设备:工控机:触控一体机:车载终端 板载:4G全网 ...

  5. 001 爬虫的基本概念以及urllib的request和parse

    1.http的请求方式: get请求 优点:比较便捷 缺点:不安全.长度有限制post请求 优点:比较安全.数据整体没有限制.可以上传文件putdelete(删除一些信息) 发送网络请求(可以带一定的 ...

  6. PostgreSQL学习笔记(二)-安装pgAdmin

    继上篇安装PostgreSQL后,我们需要安装一个PostgreSQL的图形化管理工具. pgadmin管理工具 创建Python的虚拟环境 cd /root/venv python -m venv ...

  7. c# 利用IEqualityComparer接口去除DataTable重复数据

    IEqualityComparer主要适用于定义方法以支持对象的相等比较.可以实现集合的自定义相等比较.即,您可以创建自己的相等定义,并指定此定义与接受 IEqualityComparer 接口的集合 ...

  8. python安装过程中的一些问题

    因为看到大神的教程是基于python V2.7,下载该版本且安装成功. 安装目录: https://www.python.org/download/releases/2.7/ 根据系统进行安装包下载 ...

  9. Linux中Too many open files

    1.ulimit –a open files一项就是默认的句柄数,最大为 65536 2.修改最大open files /etc/security/limits.conf文件中,加入以下配置: * s ...

  10. C++入门篇十三

    常对象: const  Person p1; 不可以调用普通成员函数,除非前面加了函数前面加了const可以调用常函数在对象之前加入const修饰 const Person p1; 常函数:void ...