摘要：Miller Rabin 算法 参考：【朝夕的ACM笔记】数论-Miller Rabin素数判定 - ~(o°ω°o) 作用在于判断素数，它有前置两个定理 费马小定理 Theorem 设 \(p\) 是一个素数，\(a\in \Z^+\) 且不是 \(p\) 的倍数，那么有 \(a^{p-1}\eq 阅读全文

摘要：Fibonacci 常见性质 参考：斐波那契数列的性质 - Milkor - 博客园 定义 \(f[1]=f[2]=1\) \(f[i]=f[i-1]+f[i-2]\) 可以构造矩阵： 和矩阵： 二者乘积为： 所以可以矩阵快速幂 \(log(n)\) 求 \(f[n]\) 性质 \(\gcd(f[i 阅读全文

摘要：最大公约数详解 一般的，设 \(a_1,a_2,...a_n\) ，是 \(n\) 个非零整数，如果存在一个非零整数 \(d\), 使得 \(d\mid a_1,d\mid a_2,...d\mid a_n\) ，那么称 \(d\) 是这 \(n\) 个数的公约数。显然可能存在多个公约数，将这些公约 阅读全文

摘要：@ 前言 彩笔作者最开始以为是真分块。。。 作者接触到这玩意是因为这道题：luogu2261 [CQOI2007]余数求和 题目简洁明了，一看就是推式子的题： 给出正整数 \(n\) 和 \(k\)，请计算 \(G(n, k) = \sum_{i = 1}^n k \bmod i\) 其中 \(k\ 阅读全文

摘要：#include<bits/stdc++.h> using namespace std; //试除法 bool is_prime(int n){ if(n<2)return false; int temp=sqrt(n); for(int i=2;i<=temp;i++){ if(n%i==0)re 阅读全文