UOJ #34. 多项式乘法
来源:互联网 发布:node.js高级编程 微盘 编辑:程序博客网 时间:2024/06/06 13:42
#34. 多项式乘法
统计这是一道模板题。
给你两个多项式,请输出乘起来后的多项式。
输入格式
第一行两个整数
第二行
第三行
输出格式
一行
样例一
input
1 21 21 2 1
output
1 4 5 2
explanation
限制与约定
时间限制:
空间限制:
下载
样例数据下载
【分析】
NTT模板,silvernebula写的果然和FFT几乎一样蛤蛤
【代码】
#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<cmath>#define pi acos(-1)#define ll long long#define fo(i,j,k) for(i=j;i<=k;i++)using namespace std;const int mxn=400005;const int mod=479*(1<<21)+1;int n,m,N,M,L;int a[mxn],b[mxn],R[mxn];inline int ksm(int a,int k){int res=1;while(k){if(k&1) res=((ll)res*a)%mod;a=((ll)a*a)%mod;k>>=1;}return res;}inline void NTT(int *a,int f){int i,j,k;fo(i,0,n-1) if(i<R[i]) swap(a[i],a[R[i]]);for(i=1;i<n;i<<=1){int wn=ksm(3,(mod-1)/(i<<1));for(j=0;j<n;j+=(i<<1)){int w=1;for(k=0;k<i;k++,w=((ll)w*wn)%mod){int x=a[j+k],y=((ll)w*a[i+j+k])%mod;a[j+k]=(x+y)%mod;a[i+j+k]=(x-y+mod)%mod;}}}if(f==-1){reverse(a+1,a+n);int inv=ksm(n,mod-2);fo(i,0,n-1) a[i]=((ll)a[i]*inv)%mod;}}int main(){int i,j;scanf("%d%d",&N,&M);fo(i,0,N) scanf("%d",&a[i]);fo(i,0,M) scanf("%d",&b[i]);int m=N+M;for(n=1;n<=m;n<<=1) L++;fo(i,0,n-1) R[i]=(R[i>>1]>>1)|((i&1)<<L-1);NTT(a,1),NTT(b,1);fo(i,0,n-1) a[i]=(ll)a[i]*b[i]%mod;NTT(a,-1);fo(i,0,m) printf("%d ",a[i]);return 0;}
阅读全文
1 0
- UOJ #34. 多项式乘法
- 【UOJ 34】多项式乘法
- [UOJ#34]多项式乘法
- UOJ 34 多项式乘法
- UOJ 34 多项式乘法
- UOJ 34 多项式乘法
- uoj#34多项式乘法
- UOJ #34 多项式乘法
- [UOJ 34]多项式乘法(FFT)
- [挖坑][uoj]多项式乘法 FFT
- 【快速傅里叶变换】[UOJ#34]多项式乘法
- 【NTT模板】[UOJ#34]多项式乘法
- 【FFT优化】[UOJ#34]多项式乘法
- 【UOJ#34】 多项式乘法(FFT && NTT)
- [FFT 模板题] UOJ #34 多项式乘法
- [UOJ#34]多项式乘法(FFT)
- [uoj 34 多项式乘法] FFT&NTT 模板
- uoj 34 多项式乘法(fft入门)
- SQL Server 触发器的使用,实现多字段自增长
- HeadFirstJava 12 图形用户接口
- 天天特惠系统秒杀优化方案
- Castle.MVC框架介绍
- C++重载输入运算符<<自定义打印类信息
- UOJ #34. 多项式乘法
- 在asp.net页面上得到Castle容器的实例
- maven整合ssm遇到的小问题2
- CNTK-window-cpu-only编译包安装
- Ip地址查询
- 第一讲 MyBatis开端
- qduoj 140 题目6 142 ycb的ACM进阶之路
- MVC结构简介
- Kali渗透测试——UPNP网关发现工具Miranda