随笔分类 - 数论
摘要：ll ai[maxn], bi[maxn]; ll mul(ll a, ll b, ll mod) { ll res = 0; while (b > 0) { if (b & 1) res = (res + a) % mod; a = (a + a) % mod; b >>= 1; } return
阅读全文
摘要：使用JAVA内置函数 /* import java.math.BigInteger; import java.io.BufferedReader; import java.io.IOException; import java.io.InputStreamReader; import java.io
阅读全文
摘要：设m是正整数，a是整数，若a模m的阶等于φ(m)，则称a为模m的一个原根。（其中φ(m)表示m的欧拉函数） 给出1个质数P，找出P最小的原根。 ll m; ll fac[maxn], cnt; ll quickPower(ll a, ll b, ll M) { ll ans = 1ll; ll ba
阅读全文
摘要：给出一个n，求1-n这n个数，同n的最小公倍数的和。例如：n = 6，1,2,3,4,5,6 同6的最小公倍数分别为6,6,6,12,30,6，加在一起 = 66。 由于结果很大，输出Mod 1000000007的结果 ai <= 1e9 公式推导如上。 最后只需要预处理出质因子，用积性函数的性质在
阅读全文
摘要：B签到成功 打表出c的存在没什么作用。 容易找出规律 答案就是c的某个次数，这个次数就是n的质因子的幂次和。 因此只需线性筛出幂次和，再快速幂就好了 #pragma warning(disable:4996) #include<iostream> #include<algorithm> #inclu
阅读全文
摘要：论文题 1e11内找素数个数。 Meissel-Lehmer 黑科技板子 #pragma warning(disable:4996) #include<iostream> #include<algorithm> #include<bitset> #include<tuple> #include<un
阅读全文
摘要：给定一个n表示a的质因数个数。 接下来n行给出质数及其指数。 按要求输出其因数，满足如下要求： 当前数是前一个数通过乘一个质数或者除以一个质数得到。 反正就是构造嘛。对于每一个i，必然要遍历前面的所有情况。其实就是升升降降。 vector<ll> ans; vector<ll> p(16); vec
阅读全文
摘要：int prime[maxn], prime_tot; int is_prime[maxn]; int mu[maxn]; void pre_calc(int lim) { mu[1] = 1; for (int i = 2; i <= lim; i++) { if (!is_prime[i]) {
阅读全文
摘要：给定整数 n，求它的所有原根。 为了减小你的输出量，给出输出参数 dd，设 nn 的所有原根有 cc 个，从小到大分别为 g_1,\ldots,g_cg1,…,gc，你只需要依次输出 int phi[maxn]; int prime[maxn]; int vis[maxn]; int rt[ma
阅读全文
摘要：非常好的数论题 结合了很多不知道的知识点 易得代码 #include<iostream> #include<unordered_map> #include<algorithm> #include<string> #include<cmath> #include<cstring> #include<v
阅读全文
摘要：第一问： 我们知道10的n次位数是n+1，2的p次等于10的log(10)2p ，所以结果就是log(10)2p+1 ，直接用库函数即可 第二问 考虑用大数快速幂对最后500位计算答案即可 此题学到的知识点 java大数取模，pow，快速幂模板 import java.util.*; import
阅读全文
摘要：ll f(ll a,ll b,ll c,ll n) ///用于求解$Sigma_{i=0}^{n} floor((a*i+b)/c)$ { ll m = (a*n+b)/c; if(n==0||m==0) return (b/c); if(n==1) return ((b/c)+((a+b)/c))
阅读全文
摘要：java初学 借鉴他人代码 求R的n次，R是小数 import java.util.*; import java.math.*; public class Main { public static void main(String[] args) { Scanner in=new Scanner(S
阅读全文
摘要：莫比乌斯函数 求单个莫比乌斯函数 O（sqrt(n)） ll getu(ll x) { ll v = 1; for (int i = 2; i * i <= x; i++) { if (x % i == 0) { v = -v, x /= i; if (x % i == 0) return 0; }
阅读全文
摘要：洛谷P1495 https://www.luogu.com.cn/problem/P1495 #include<iostream> #include<string> #include<cmath> #include<cstring> #include<vector> #include<map> #i
阅读全文
摘要：#include<iostream> #include<string> #include<cmath> #include<cstring> #include<vector> #include<map> #include<set> #include<algorithm> #include<queue>
阅读全文
摘要：k个数的公约数的含义是这k个数中均含有某个因数， 因此我们只需把所有数的因数求出来，发现有k个数均含有某个因数，那么这个数必然是这k个数的公约数 因此此题的思路就是求出每个数的因子，并存储每个因子出现的次数。最后搜索这个数组就行 题目还保证了单调性 #include<iostream> #inclu
阅读全文
摘要：求因子数一定的最小数(反素数) #include<iostream> #include<string> #include<cmath> #include<cstring> #include<vector> #include<map> #include<set> #include<algorithm>
阅读全文
摘要：算不出的等式 BJOI2012 看到这题 真没什么办法 无奈看题解 1.注意到p/q 联想到斜率 2.注意到 [ ] 联想到整点 注意到k在变化，构造一次函数 f(x)=p/q*x ,g(x)=q/p*x 收到【】 的影响，y值即为f(x)下取整后的值，即垂线上整点的个数 又考虑到p==q时 需特判
阅读全文
摘要：符号说明 a|b a整除b (a,b) a与b的最大公因数 [a,b] a与b的最小公倍数 pα||a pα|a但pα+1∤a a≡b(mod m) a与b对模m同余 a-1 (mod m） a对模m的数论倒数 性质1 如果a|b，那么(-a)|b,反过来也成立 性质2 如果a|b,b|c,那么a|
阅读全文
浙公网安备 33010602011771号