多项式除法，多项式取模

最新推荐文章于 2022-04-10 13:07:29 发布
最新推荐文章于 2022-04-10 13:07:29 发布
阅读量1.3k 收藏
点赞数 1
分类专栏： 多项式及生成函数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45483201/article/details/113685002
版权

多项式除法

给定一个 n n n次多项式 F ( x ) F(x) F(x)和 m m m次多项式 G ( x ) G(x) G(x)，要求 R ( x ) , Q ( x ) R(x), Q(x) R(x),Q(x)，满足 F ( x ) = R ( x ) G ( x ) + Q ( x ) F(x) = R(x)G(x) + Q(x) F(x)=R(x)G(x)+Q(x)。

R ( x ) R(x) R(x)是一个 n − m n - m n−m阶多项式， Q ( x ) Q(x) Q(x)是一个小于 m m m阶的多项式。
有 F ( x ) ≡ R ( x ) G ( x ) + Q ( x ) ( m o d x n + 1 ) F ( 1 x ) ≡ R ( 1 x ) G ( 1 x ) + Q ( 1 x ) ( m o d x n + 1 ) 同 时 乘 上 一 个 x n , F r e v ( x ) ≡ ( x m R r e v ( x ) ) ( x n − m G r e v ( x ) ) + x n − d e g Q Q r e v ( x ) ( m o d x n + 1 ) F r e v ( x ) ≡ R r e v ( x ) G r e v ( x ) + Q r e v ( x ) x n − d e g Q ( m o d x n + 1 ) 有 d e g Q < m , n − d e g Q > = n − m + 1 , 所 以 有 F r e v ( x ) ≡ R r e v ( x ) G r e v ( x ) ( m o d x n − m + 1 ) 有F(x) \equiv R(x) G(x) + Q(x) \pmod{x ^ {n + 1}}\\ F(\frac{1}{x}) \equiv R(\frac{1}{x})G(\frac{1}{x}) + Q(\frac{1}{x}) \pmod {x ^{n + 1}}\\ 同时乘上一个x ^ n, F^{rev}(x) \equiv \left(x ^m R ^{rev}(x)\right) \left(x ^{n - m}G ^{rev}(x)\right) + x ^{n - deg_Q} Q ^{rev} (x) \pmod{x ^{n + 1}} \\ F^{rev}(x) \equiv R^{rev}(x) G^{rev}(x) + Q^{rev}(x) x ^{n - deg_Q}\pmod{x ^{n + 1}}\\ 有deg_Q < m, n - deg_Q >= n - m + 1,所以有F^{rev}(x) \equiv R^{rev}(x) G^{rev}(x) \pmod{x ^{n - m + 1}}\\ 有F(x)≡R(x)G(x)+Q(x)(modxn+1)F(x1​)≡R(x1​)G(x1​)+Q(x1​)(modxn+1)同时乘上一个xn,Frev(x)≡(xmRrev(x))(xn−mGrev(x))+xn−degQ​Qrev(x)(modxn+1)Frev(x)≡Rrev(x)Grev(x)+Qrev(x)xn−degQ​(modxn+1)有degQ​<m,n−degQ​>=n−m+1,所以有Frev(x)≡Rrev(x)Grev(x)(modxn−m+1)
只要多项式求逆，即可得到 R ( x ) R(x) R(x)，然后代入原式求得 Q ( x ) Q(x) Q(x)。

#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 = 6e5 + 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作为中间转移数组,
  如果要用到polyinv，得提前调用get_inv(n)先预先得到我们想要得到的逆元范围。
*/

int f[N], fr[N], g[N], gr[N], rr[N], n, m;

int main() {
  // freopen("in.txt", "r", stdin);
  // freopen("out.txt", "w", stdout);
  // ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
  scanf("%d %d", &n, &m);
  for (int i = 0; i <= n; i++) {
    scanf("%d", &f[i]);
    fr[n - i] = f[i];
  }
  for (int i = 0; i <= m; i++) {
    scanf("%d", &g[i]);
    gr[m - i] = g[i];
  }
  for (int i = n - m + 1; i <= n; i++) {
    fr[i] = gr[i] = 0;
  }
  polyinv(gr, b, n - m + 1);
  for (int i = 0; i < n - m + 1; i++) {
    gr[i] = b[i];
    b[i] = 0;
  }
  int lim = 1;
  while (lim < 2 * (n - m + 1)) {
    lim <<= 1;
  }
  get_r(lim);
  NTT(fr, lim, 1);
  NTT(gr, lim, 1);
  for (int i = 0; i < lim; i++) {
    fr[i] = 1ll * fr[i] * gr[i] % mod;
  }
  NTT(fr, lim, -1);
  for (int i = 0; i <= n - m; i++) {
    rr[i] = fr[n - m - i];
  }
  for (int i = 0; rr[i]; i++) {
    printf("%d ", rr[i]);
  }
  puts("");
  lim = 1;
  while (lim <= 2 * n) {
    lim <<= 1;
  }
  get_r(lim);
  NTT(rr, lim, 1);
  NTT(g, lim, 1);
  for (int i = 0; i < lim; i++) {
    g[i] = 1ll * g[i] * rr[i] % mod;
  }
  NTT(g, lim, -1);
  for (int i = 0; i < m; i++) {
    f[i] = (f[i] - g[i] + mod) % mod;
  }
  for (int i = 0; i < m; i++) {
    printf("%d ", f[i]);
  }
  puts("");
  return 0;
}
_lifehappy_
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
多项式mod的运算方法
yaongtime的专栏
12-08 2万+
伽罗华域(Galois Field，GF)中关于多项式的mod运算过程，也不知道这样称合适不。这里主要是以构造GF(23)域为例，说明mod的运算过程， 假设本原多项式为p(x)=x3+x+1，α定义为 p(x)= 0的根，即 α3＋α＋1 = 0， GF(23)中的元素可计算如下： 0 mod(α3＋α＋1) = 0 α0 mod(α3
多项式mod的运算（例题）
OOOOOOAAAAA的博客
07-19 3702
https://blog.csdn.net/yaongtime/article/details/17200401 根据上述作者方法此处作一例题，题目对应https://www.bilibili.com/video/BV1zt411S7r8 信息安全数学基础——哈尔滨工业大学P27题目
多项式求逆，多项式取模，多项式开方 学习笔记
热门推荐
KsCla
02-23 1万+
前言 还记得上个学期tututu跟我提过多项式的很多操作，还有一些优化常数的奇技淫巧，然而那个时候我一脸懵逼。最近几天无所事事，去洛谷做比赛又整天被吊着打，闲暇之余就想着学一下多项式的几个基本操作。 其实一开始我是想学CZT的，根据myy的论文它能把BZOJ3992那题优化到O(mlog(m)+mlog(n))O(mlog⁡(m)+mlog⁡(n))O(m\log(m)+m\log(n))。然...
多项式取模
gigo_64 sink in my world
10-08 2078
又是一个毒瘤操作。 首先我们可以干这么一件很有趣的事情。 把一个多项式的所有系数reverse了。 然后再把代进去，然后再乘上一个，又变回原来的多项式了。 数学讲就是 所以它转回去了。 切入正题：给定n次多项式A，m次多项式B，求的一个D和R，使得 ，且R的最高次小于m。dzyo说可以证明DR唯一。 所以我们可以这么干。 两边同乘x^n 所以...
常系数齐次递推 / 多项式取模
一位蒟蒻的小博客
05-01 612
例题 给定G(次数界3w左右)，求F(x)=11−G(x)F(x)=\frac 1 {1-G(x)}F(x)=1−G(x)1​的第n项，保证G(x)[x0]=0G(x)[x^0]=0G(x)[x0]=0。 n在1e9以上。 有毛病 毒瘤 毒瘤 毒瘤 式子就是F(x)=1+F(x)G(x)F(x)=1+F(x)G(x)F(x)=1+F(x)G(x) 这是什么意思呢 其实就相当于一个齐次线性递推。 f...
2.带余除法的c语言程序
06-11
在模运算中，我们通常对一个数进行除法后取模，以获取该数在模意义下的等价类。例如，7 mod 3 = 1，意味着7在模3的情况下与1相同。模运算在密码学、数据编码和计算机算法设计中都有重要应用。 总之，带余除法是...
CRC16校验模二除法计算器
10-16
模二除法是实现CRC校验的基础算法，它不同于我们平常所学的十进制除法，因为在这个系统中，除数和被除数都是二进制数，并且遵循模二运算规则，即两个数相加如果结果超过1，则结果为结果对2取模。 CRC16校验指的是...
组合数学- 组合数取模- 卢卡斯定理与扩展卢卡斯定理.rar
09-16
在计算上，卢卡斯定理和扩展卢卡斯定理的优势在于它们避免了直接计算大数的除法和余数运算，这对于计算机程序来说是非常有益的。在实际应用中，如在密码学中的快速幂运算、计算有限域上的多项式系数等，这两个定理都...
CRC校验.pdf
03-11
发送方首先将信息字段（K位）与生成多项式相乘，然后取模（即进行多项式除法），余数即为R位的校验字段，最终形成N=K+R位的CRC码字。 在CRC码集选择时，需要确保码字与一个R次生成多项式g(x)的关系满足V(x)=A(x)g(x...
BCH.rar_BCH码_BCH码译码_BCH编译码_bch译码_单片机 编译码
09-23
BCH码的核心是基于伽罗华域（Galois Field）的数学理论，特别是有限域上的多项式除法和剩余类环。它定义了一个特定的生成多项式，所有合法的码字都是这个生成多项式的倍数。 **BCH码的构造** 1. **生成多项式选择*...
多项式相乘求模运算，模多项式可以改变
12-11
实现多项式相乘求模运算，并且模多项式是可以改变的，方便，简单
多项式模2运算及求逆元
weixin_34152820的博客
10-13 2934
在GF(28)域上，多项式相加相减结果相同，均为异或操作 x3+x2+1 对应的二进制为 1101 用整数表示就是 13 x2+x+1 对应的二进制为 0111　　用整数表示就是 7 x3+x2+1 + x2+x+1 =x3+2x2+x+2 = x3+x 等同于1101⊕0111 = 1010 用整数表示就是 10 既是 1...
【模板】多项式除法
远行客
07-29 5946
传送门:洛谷：【模板】多项式除法 题意 给定一个 nnn 次多项式 F(x)F(x)F(x) 和一个 mmm 次多项式 G(x)G(x)G(x) ，请求出多项式 Q(x)Q(x)Q(x) , R(x)R(x)R(x) ，满足以下条件： Q(x)Q(x)Q(x) 次数为 n−mn−mn-m， R(x)R(x)R(x) 次数小于 mmm F(x)=Q(x)×G(x)+R(x)F(x)=...
多项式的求逆、取模和多点求值学习小记
Felix-Lee的博客
01-19 1562
最近学习了多项式的求逆、取模和多点求值，这些方法能够解决很多多项式问题。 这三个操作是环环相扣的，很有趣，学完后不妨记录一下。 多项式求逆 给出一个次数界为 nnn 的多项式 A(x)A(x)A(x) ，需要求 B(x)B(x)B(x) 满足： A(x)B(x)≡1(mod&nbsp;xn)A(x)B(x)\equiv1(mod\ x^n)A(x)B(x)≡1(mod&nbsp;x...
多项式求逆与多项式除法/取模
Never give in.
10-17 6108
多项式求逆 多项式求逆是多项式模块中的一个重要操作（“操作”这个词看出如今多项式题是多么的工业化，犹如毒瘤8操作LCT）,在做生成函数/多项式除法、多项式取模/多项式多点求值中均有应用 对于一个n次多项式F(x)F(x)F(x)，我们希望求出一个m-1次多项式G(x)G(x)G(x)，满足F(x)G(x)≡1(modxm)F(x)G(x)\equiv 1\pmod {x^m}F(x)G(x)≡1(...
多项式除法
qq_60899598的博客
04-10 1277
定义： 对于一个给定的多项式： 其中： A(X):被除数多项式，次数界为n B(X):除数多项式，次数界为m(m<n) D(X):商多项式，次数界为n-m+1 R(X):余数多项式，次数界小于m(余数要小于除数，对于多项式来说就是余数的次数界要小于除数的次数界) 逆元 A(X)与互为逆元 mod :模就是去除掉多项式里n次的项 取模 了解了上述的知识后，我们就可以开始对多项是取模了 我们要求的有两个：商多项式(D(X)),余数多项式(R(X)) 我们先考虑求商，.
aes加密算法中的多项式求模过程
weixin_38161055的博客
06-09 3319
在写AES实验报告时，对多项式求模那一预备知识不明白，上网看到了一个文档就明白了，下面是截图。查考文档地址：https://wenku.baidu.com/view/cd0a0146b307e87101f6962a.html...
多项式基础II 多项式求逆 多项式取模 多项式开方 多项式牛顿迭代 多项式ln 多项式exp
myjs999的博客
12-05 514
用牛顿迭代。
matlab计算多项式除法
最新发布
06-12
在Matlab中，可以使用`deconv`函数来进行多项式除法的计算。 假设要计算多项式 A(x) 除以多项式 B(x)，可以将它们表示为向量的形式，其中向量的每个元素表示一个项的系数，例如： ``` A(x) = 3x^3 + 2x^2 + x + 1 B(x) = x^2 + 2x + 1 ``` 可以表示为向量： ``` A = [3 2 1 1] B = [1 2 1] ``` 然后，使用`deconv`函数进行多项式除法的计算： ``` [Q, R] = deconv(A, B) ``` 其中，`Q`表示商式，`R`表示余式。在上面的例子中，`Q`为`[3 0 -5]`，`R`为`[8 -13]`。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值