(CCPC 2020 网络选拔赛)HDU 6900 Residual Polynomial(分治 + NTT)

Residual Polynomial

写出所有的fi(x)f_i(x)fi(x)出来,fi,jf_{i, j}fi,j表示fi(x)f_i(x)fi(x)的第jjj项系数

{f1,0f1,1f1,2…f1,n−1f1,nf2,0f2,1f2,2…f2,n−1f2,nf3,0f3,1f3,2…f3,n−1f3,n⋮⋮⋮⋱⋮⋮fn−1,0fn−1,1fn−1,2…fn−1,n−1fn−1,nfn,0fn,1fn,2…fn,n−1fn,n}\left\{ \begin{matrix} f_{1, 0} & f_{1, 1} & f_{1, 2} & \dots & f_{1, n - 1} & f_{1, n}\\ f_{2, 0} & f_{2, 1} & f_{2, 2} & \dots & f_{2, n - 1} & f_{2, n}\\ f_{3, 0} & f_{3, 1} & f_{3, 2} & \dots & f_{3, n - 1} & f_{3, n}\\ \vdots & \vdots & \vdots & \ddots & \vdots & \vdots\\ f_{n - 1, 0} & f_{n - 1, 1} & f_{n - 1, 2} & \dots & f_{n - 1, n - 1} & f_{n - 1, n}\\ f_{n, 0} & f_{n, 1} & f_{n, 2} & \dots & f_{n, n - 1} & f_{n, n}\\ \end{matrix} \right\} f1,0f2,0f3,0fn1,0fn,0f1,1f2,1f3,1fn1,1fn,1f1,2f2,2f3,2fn1,2fn,2f1,n1f2,n1f3,n1fn1,n1fn,n1f1,nf2,nf3,nfn1,nfn,n
第一行是a0,a1,a2,…,an−1,ana_0, a_1, a_2, \dots, a_{n - 1}, a_{n}a0,a1,a2,,an1,an

考虑下一行如何从上一行转移,分两种:

一、fi,j×ci+1=fi+1,jf_{i, j} \times c_{i + 1} = f_{i + 1, j}fi,j×ci+1=fi+1,j

二、fi,j×j×bi+1=fi+1,j−1f_{i, j} \times j \times b_{i + 1} = f_{i + 1, j - 1}fi,j×j×bi+1=fi+1,j1

可以想想成每个点有一条向下的连边value=cvalue = cvalue=c,同时有一条向左下的连边value=b×jvalue = b \times jvalue=b×j

我们从f1,if_{1, i}f1,i出发,走向fn,jf_{n, j}fn,j,显然选择了i−ji- jijb×jb \times jb×j,所以还有n−1−(i−j)n - 1 - (i - j)n1(ij)ccc

先不考虑jjj,只算b,cb, cb,c的贡献,设选bbb的生成函数为F(x)F(x)F(x),对于某个点有F(x)=c+bxF(x) = c + bxF(x)=c+bx,因为上面的转移是具有区间性的。

所以F(l,r,x)=F(l,mid,x)∗F(mid,r,x)F(l, r, x) = F(l, mid, x) * F(mid, r, x)F(l,r,x)=F(l,mid,x)F(mid,r,x),也就是在左区间选几个,右区间选几个,然后组合一下,这个组合刚好满足卷积的形式。

这一步就可以递归O(nlog⁡nlog⁡n)O(n \log n \log n)O(nlognlogn)求解了。

接下来考虑jjj的影响:

fn,k=∑i−k=jF(j)f1,ii!k!f_{n, k} = \sum\limits_{i - k = j} F(j)f_{1, i} \frac{i!}{k!}fn,k=ik=jF(j)f1,ik!i!

这里可以简单理解一下,每次只能向左向下移动,k是我们最后在的列,只能从初始列i≥ki \geq kik的列转移过来。

(k!×fn,k)=∑i−k=jF(j)(f1,i×i!)\left(k! \times f_{n, k}\right) = \sum\limits_{i - k = j} F(j) \left(f_{1, i} \times i!\right)(k!×fn,k)=ik=jF(j)(f1,i×i!)

设H(k)=k!×fn,k,G(i)=f1,i×i!设H(k) = k! \times f_{n, k}, G(i) = f_{1, i} \times i !H(k)=k!×fn,k,G(i)=f1,i×i!考虑把G翻转,G(n−i)=G(i)考虑把G翻转,G(n - i) = G(i)GG(ni)=G(i)

H(k)=∑i+j=n−kF(j)G(i)H(k) = \sum_{i + j = n - k} F(j) G(i)H(k)=i+j=nkF(j)G(i),于是再做一次卷积即可求得答案。

#include <bits/stdc++.h>using namespace std;const int mod = 998244353, inv2 = mod + 1 >> 1;namespace Quadratic_residue {struct Complex {int r, i;Complex(int _r = 0, int _i = 0) : r(_r), i(_i) {}};int I2;Complex operator * (const Complex &a, Complex &b) {return Complex((1ll * a.r * b.r % mod  + 1ll * a.i * b.i % mod * I2 % mod) % mod, (1ll * a.r * b.i % mod + 1ll * a.i * b.r % mod) % mod);}Complex quick_pow(Complex a, int n) {Complex ans = Complex(1, 0);while (n) {if (n & 1) {ans = ans * a;}a = a * a;n >>= 1;}return ans;}int get_residue(int n) {mt19937 e(233);if (n == 0) {return 0;}if(quick_pow(n, (mod - 1) >> 1).r == mod - 1) {return -1;}uniform_int_distribution<int> r(0, mod - 1);int a = r(e);while(quick_pow((1ll * a * a % mod - n + mod) % mod, (mod - 1) >> 1).r == 1) {a = r(e);}I2 = (1ll * a * a % mod - n + mod) % mod;int x = quick_pow(Complex(a, 1), (mod + 1) >> 1).r, y = mod - x;if(x > y) swap(x, y);return x;}
}const int N = 1e6 + 10;int r[N], inv[N], b[N], c[N], d[N], e[N], t[N];int quick_pow(int a, int n) {int ans = 1;while (n) {if (n & 1) {ans = 1ll * a * ans % mod;}a = 1ll * a * a % mod;n >>= 1;}return ans;
}void get_r(int lim) {for (int i = 0; i < lim; i++) {r[i] = (i & 1) * (lim >> 1) + (r[i >> 1] >> 1);}
}void get_inv(int n) {inv[1] = 1;for (int i = 2; i <= n; i++) {inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod;}
}void NTT(int *f, int lim, int rev) {for (int i = 0; i < lim; i++) {if (i < r[i]) {swap(f[i], f[r[i]]);}}for (int mid = 1; mid < lim; mid <<= 1) {int wn = quick_pow(3, (mod - 1) / (mid << 1));for (int len = mid << 1, cur = 0; cur < lim; cur += len) {int w = 1;for (int k = 0; k < mid; k++, w = 1ll * w * wn % mod) {int x = f[cur + k], y = 1ll * w * f[cur + mid + k] % mod;f[cur + k] = (x + y) % mod, f[cur + mid + k] = (x - y + mod) % mod;}}}if (rev == -1) {int inv = quick_pow(lim, mod - 2);reverse(f + 1, f + lim);for (int i = 0; i < lim; i++) {f[i] = 1ll * f[i] * inv % mod;}}
}void polyinv(int *f, int *g, int n) {if (n == 1) {g[0] = quick_pow(f[0], mod - 2);return ;}polyinv(f, g, n + 1 >> 1);for (int i = 0; i < n; i++) {t[i] = f[i];}int lim = 1;while (lim < 2 * n) {lim <<= 1;}get_r(lim);NTT(t, lim, 1);NTT(g, lim, 1);for (int i = 0; i < lim; i++) {int cur = (2 - 1ll * g[i] * t[i] % mod + mod) % mod;g[i] = 1ll * g[i] * cur % mod;t[i] = 0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}
}void polysqrt(int *f, int *g, int n) {if (n == 1) {g[0] = Quadratic_residue::get_residue(f[0]);return ;}polysqrt(f, g, n + 1 >> 1);polyinv(g, b, n);int lim = 1;while (lim < 2 * n) {lim <<= 1;}get_r(lim);for (int i = 0; i < n; i++) {t[i] = f[i];}NTT(g, lim, 1);NTT(b, lim, 1);NTT(t, lim, 1);for (int i = 0; i < lim; i++) {g[i] = (1ll * inv2 * g[i] % mod + 1ll * inv2 * b[i] % mod * t[i] % mod) % mod;b[i] = t[i] = 0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}
}void derivative(int *a, int *b, int n) {for (int i = 0; i < n; i++) {b[i] = 1ll * a[i + 1] * (i + 1) % mod;}
}void integrate(int *a, int n) {for (int i = n - 1; i >= 1; i--) {a[i] = 1ll * a[i - 1] * inv[i] % mod;}a[0] = 0;
}void polyln(int *f, int *g, int n) {polyinv(f, b, n);derivative(f, g, n);int lim = 1;while (lim < 2 * n) {lim <<= 1;}get_r(lim);NTT(g, lim, 1);NTT(b, lim, 1);for (int i = 0; i < lim; i++) {g[i] = 1ll * g[i] * b[i] % mod;b[i] = 0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}integrate(g, n);
}void polyexp(int *f, int *g, int n) {if (n == 1) {g[0] = 1;return ;}polyexp(f, g, n + 1 >> 1);int lim = 1;while (lim < 2 * n) {lim <<= 1;}polyln(g, d, n);for (int i = 0; i < n; i++) {t[i] = (f[i] - d[i] + mod) % mod;}t[0] = (t[0] + 1) % mod;get_r(lim);NTT(g, lim, 1);NTT(t, lim, 1);for (int i = 0; i < lim; i++) {g[i] = 1ll * g[i] * t[i] % mod;t[i] = d[i] =  0;}NTT(g, lim, -1);for (int i = n; i < lim; i++) {g[i] = 0;}
}/*b存放多项式逆,c存放多项式开根,d存放多项式对数ln,e存放多项式指数exp,t作为中间转移数组,如果要用到polyln,得提前调用get_inv(n)先预先得到我们想要得到的逆元范围。
*/int B[N], C[N], f[N], fac[N], inv1[N], n;vector<int> F[N];int A1[N], B1[N];void init() {fac[0] = 1;for (int i = 1; i < N; i++) {fac[i] = 1ll * fac[i - 1] * i % mod;}inv1[N - 1] = quick_pow(fac[N - 1], mod - 2);for (int i = N - 2; i >= 0; i--) {inv1[i] = 1ll * inv1[i + 1] * (i + 1) % mod;}
}void merge(int rt) {int ls = rt << 1, rs = rt << 1 | 1;int n = F[ls].size(), m = F[rs].size();F[rt].resize(n + m);for (int i = 0; i < n + m; i++) {F[rt][i] = 0;}if (n <= 50 && m <= 50) {for (int i = 0; i < n; i++) {for (int j = 0; j < m; j++) {F[rt][i + j] = (F[rt][i + j] + 1ll * F[ls][i] * F[rs][j] % mod) % mod;}}return ;}for (int i = 0; i < n; i++) {A1[i] = F[ls][i];}for (int i = 0; i < m; i++) {B1[i] = F[rs][i];}int lim = 1;while (lim < n + m) {lim <<= 1;}get_r(lim);NTT(A1, lim, 1);NTT(B1, lim, 1);for (int i = 0; i < lim; i++) {A1[i] = 1ll * A1[i] * B1[i] % mod;}NTT(A1, lim, -1);for (int i = 0; i < n + m; i++) {F[rt][i] = A1[i];}for (int i = 0; i < lim; i++) {A1[i] = B1[i] = 0;}
}void divide(int rt, int l, int r) {if (l == r) {F[rt].push_back(C[l]);F[rt].push_back(B[l]);return ;}int mid = l + r >> 1;divide(rt << 1, l, mid), divide(rt << 1 | 1, mid + 1, r);merge(rt);F[rt << 1].clear(), F[rt << 1 | 1].clear();
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);int T;init();scanf("%d", &T);while (T--) {scanf("%d", &n);for (int i = 0; i <= n; i++) {scanf("%d", &f[i]);f[i] = 1ll * f[i] * fac[i] % mod;}for (int i = 0; i <= n - 2; i++) {scanf("%d", &B[i]);}for (int i = 0; i <= n - 2; i++) {scanf("%d", &C[i]);}divide(1, 0, n - 2);int lim = 1;while (lim <= 2 * n) {lim <<= 1;}get_r(lim);for (int i = 0; i <= n; i++) {A1[i] = f[n - i];}for (int i = 0; i < F[1].size(); i++) {B1[i] = F[1][i];}F[1].clear();NTT(A1, lim, 1);NTT(B1, lim, 1);for (int i = 0; i < lim; i++) {A1[i] = 1ll * A1[i] * B1[i] % mod;}NTT(A1, lim, -1);for (int i = 0; i <= n; i++) {printf("%lld%c", 1ll * inv1[i] * A1[n - i] % mod, i == n ? '\n' : ' ');}for (int i = 0; i < lim; i++) {A1[i] = B1[i] = 0;}}return 0;
}

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

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

相关文章

使用Elastic APM监控你的.NET Core应用

前言在应用实际的运维过程中&#xff0c;我们需要更多的日志和监控来让我们对自己的应用程序的运行状况有一个全方位的了解。然而对于大部分开发者而言&#xff0c;平时大家所关注的更多的是如何更优雅的实现业务&#xff0c;或者是如何让应用的响应速度更快等等与编码相关的技…

HDU 6755 Fibonacci Sum(二次剩余 + 二项式展开)

Fibonacci Sum 斐波那契通项有an15((152)n−(1−52)n)(15)k∑i0n((152)ic−(1−52)ic)kA152,B1−52(15)k∑i0n∑j0k(−1)k−jCkjAicjBic(k−j)(15)k∑j0k(−1)k−jCkj(∑i0nAcjiBc(k−j)i)斐波那契通项有a_n \frac{1}{\sqrt 5}\left((\frac{1 \sqrt 5}{2}) ^ n - (\frac{1 - …

Java 8 新的时间处理API

一&#xff1a;时间日期API的演进&#xff0c;及存在的问题 JDK 1.0 时期&#xff1a; 对于日期和时间的支持只能依赖于java.util.Date类。它的最小精度是毫秒起始年份为1900年&#xff0c;起始月份为0。20180822表示为new Date (118,7,22)返回值使用JVM默认时区&#xff1a;…

ASP.NET Core结合Nacos来完成配置管理和服务发现

前言今年4月份的时候&#xff0c;和平台组的同事一起调研了一下Nacos&#xff0c;也就在那个时候写了.net core版本的非官方版的SDK。虽然公司内部由于某些原因最后没有真正的用起来&#xff0c;但很多人还是挺看好的。在和镇汐大大沟通后&#xff0c;决定写一篇博客简单介绍一…

2020 ICPC 济南 F. Gcd Product

Gcd Product Cm∑i1mAgcd⁡(i,m)Bgcd⁡(k1−i,m)∑d1∣mAd1∑d2∣mBd2∑i1m([gcd⁡(id1,md1)1][d1∣i])([gcd⁡(m1−id2,md2)1][d2∣m1−i])∑d1∣mAd1∑d2∣mBd2∑k1∣md1μ(k1)∑k2∣md2μ(k2)∑i1m([d1∣i][k1∣id1])([d2∣m1−i][k2∣m1−id2])T1d1k1,T2d2k2∑T1∣m∑d1∣T…

Java 时间处理

时区、冬令时和夏令时、时间戳 时间戳 距离一个标准参照时间经过的秒数&#xff08;毫秒数&#xff09; 有两个常用参照时间&#xff1a; 1970-01-01 00:00:00 应用最广泛的时间戳参照点2001-01-01 00:00:00 常被苹果系统使用 注意&#xff1a;以上时间节点皆采用UTC的标准时…

试试这个Excel知识测验,得分超过80分算你赢

大家可能都知道&#xff0c;全世界使用Excel的用户超过了10亿。Excel的知识真所谓是博大精深&#xff0c;并且还很有趣味。我最近编写了一个Excel小工具&#xff0c;可以让大家可以在Excel里面进行各种知识小测验&#xff0c;并且与全世界的高手一比高低。这个小工具&#xff0…

SimpleDateFormat与线程安全

SimpleDateFormat不是线程安全的。 SimpleDateFormat(下面简称sdf)类内部有一个Calendar对象引用&#xff0c;它用来储存和这个sdf相关的日期信息&#xff0c;例如sdf.parse(dateStr)&#xff0c;sdf.format(date) 诸如此类的方法参数传入的日期相关String, Date等等&#xff…

几道偏序问题(数据结构)

P3157 [CQOI2011]动态逆序对 #include <bits/stdc.h>using namespace std;typedef long long ll;const int N 1e5 10;int root[N], ls[N << 8], rs[N << 8], sum[N << 8], cnt;int n, m, pos[N];inline int lowbit(int x) {return x & (-x); }v…

自学架构设计?帮你总结了 4 个方法

从编程思维到架构思维的升级&#xff0c;是工作 3、5 年的程序员遇到的第一个槛&#xff0c;特别是当你准备晋升考核时。我有个哥们&#xff0c;技术和业务都很不错&#xff0c;腾讯 T2.3 升 T3.1&#xff0c;就卡在了架构设计这部分。架构这个事儿&#xff0c;不像算法和代码&…

如何在东八区的计算机上获取美国时间

既可以用旧API&#xff08;JDK8之前&#xff09;&#xff0c;也可以使用新API。以下用旧API为例&#xff1a; 在Java语言中&#xff0c;可以通过java.util.Calendar类取得一个本地时间或者指定时区的时间实例&#xff0c;如下&#xff1a; // 取得本地时间&#xff1a; Calen…

Wannafly挑战赛24 无限手套(生成函数)

无限手套 每种宝石的生成函数为∑n≥0xn(ain2bin1)对其进行化简∑n≥0xn∑n≥0binxn∑n≥0ain2xn11−xbix(1−x)2aix(1x)(1−x)3最后答案∏i1m((ai−bi1)x2(aibi−2)x1)(1−x)3m每种宝石的生成函数为\sum_{n \geq 0} x ^ n(a_i n ^ 2 b_i n 1)\\ 对其进行化简\sum_{n \geq 0}…

.NET Core 3.0及ASP.NET Core 3.0 前瞻

前几天微软发布了 .NET Core 3.0 Preview 9 &#xff0c;这是.NET Core 3.0 最后一个预览版。[翻译] .NET Core 3.0 Preview 9 发布.NET Core 3.0 正式发布将在.NET Conf 上发布&#xff0c;.NET Conf 时间是9月23日至25日。Visual Studio 2019 16.3预览版3和Visual Studio for…

有了Unicode为啥还需要UTF-8

有了Unicode为啥还需要UTF-8 要回答这个问题&#xff0c;需要吃透“编码”的概念&#xff0c;刚好看到大神阮一峰写的文章&#xff1a;字符编码笔记&#xff1a;ASCII&#xff0c;Unicode 和 UTF-8 抄录如下&#xff0c;便于查找&#xff1a; 一、ASCII 码 我们知道&#xf…

2019ICPC西安邀请赛 E. Tree(树剖 + 线段树)

Tree 给定一棵树&#xff0c;节点有点权&#xff0c;然后有三种操作&#xff1a; 一、修改1−>s1->s1−>s的路径上的点权与ttt进行按位或。 二、修改1−>s1->s1−>s的路径上的点权与ttt进行按位与。 三、查询1−>s1->s1−>s的路径上的点权异或和…

升职却不加薪,为什么我还觉得老板说的挺有道理

前几天晚上&#xff0c;DevOps交流群里&#xff0c;有人抛出这样一个观点&#xff1a;如果有人来找我加薪&#xff0c;我一定告诉他我要给他升职&#xff0c;因为升职是免费的&#xff0c;加薪可是真的要花钱。但是我也会许诺他&#xff0c;如果他能把那个团队搞好&#xff0c;…

GBK、GB2312、GB18030

GBK、GB2312都是编码方式&#xff0c;都是双字节字符集&#xff0c;GB18030编码是一二四字节变长编码 GB&#xff1a;国标&#xff0c;国家标准 BIG5&#xff1a; 支持繁体中文 GB2312&#xff1a; 支持简体中文&#xff08;GBK的子集&#xff09; 共收录 6763 个汉字 GBK&…

2019 ICPC 南京 F. Paper Grading(字典树dfs序上树套树)

Paper Grading 题意&#xff1a;给定nnn个字符串&#xff0c;有两种操作&#xff1a; 一、给定i,ji, ji,j&#xff0c;交换第iii个跟第jjj个字符串。 二、给定 str &#xff0c;k,l,rk, l, rk,l,r&#xff0c;问你在区间[l,r][l, r][l,r]中的字符&#xff0c;与 str 至少有k…

细节之中自有天地,整洁成就卓越代码

溪源 | 长沙.NET技术社区开篇我们总是很容易就能写出满足某个特定功能的代码&#xff0c;却很难写出优雅代码。又最欣赏那些优雅的代码&#xff0c;因为优雅代码更能体现一个开发者的积累。就像写一篇散文&#xff0c;有的就像初学者不得其门而入&#xff0c;遣词造句都非常困难…

#2693. jzptab

jzptab ∑i1n∑j1mlcm(i,j)∑i1n∑j1mijgcd⁡(i,j)∑d1nd∑i1nd∑j1mdij[gcd⁡(i,j)1]∑d1nd∑k1ndk2μ(k)∑i1nkdi∑j1mkdjTkd,f(n)∑i1ni∑T1nf(nT)f(mT)(T∑k∣Tμ(k)k)设g(n)n∑d∣nμ(d)d先令g(n)g(n)ng(1)1,g(p)μ(1)μ(p)p1−p,g(pk,k≥2)1−p同时是积性函数&#xff0c…