传送门

人傻常数大.jpg

因为求逆的时候没清零结果调了几个小时……

前置芝士

多项式除法,多项式求逆

什么?你不会?左转你谷模板区,包教包会

题解

首先我们要知道一个结论$$f(x_0)\equiv f(x)\pmod{(x-x_0)}$$

其中\(x_0\)为一个常量,\(f(x_0)\)也为一个常量

证明如下,设\(f(x)=g(x)(x-x_0)+A\),也就是说\(A\)是\(f(x)\)对\((x-x_0)\)这个多项式取模之后的结果

因为\((x-x_0)\)的最高次项为\(1\),所以\(A\)的最高次项为\(0\),也就是说\(A\)是一个常数,即\(f(x)\equiv A\pmod{(x-x_0)}\)

我们把\(x_0\)代入上式,得\(f(x_0)=g(x_0)(x_0-x_0)+A\),同理可得\(f(x_0)\equiv A\pmod{(x-x_0)}\)

于是我们知道上式成立

这有毛用啊\(O(n\log n)\)多项式取模还没我暴力快

乍一看的确没啥卵用,但是考虑取模的过程是否能优化呢?

答案是可以的,我们考虑分治。设当前分治区间为\([l,r]\),令\(P_0=\prod_{i=l}^{mid}(x-x_i)\),\(P_1=\prod_{i=mid+1}^r (x-x_i)\),当前已经算出了\(A\equiv f(x)\pmod{\prod_{i=l}^r(x-x_i)}\),那么只要分别用\(A\)对\(P_0\)和\(P_1\)取模,然后继续递归下去就行了。取模之后\(A(x)\)的最高次项的次数变为原来的一半,问题规模也就变为原来的一半。继续递归下去就行了

时间复杂度为\(O(n\log^2n)\)

upd:改了改代码,常数应该会小一点,比方说分治到某个时候暴力秦九韶展开

  1. //minamoto
  2. #include<bits/stdc++.h>
  3. #define R register
  4. #define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
  5. #define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
  6. #define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
  7. using namespace std;
  8. char buf[1<<21],*p1=buf,*p2=buf;
  9. inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}
  10. int read(){
  11. R int res,f=1;R char ch;
  12. while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1);
  13. for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0');
  14. return res*f;
  15. }
  16. char sr[1<<21],z[20];int C=-1,Z=0;
  17. inline void Ot(){fwrite(sr,1,C+1,stdout),C=-1;}
  18. void print(R int x){
  19. if(C>1<<20)Ot();if(x<0)sr[++C]='-',x=-x;
  20. while(z[++Z]=x%10+48,x/=10);
  21. while(sr[++C]=z[Z],--Z);sr[++C]='\n';
  22. }
  23. const int N=(1<<17)+5,P=998244353;
  24. inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;}
  25. inline int dec(R int x,R int y){return x-y<0?x-y+P:x-y;}
  26. inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;}
  27. int ksm(R int x,R int y){
  28. R int res=1;
  29. for(;y;y>>=1,x=mul(x,x))(y&1)?res=mul(res,x):0;
  30. return res;
  31. }
  32. int r[19][N],w[2][N],lg[N],inv[19];
  33. void Pre(){
  34. fp(d,1,17){
  35. fp(i,1,(1<<d)-1)r[d][i]=(r[d][i>>1]>>1)|((i&1)<<(d-1));
  36. lg[1<<d]=d,inv[d]=ksm(1<<d,P-2);
  37. }
  38. for(R int t=(P-1)>>1,i=1,x,y;i<131072;i<<=1,t>>=1){
  39. x=ksm(3,t),y=ksm(332748118,t),w[0][i]=w[1][i]=1;
  40. fp(k,1,i-1)
  41. w[1][k+i]=mul(w[1][k+i-1],x),
  42. w[0][k+i]=mul(w[0][k+i-1],y);
  43. }
  44. }
  45. int lim,d,n,m;
  46. inline void init(R int len){lim=1,d=0;while(lim<len)lim<<=1,++d;}
  47. void NTT(int *A,int ty){
  48. fp(i,0,lim-1)if(i<r[d][i])swap(A[i],A[r[d][i]]);
  49. for(R int mid=1;mid<lim;mid<<=1)
  50. for(R int j=0,t;j<lim;j+=(mid<<1))
  51. fp(k,0,mid-1)
  52. A[j+k+mid]=dec(A[j+k],t=mul(w[ty][mid+k],A[j+k+mid])),
  53. A[j+k]=add(A[j+k],t);
  54. if(!ty)fp(i,0,lim-1)A[i]=mul(A[i],inv[d]);
  55. }
  56. void Inv(int *a,int *b,int len){
  57. if(len==1)return b[0]=ksm(a[0],P-2),void();
  58. Inv(a,b,len>>1),lim=(len<<1),d=lg[lim];
  59. static int A[N],B[N];
  60. fp(i,0,len-1)A[i]=a[i],B[i]=b[i];fp(i,len,lim-1)A[i]=B[i]=0;
  61. NTT(A,1),NTT(B,1);
  62. fp(i,0,lim-1)A[i]=mul(A[i],mul(B[i],B[i]));
  63. NTT(A,0);
  64. fp(i,0,len-1)b[i]=dec(add(b[i],b[i]),A[i]);
  65. fp(i,len,lim-1)b[i]=0;
  66. }
  67. struct node{
  68. node *lc,*rc;vector<int>vec;int deg;
  69. void Mod(const int *a,int *r,int n){
  70. static int A[N],B[N],D[N];
  71. int len=1;while(len<=n-deg)len<<=1;
  72. fp(i,0,n)A[i]=a[n-i];fp(i,0,deg)B[i]=vec[deg-i];
  73. fp(i,n-deg+1,len-1)B[i]=0;
  74. Inv(B,D,len);
  75. lim=(len<<1),d=lg[lim];
  76. fp(i,n-deg+1,lim-1)A[i]=D[i]=0;
  77. NTT(A,1),NTT(D,1);
  78. fp(i,0,lim-1)A[i]=mul(A[i],D[i]);
  79. NTT(A,0);
  80. reverse(A,A+n-deg+1);
  81. init(n+1);
  82. fp(i,n-deg+1,lim-1)A[i]=0;
  83. fp(i,0,deg)B[i]=vec[i];fp(i,deg+1,lim-1)B[i]=0;
  84. NTT(A,1),NTT(B,1);
  85. fp(i,0,lim-1)A[i]=mul(A[i],B[i]);
  86. NTT(A,0);
  87. fp(i,0,deg-1)r[i]=dec(a[i],A[i]);
  88. }
  89. void Mul(){
  90. static int A[N],B[N];deg=lc->deg+rc->deg,vec.resize(deg+1),init(deg+1);
  91. fp(i,0,lc->deg)A[i]=lc->vec[i];fp(i,lc->deg+1,lim-1)A[i]=0;
  92. fp(i,0,rc->deg)B[i]=rc->vec[i];fp(i,rc->deg+1,lim-1)B[i]=0;
  93. NTT(A,1),NTT(B,1);
  94. fp(i,0,lim-1)A[i]=mul(A[i],B[i]);
  95. NTT(A,0);
  96. fp(i,0,deg)vec[i]=A[i];
  97. }
  98. }pool[N],*rt;
  99. int A[N],a[N],tot;
  100. inline node* newnode(){return &pool[tot++];}
  101. void solve(node* &p,int l,int r){
  102. p=newnode();
  103. if(l==r)return p->deg=1,p->vec.resize(2),p->vec[0]=P-a[l],p->vec[1]=1,void();
  104. int mid=(l+r)>>1;
  105. solve(p->lc,l,mid),solve(p->rc,mid+1,r);
  106. p->Mul();
  107. }
  108. int b[25];
  109. void calc(node* p,int l,int r,const int *A){
  110. if(r-l<=512){
  111. fp(i,l,r){
  112. int x=a[i],c1,c2,c3,c4,now=A[r-l];
  113. b[0]=1;fp(j,1,16)b[j]=mul(b[j-1],x);
  114. for(R int j=r-l-1;j-15>=0;j-=16){
  115. c1=(1ll*now*b[16]+1ll*A[j]*b[15]+1ll*A[j-1]*b[14]+1ll*A[j-2]*b[13])%P,
  116. c2=(1ll*A[j-3]*b[12]+1ll*A[j-4]*b[11]+1ll*A[j-5]*b[10]+1ll*A[j-6]*b[9])%P,
  117. c3=(1ll*A[j-7]*b[8]+1ll*A[j-8]*b[7]+1ll*A[j-9]*b[6]+1ll*A[j-10]*b[5])%P,
  118. c4=(1ll*A[j-11]*b[4]+1ll*A[j-12]*b[3]+1ll*A[j-13]*b[2]+1ll*A[j-14]*b[1])%P,
  119. now=(0ll+c1+c2+c3+c4+A[j-15])%P;
  120. }
  121. fd(j,(r-l)%16-1,0)now=(1ll*now*x+A[j])%P;
  122. print(now);
  123. }
  124. return;
  125. }
  126. int mid=(l+r)>>1,b[p->deg+1];
  127. p->lc->Mod(A,b,p->deg-1),calc(p->lc,l,mid,b);
  128. p->rc->Mod(A,b,p->deg-1),calc(p->rc,mid+1,r,b);
  129. }
  130. int main(){
  131. // freopen("testdata.in","r",stdin);
  132. n=read(),m=read();if(!m)return 0;
  133. Pre();
  134. fp(i,0,n)A[i]=read();
  135. fp(i,1,m)a[i]=read();
  136. solve(rt,1,m);
  137. if(n>=m)rt->Mod(A,A,n);
  138. calc(rt,1,m,A);
  139. return Ot(),0;
  140. }

洛谷P5050 【模板】多项式多点求值的更多相关文章

  1. 洛谷P5282 【模板】快速阶乘算法(多项式多点求值+MTT)

    题面 传送门 前置芝士 \(MTT\),多项式多点求值 题解 这题法老当初好像讲过--而且他还说这种题目如果模数已经给定可以直接分段打表艹过去 以下是题解 我们设 \[F(x)=\prod_{i=0} ...

  2. luogu P5667 拉格朗日插值2 拉格朗日插值 多项式多点求值 NTT

    LINK:P5667 拉格朗日插值2 给出了n个连续的取值的自变量的点值 求 f(m+1),f(m+2),...f(m+n). 如果我们直接把f这个函数给插值出来就变成了了多项式多点求值 这个难度好像 ...

  3. 洛谷.3803.[模板]多项式乘法(FFT)

    题目链接:洛谷.LOJ. FFT相关:快速傅里叶变换(FFT)详解.FFT总结.从多项式乘法到快速傅里叶变换. 5.4 又看了一遍,这个也不错. 2019.3.7 叕看了一遍,推荐这个. #inclu ...

  4. 洛谷.3803.[模板]多项式乘法(NTT)

    题目链接:洛谷.LOJ. 为什么和那些差那么多啊.. 在这里记一下原根 Definition 阶 若\(a,p\)互质,且\(p>1\),我们称使\(a^n\equiv 1\ (mod\ p)\ ...

  5. 【洛谷P5050】 【模板】多项式多点求值

    code: #include <bits/stdc++.h> #define ll long long #define ull unsigned long long #define set ...

  6. 洛谷.4238.[模板]多项式求逆(NTT)

    题目链接 设多项式\(f(x)\)在模\(x^n\)下的逆元为\(g(x)\) \[f(x)g(x)\equiv 1\ (mod\ x^n)\] \[f(x)g(x)-1\equiv 0\ (mod\ ...

  7. 洛谷 P4238 [模板] 多项式求逆

    题目:https://www.luogu.org/problemnew/show/P4238 看博客:https://www.cnblogs.com/xiefengze1/p/9107752.html ...

  8. 洛谷.4512.[模板]多项式除法(NTT)

    题目链接 多项式除法 & 取模 很神奇,记录一下. 只是主要部分,更详细的和其它内容看这吧. 给定一个\(n\)次多项式\(A(x)\)和\(m\)次多项式\(D(x)\),求\(deg(Q) ...

  9. 洛谷 P4512 [模板] 多项式除法

    题目:https://www.luogu.org/problemnew/show/P4512 看博客:https://www.cnblogs.com/owenyu/p/6724611.html htt ...

随机推荐

  1. 论immutable不可变性

    什么叫immutable和mutable?简单来讲,一个immutable的对象一旦被创建好,它的状态将不会改变.反过来,如果一个类的实例是immutable的,那么我们把这个类也称作immutabl ...

  2. Unity发布安卓Splash Image适应手机、平板

    问题描述: Unity发布安卓应用时,遇到只能添加一张SplashImage而需要适配手机及平板不同分辨率处理 解决方法: http://blog.csdn.net/aiyan1111111/arti ...

  3. Codeforces 76D 位运算

    题意:给你两个数x 和 y, x = a + b, y = a XOR b,问有没有合法的a和b满足这个等式? 思路:有恒等式: a + b = ((a & b) << 1) + ...

  4. C语言中static修饰符的意义

    在C语言中,static通常有2种含义:1)定义变量的生命周期:2)定义变量或者函数的作用域. 变量的生命周期是指,相对于程序运行的进程生命周期,变量存在的时间段.变量的生命周期由变量的存储类型(位置 ...

  5. Python实现常见算法[3]——汉罗塔递归

    #!/usr/bin/python # define three list var. z1 = [1,2,3,4,5,6,7,"1st zhu"] z2 = ["2st ...

  6. SqlSugar Asp.Net 高性能ORM框架

    SqlSugar从去年到现在已经一年了,版本从1.0升到了现在的2.4.1 ,这是一个稳定版本 ,有数家公司已经项目上线,在这里我将SqlSugar的功能重新整理成一篇新的贴子,希望大家喜欢. 公司团 ...

  7. layer使用注意事项

    ajax一定要设置为异步

  8. IntelliJ IDEA开发golang环境配置

    IntelliJ IDEA开发golang环境配置 首先把GO安装好...(自行安装,附上一篇我之前写的MAC安装GO) 安装IntelliJ IDEA,下载地址: https://www.jetbr ...

  9. Spring Boot☞ 多数据源配置(二):Spring-data-jpa

    效果图: 代码区: package com.wls.integrateplugs.datasource; import org.springframework.beans.factory.annota ...

  10. Linux Mint 17使用配置

    => 调亮度: 方法一: 每次开机或注销登录之后需要重新设置 > /sys/class/backlight/intel_backlight/brightness #根据自己需要调值,如20 ...