生成函数好题!

搬一手铃悬的题解(侵删)

现在只需要考虑怎么求出g和逆变换即可,其实也就是对函数F(x)求F(x+1)和F(x-1)。

直接二项式定理展开发现是个卷积的形式,大力NTT即可。

#include<bits/stdc++.h>
#define N 440000
#define eps 1e-7
#define inf 1e9+7
#define db double
#define ll long long
#define ldb long double
using namespace std;
inline int read()
{
char ch=0;
int x=0,flag=1;
while(!isdigit(ch)){ch=getchar();if(ch=='-')flag=-1;}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*flag;
}
const int d=3,mo=998244353;
int ksm(int x,int k)
{
int ans=1;
while(k)
{
if(k&1)ans=1ll*ans*x%mo;
k>>=1;x=1ll*x*x%mo;
}
return ans;
}
int rev[N];
void ntt(int *f,int n,int flag)
{
for(int i=0;i<n;i++)
{
rev[i]=(rev[i>>1]>>1)+(i&1)*(n>>1);
if(i<rev[i])swap(f[i],f[rev[i]]);
}
for(int k=2,kk=1;k<=n;k<<=1,kk<<=1)
{
int wn=ksm(d,(mo-1)/k);
if(flag==-1)wn=ksm(wn,mo-2);
for(int i=0;i<n;i+=k)
for(int j=0,w=1;j<kk;j++,w=1ll*w*wn%mo)
{
int t=1ll*w*f[i+j+kk]%mo;
f[i+j+kk]=(f[i+j]-t)%mo;
f[i+j]=(f[i+j]+t)%mo;
}
}
if(flag==-1)
{
int k=ksm(n,mo-2);
for(int i=0;i<n;i++)f[i]=1ll*f[i]*k%mo;
}
}
int a[N],b[N];
void mul(int len)
{
ntt(a,len,+1);ntt(b,len,+1);
for(int i=0;i<len;i++)a[i]=1ll*a[i]*b[i]%mo;
ntt(a,len,-1);
}
int n,m,len,f[N],g[N],fac[N],vac[N];
int main()
{
n=read();ll t;cin>>t;m=(t%(mo-1));len=1;
while(len<2*(n+1))len<<=1;
for(int i=0;i<=n;i++)f[i]=read();
fac[0]=vac[0]=1;
for(int i=1;i<=len;i++)fac[i]=1ll*fac[i-1]*i%mo;
vac[len]=ksm(fac[len],mo-2);
for(int i=len-1;i>=1;i--)vac[i]=1ll*vac[i+1]*(i+1)%mo;
//get g(x)=f(x+1)
for(int i=0;i<=n;i++)a[i]=1ll*f[i]*fac[i]%mo,b[i]=vac[i];
for(int i=n+1;i<len;i++)a[i]=b[i]=0;
reverse(a,a+n+1);mul(len);
for(int i=0;i<=n;i++)g[i]=1ll*vac[i]*a[n-i]%mo;
//solve get g*(x)
for(int i=0;i<=n;i++)g[i]=1ll*ksm(ksm(i+1,m),mo-2)*g[i]%mo;
//get f*(x)=g(x-1)
for(int i=0;i<=n;i++)a[i]=1ll*g[i]*fac[i]%mo,b[i]=1ll*ksm(-1,i)*vac[i]%mo;
for(int i=n+1;i<len;i++)a[i]=b[i]=0;
reverse(a,a+n+1);mul(len);
for(int i=0;i<=n;i++)f[i]=1ll*vac[i]*a[n-i]%mo;
//print f(x)
for(int i=0;i<=n;i++)printf("%d ",(f[i]%mo+mo)%mo);
return 0;
}

CF923E Perpetual Subtraction的更多相关文章

  1. 【CF932E】Perpetual Subtraction(NTT,线性代数)

    [CF932E]Perpetual Subtraction(NTT,线性代数) 题面 洛谷 CF 题解 设\(f_{i,j}\)表示\(i\)轮之后这个数恰好为\(j\)的概率. 得到转移:\(\di ...

  2. Codeforces 947E Perpetual Subtraction (线性代数、矩阵对角化、DP)

    手动博客搬家: 本文发表于20181212 09:37:21, 原地址https://blog.csdn.net/suncongbo/article/details/84962727 呜啊怎么又是数学 ...

  3. Codeforces 923E - Perpetual Subtraction(微积分+生成函数+推式子+二项式反演+NTT)

    Codeforces 题目传送门 & 洛谷题目传送门 神仙题 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 首先考虑最朴素的 \(dp\),设 \(dp_{z,i}\) 表示经 ...

  4. Solution -「CF 923E」Perpetual Subtraction

    \(\mathcal{Description}\)   Link.   有一个整数 \(x\in[0,n]\),初始时以 \(p_i\) 的概率取值 \(i\).进行 \(m\) 轮变换,每次均匀随机 ...

  5. ZJOI2018游记Round1

    广告 ZJOI2018Round2游记 All Falls Down 非常感谢学弟学妹们捧场游记虽然这是一篇假游记 ZJOI Round1今天正式落下帷幕.在这过去的三天里遇到了很多朋友,见识了很多有 ...

  6. PAT 解题报告 1050. String Subtraction (20)

    1050. String Subtraction (20) Given two strings S1 and S2, S = S1 - S2 is defined to be the remainin ...

  7. [leetcode-592-Fraction Addition and Subtraction]

    Given a string representing an expression of fraction addition and subtraction, you need to return t ...

  8. [LeetCode] Fraction Addition and Subtraction 分数加减法

    Given a string representing an expression of fraction addition and subtraction, you need to return t ...

  9. [Swift]LeetCode592. 分数加减运算 | Fraction Addition and Subtraction

    Given a string representing an expression of fraction addition and subtraction, you need to return t ...

随机推荐

  1. Go 初体验 - 并发与锁.3 - 竞态

    竞态,就是多个协程同时访问临界区,由并发而产生的数据不同步的状态. 这个说的有点low,没办法,我就是这么表达的,官方的请度娘. 先上代码: 输出: 为何不是1000?就是因为竞态,发生竞态后,最终的 ...

  2. UEP-添加

    添加数据时候自动更新时间,注意添加数据时,要设置最大主键 时间的工具类:CommonUtil 人的工具类: ContextUtil ScmCompanyAccount scmCompanyAccoun ...

  3. Redis哨兵模式(sentinel)部署记录(主从复制、读写分离、主从切换)

    部署环境: CentOS7.5  192.168.94.11 (master) 192.168.94.22 (slave0) 192.168.94.33 (slave1) 192.168.94.44 ...

  4. Python批量合并处理B站视频

    最近想学习后端,又不想花钱,怎么办呢?于是在手机端B站(哔哩哔哩)上面找到了满意的免费视频教程,但是手机端看起来很不方便啊.于是,我通过在手机端缓存下来后,导入到了电脑端,但是我后面了发现两个问题: ...

  5. matlab 表示一个导数的函数

    对此式在matlab做表示:

  6. Tomcat笔记 #01# WEB应用管理工具简介

    索引 查看JVM以及SERVLET/接口的情况 动态管理WEB应用 Tomcat自带了一个基于网页的web应用管理工具,可以帮助我们监控&管理部署上去的WEB APP,特别方便!恰好之前碰到的 ...

  7. Java线程基础(二)

    今天上午考完了计算机二级,也算卸掉了一个大包袱吧,希望能过!(其实也就考着玩的,不来点考试就要发霉了) 好了,趁着难得的考后休息时间我就接着上一次没写完的继续更新吧. 上一篇文章——>Java核 ...

  8. Java线程基础(一)

    说在前面,经过一段学习过后,自己发觉线程在Java中占有举足轻重的地位,总觉得如此复杂的线程知识点一定要好好理清才好消化,因而有了这篇文章. 但因鄙人资历尚浅,如有遗漏错误之处还请广大网友不吝赐教. ...

  9. image management in kubernet

    Image How can I edit an existing docker image metadata? docker-copyedit Registry Disk kubevirtis a g ...

  10. SQL语句基本语法总结

    SQL语句基本语法 表的创建.修改.删除: 表中数据的增加.修改.删除.查看: SQL的 语法.数据类型.约束.删除表中数据(2种方法).语句查询的执行顺序: 查询 单列.多列.全部 子查询 单行子查 ...