数论
Fizzmy
这个作者很懒,什么都没留下…
展开
-
Lucas定理的应用
题意:有n盏灯环形排列,顺时针依次标号为1⋯n。初始时刻为0,初始时刻第i盏灯的亮灭a[i]给定,0表示灭,1表示亮。下一时刻每盏灯的亮灭取决于当前时刻这盏灯与顺时针方向下一盏灯的亮灭。若两盏灯状态相同,则下一时刻该灯灭,否则该灯亮。求时刻t第k盏灯的状态。 (n,t,k≤1e7) Solution: 异或可以看成不进位的二进制加法,于是我们可以先做出一个图找找规律: 可见,对于第i个数原创 2017-10-22 15:24:27 · 758 阅读 · 0 评论 -
BZOJ4869: [Shoi2017]相逢是问候-线段树+数论
传送门题意:给出一个序列a,两种操作1.将[l,r]这段区间所有的数aiaia_i换为caicaic^{a_i}2.求[l,r]这段区间的和,对p取模1≤n≤50000; 1≤m≤50000; 1≤p≤100000000; 0<c<p; 0≤ai<p1≤n≤50000; 1≤m≤50000; 1≤p≤...原创 2018-03-14 19:16:03 · 367 阅读 · 0 评论 -
Codeforces 906D:Power Tower-扩展欧拉定理
传送门题意:给你n个数w和一个数p,q个询问,每次询问一个区间[l,r],求w_l^{(w_{l+1}^{(w_{l+2}…^{w_r})})}\mod pSolution:这道题类似于BZOJ3884但是稍有一点不同:根据扩展欧拉定理:ab≡{ab%ϕ(p) &原创 2018-02-04 14:59:38 · 856 阅读 · 0 评论 -
BZOJ2987-earthquake
权限题。题意:给定a,b,c,求满足方程Ax+By=CAx+By的非负整数解A,B=109,C=min(A,B)∗109A,BSolution:类欧几里得算法模板题没学过的点击这里变形一下这个式子:y=C−AxBy但是这个式子如果直接套用类欧几里得是包含负数的,我们需要想办法帮他变变形:考虑这个式子的几何意义:线段y=C−AxBy=\frac {C-Ax} B和原创 2018-01-10 21:53:19 · 330 阅读 · 0 评论 -
BZOJ 3930 [CQOI2015]选数-数论+递推
传送门题意:从区间[L,H](L和H为整数)中选取N个整数,求N个整数最大公约数刚好为K的选取方案有多少个,答案模1e9+7 1≤N,K≤10^9,1≤L≤H≤10^9,H-L≤10^5Solution:并不是非常理解莫反怎么做= =注意到H-L≤10^5 尝试在这上面做文章 首先简化题意: [⌈Lk⌉,⌊Hk⌋][\lceil\frac L k\rceil,\原创 2018-01-18 21:25:36 · 182 阅读 · 0 评论 -
BZOJ2154+BZOJ2693 Crash的数字表格&jzptab-莫比乌斯反演
传送门 BZOJ2154 题意:求∑ni=1∑mj=1lcm(i,j)∑i=1n∑j=1mlcm(i,j)\sum_{i=1}^n\sum_{j=1}^m lcm(i,j)Solution:看到此题首先想到这步:∑i=1n∑j=1nlcm(i,j)=∑i=1n∑j=1nijgcd(i,j)∑i=1n∑j=1nlcm(i,j)=∑i=1n∑j=1nijgcd(i,j)\sum...原创 2018-01-09 20:58:41 · 324 阅读 · 0 评论 -
BZOJ2142 礼物-扩展lucas
传送门 题意:给定n及m个数a[1…m],再给定一个模数P,求∑ni=1Cain−∑i−1j=1aj%P\sum_{i=1}^nC_{n-\sum_{j=1}^{i-1}a_j}^{a_i}\%P(1≤n≤109,1≤m≤51≤n≤10_9,1≤m≤5,设P=pc11∗pc22∗pc33∗…∗pcttP=p_1^{c_1}*p_2^{c_2}*p_3^{c_3}*…*p_t^{c_t原创 2018-01-18 13:55:35 · 243 阅读 · 0 评论 -
BZOJ3884 上帝与集合的正确用法-扩展欧拉定理
传送门 题意:求2(2(2...))modp2^{(2^{(2...)})}\mod p的值,多组询问。p≤107p≤10^7。Solution:根据扩展欧拉定理,我们知道:当b≥φ(p)b≥\varphi(p)时,ab≡ab%φ(p)+φ(p)modpa^b≡a^{b\%\varphi(p)+\varphi(p) }\mod p此公式的具体证明可以百度扩展欧拉定理我们假设f原创 2018-01-09 16:59:14 · 725 阅读 · 0 评论 -
BZOJ 4596 黑暗前的幻想乡-矩阵树定理+容斥原理
传送门题意:幽香要修建幻想乡的公路。幻想乡有 N 个城市,之间原来没有任何路。幽香向选民承诺要减税,所以她打算只修 N- 1 条路将这些城市连接起来。但是幻想乡有正好 N- 1 个建筑公司,她打算让每个建筑公司都负责一条路来修。 每个建筑公司都告诉了幽香自己有能力负责修建的路是哪些城市之间的。所以幽香打算选择 N-1 条能够连接幻想乡所有城市的边,然后每条边都交给一个能够负责该边的建筑公原创 2018-01-17 19:09:43 · 224 阅读 · 0 评论 -
Codeforces 839D Winter is here-容斥
传送门题意:给出一个n个数的数列ai" role="presentation">aiaia_i,考虑其中所有gcd大于1的集合,每个集合的贡献是gcd乘上集合的大小,求这个数列的总贡献%1e9+7(1 ≤ n ≤ 200000,1 ≤ ai原创 2018-01-25 08:07:07 · 247 阅读 · 0 评论 -
基础数论入门
(一)定理和性质一、裴蜀定理如果 a,b∈Na,b∈Na,b ∈ N , (a,b)=d(a,b)=d(a,b) = d 那么一定存在 x,yx,yx,y 使得 d|(a∗x+b∗y)d|(a∗x+b∗y)d | (a*x + b*y) 证明:非常简单,鉴于可能有数论刚入门的OIer所以这里简单证明一下: 因为(a,b)=d(a,b)=d(a,b) = d 所以我们就可以假设...原创 2017-08-12 17:02:49 · 6566 阅读 · 0 评论 -
bzoj1101[POI2007]Zap-莫比乌斯反演
bzoj 1101[POI2007]Zap-莫比乌斯反演题意:T组数据,每组给出m,n,d 求∑ni=1∑nj=1[gcd(i,j)==d](T≤50000,n,m,d<=50000)\sum_{i=1}^{n}\sum_{j=1}^{n} [gcd(i,j)==d](T≤50000,n,m,d<=50000)Solution:首先大家要记住一些结论:μ∗1=[n==1]μ*1=[n==1]类似地,原创 2017-12-04 16:59:54 · 321 阅读 · 0 评论 -
BZOJ4833: [Lydsy1704月赛]最小公倍佩尔数-数论
传送门题意:令(1+2–√)n=e(n)+f(n)∗2–√(1+2)n=e(n)+f(n)∗2(1+\sqrt2)^n=e(n)+f(n)*\sqrt2,其中e(n),f(n)e(n),f(n)e(n),f(n)都是整数。令g(n)g(n)g(n)表示f(1),f(2)…f(n)f(1),f(2)…f(n)f(1),f(2)…f(n)的最小公倍数.给定两个正整数n和p,其中p是...原创 2018-04-30 21:17:19 · 524 阅读 · 0 评论