Luogu 4721 【模板】分治 FFT

还不会这题的多项式求逆的算法。

发现每一项都是一个卷积的形式,那么我们可以使用$NTT$来加速,直接做是$O(n^2logn)$的,我们考虑如何加速转移。

可以采用$cdq$分治的思想,对于区间$[l, r]$中的数,先计算出$[l, mid]$中的数对$[mid + 1, r]$中的数的贡献,然后直接累加到右边去。

容易发现,这样子每一次需要用向量$[l,l + 1, l +  2, \dots, mid]$卷上$g$中$[1, 2, \dots, r - l]$。

时间复杂度$O(nlog^2n)$,感觉这东西跑得并不慢鸭。

Code:

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;const int N = 3e5 + 5;
const ll P = 998244353LL;int n, lim, pos[N];
ll f[N], g[N], a[N], b[N];template <typename T>
inline void read(T &X) {X = 0; char ch = 0; T op = 1;for (; ch > '9'|| ch < '0'; ch = getchar())if (ch == '-') op = -1;for (; ch >= '0' && ch <= '9'; ch = getchar())X = (X << 3) + (X << 1) + ch - 48;X *= op;
}template <typename T>
inline void swap(T &x, T &y) {T t = x; x = y; y = t;
}inline ll fpow(ll x, ll y) {ll res = 1LL;for (; y > 0; y >>= 1) {if (y & 1) res = res * x % P;x = x * x % P;}return res;
}inline void prework(int len) {int l = 0;for (lim = 1; lim <= len; lim <<= 1, ++l);for (int i = 0; i < lim; i++)pos[i] = (pos[i >> 1] >> 1) | ((i & 1) << (l - 1));
}inline void ntt(ll *c, int opt) {for (int i = 0; i < lim; i++)if (i < pos[i]) swap(c[i], c[pos[i]]);for (int i = 1; i < lim; i <<= 1) {ll wn = fpow(3, (P - 1) / (i << 1));if (opt == -1) wn = fpow(wn, P - 2);for (int len = i << 1, j = 0; j < lim; j += len) {ll w = 1;for (int k = 0; k < i; k++, w = w * wn % P) {ll x = c[j + k], y = c[j + k + i] * w % P;c[j + k] = (x + y) % P, c[j + k + i] =(x - y + P) % P;}}}if (opt == -1) {ll inv = fpow(lim, P - 2);for (int i = 0; i < lim; i++) c[i] = c[i] * inv % P;}
}void solve(int l, int r) {if (l == r) {a[l] = (a[l] + b[l]) % P;return;}int mid = ((l + r) >> 1);solve(l, mid);prework(r - l + 1);for (int i = 0; i < lim; i++) g[i] = f[i] = 0;for (int i = l; i <= mid; i++) f[i - l] = a[i];for (int i = 1; i <= r - l; i++) g[i - 1] = b[i];ntt(f, 1), ntt(g, 1);for (int i = 0; i < lim; i++) f[i] = f[i] * g[i] % P;ntt(f, -1);for (int i = mid + 1; i <= r; i++) a[i] = (a[i] + f[i - l - 1]) % P;solve(mid + 1, r);
}int main() {read(n); n--;for (int i = 1; i <= n; i++) read(b[i]);a[0] = 1;solve(1, n);for (int i = 0; i <= n; i++)printf("%lld%c", a[i], i == n ? '\n' : ' ');return 0;    
}
View Code

 

转载于:https://www.cnblogs.com/CzxingcHen/p/10197696.html

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

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

相关文章

根文件系统构建(Buildroot 方式)

目录何为buildroot&#xff1f;buildroot 简介buildroot 下载buildroot 构建根文件系统配置buildroot编译buildrootbuildroot 根文件系统测试buildroot 第三方软件和库的配置buildroot 下的busybox 配置busybox 配置busybox 中文字符的支持编译busybox根文件系统测试前面我们学…

Kernel panic -not syncing:Attempted to kill init!

DELL R210II安装RHEL后启动系统正常&#xff0c;关机再装磁盘阵列卡接好磁盘阵列&#xff0c;开机后CTRLR进入界面配置RAID5。启动后出现错误&#xff1a; Kernel panic -not syncing:Attempted to kill init! 解决过程&#xff1a;1、强行关机&#xff0c;再重启进入系统选择界…

根文件系统构建(Yocto方式)

本章有独立的教程&#xff0c;请参考《【正点原子】I.MX6U Yocto 构建根文件系统V1.0.pdf》&#xff01;

Jquery weui picker 支持label和value

万年没更新了. 最近用jquery weui. 在使用picker时需要一些问题. 就是让picker 显示label, 但是取值的时候取value用于存储. 官网例子如下 Jquery-weui 官网 : 点这里 写这篇文章时&#xff0c;用的版本是1.2.1版本 $("#picker-name").picker({title: "请选择您…

Linux—任务计划

1、周期性任务计划&#xff1a;管理系统的过程中&#xff0c;我们并不是每次都直接执行脚本或者命令&#xff0c;有时候需要让脚本、命令以及系统等在指定的时间按照我们的意愿执行我们设定好的脚本或者命令&#xff0c;这时候就需要使用到Linux任务计划的功能了&#xff0c;但…

计算机网络第1章(概述)

B站视频&#xff1a;计算机网络微课堂&#xff08;有字幕无背景音乐版&#xff09; 网址&#xff1a;https://www.bilibili.com/video/BV1c4411d7jb?p61 目录图示说明1.1、计算机网络在信息时代的作用我国互联网发展状况1.2、因特网概述1、网络、互连网&#xff08;互联网&…

django 中使用 channels 实现websocket

.... 转载于:https://www.cnblogs.com/zxmbky/p/10203009.html

HDU 1754线段树

第一个自己动手写的线段树&#xff0c;1Y还是有点小激动哈&#xff08;虽然是模版题&#xff09; 1 #include<cstdio>2 #include<cstring>3 #include<algorithm>4 using namespace std;5 const int SIZE200005;6 const int INF1000000;7 int maxv[SIZE<&l…

Linux命令(七)Linux用户管理和修改文件权限

1. 用户管理 1.1 创建用户/设置密码/删除用户 &#xff08;-m很重要&#xff0c;自动添加用户家目录&#xff09; 创建用户组dev, 给用户组dev新建xiaoqin用户&#xff0c;给新用户设置密码&#xff01; 1.2 查看用户信息 1.3 设置用户主组/附加组 &#xff08;常使用修改用户附…

元模型驱动

感谢徐昊给我们传授这么好的思想。虽然我现在还没有全部弄明白——————————在网上查到了元数据的解释是关于数据的数据。那么元模型就是关于模型的模型。抽象的过程就是形成元模型的过程。建模的过程就是抽象的过程。我们要形成object的模型&#xff0c;所以建立了的是…

2月末周全球域名解析商Top15:万网DNSPod份额均上涨

IDC评述网&#xff08;idcps.com&#xff09;03月11日报道&#xff1a;据域名统计机构WebHosting.info最新数据&#xff0c;在2月末周&#xff08;截至3月3日&#xff09;期间&#xff0c;我国仍是两名域名解析服务商跻身进全球前十五强。其中&#xff0c;中国万网排名第9&…

IHttpHandler与IHttpHandlerFactory的区别一例

在进行访问请求拦截处理时&#xff0c;处理类可以继承自IHttpHandler也可以继承自IHttpHandlerFactory&#xff0c;但他们之间有什么区别呢&#xff1f;近日开发一个需要拦截访问处理的项目&#xff0c;由于需要使用到SESSION&#xff0c;在这里发出了求助贴子&#xff08;积木…

ubuntu - 安装hive

粗略步骤&#xff1a; 详细参考&#xff1a;https://www.2cto.com/net/201804/735478.html 环境&#xff1a;ubunut jdk hadoop mysql 一、下载hive 二、解压&#xff08;到指定目录&#xff09; 三、配置 /etc/profile [环境变量] - 执行&#xff1a; source /etc/profil…

今天拿到了同事给我买的《java与模式》

自从有了设计模式的概念后一直想买本模式方面的书&#xff0c;所以买了《C#设计模式》&#xff0c;但是读了几个月的时间发现自己还只是看了简单工厂&#xff0c;工厂&#xff0c;抽象工厂和单件模式&#xff0c;已经被此书搞得晕头转向了&#xff08;书也烂&#xff0c;自己水…

计算机网络第4章(网络层)

B站视频&#xff1a;计算机网络微课堂&#xff08;有字幕无背景音乐版&#xff09; 网址&#xff1a;https://www.bilibili.com/video/BV1c4411d7jb?p61 目录4.1、网络层概述简介总结4.2、网络层提供的两种服务面向连接的虚电路服务无连接的数据报服务虚电路服务与数据报服务的…

Struts2中action接受参数方法

Struts2中Action接收参数的方法主要有以下三种&#xff1a;1.使用Action的属性接收参数&#xff1a; a.定义&#xff1a;在Action类中定义属性&#xff0c;创建get和set方法&#xff1b; b.接收&#xff1a;通过属性接收参数&#xff0c;如&#xff1a;userName&#xff…

Python __slots__ 作用

参考&#xff1a;https://blog.csdn.net/u010733398/article/details/52803643 https://blog.csdn.net/sxingming/article/details/52892640 转载于:https://www.cnblogs.com/xiaoshayu520ly/p/10207390.html

arm's multiply(arm的乘)之三

路转峰回发现编译器的优化比较复杂&#xff0c;在代码结构不同的情况下做的优化也完全不同。又做了一个实验&#xff0c;这次evc很笨没有做优化int array[1000];int res;for(int p10;p1<1000;p1)for(int p20;p2<1000;p2)for(int p30;p3<1000;p3)res array[1000]*105;…

CSS 实现背景半透明

IE过渡滤镜 CSS3 rgba 即可完美实现。 具体实现代码如下&#xff1a; .transparent {background:rgba(0, 0, 0, 0.3);filter:progid:DXImageTransform.Microsoft.Gradient(GradientType0, StartColorStr#4c000000, EndColorStr#4c000000); } :root .transparent {filter:progi…

计算机网络第5章(传输层)

B站视频&#xff1a;计算机网络微课堂&#xff08;有字幕无背景音乐版&#xff09; 网址&#xff1a;https://www.bilibili.com/video/BV1c4411d7jb?p61 目录5.1、运输层概述概念总结5.2、运输层端口号、复用与分用的概念为什么用端口号发送方的复用和接收方的分用TCP/IP体系的…