数论
文章平均质量分 58
RustSword_RWZ
OI蒟蒻 嘴巴型选手233
展开
-
质因数分解QwQ 2043洛谷oj
#include<bits/stdc++.h>using namespace std;#define MN 100005int ans[MN],n,pp[MN],tot;bool ispp[MN];void OLS(){ memset(ispp,1,sizeof(ispp)); for(int i=2;i<=n;i++) { if(ispp[i])原创 2017-05-11 16:24:53 · 494 阅读 · 0 评论 -
NOI2012 随机数生成器 一道noi水题 error:爆int
题目描述栋栋最近迷上了随机算法,而随机数是生成随机算法的基础。栋栋准备使用线性同余法(Linear Congruential Method)来生成一个随机数列,这种方法需要设置四个非负整数参数m,a,c,X[0],按照下面的公式生成出一系列随机数{Xn}: X[n+1]=(aX[n]+c) mod m其中mod m表示前面的数除以m的余数。从这个式子可以看出,这原创 2017-07-24 22:56:10 · 279 阅读 · 0 评论 -
hankson的BT(bio-tech 变态)题 关于for循环边界
题目描述Hanks 博士是 BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫 Hankson。现在,刚刚放学回家的 Hankson 正在思考一个有趣的问题。今天在课堂上,老师讲解了如何求两个正整数 c1 和 c2 的最大公约数和最小公倍数。现在 Hankson 认为自己已经熟练地掌握了这些知识,他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”,这个问题是这样的:已知正整原创 2017-08-21 17:46:23 · 1323 阅读 · 0 评论 -
转自mbxc816 欧拉函数递推
欧拉函数的定义:E(k)=([1,n-1]中与n互质的整数个数). 因为任意正整数都可以唯一表示成如下形式: k=p1^a1*p2^a2*……*pi^ai;(即分解质因数形式)可以推出:E(k)=(p1-1)(p2-1)……(pi-1)*(p1^(a1-1))(p2^(a2-1))……(pi^(ai-1)) =k*(p1-1)(p2-1)……转载 2017-09-01 20:38:25 · 383 阅读 · 0 评论 -
乘法逆元
#include<bits/stdc++.h>using namespace std;int inv[30000100],n,mod;int read(){ char ch=getchar();int k=0; while(ch<'0'||ch>'9') ch=getchar(); while(ch>='0'&&ch<='9') k*=10,k+=ch-'0',ch=原创 2017-09-01 20:39:20 · 357 阅读 · 0 评论 -
关于莫比乌斯反演与除法枚举
背景某只蒟蒻叫做LittlePrincess,有一天他学会了莫比乌斯反演,但是他做题总发现自己的反演和别人的暴力是一个世界复杂度!蒟蒻十分伤心,于是钻研dalao博客。dalao:剩下的就是个n−−√n\sqrt{n}的除法分块了。这令蒟蒻十分蒙蔽,决心一探究竟——终于在万能的群友的帮助下….例题:YY的GCD我们经过GCD的历练发现我们将∑i∑j(gcd(i,j)==d)∑i∑j...原创 2018-02-22 20:34:15 · 210 阅读 · 0 评论 -
SDOI2014 数表
题目描述有一张N*m的数表,其第i行第j列(1 < =i < =n,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和。给定a,计算数表中不大于a的数之和。输入输出格式输入格式: 输入包含多组数据。 输入的第一行一个整数Q表示测试点内的数据组数,接下来Q行,每行三个整数n,m,a(|a| < =10^9)描述一组数据。输出格式: 对每...原创 2018-03-05 23:02:05 · 254 阅读 · 0 评论 -
对于莫比乌斯相关函数总结
φ(x)φ(x)φ(x)若n=pkn=pkn=p^k则φ(n)=pk−pk−1=(p−1)pk−1φ(n)=pk−pk−1=(p−1)pk−1φ(n)=p^k-p^{k-1}=(p-1)p^{k-1}积性函数 若n=pk11∗pk22∗...∗pkmmn=p1k1∗p2k2∗...∗pmkmn=p_1^{k_1}*p_2^{k_2}*...*p_m^{k_m}则有φ(n)=n∗(1−1p...原创 2018-03-31 15:41:32 · 334 阅读 · 0 评论