![](https://img-blog.csdnimg.cn/20201014180756923.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
文章平均质量分 63
Vectorxj
这个作者很懒,什么都没留下…
展开
-
[数论][莫队][莫比乌斯反演] hdu 4676 Sum Of Gcd
DescriptionDescription给定一个nn排列,qq次询问,每次询问∑L≤i<j≤Rgcd(ai,aj)\sum_{L\le i\lt j\le R}gcd(a_i,a_j)。SolutionSolution首先先推一推柿子:∑L≤i<j≤R(ai,aj)==∑L≤i<j≤R∑d|(ai,aj)d∑k|(aid,ajd)μ(k)∑L≤i<j≤R∑k|(ai,aj)∑d|kμ(kd)d原创 2017-10-01 17:18:02 · 426 阅读 · 0 评论 -
[数论][莫比乌斯反演][杜教筛] BZOJ 3512: DZY Loves Math IV
DescriptionDescription求∑i=1n∑j=1mφ(ij)\sum_{i=1}^n\sum_{j=1}^m\varphi(ij)SolutionSolution好强的题 设S(n,m)=∑i=1mφ(ni)S(n,m) =\sum_{i=1}^m\varphi(ni) 当|μ(n)|=1|\mu(n)|=1时φ(ni)=φ(i)∑d|(n,k)φ(nd)\begin{eqnar原创 2017-12-20 20:18:26 · 595 阅读 · 0 评论 -
[二次剩余]求解二次剩余
DescriptionDescription求解x2≡n(modp)x^2 \equiv n \pmod p。pp是一个奇质数。SolutionSolution由费马小定理np−1≡1(modp)n^{p-1}\equiv 1 \pmod p所以np−12≡±1(modp)n^{{p-1}\over 2}\equiv \pm1 \pmod p由欧拉准则(np)≡np−12(modp原创 2017-08-01 10:44:42 · 2383 阅读 · 0 评论 -
[子集DP][Lucas定理] BZOJ 4903: [Ctsc2017]吉夫特
SolutionSolution原来可以把简单的题面写的那么长~ 根据Lucas定理,(nm)≡(⌊n2⌋⌊m2⌋)(nmod2mmod2)(mod2){n\choose m}\equiv {{\lfloor{n\over2}\rfloor}\choose{\lfloor{m\over2}\rfloor}}{n\bmod2\choose m\bmod 2}\pmod2有(00)=1,(01)原创 2018-01-15 16:32:13 · 504 阅读 · 0 评论 -
[原根][数论][二项式定理] LOJ #6247. 九个太阳
SolutionSolution原题 wans====GP−1k∑i=0n(ni)[k∣i]∑i=0n(ni)1k∑j=0k−1wij1k∑j=0k−1(wj+1)n\begin{eqnarray}w & = & G^{P-1\over k}\\ans & = & \sum_{i=0}^n{n\choose i}[k\mid i]\\ & = & \sum_{i=0}^n{n\cho原创 2018-01-24 14:42:55 · 443 阅读 · 0 评论 -
[FFT][离散对数][原根] LOJ #6156. A * B Problem
SolutionSolutionSolution刚开始看错题了。。。其实就是一个FFT裸题 如果p∤a,p∤bp∤a,p∤bp\nmid a,p\nmid bab≡k(modp)ab≡k(modp)ab\equiv k\pmod p是与logGa+logGb≡logGk(modφ(p))logGa+logGb≡logGk(modφ(p))\log_Ga+\log_Gb\equiv \l...原创 2018-01-18 10:49:23 · 680 阅读 · 0 评论 -
[杜教筛][莫比乌斯反演] LOJ #6229. 这是一道简单的数学题
SolutionSolution求∑i=1n∑j=1iij(i,j)2\sum_{i=1}^n\sum_{j=1}^i{ij\over(i,j)^2}可能算是裸题? 相当于求12∑d=1n(d[d=1]+d2φ(d))⌊nd⌋{1\over2}\sum_{d=1}^n\big(d[d=1]+d^2\varphi(d)\big)\lfloor{n\over d}\rfloor 关于idn⋅φ原创 2018-01-18 12:36:05 · 758 阅读 · 0 评论 -
[数论][构造][离散对数] Codeforces 913 G. Power Substring
SolutionSolution设aa长度为nn。 可以尝试构造一个数bb,使得a⋅10m+ba\cdot10^m+b成为其后缀。 就等价于x=a⋅10m+b≡2k(mod10n+m)x=a\cdot10^m+b\equiv2^k\pmod {10^{n+m}} 构造k≥n+mk\ge n+m,那要满足这样的kk就等价于2n+m∣x∧5∤x2^{n+m}\mid x\wedge5\nmi原创 2018-01-10 12:59:28 · 636 阅读 · 0 评论 -
[Bluestein's Algorithm][DFT] CodeChef REALSET
SolutionSolutionSolution由F(b)≠0,F(a∗b)=0F(b)≠0,F(a∗b)=0\mathcal{F}(\textbf{b}) \neq \textbf{0},\mathcal{F}(\textbf{a}*\textbf{b})=0可以得到F(a)F(a)\mathcal{F}(\textbf{a})至少有一项为000。 考虑aa\textbf{a}向量的DFT...原创 2018-03-01 21:42:56 · 628 阅读 · 0 评论 -
[数论][莫比乌斯反演] 51Nod 1222 最小公倍数计数
DescriptionDescriptionf(n)=∑i∑j[lcm(i,j)=n]f(n)=\sum_{i}\sum_{j}[lcm(i,j)=n]求∑n=abf(n)\sum_{n=a}^bf(n)SolutionSolutionS(n)=∑k=1n√μ(k)∑a∑b∑d[abd≤nk2]\begin{eqnarray}S(n)&=&\sum_{k=1}^{\sqrt{n}}\mu(k)\s原创 2017-12-22 13:21:03 · 316 阅读 · 0 评论 -
[数论][莫比乌斯反演][杜教筛] 51Nod 1220 约数之和
DescriptionDescription求∑i=1n∑i=1nd(ij)\sum_{i=1}^n\sum_{i=1}^nd(ij)SolutionSolution∑i=1n∑j=1nd(ij)A(n)B(n)S(n)=======∑i=1n∑j=1n∑a∣i∑b∣j[(a,b)=1]ajb∑k=1nμ(k)kA(⌊nk⌋)B(⌊nk⌋)∑a=1na⌊na⌋∑b=1n∑j=1⌊nb⌋j∑i=1n(原创 2017-12-21 17:16:24 · 323 阅读 · 0 评论 -
[数论][莫比乌斯反演] BZOJ 4816: 数字表格
DescriptionDescription求∏i=1n∏i=1mF(i,j)\prod_{i=1}^n\prod_{i=1}^mF_{(i,j)} 1≤n,m≤1061\le n,m\le10^6。SolutionSolution推一下柿子: ∏i=1n∏i=1mF(i,j)====∏d=1nF∑ni=1∑mj=1[(i,j)=d]d∏d=1nF∑⌊nd⌋k=1μ(k)⌊ndk⌋⌊mdk⌋d∏原创 2017-10-01 19:52:12 · 311 阅读 · 0 评论 -
[数论][二项式定理][矩阵乘法] BZOJ 3328: PYXFIB
DescriptionDescription求∑i=0⌊nk⌋(nik)Fik\sum_{i=0}^{\lfloor{n\over k}\rfloor}{{n\choose {ik}}F_{ik}}1≤n≤1018,1≤k≤2×104,p1\le n\le10^{18},1\le k\le2\times10^4,p为质数且pmodk=1p\bmod k = 1SolutionSolution设gg为原创 2017-10-03 13:04:58 · 1241 阅读 · 0 评论 -
[数论]Codeforces 757E. Bash Plays with Functions
DescriptionDescriptionf0(n)fr+1==∑uv=n[(u,v)=1]1∗fr\begin{eqnarray}f_0(n)&=&\sum_{uv=n}[(u,v)=1]\\f_{r+1}&=&1*f_r\\\end{eqnarray}求fr(n)f_r(n)。SolutionSolution因为nn的所有质因子之间对f0(n)f_0(n)贡献是独立的,所以f0f_0是个原创 2017-10-21 08:37:40 · 373 阅读 · 0 评论 -
BZOJ 3239 [BSGS]
DescriptionDescription给定A,BA,B,求最小的xx,满足Ax≡B(modP)A^x \equiv B \pmod PSolutionSolutionBSGSBSGS裸题#include <bits/stdc++.h>using namespace std;typedef long long ll;int A, B, P, Ans;map<int, int> ap;inli原创 2017-09-05 19:52:34 · 533 阅读 · 0 评论 -
[数论] 51nod 1365 Fib(N) mod Fib(K)
DescriptionDescriptionFib(N)Fib(N)表示斐波那契数列的第NN项(F(0)=0,F(1)=1)(F(0) = 0, F(1) = 1),给出NN和KK,求Fib(N)modFib(K)Fib(N) \bmod Fib(K)。由于结果太大,输出Mod 1000000007Mod~1000000007的结果。SolutionSolution先贴几个公式。F−n=(−1)n−原创 2017-09-29 15:21:54 · 830 阅读 · 0 评论 -
BZOJ 4916 [莫比乌斯反演][杜教筛]
DescriptionDescription求∑i=1Nμ(i2)\sum_{i=1}^N{\mu (i^2)}以及∑i=1Nφ(i2)\sum_{i=1}^N{\varphi (i^2)}SolutionSolution第一问puts("1"); 第二问考虑杜教筛,首先发现式子求的就是S(N)=∑i=1Nφ(i)∗iS(N) = \sum_{i=1}^N{\varphi (i)*i}我们令f(x原创 2017-06-03 20:28:33 · 534 阅读 · 0 评论 -
[数论][二次剩余][BSGS] CodeChef FN
DescriptionDescription设FnF_n是斐波那契数列第nn项,给定FnmodPF_n\bmod P和PP求最小的nn。 保证55在膜PP意义下有二次剩余。SolutionSolution经过一系列xjb推会发现FnmodP{F_n \bmod P}是循环的。所以答案小于PP。 首先用CipollaCipolla算法求出55的二次剩余。 考虑FnF_n的通项公式。设x=5√x=原创 2017-09-05 20:02:27 · 578 阅读 · 0 评论 -
[数论][莫比乌斯反演][杜教筛] 51Nod 1238 最小公倍数之和 V3
DescriptionDescription求∑i=1n∑j=1nlcm(i,j)\sum_{i=1}^n\sum_{j=1}^nlcm(i,j) 推柿子:A(n)Ans========∑i=1nlcm(i,n)n∑d∣n∑i=1nid[(nd,id)=1]n∑d∣n∑i=1di[(d,i)=1]n∑d∣n[d=1]+dφ(d)2n2+12∑d∣ndφ(d)2∑i=1nA(i)−n(n+1)2∑i原创 2017-12-21 15:50:15 · 455 阅读 · 0 评论 -
[数论][Lucas定理] 51Nod 1778 小Q的集合
SolutionSolutionSolution通过枚举⌊im⌋⌊im⌋\lfloor{i\over m}\rfloor和imodmimodmi\bmod m。由Lucas定理是可以独立计算贡献的。 ans====∑i=0n(ni)(ik−(n−i)k)2∑i=0n(⌊nm⌋⌊im⌋)(nmodmimodm)(ik−(n−i)k)2∑i=0⌊nm⌋(⌊nm⌋i)∑j=0nmodm(nmodm...原创 2018-02-25 11:37:08 · 319 阅读 · 0 评论