[LOJ 6485]LJJ学二项式定理(单位根反演)
也许更好的阅读体验
\(\mathcal{Description}\)
\(T\)组询问,每次给\(n,s,a_0,a_1,a_2,a_3\)求
\(\begin{aligned}\left(\sum ^{n}_{i=0}\begin{pmatrix} n \\ i \end{pmatrix}\cdot s^{i}\cdot a_{i\ mod\ 4}\right)mod\ 998244353\end{aligned}\)
\(\mathcal{Solution}\)
这道题要用单位根反演
\(\begin{aligned}\left[ n\ |\ k\right] =\dfrac {1}{n}\sum ^{n-1}_{i=0}w^{ik}_{n}\end{aligned}\)
我们先打表,或者查表找出在模\(998244353\)下的原根\(\omega_4^0\)并处理好\(\omega_4^1,\omega_4^2,\omega_4^3\)
考虑对每个\(a_i\)单独计算其答案
以\(a_0\)为例
\(\begin{aligned}
ans_{a_0}&=\frac{1}{4}a_0\sum_{i=0}^n [4\ |\ i]{n\choose i}s^i\\
&=\frac{1}{4}a_0\sum_{i=0}^n{n\choose i}s^i\sum_{j=0}^3 (\omega_4^{j})^i\\
&=\frac{1}{4}a_0\sum_{j=0}^3\sum_{i=0}^n {n\choose i}s^i(\omega_4^j)^i\\
&=\frac{1}{4}a_0\sum_{j=0}^3(s\omega_4^j+1)^n
\end{aligned}\)
对于其他\(a_i\),我们可以将其写成\([4\ |\ i+4-k]\)的形式,并提一个\(\omega^x\)出来
下面完整的推一遍
\(\begin{aligned}ans&=\sum ^{3}_{k=0}a_{k}\sum ^{n}_{i=0}\left[ 4\ |\ i+4-k\right] s^{i}\cdot \begin{pmatrix} n \\ i \end{pmatrix}\\
&=\sum ^{3}_{k=0}\dfrac {1}{4}a_{k}\sum ^{n}_{i=0}\sum ^{3}_{j=0}\omega^{j\cdot\left( i+4-k\right) }_{4}\cdot s^{i}\cdot \begin{pmatrix} n \\ i \end{pmatrix}\\
&=\sum ^{3}_{k=0}\dfrac {1}{4}a_{k}\sum ^{n}_{i=0}\sum ^{3}_{j=0}\omega^{\left(4j-jk\right) }_{4}\cdot \omega^{j^i}_4\cdot s^{i}\cdot \begin{pmatrix} n \\ i \end{pmatrix}\\
&=\sum ^{3}_{k=0}\dfrac {1}{4}a_{k}\sum ^{3}_{j=0}\omega^{\left(4j-jk\right) }_{4}\sum ^{n}_{i=0} \omega^{j^i}\cdot s^{i}\cdot \begin{pmatrix} n \\ i \end{pmatrix}\\
&=\sum ^{3}_{k=0}\dfrac {1}{4}a_{k}\sum ^{3}_{j=0}\omega^{\left(4j-jk\right) }_{4}\left(s\cdot \omega^j_4+1\right)^n \\
\end{aligned}\)
拿出来,再写一遍
\(\begin{aligned}ans=\sum ^{3}_{i=0}\dfrac {1}{4}a_{i}\sum ^{3}_{j=0}\omega^{\left(4j-ij\right) }_{4}\left(s\cdot \omega^j_4+1\right)^n \\
\end{aligned}\)
\(\mathcal{Code}\)
#include <cstdio>
#define ll long long
using namespace std;
const int mod = 998244353;
const int w [] = {1,911660635,998244352,86583718};
ll T,n,s,ans,res,inv;
ll a[4];
ll ksm (ll a,ll b)
{
a%=mod;
ll s=1;
for (;b;b>>=1,a=a*a%mod)
if (b&1) s=s*a%mod;
return s;
}
int main ()
{
scanf("%lld",&T);
inv=ksm(4,mod-2);
while (T--){
scanf("%lld%lld%lld%lld%lld%lld",&n,&s,&a[0],&a[1],&a[2],&a[3]);
ans=0;
for (int i=0;i<=3;++i){
res=0;
for (int j=0;j<=3;++j)
res=(res+w[(4*j-i*j)%4]*ksm(s*w[j]+1,n)%mod)%mod;
ans=(ans+a[i]*res%mod)%mod;
}
printf("%lld\n",ans*inv%mod);
}
return 0;
}
如有哪里讲得不是很明白或是有错误,欢迎指正
如您喜欢的话不妨点个赞收藏一下吧
[LOJ 6485]LJJ学二项式定理(单位根反演)的更多相关文章
- loj 6485 LJJ学二项式定理 —— 单位根反演
题目:https://loj.ac/problem/6485 先把 \( a_{i mod 4} \) 处理掉,其实就是 \( \sum\limits_{i=0}^{3} a_{i} \sum\lim ...
- LOJ 6485 LJJ 学二项式定理——单位根反演
题目:https://loj.ac/problem/6485 \( \sum\limits_{k=0}^{3}\sum\limits_{i=0}^{n}C_{n}^{i}s^{i}a_{k}[4|(i ...
- loj #6485. LJJ 学二项式定理 单位根反演
新学的黑科技,感觉好nb ~ #include <bits/stdc++.h> #define ll long long #define setIO(s) freopen(s". ...
- loj#6485. LJJ 学二项式定理(单位根反演)
题面 传送门 题解 首先你要知道一个叫做单位根反演的东西 \[{1\over k}\sum_{i=0}^{k-1}\omega^{in}_k=[k|n]\] 直接用等比数列求和就可以证明了 而且在模\ ...
- LOJ #6485 LJJ 学二项式定理
QwQ LOJ #6485 题意 求题面中那个算式 题解 墙上暴利 设$ f(x)=(sx+1)^n$ 假设求出了生成函数$ f$的各项系数显然可以算出答案 因为模$ 4$的缘故只要对于每个余数算出次 ...
- loj #6485. LJJ 学二项式定理 (模板qwq)
$ \color{#0066ff}{ 题目描述 }$ LJJ 学完了二项式定理,发现这太简单了,于是他将二项式定理等号右边的式子修改了一下,代入了一定的值,并算出了答案. 但人口算毕竟会失误,他请来了 ...
- LOJ 6485 LJJ学多项式
前言 蒟蒻代码惨遭卡常,根本跑不过 前置芝士--单位根反演 单位根有这样的性质: \[ \frac{1}{n}\sum_{i=0}^{n-1}\omega_{n}^{ki}=\left[n|k\rig ...
- 【LOJ#6485】LJJ 学二项式定理(单位根反演)
[LOJ#6485]LJJ 学二项式定理(单位根反演) 题面 LOJ 题解 显然对于\(a0,a1,a2,a3\)分开算答案. 这里以\(a0\)为例 \[\begin{aligned} Ans&am ...
- LOJ6485 LJJ 学二项式定理 解题报告
LJJ 学二项式定理 题意 \(T\)组数据,每组给定\(n,s,a_0,a_1,a_2,a_3\),求 \[ \sum_{i=0}^n \binom{n}{i}s^ia_{i\bmod 4} \] ...
随机推荐
- Spring Boot 教程系列学习
Spring Boot基础教程1-Spring Tool Suite工具的安装 Spring Boot基础教程2-RESTful API简单项目的快速搭建 Spring Boot基础教程3-配置文件详 ...
- 使用adb 命令(atrace)抓起systrace的方法。【转】
本文转载自:https://www.cnblogs.com/liuliu-word/p/9963017.html adb shell atrace -c -b 10240 --async_start ...
- java.lang.UnsupportedClassVersionError: com/mysql/cj/jdbc/Driver : Unsupported major.minor version 52.0 (unable to load class [com.mysql.cj.jdbc.Driver])
原因: com/mysql/cj/jdbc/Driver是6.0版本的驱动,兼容JDK8环境,不兼容JDK7环境,在基于jdk7的tomcat中编译运行会出错,在基于jdk8的tomcat中编译运行则 ...
- win7 激活码 秘钥
019.06最新windows7旗舰版系统激活码: 目前市面上的win7旗舰版激活码大部分都已经过期或失效了,下面来分享一些最新的. win7旗舰版激活密钥: BG2KW-D62DF-P4HY6-6J ...
- http://stblog.baidu-tech.com/?p=1684) coredump调试记录 - PHP篇 原创: 扶墙 贝壳产品技术 今天
http://stblog.baidu-tech.com/?p=1684) coredump调试记录 - PHP篇 原创: 扶墙 贝壳产品技术 今天
- 从0开始学爬虫2之json的介绍和使用
从0开始学爬虫2之json的介绍和使用 Json 一种轻量级的数据交换格式,通用,跨平台 键值对的集合,值的有序列表 类似于python中的dict Json中的键值如果是字符串一定要用双引号 jso ...
- 我是如何同时拿到阿里和腾讯offer的 【转载】
前言 三月真是一个忙碌的季节,刚刚开学就需要准备各种面试和笔试(鄙视).幸运的是,在长达一个月的面试内推季之后,终于同时拿到了阿里和腾讯的offer,还是挺开心的.突而想起久未更新的博客,就冒昧学一学 ...
- Qt编写气体安全管理系统12-设备双击
一.前言 在编写这个项目的过程中,有个得到客户夸赞的小功能就是,设备按钮双击,在离线的时候是双击重连设备,在线的时候是双击弹出具体详情界面,回控设备,参数设置等.在modbus设备通信过程中,设定了超 ...
- 【Leetcode_easy】874. Walking Robot Simulation
problem 874. Walking Robot Simulation solution1: 思路:1)如何表示移动的方向以及移动的位置坐标; 2)障碍物坐标如何检查;3)求解的是最大距离; cl ...
- python logging模块日志输出
import logging logger = logging.getLogger(__name__) logger.setLevel(level = logging.INFO) handler = ...