【BZOJ3328】PYXFIB【矩阵快速幂】【单位根反演】【二项式定理】

传送门

题意:

∑i=0⌊nk⌋(nik)Fik\sum_{i=0}^{\lfloor\frac nk\rfloor}\binom n{ik}F_{ik}i=0kn(ikn)Fik

FFF为斐波拉契数列

n≤1e18,k≤2e4,p≤1e9n\leq 1e18,k\leq2e4,p\leq1e9n1e18,k2e4,p1e9且为质数且模kkk111


显然就是求

∑i=0n[k∣i](ni)Fi\sum_{i=0}^n[k\mid i]\binom niF_ii=0n[ki](in)Fi

显然FFF可以通过矩阵表示出来

A=(1110),S=(10)A=\left(\begin{matrix}1&1\\1&0\end{matrix}\right),S=\left(\begin{matrix}1\\0\end{matrix}\right)A=(1110),S=(10)

那么只需要求出

∑i=0n[k∣i](ni)AiS\sum_{i=0}^n[k\mid i]\binom niA^iSi=0n[ki](in)AiS

相当于求

∑i=0n[k∣i](ni)Ai\sum_{i=0}^n[k\mid i]\binom niA^ii=0n[ki](in)Ai

然后日常把这个幻视成二项式定理害得我调了两个小时

显然可以单位根反演

1k∑j=0k−1∑i=0n(ni)Aiωij\frac 1k\sum_{j=0}^{k-1}\sum_{i=0}^n\binom niA^i\omega^{ij}k1j=0k1i=0n(in)Aiωij

=1k∑j=0k−1∑i=0n(ni)(Aωj)i=\frac 1k\sum_{j=0}^{k-1}\sum_{i=0}^n\binom ni(A\omega^j)^i=k1j=0k1i=0n(in)(Aωj)i

二项式定理走一波

=1k∑j=0k−1(Aωj+I)n=\frac 1k\sum_{j=0}^{k-1}(A\omega^j+I)^n=k1j=0k1(Aωj+I)n

然后就可以算了

因为p%k=1p\%k=1p%k=1,所以可以找一个原根算出ω\omegaω

复杂度O(klog⁡n)O(k\log n)O(klogn)

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
using namespace std;
typedef long long ll;
ll n;
int k,MOD;
inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;}
inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;}
inline int qpow(int a,int b)
{int ans=1;while (b){if (b&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD;b>>=1;}return ans;
}
int g,Wn;
inline bool check(int x)
{int n=x;for (int i=2;i*i<=x;i++)if (x%i==0){if (qpow(g,n/i)==1) return false;while (x%i==0) x/=i;}if (x>1&&qpow(g,n/x)==1) return false;return true;
}
struct mat
{int e[2][2];mat(const int& x=0){e[0][0]=e[1][1]=x;e[0][1]=e[1][0]=0;}int* operator [](int i){return e[i];}	
};
inline int det(mat m){return dec((ll)m[0][0]*m[1][1]%MOD,(ll)m[0][1]*m[1][0]%MOD);}
inline mat operator *(mat a,mat b)
{mat ans;ans[0][0]=add((ll)a[0][0]*b[0][0]%MOD,(ll)a[0][1]*b[1][0]%MOD);ans[0][1]=add((ll)a[0][0]*b[0][1]%MOD,(ll)a[0][1]*b[1][1]%MOD);ans[1][0]=add((ll)a[1][0]*b[0][0]%MOD,(ll)a[1][1]*b[1][0]%MOD);ans[1][1]=add((ll)a[1][0]*b[0][1]%MOD,(ll)a[1][1]*b[1][1]%MOD);	return ans;
}
inline mat operator +(mat a,mat b)
{mat ans;ans[0][0]=add(a[0][0],b[0][0]);ans[0][1]=add(a[0][1],b[0][1]);ans[1][0]=add(a[1][0],b[1][0]);ans[1][1]=add(a[1][1],b[1][1]);	return ans;
}
inline mat operator *(mat a,const int& t)
{mat ans;ans[0][0]=(ll)a[0][0]*t%MOD;ans[0][1]=(ll)a[0][1]*t%MOD;ans[1][0]=(ll)a[1][0]*t%MOD;ans[1][1]=(ll)a[1][1]*t%MOD;	return ans;	
}
inline mat qpow(mat a,ll p)
{mat ans(1);while (p){if (p&1) ans=ans*a;a=a*a;p>>=1;}return ans;
}
int main()
{int T;scanf("%d",&T);while (T--){scanf("%lld%d%d",&n,&k,&MOD);for(g=2;!check(MOD-1);g++);Wn=qpow(g,(MOD-1)/k);mat T,sum;T[0][0]=T[0][1]=T[1][0]=1;T[1][1]=0;for (int i=0;i<k;i++,T=T*Wn) sum=sum+qpow(T+mat(1),n);sum=sum*qpow(k,MOD-2);printf("%d\n",sum[0][0]);}return 0;
}

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

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

相关文章

.NET分布式框架 | Orleans 知多少

引言公司物联网项目集成Orleans以支持高并发的分布式业务&#xff0c;对于Orleans也是第一次接触&#xff0c;本文就分享下个人对Orleans的理解。这里先抛出自己的观点&#xff1a;Orleans 是一个支持有状态云生应用/服务水平伸缩的基于Virtual Actor 模型的.NET分布式框架。下…

nowcoder 牛牛的最大兴趣组 质因子 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先nnn很小的话可以暴力连边&#xff0c;让后染个色求一个颜色最多的即可。但是这个题显然不行&#xff0c;由于是三次方&#xff0c;所以考虑质因子入手。 首先很容易就能想到将所有的数…

ASP.NET Core on K8S学习初探(3)部署API到K8S

“ 终于可以部署ASP.NET Core到K8S中了...”在上一篇《基本概念快速一览》中&#xff0c;我们把基本的一些概念快速地简单地不求甚解地过了一下&#xff0c;本篇开始我们会将ASP.NET Core WebAPI部署到K8S&#xff0c;从而结束初探的旅程。01—准备一个WebAPI这里准备一个空的A…

K-D Tree学习笔记

引入 K-D Tree 是一种处理高维空间的数据结构。 支持O(nk−1k)O(n^{\frac {k-1}k})O(nkk−1​)查询给定超矩形内的点的信息&#xff0c; kkk 为维数。 可以用替罪羊树的思想实现动态插入。 但其实用的最多的是在错误的复杂度内查询奇奇怪怪的点信息。 构建 K-D Tree 是平…

Codeforces Round #633 B. Edge Weight Assignment 结论题 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先考虑最小值&#xff0c;如果从一个叶子结点出发到任意叶子的距离都为偶数&#xff0c;那么只需要一个值就可以满足条件。如果有奇数的&#xff0c;考虑111 ^ 222 ^ 303030&#xff0c;…

ASP.NET Core Web API中使用Swagger

本节导航Swagger介绍在ASP.NET CORE 中的使用swagger在软件开发中,管理和测试API是一件重要而富有挑战性的工作。在我之前的文章《研发团队,请管好你的API文档》 也专门阐述了通过文档管理工具,来保证API文档和代码的一致性,这样更加有助于团队的协作。以往我们总是通过第三方平…

【ZJOI2015】幻想乡战略游戏【点分树】【带权重心】

题意&#xff1a;nnn个点带边权的树&#xff0c;动态修改点权viv_ivi​&#xff0c;最小化 钦定一个点xxx 后 ∑idist(x,i)∗vi\sum\limits_{i} dist(x,i)*v_ii∑​dist(x,i)∗vi​的值。 n,q≤105n,q \leq10^5n,q≤105&#xff0c;度数不超过202020 限制度数的树上的一些诡异…

P4559 [JSOI2018]列队 主席树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个学生以及其位置&#xff0c;mmm个询问&#xff0c;每次询问[l,r][l,r][l,r]的学生跑到[k,kr−l][k,kr-l][k,kr−l]的位置的最小总的移动距离。当然不同的学生不能站在同一位置。 n,m≤5e5,1≤ai,k≤…

C# 8.x 先睹为快

尽管 C# 8.0 还需要几个月的时间才能发布&#xff0c;但是 C# 8.x 的特性计划已经开始了。该计划中的特性有一些是新提出的&#xff0c;而另外一些则是之前在 C# 8 中考虑过的。和往常一样&#xff0c;这个特性列表仍可能会发生变化。如之前的一个报告所讨论的那样&#xff0c;…

【WC2014】紫荆花之恋【替罪羊思想】【动态点分树】【替罪羊树】

传送门 题意&#xff1a;有一棵开始时没有结点的树&#xff0c;nnn次询问&#xff0c;每次新加一点并给定父结点、到父亲的距离、参数rir_iri​,并询问满足dist(u,v)≤rurvdist(u,v)\leq r_ur_vdist(u,v)≤ru​rv​的点对(u,v)(u,v)(u,v)的对数。 n≤105n\leq 10^5n≤105&…

NuGet是什么?理解与使用

如果你了解python&#xff0c;那么它类似pip。如果你了解nodejs&#xff0c;那么它类似npm。如果你了解ruby&#xff0c;那么它类似gem。对&#xff0c;它就是一个包&#xff08;package&#xff09;管理平台&#xff0c;确切的说是 .net平台的包管理工具&#xff0c;它提供了一…

P2468 [SDOI2010]粟粟的书架 主席树 + 二分 + 二维前缀和

传送门 文章目录题意&#xff1a;题意&#xff1a;题意&#xff1a; 给你一个n∗mn*mn∗m的矩阵&#xff0c;每次询问一个矩形&#xff0c;左上角是(x1,y1)(x_1,y_1)(x1​,y1​)&#xff0c;右上角是(x2,y2)(x_2,y_2)(x2​,y2​)&#xff0c;和一个hhh&#xff0c;问是否能选出…

带花树学习笔记

引入 带花树解决的是一般图的最大匹配问题。 学习此算法前你需要掌握二分图匹配的匈牙利算法&#xff0c;至少需要知道它的思想。 理论 众所周知&#xff0c;二分图匹配的思想是不断地找增广路。 严谨地讲&#xff0c;是找到了一条简单路径&#xff0c;它的起点和终点都未…

你竟然没用 Nuget 构建项目?

想写这篇文章缘起于公众号文章里的一条留言&#xff1a;对于一个现代化的开发平台&#xff0c;建立一种让开发者创建&#xff0c;分享与使用可复用代码的机制是十分必要的。这种“可复用代码”被打包后的文件通常被称作“包”&#xff08;package&#xff09;&#xff0c;对于.…

Codeforces Round #339 (Div. 1) D. Kingdom and its Cities 虚树 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤1e5,∑k≤1e5,q≤1e5n\le1e5,\sum k\le1e5,q\le1e5n≤1e5,∑k≤1e5,q≤1e5。 思路&#xff1a; 经过分析&#xff0c;外敌占领的城市只有可能是两点的lcalcalca或者两点之间任意一点&#xff0c;是哪个点…

小白开学Asp.Net Core《四》 —— 使用AspectCore-Framework

小白开学Asp.Net Core《四》—— 使用AspectCore-Framework一、AspectCore-Framework说AspectCore-Framework不得不先谈谈的AOP&#xff0c;AOP&#xff1a;在软件业&#xff0c;AOP为Aspect Oriented Programming的缩写&#xff0c;意为&#xff1a;面向切面编程&#xff0c;通…

【WC2016】挑战NPC 【带花树】【建图】

传送门 题意&#xff1a;有nnn个球和mmm个筐&#xff0c;每个筐最多放333个球&#xff0c;每个球只能放入特定的一些筐中&#xff0c;在题中给出。构造一种放球的方案使得nnn个球都被放在某个筐中且 球的个数不超过111 的筐的数量尽量大。 m≤100,n≤3mm\leq 100,n\leq 3mm≤1…

P4103 [HEOI2014]大工程 虚树 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 简化一下题意就是求树上给定点集中每两个点之间的距离之和&#xff0c;相距最远的点&#xff0c;相距最近的点。 对于距离和我们统计一下边的贡献就好啦&#xff0c;边两头的sizesizesize乘…

【翻译】Keras.NET简介 - 高级神经网络API in C#

用C#&#xff0c;远离996Keras.NETKeras.NET是一个高级神经网络API&#xff0c;它使用C#编写&#xff0c;并带有Python绑定&#xff0c;可以在Tensorflow、CNTK或Theano上运行。其关注点是实现快速实验。因为做好研究的关键是&#xff1a;能在尽可能短的时间内从一个想法发展出…

TIOBE 6月排行:C# 以微弱的优势超过了 Visual Basic .NET 的排名,再次进入 TOP 5

月经贴警告……TIOBE 编程语言排行榜 7 月更新已公布&#xff0c;排名前十的分别是&#xff1a;Java, C, Python, C, C#, Visual Basic .NET, JavaScript, PHP, SQL 和汇编语言。和上个月的不同之处主要是 C# 以微弱的优势超过了 Visual Basic .NET 的排名&#xff0c;再次进入…