原文链接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. LeetCode-11. 盛最多水的容器

    给定 n 个非负整数 a1,a2,...,an,每个数代表坐标中的一个点 (i, ai) .在坐标内画 n 条垂直线,垂直线 i 的两个端点分别为 (i, ai) 和 (i, 0).找出其中的两条线, ...

  2. MySql分页查询慢|这里告诉你答案

    一.背景 我们在开发的过程中使用分页是不可避免的,通常情况下我们的做法是使用limit加偏移量:select * from table where column=xxx order by xxx li ...

  3. 爬虫框架Scrapy 之(二) --- scrapy文件

    框架简介 核心部分: 引擎.下载器.调度器 自定义部分: spider(自己建的爬虫文件).管道(pipelines.py) 目录结构 firstSpider firstSpider spiders ...

  4. mp的猜猜看

    ~~~~|yjb1072452141---dc9339b4c33103abc4919375203e7a24|A1482583628---0142e0b6090b9b2838328445a79cd1b8 ...

  5. awk删除重复文件

    #!/bin/bash #查找并删除重复文件,每个文件只保留1份 ls -LS --time-style=long-iso | awk 'BEGIN { getline; getline; name1 ...

  6. 树莓派设置固定IP地址

    vi /etc/dhcpcd.conf # 使用 vi 编辑文件,增加下列配置项 # 指定接口 eth0 interface eth0 # 指定静态IP,/24表示子网掩码为 255.255.255. ...

  7. Aras 发布Web Services

    https://blog.csdn.net/plm888/article/details/10890173

  8. How to delete VSTS Project

    Buiding is so easy , Where is deleting ? C:\Program Files (x86)\Microsoft Visual Studio 12.0\Common7 ...

  9. jquery实现点击页面空白处,弹框消失

    要求:点击1,弹框2显示,点击空白处,弹框2消失 $("#AddDevices"):按钮1 $(".addDeviceBox")弹框2 //点击添加设备弹框 $ ...

  10. python numpy 间的的数据变算公式

    import numpy as np a = np.arange(100) print(np.sum(a))#求和 print(np.mean(a))#平均值 print(np.max(a))#最大值 ...