#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
struct cp{double x,y;};
int n1,n2,n,m;
double pi=acos(-);
cp a[],b[],cur[];
cp operator *(cp x,cp y){return (cp){x.x*y.x-x.y*y.y,x.x*y.y+x.y*y.x};}
cp operator +(cp x,cp y){return (cp){x.x+y.x,x.y+y.y};}
cp operator -(cp x,cp y){return (cp){x.x-y.x,x.y-y.y};}
void fft(cp *a,int n,int fl)
{
for (int i=n>>,j=;j<n;j++)
{
if (i<j) swap(a[i],a[j]);
int k=n>>;
for (;k&i;i^=k,k>>=);
i^=k;
}
for (int m=;m<=n;m<<=)
{
cp w=(cp){cos(*pi*fl/m),sin(*pi*fl/m)};
cur[]=(cp){,};
for (int i=;i<m;i++) cur[i]=cur[i-]*w;
for (int i=;i<n;i+=m)
for (int j=i;j<i+(m>>);j++)
{
cp u=a[j],v=a[j+(m>>)]*cur[j-i];
a[j]=u+v;
a[j+(m>>)]=u-v;
}
}
}
int main()
{
scanf("%d%d",&n1,&n2);n1++;n2++;
for (int i=;i<n1;i++) scanf("%lf",&a[i].x);
for (int i=;i<n2;i++) scanf("%lf",&b[i].x);
n=max(n1,n2);
m=;while (m<=n*) m*=;
fft(a,m,);fft(b,m,);
for (int i=;i<=m;i++) a[i]=a[i]*b[i];
fft(a,m,-);
for (int i=;i<n1+n2-;i++) printf("%d ",(int)(a[i].x/m+0.5));
}

最新文章

  1. solr服务(搜索服务)
  2. mysql 表被锁时,需要执行的命令
  3. python查找并删除相同文件-UNIQ File-wxPython版本
  4. Docker之功能汇总
  5. NavMesh系统动态碰撞的探讨
  6. 配置ssh免密码连接
  7. Python: Lambda Functions
  8. [ActionScript 3.0] AS3 时间日期格式化DateTimeFormatter类的运用
  9. rabbitmq 消息持久化
  10. thinkphp3.2自定义常量
  11. FineUI经典项目展示(1)生产在线管理系统
  12. 【BZOJ4873】[六省联考2017]寿司餐厅(网络流)
  13. noip模拟【service】
  14. POJ-2386.Lakecounting(DFS求连通块)
  15. strcpy、memcpy和memset的区别
  16. smartgit 过期
  17. 主席树 或者 离散化+分块 BZOJ 4636
  18. python 日志的配置,python对日志封装成类,日志的调用
  19. DNS原理浅析
  20. 2 分支语句——《Swift3.0 从入门到出家》

热门文章

  1. 【洛谷2617_BZOJ1901】Dynamic Rankings(树套树)
  2. 【知识总结】多项式全家桶(三)(任意模数NTT)
  3. ACM_括号匹配
  4. magento getUrl函数跳转admin模块问题
  5. Oracle数据库数据显示乱码问题解决方法。
  6. Java 8 (6) Stream 流 - 并行数据处理与性能
  7. 在Azure Ubunt Server 14.04虚机中使用Deep-Visualization-Toolbox
  8. C#和Java在语法上的差异(原创,持续更新中)
  9. Java编程思想读书笔记_第8章
  10. mvc使用linq to sql进行sum统计遇到查询为null的问题