P6800-[模板]Chirp Z-Transform【NTT】

正题

题目链接:https://www.luogu.com.cn/problem/P6800


题目大意

给出一个nnn此多项式PPP,对于k∈[0,m−1]k\in[0,m-1]k[0,m1]所有的求P(ck)P(c^k)P(ck)
输出答案对998244353998244353998244353取模
1≤n,m≤1061\leq n,m\leq 10^61n,m106


解题思路

g(n)=∑i=0n−1aici×ng(n)=\sum_{i=0}^{n-1}a_ic^{i\times n}g(n)=i=0n1aici×n
然后根据i×n=(i+n2)−(i2)−(n2)i\times n=\binom{i+n}{2}-\binom{i}{2}-\binom{n}{2}i×n=(2i+n)(2i)(2n)
g(n)=c−(n2)∑i=0n−1aic(i+n2)c−(i2)g(n)=c^{-\binom{n}{2}}\sum_{i=0}^{n-1}a_ic^{\binom{i+n}{2}}c^{-\binom{i}{2}}g(n)=c(2n)i=0n1aic(2i+n)c(2i)
然后这是一个反着卷积的形式,直接上NTT就好了

时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=4e6+10,P=998244353;
ll n,m,c,a[N],r[N],F[N],G[N];
ll power(ll x,ll b){ll ans=1;b%=P-1;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
ll C(ll n)
{return n*(n-1)/2;}
void NTT(ll *f,ll n,ll op){for(ll i=0;i<n;i++)if(i<r[i])swap(f[i],f[r[i]]);for(ll p=2;p<=n;p<<=1){ll tmp=power(3,(P+1)/p),len=p>>1;if(op==-1)tmp=power(tmp,P-2);for(ll k=0;k<n;k+=p){ll buf=1;for(ll i=k;i<k+len;i++){ll tt=f[i+len]*buf%P;f[i+len]=(f[i]-tt+P)%P;f[i]=(f[i]+tt)%P;buf=buf*tmp%P;}}}if(op==-1){ll invn=power(n,P-2);for(ll i=0;i<n;i++)f[i]=f[i]*invn%P;}return;
}
signed main()
{scanf("%lld%lld%lld",&n,&c,&m);ll inv=power(c,P-2);for(ll i=0;i<n;i++)scanf("%lld",&a[i]);for(ll i=0;i<n+m;i++)F[i]=power(c,C(n+m-i-1));for(ll i=0;i<n;i++)G[i]=a[i]*power(inv,C(i))%P;ll len=1;while(len<n+m)len<<=1;for(ll i=0;i<len;i++)r[i]=(r[i>>1]>>1)|((i&1)?(len>>1):0);NTT(F,len,1);NTT(G,len,1);for(ll i=0;i<len;i++)F[i]=F[i]*G[i]%P;NTT(F,len,-1);for(ll i=n+m-1;i>=n;i--)printf("%lld ",F[i]*power(inv,C(n+m-i-1))%P);return 0;
}

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

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

相关文章

AtCoder Grand Contest 021 D - Reversed LCS(区间dp)

D - Reversed LCS 繁凡さん 设 f[l,r,k]f [ l , r , k ]f[l,r,k] 表示区间 [l,r][ l , r ][l,r] 中修改 kkk 次能得到的最长回文子序列的长度&#xff0c;直接区间DP转移即可。 字符串的 最长回文子序列&#xff08;lps&#xff09; 长度等于其自身与反转的 最长公共子序列&am…

牛客题霸 [滑动窗口的最大值] C++题解/答案

牛客题霸 [滑动窗口的最大值] C题解/答案 题目描述 给定一个数组和滑动窗口的大小&#xff0c;找出所有滑动窗口里数值的最大值。例如&#xff0c;如果输入数组{2,3,4,2,6,2,5,1}及滑动窗口的大小3&#xff0c;那么一共存在6个滑动窗口&#xff0c;他们的最大值分别为{4,4,6,…

基于.NET Standard的分布式自增ID算法--Snowflake

概述本篇文章主要讲述分布式ID生成算法中最出名的Snowflake算法。搞.NET开发的&#xff0c;数据库主键最常见的就是int类型的自增主键和GUID类型的uniqueidentifier。那么为何还要引入snowflake呢&#xff1f;INT自增主键自增主键是解决主键生成的最简单方案&#xff0c;它有如…

水题杂做(

前言 平时闲的无聊写的或者口嗨的水题就不水博客了。 CF1545A AquaMoon and Strange Sort 题意 给出nnn个数字的一个序列&#xff0c;每次可以交换相邻的两个数字&#xff0c;求是否能够将序列从小到大排序且每个数字被交换次数为偶数。 n≤105n\leq 10^5n≤105 解题思路 …

2021“MINIEYE杯”中国大学生算法设计超级联赛(4)Display Substring(后缀数组+二分)

Display Substring #include<bits/stdc.h> using namespace std; typedef long long ll; // sa[i]: 排名是i位的是第几个后缀 // rk[i]: 第i个后缀的排名是多少 // height[i]: sa[i]与sa[i-1] const int N100010; char s[N]; int rk[N],sa[N],cnt[N],height[N]; int x[N]…

线性筛素数的实现与证明

情境 大家应该都知道用nlogn的时间复杂度筛出[1,n]的所有素数&#xff0c;但是当n的范围较大时&#xff0c;这个方法就不奏效了 今天我们谈谈素数的线性筛法&#xff0c;也就是用On的时间复杂度筛出[1,n]的所有素数 解析 nlog的算法之所以会慢&#xff0c;是因为它进行了很多…

牛客题霸 [在转动过的有序数组中寻找目标值] C++题解/答案

牛客题霸 [在转动过的有序数组中寻找目标值] C题解/答案 题目描述 给出一个转动过的有序数组&#xff0c;你事先不知道该数组转动了多少 (例如,0 1 2 4 5 6 7可能变为4 5 6 7 0 1 2). 在数组中搜索给出的目标值&#xff0c;如果能在数组中找到&#xff0c;返回它的索引&#…

领域驱动设计,让程序员心中有码(二)

引子&#xff0c;软件工程没有银弹上一篇博文领域驱动设计&#xff0c;让程序员心中有码&#xff0c;抛出了一个问题&#xff0c;领域驱动设计真的是万能的良方吗&#xff1f;对于这个问题&#xff0c;大家的答案无疑是一致的&#xff0c;作为一种非常受软件行业欢迎的软件思想…

bzoj#2407-探险【最短路,二进制分组】

正题 题目链接:https://darkbzoj.tk/problem/2407 题目大意 nnn个点的一张无向图&#xff08;但是正反权值不同&#xff09;&#xff0c;求一个从111出发回到111且不经过重复边的最短路径。 1≤n≤10000,1≤m≤21051\leq n\leq 10000,1\leq m\leq 2\times 10^51≤n≤10000,1≤…

The 2020 ICPC Asia Macau Regional Contest A. Accelerator(分治+NTT)

A. Accelerator 实质上就是求n个多项式相乘 (a1x)(a2x)…(anx)(a_1x)(a_2x)\dots(a_nx) (a1​x)(a2​x)…(an​x) 对于[xm][x^m][xm]的系数即是从nnn个选出n−mn-mn−m的乘积之和。 每次从nnn个选出n−mn-mn−m的概率显然为 (nm)n!\frac{\binom{n}{m}}{n!} n!(mn​)​ 然后期…

牛客题霸 [栈和排序] C++题解/答案

牛客题霸 [栈和排序] C题解/答案 题目描述 给你一个1->n的排列和一个栈&#xff0c;入栈顺序给定 你要在不打乱入栈顺序的情况下&#xff0c;对数组进行从大到小排序 当无法完全排序时&#xff0c;请输出字典序最大的出栈序列 题解&#xff1a; 栈的性质是先进后出 当栈…

邮局-[IOI2000](四边形不等式)

概要 四边形不等式的核心在于缩小最优转移的范围 题目描述 传送门 解析 这道题说是不等式&#xff0c;但其实也可以感性理解 &#xff08;其实就是不想证明&#xff09; 定义pl[i][k]: i到n的村庄建造k座邮局时&#xff0c;第一座管辖的范围是i-pl[i][k] (也就是最优决策…

.NET Core实战项目之CMS 第九章 设计篇-白话架构设计

前面两篇文章给大家介绍了我们实战的CMS系统的数据库设计&#xff0c;源码也已经上传到服务器上了。今天我们就好聊聊架构设计&#xff0c;在开始之前先给大家分享一下这几天我一直在听的《从零开始学架构》里面关于架构设计的定义以及架构设计的三大原则&#xff0c;希望能对大…

P6847-[CEOI2019]Magic Tree【dp,线段树合并】

正题 题目链接:https://www.luogu.com.cn/problem/P6847 题目大意 nnn个点的一棵树上&#xff0c;每个时刻可以割掉一些边&#xff0c;一些节点上有果实表示如果在did_idi​时刻这个点恰好不与111联通&#xff0c;那么就可以获得wiw_iwi​的价值。 1≤n,k≤1051\leq n,k\leq 1…

P3899 [湖南集训]更为厉害(线段树合并、长链剖分、二维数点)

P3899 [湖南集训]更为厉害 若 deepb<deepa\text{deep}_b<\text{deep}_adeepb​<deepa​&#xff1a;c 在点 a 的子树中&#xff0c;根据乘法原理计算答案为 min⁡(deepa,k)∗(sza−1)\min(\text{deep}_a,k)∗(\text{sz}_a−1)min(deepa​,k)∗(sza​−1)若 deepb>…

树链剖分-学习报告

概要 &#xff08;这是笔者刚学完趁热写的博客 由于也是初学&#xff0c;所以可能会对初学读者比较友好吧 毕竟你不明白的地方我可能也迷惑过&#xff09; 学完感觉树链剖分的设计真的很巧妙 其主要用途是压缩树上路径 又是一个n到log的优化 原理 先引定义&#xff1a; 设…

牛客题霸 [拼接所有的字符串产生字典序最小的字符串] C++题解/答案

牛客题霸 [拼接所有的字符串产生字典序最小的字符串] C题解/答案 题目描述 给定一个字符串的数组strs&#xff0c;请找到一种拼接顺序&#xff0c;使得所有的字符串拼接起来组成的字符串是所有可能性中字典序最小的&#xff0c;并返回这个字符串。 题解&#xff1a; 用sort…

今日头条Marketing API小工具(.Net Core版本)

前言由于工作原因&#xff0c;需要用到今日头条的Marketing API做一些广告投放的定制化开发。然后看现在网上也没多少关于头条Marketing API的文章&#xff0c;于是便就有了该篇文章。头条Marketing API主页地址&#xff1a;https://ad.toutiao.com/openapi/index.html。头条Ma…

P6880-[JOI 2020 Final]オリンピックバス【最短路】

正题 题目链接:https://www.luogu.com.cn/problem/P6880 题目大意 给出nnn个点mmm条边的有向图&#xff0c;边有边权和一个翻转权值。 翻转至多一条边使得1−>n−>11->n->11−>n−>1往返的权值加上翻转权值最小。 1≤n≤200,1≤m≤51041\leq n\leq 200,1\…

The 2020 ICPC Asia Macau Regional Contest J. Jewel Grab(数颜色+链表)

J. Jewel Grab Tartarus _Wallace_ 转化询问&#xff1a;对于一个询问 [s,k]&#xff0c;找到一个最长的区间 [s,t]&#xff0c;满足区间中出现次数超过一次的元素&#xff0c;的出现次数减一&#xff0c;的和&#xff0c;不超过 k。 对于该区间[s,t] 区间数颜色&#xff1a;…