随笔分类 - 数论，数学
摘要："题目链接" "高斯消元详解"
阅读全文
摘要："题目链接" 扩展中国剩余定理： "1(直观的)" 、 "2(详细证明)" 。 [Upd：] https://www.luogu.org/problemnew/solution/P4774 cpp include include define gc() getchar() typedef long
阅读全文
摘要："题目链接" "扩展Lucas" cpp //求C_n^k%m include typedef long long LL; LL FP(LL x,LL k,LL p) { LL t=1ll; for(; k; k =1,x=x x%p) if(k&1) t=t x%p; return t; } vo
阅读全文
摘要：题目链接 Lucas定理 日常水题...sublime和C++字体死活不同步怎么办... 复杂度： 如果能$O(p)$预处理，复杂度为$O(\log_pn)$，否则要快速幂为$O(\log_pn*\log p)$。 当$p$不是质数时，用扩展Lucas，要加CRT，复杂度看起来比较麻烦。。 //想错
阅读全文
摘要："题目链接" 颓了一天 写个模板吧。。 Chinese_Remainder_Theorem: "MashiroSky" 、 "远航之曲" cpp include include define gc() getchar() typedef long long LL; const int N=13; L
阅读全文
摘要：题目链接 Pollard_Rho:http://blog.csdn.net/thy_asdf/article/details/51347390 https://zhuanlan.zhihu.com/p/267884783 复杂度：\(O(n^{\frac14})\)，即可以跑$long\ long$
阅读全文
摘要：题目链接： "洛谷" 、 "BZOJ2179" cpp //将乘数拆成 a0 10^n + a1 10^(n 1) + ... + a_n 1的形式 //可以发现多项式乘法就模拟了竖式乘法 所以用FFT即可 注意处理进位 //n位 n位最多就只有2n位了 //putchar的速度。。还是快的 inc
阅读全文
摘要："题目链接" $Description$ 求$A_0,A_1,A_2,\cdots,A_{n 1}$，满足 $$A_0 1^0+A_1 1^1+\ldots+A_{n 1} 1^{n 1}\equiv B "1" $$ $$A_0 2^0+A_1 2^1+\ldots+A_{n 1} 2^{n 1}
阅读全文
摘要："题目链接" 显然我们需要使每个i满足$$( ∑_{j} X[j] A[i][j] ) mod\ 2 = B[i]$$ 求这个方程自由元Xi的个数ans，那么方案数便是$2^{ans}$ %2可以用^代替，不难看出 B[i]=st[i]^ed[i] 如果X[j]=1，假设j会影响i，那么X[j] A
阅读全文
摘要：题目链接： "洛谷" 、 "LOJ" . FFT相关： "快速傅里叶变换(FFT)详解" 、 "FFT总结" 、 "从多项式乘法到快速傅里叶变换" . 5.4 又看了一遍， "这个" 也不错。 2019.3.7 叕看了一遍，推荐 "这个" 。 cpp include include include
阅读全文
摘要："题目链接" cpp include include include const int N=105; const double eps=1e 10; int n; inline bool bigger(double a,double b) {return std::fabs(a) std::fab
阅读全文
摘要："https://www.zybuluo.com/SovietPower/note/1016329" "参考原文"
阅读全文
摘要："题目链接"
阅读全文
摘要："题目链接" $Description$ $n\leq 10^{10}$，求 $$\sum_{i=1}^n\sum_{j=1}^ngcd(i,j)\ mod\ (1e9+7)$$ $Solution$ 首先 $$\sum_{i=1}^n\sum_{j=1}^ngcd(i,j)=\sum_{d=1}^
阅读全文
摘要：$Description$ 给定p， $Solution$ 欧拉定理：$若(a,p)=1$,则$a^b\equiv a^{b\%\varphi(p)}(mod\ p)$. 扩展欧拉定理：$a^b\equiv a^{b\%\varphi(p)+\varphi(p)}(mod\ p)$ (a为任意整数，
阅读全文
摘要：\(Description\) \(n^2-3n+2=\sum_{d|n}f(d)\) 求$\sum_{i=1}nf(i)\ mod \ 109+7$. \(Solution\) 参考. 设$g(n)=n^2-3n+2$，那么$f*I=g$。 如果狄利克雷卷积中左边的一个函数是待求前缀和的，而另外两
阅读全文
摘要：$Description$ $$g(i)=\sum_{i_1|i}\sum_{i_2|i_1}\sum_{i_3|i_2}\cdots\sum_{i_k|i_{k 1}}f(i_k)\ mod\ 1000000007$$ 给出$n,k,f[1\sim n]$，求$g[1\sim n]$. $Solu
阅读全文
摘要："题目链接" map： cpp //杜教筛 include include typedef long long LL; const int N=5e6; int mu[N+3],P[N+3],cnt; bool Not_P[N+3]; std::map sum; //std::map::iterat
阅读全文
摘要："题目链接" cpp / http://www.spoj.com/problems/VLATTICE/ 题意：求一个n n n的晶体，有多少点可以在(0,0,0)处可以直接看到。 同BZOJ.2301 题目即要求gcd(i,j,k)=1的(i,j,k)数对个数，1 include define gc
阅读全文
摘要：~~ [Update] 我好像现在都看不懂我当时在写什么了= =~~ $Description$ 求$\sum_{i=a}^b\sum_{j=c}^d[(i,j)=k]$ $Solution$ 首先是把下界作为1.可以化为求 $$\sum_{i=1}^{\lfloor\frac{N}{k}\rflo
阅读全文
浙公网安备 33010602011771号