51nod 数论专题 (按分值排序)（持续更新）

最新推荐文章于 2023-08-09 17:02:32 发布
VIP文章 最新推荐文章于 2023-08-09 17:02:32 发布
阅读量489 收藏 2
点赞数
分类专栏： 51nod
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yczhaoxun/article/details/81842034
版权

#查找某篇题解：ctrl+f

#1 .51nod 1011

gcd,不解释

代码：

int gcd(a, b){
	return b ? gcd(b, a%b) : a;
}

#2 .51nod 1135

求最小原根

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#define mem(a,b) memset(a, b, sizeof(a))
#define LL long long

using namespace std;

const int maxn = 1000005;

int yue[maxn], tot, cnt;
int v[maxn], prime[maxn]; 

void is_prime(){
	mem(v,1);
	for(int i = 2; i <= maxn; i++){
		if(v[i]){
			prime[++cnt] = i;
			for(int j = i; j <= maxn; j += i){
				v[j] = 0;
			}
		}
	}
}

void div(int x){
	tot = 0;
	int t = (int)sqrt(1.0*x);
	for(int i = 1; prime[i] <= t; i++){
		if(x % prime[i] == 0){
			yue[++tot] = prime[i];
			while(x % prime[i] == 0) x /= prime[i];
		}
	}
	if(x > 1)
		yue[++tot] = x;
}

LL Pow(LL a, LL b, LL m){
	LL ans = 1;
	a %= m;
	while(b){
		if(b & 1)
			ans = ans * a % m;
		a = a * a % m;
		b >>= 1;
	}
	return ans;
}


int main(){
	int p;
	is_prime();
	while(scanf("%d", &p) == 1 && p){
		div(p-1);
		for(int i = 2; i <= p-1; i++){
			bool flag = 1;
			for(int j = 1; j <= tot; j++){
				int t = (p-1)/yue[j];
				if(Pow((LL)i, (LL)t, (LL)p) == 1){
					flag = 0;
					break;
				}
			}
			if(flag){
				printf("%d\n", i);
				break;
			}
		}
	}
return 0;
}

#3 .51nod 1046

快速幂

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>

#define LL long long

using namespace std;

LL Pow(LL a, LL b, LL p){
	LL ans = 1LL;
	a %= p;
	while(b){
		if(b & 1)
			ans = ans * a % p;
		a = a * a % p;
		b >>= 1;
	}
	return ans;
}

int main(){
	LL a, b, c;
	scanf("%lld%lld%lld", &a, &b, &c);
	printf("%lld\n", Pow(a, b, c));
return 0;
}

#4. 51nod 1073

约瑟夫环

代码：

#include <iostream>

using namespace std;
int main(){
	int n, k;
	cin>>n>>k;
	int s = 0;
	for(int i = 1; i <= n; i++){
		s = (s + k) % i;
	}
	cout<<s+1<<endl;
return 0;
}

#5. 51nod 1256

乘法逆元（因为N不一定是素数，所以用exgcd求解）

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>

#define LL long long

using namespace std;

void exgcd(LL a, LL b, LL &x, LL &y){
	if(b == 0){
		x = 1, y = 0;
		return;
	}
	exgcd(b, a%b, x, y);
	LL tmp = x;
	x = y;
	y = tmp - y * (a / b);
}

int main(){
	LL a, b, x, y;
	scanf("%lld%lld", &a, &b);
	exgcd(a, b, x, y);
	x = (x % b + b) % b;
	printf("%lld\n",x);
return 0;
}

#6. 51nod 1136

求欧拉函数值

代码：

#include <iostream>
#include <cmath>

using namespace std;

int phi(int n){
	int ans = n;
	for(int i = 2; i <= sqrt(n); i++){
		if(n % i == 0){
			ans = ans / i * (i-1);
			while(n % i == 0) n /= i;
		}
	}
	if(n > 1)	ans = ans / n * (n-1);
	return ans;
}

int main(){
	int n;
	cin>>n;
	cout<<phi(n)<<endl;
return 0;
}

#7 .51nod 1079

中国剩余定理

推荐博客：https://blog.csdn.net/acdreamers/article/details/8050018

https://www.cnblogs.com/MashiroSky/p/5918158.html

代码：

#include <iostream>
#include <cstdio>
#include <cmath>
#define LL long long

LL exgcd(LL a, LL b, LL &x, LL &y){
	if(!b){
		x = 1, y = 0;
		return a;
	} 
	int d = exgcd(b, a%b, x, y);
	int tmp = x;
	x = y;
	y = tmp - (a/b)*y;
	return d;
}

LL CSD(int m[], int a[], int n){
	LL x  = 0, y, d;
	LL M = 1;
	LL ans = 0;
	for(int i = 1; i <= n; i++){
		M *= m[i];
	}
	for(int i = 1; i <= n; i++){
		LL Mi = M / m[i];
		d = exgcd(Mi, m[i], x, y);
		x = x/d;
		ans = (ans + a[i]*Mi*x) % M;
	}
	return (ans + M) % M;
}

int main(){
	int n, m[15], a[15];
	scanf("%d", &n);
	for(int i = 1; i <= n; i++){
		scanf("%d %d", &m[i], &a[i]);
	}
	printf("%lld\n", CSD(m, a, n));
return 0;
}

#8. 51nod 1240
##莫比乌斯(mobius)函数
代码：

//求单个数的mobius值
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

int c[1000005];

int main(){
	int n, ans = 0, m = 0;
	scanf("%d", &n);
	int nn = n;
	for(int i = 2; i <= sqrt(nn); i++){
		if(n % i == 0){
			c[++m] = 0;
			while(n % i == 0){
				n /= i;
				c[m]++;
			}
		}
	}	
	if(n > 1){
		c[++m] = 1;
	}
	for(int i = 1; i <= m; i++){
		if(c[i] >= 2){
			printf("0\n");
			return 0;
		}
	}
	if(m & 1){
		printf("-1\n");
	}
	else{
		printf("1\n");
	}
return 0;
}

//求2到n的mobius值
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

int miu[1000005];
int v[1000005];
int n;

int main(){
	scanf("%d", &n);
	for(int i = 2; i <= n; i++) miu[i] = 1, v[i] = 0;
	for(int i = 2; i <= n; i++){
		if(v[i]) continue;
		miu[i] = -1;
		for(int j = 2*i; j <= n; j += i){
			v[j] = 1;
			if( (j / i) % i == 0 ) miu[j] = 0;
			else miu[j] *= -1;
		}
	}
	for(int i = 2; i <= n; i++){
		printf("%d ", miu[i]);
	}
return 0;
}

#9. 51nod 1106

质数检测, 水题

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;

int n, t;

bool is_prime(int x){
	for(int i = 2; i <= sqrt(x); i++){
		if(x % i == 0){
			return false;
		}
	}
	return true;
}

int main(){
	scanf("%d", &n);
	for(int i = 1; i <= n; i++){
		scanf("%d", &t);
		if(is_prime(t))
			printf("Yes\n");
		else printf("No\n");
	}
return 0;
}

#10. 51nod 1012
##lcm 算了还是不一个一个做了，这种水题。。。还是二分一下适合我的题。

#include <iostream>
#include <cstdio>
#include <cmath>
#define LL long long

using namespace std;

LL gcd(LL a, LL b){
	return b ? gcd(b, a%b) : a;
}



int main(){
	LL a, b;
	scanf("%lld %lld", &a, &b);
	printf("%lld\n", a/gcd(a ,b)*b);
return 0;
}

#11. 51nod 1181
题意：如果一个质数，在质数列表中的编号也是质数，那么就称之为质数中的质数。例如：3 5分别是排第2和第3的质数，所以他们是质数中的质数。现在给出一个数N，求>=N的最小的质数中的质数是多少？
题解：N范围是1e6，考虑用线性筛筛出1e7以下的素数，再把下标是素数的另存一个数组，二分查找答案即可。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;

const int maxn = (int)(1e7+100);


int p[maxn], v[maxn];
int pp[maxn];
int m, t;

void pre(){
	memset(v, 0, sizeof(v));
	for(int i = 2; i <= maxn; i++){
		if(v[i] == 0){
			v[i] = i;
			p[++m] = i;
		}
		for(int j = 1; j <= m; j++){
			if(v[i] < p[j] || p[j] > maxn/i) break;
			v[i*p[j]] = p[j];
		}
	}
	for(int i = 1; i <= m; i++){
		if(v[i] == i){
			pp[++t] = p[i];
			//printf("%d ", i);
		}
	}
}

int main(){
	pre();
	int n;
	scanf("%d", &n);
	int x = lower_bound(pp+1, pp+t+1, n) - pp;
	printf("%d\n", pp[x]);
return 0;
}

#12. 51nod 1080
题意：给出一个整数N，将N表示为2个整数i与j的平方之和（i <= j)，如果有多种表示，按照i的递增序输出。
例如：N = 130，130 = 3^2 + 11^2 = 7^2 + 92（注：32 + 112同112 + 3^2算1种）
题解：i最大sqrt(n/2),暴力枚举即可。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;

int v[40005];

int main(){
	int n, t, flag = 1;
	memset(v, 0, sizeof(v));
	scanf("%d", &n);
	t = (int)sqrt(n/2);
	for(int i = 0; i <= t; i++){
		int x = n - i*i;
		int y = (int)sqrt(x);
		if(y*y == x && !v[y]){
			printf("%d %d\n", i, y);
			v[i] = 1;
			flag = 0;
		}
	}
	if(flag){
		printf("No Solution\n");
	}
return 0;
}

#13. 51nod 1010
题意： K的因子中只包含2 3 5。满足条件的前10个数是：2,3,4,5,6,8,9,10,12,15。
所有这样的K组成了一个序列S，现在给出一个数n，求S中 >= 给定数的最小的数。
例如：n = 13，S中 >= 13的最小的数是15，所以输出15。
题解： 预处理出1e18内的这样的数，二分查找即可。

#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <algorithm>
#define LL long long
using namespace std;

const LL maxn = (LL)1e18+100LL;

int T;
LL n;
LL num[1000005], tot;

void pre(){
	tot = 0;
	for(LL i = 1; i <= maxn; i *= 2LL){
		for(LL j = 1; i*j <= maxn; j *= 3LL){
			for(LL k = 1; i*j*k <= maxn; k *= 5LL){
				if(i*j*k != 1)
					num[++tot] = i*j*k;
			}
		}
	}
	sort(num+1, num+tot+1);
}


int main(){
	pre();
	scanf("%d", &T);
	while(T--){
		scanf("%lld", &n);
		int x = lower_bound(num+1, num+1+tot, n) - num;
		printf("%lld\n", num[x]);
	}
return 0;
}

#14. 51nod 1126
**题意：**有一个序列是这样定义的：f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.
给出A，B和N，求f(n)的值。
题解：：矩阵快速幂，转移矩阵为 [ 0 B 1 A ] \begin{bmatrix} 0&amp;B\\1&amp;A\\ \end{bmatrix} [01​BA​]，不过要注意负数要加mod再mod。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>

using namespace std;

const int mod = 7;

int A, B, n;

void mul(int f[2], int a[2][2]){
	int c[2];
	memset(c, 0, sizeof(c));
	for(int i = 0; i < 2; i++){
		for(int j = 0; j < 2; j++){
			c[i] = ( c[i] + (long long)f[j] * a[j][i] ) % mod;
		}
	}
	memcpy(f, c, sizeof(c));
}

void mulself(int a[2][2]){
	int c[2][2];
	memset(c, 0, sizeof(c));
	for(int i = 0; i < 2; i++){
		for(int j = 0; j < 2; j++){
			for(int k = 0; k < 2; k++){
				c[i][j] = ( c[i][j] + (long long)a[i][k] * a[k][j] ) % mod;
			}
		}
	}
	memcpy(a, c, sizeof(c));
}

int main(){
	scanf("%d%d%d",&A,&B,&n);
	int f[2] = {1, 1};
	int a[2][2];
	a[0][0]= 0, a[1][0] = 1;
	a[0][1] = B, a[1][1] = A;
	n -= 1;
	for(; n; n >>= 1){
		if(n & 1) mul(f, a);
		mulself(a);
	}
	printf("%d\n", (f[0]+mod) % mod);
return 0;
}

#15. 51nod 1014
题意： X*X mod P = A，其中P为质数。给出P和A，求<=P的所有X。
**题解：**枚举x即可。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>

#define LL long long

int main(){
	int p, a, flag = 0;
	scanf("%d%d", &p, &a);
	for(int i = 0; i <= p; i++){
		LL x = (LL)i*i;
		if(x % p == a){
			printf("%d ", i);
			flag = 1;
		}
	}
	if(!flag)
		printf("No Solution\n");
return 0;
}

#16. 51nod 1352
**题意：**给出N个固定集合{1，N},{2,N-1},{3,N-2},…,{N-1,2},{N,1}.求出有多少个集合满足：第一个元素是A的倍数且第二个元素是B的倍数。提示：对于第二组测试数据，集合分别是：{1,10},{2,9},{3,8},{4,7},{5,6},{6,5},{7,4},{8,3},{9,2},{10,1}.满足条件的是第2个和第8个。
题解：
1.设满足条件的是 ( t 1 , t 2 ) (t1, t2) (t1,t2), 即 ( A x , B y ) (Ax, By) (Ax,By),所以有 A x + B y = n + 1 Ax+By=n+1 Ax+By=n+1,转化为求这个方程的满足 1 &lt; = A x &lt; = n , 1 &lt; = B y &lt; = n 1&lt;=Ax&lt;=n, 1&lt;=By&lt;=n 1<=Ax<=n,1<=By<=n的解的个数。
2.先用扩展欧几里得求出 A x + B y = g c d ( A , B ) Ax+By=gcd(A, B) Ax+By=gcd(A,B)的一个解 x 0 x_0 x0​,那么原方程的一个解就是 x 0 ∗ n + 1 d x_0*\frac {n+1}d x0​∗dn+1​.
3.公式变形： A x + B y = n + 1 Ax+By=n+1

最低0.47元/天 解锁文章
elorole
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
最长高地（51Nod-2509）.rar
09-16
最长高地（51Nod-2509）.rar
51nod的题目标答。
06-07
51nod的题目标答，大家喜欢的可以下载看看，但记住不要复制！！！
NOD 32离线更新方法和账号
10-08
NOD 32离线更新方法和账号,NOD 32离线更新方法和账号,
最长配对（51Nod-2494）.rar
09-16
最长配对（51Nod-2494）.rar
51nod:51nod.com的代码
05-03
51点 51nod.com的代码
51nod初学者首页1级题题解（5）
m0_57362922的博客
09-03 677
老师出完题后，全班同学都在埋头计算，小高斯却很快算出答案等于 5050。1～100正好可以分成这样的50对数，每对数的和都相等。输入一个数n，请你输出n-10的值。给一个正整数n，计算它最多能被2的多少次幂整除。i是计数器,每次加1，当不能被2整除时，就输出。（2）将这个数对10取余，得到的就是个位的数字。给定一个整数n，求它的阶乘，n小于等于10。输入一个数n，请你输出 n 乘以 3 的值。>>1和/2相等，输出时外面的括号不能去掉。输入一个数n，请你输出n+2的值。也可以这样，只是效率低了一点点。
51nod动态规划教程
旺旺的碎冰冰~的博客
03-05 503
51nod动态规划教程 最后，我们来提供输入输出数据，由你来写一段程序，实现这个算法，只有写出了正确的程序，才能继续后面的课程。 输入 第1行：字符串A 第2行：字符串B (A,B的长度 &amp;lt;= 1000) 输出 输出最长的子序列，如果有多个，随意输出1个。 输入示例 abcicba abdkscab 输出示例 ...
51nod 1014 暴力，细节，模拟
HPU WIN
08-24 151
题意就不在赘述了，唯一一个需要注意的一点是对于for循环中的x要用long long，虽然给予的条件是x&lt;=1e6,但是你在模拟时要用x*x,这个结果int 是无法保存到因为是1e12的量级，所以就算你是用longlong sum=x*x，也是不行的，因为计算器是先计算出x*x然后按照你给定的类型进行保存，然后再输给sum，所以结果会爆，这一点值得写篇博客来提醒自己。 #include&...
51nod - 1494 选举拉票（Java & JS & Python）
最新发布
伏城之外的博客
08-09 1189
0 ≤ bi ≤ 10^4），表示第 i 个选民选的是第 ai 号候选人，想要让他选择自己就要花 bi 的钱。你是 0 号候选人（所以，如果一个选民选你的话 ai 就是 0，这个时候 bi 也肯定是 0 ）。你已经知道每一个人要选的人是谁，以及要花多少钱才能让这个人选你。现在你想要花最少的钱使得你当上县长。你当选的条件是你的票数比任何一个其它候选人的多（严格的多，不能和他们中最多的相等）。第一行有一个整数 n（1 ≤ n ≤ 10^5），表示这个县的选民数目。输出一个整数表示花费的最少的钱。
51nod1517 递推、平方根与误差
ZCH1901的博客
08-18 675
1517 递推、平方根与误差 Noder遇到这样一个递推式： F(n)={1(n<0)F(n−1)+F(n−a−−√)(0<=n) 这个递推中存在2个未知量：n和a，给出这两个变量的值，便可以计算F(n)的值。 在此之前，Noder学习了一些简单的线性递推知识，他希望利用已经掌握的知识，来解决这个问题。经过一段时间思考，找到了这样一个思路： 定义G(n)=G(n-x)+G(n-y)（x，y为正整数），并且F(n)=G(x*n)，然后就转为了他所熟悉的线性递推，FFT、多项式求Mod等.
51Nod1018 排序（C语言）
sinat_39591298的博客
08-17 326
给出N个整数，对着N个整数进行排序 Input 第1行：整数的数量N（1 <= N <= 50000) 第2 - N + 1行：待排序的整数（-10^9  Output 共n行，按照递增序输出排序好的数据。 Input示例 5 5 4 3 2 1 Output示例 1 2 3 4 5 C语言AC代码 #include
51nod初学者首页1级题题解（3）
m0_57362922的博客
08-16 480
现在分别给出上述四种所购买的斤数（均不超过 20），请你编写程序帮助售货员阿姨计算并依次输出顾客购买四种水果需要的钱数及总钱数。由于小瓜非常挑剔，他只吃重大于等于m克的龙虾，于是他打算把所有重大于等于m克的龙虾买回家。给你两个不同的字符串，如果第一个字符串的字典序小于第二个字符串，则输出YES，如果第一个字符串的字典序大于第二个字符串，则输出NO。给你两个字符串，一个母串，一个子串，请你找出子串第一次在母串中出现的位置。遍历字符串，当字符串的某一个字符等于x时，输出i，否则在循环下面输出-1....
51nod 2658 最多分成多少块 V2 思维
pxlsdz的博客
08-13 399
2658 最多分成多少块 V2 1.0 秒 131,072.0 KB 10 分 2级题 小b有个长度为n的数组a，她想将这个数组排序。 然而小b很懒，她觉得对整个数组排序太累了，因此她请你将a分成一些块，使得她只需要对每一块分别排序，就能将整个数组排序。 请问你最多能把a分成多少块。 保证a为0...n-1的一个排列。 样例解释： 将a分成2块或者更多块，都无法...
约瑟夫环（51nod）
nowting的博客
08-27 421
先理解下什么是约瑟夫问题：https://blog.csdn.net/tingyun_say/article/details/52343897 以及解决方法：https://blog.csdn.net/wusuopubupt/article/details/18214999 题目：n个人坐成一个环（编号0~n-1），从第1个人开始报数，数到q的人出列，后面的人重新从1开始报数。问最后剩下的人的...
51nod 1240 莫比乌斯函数
AK_HuangYC的博客
10-29 1153
莫比乌斯函数，由德国数学家和天文学家莫比乌斯提出。梅滕斯(Mertens)首先使用μ(n)（miu(n)）作为莫比乌斯函数的记号。（据说，高斯(Gauss)比莫比乌斯早三十年就曾考虑过这个函数）。 具体定义如下： 如果一个数包含平方因子，那么miu(n) = 0。例如：miu(4), miu(12), miu(18) = 0。 如果一个数不包含平方因子，并且有k个不同的质因子，那么miu
码力恢复——51NOD基础篇
辞树
02-28 383
1000 A + B #include<bits/stdc++.h> using namespace std; #define clr(a) memset(a, 0, sizeof(a)) #define line cout << "--------" << endl typedef long long ll; const int maxn = 1e5 + 10; int main(){ ll a, b; cin >> a >
51nod 2583 数论只会Gcd
t56346的博客
08-20 229
题目描述： 有这样一段两两求最大公约数的程序CoGcd， int Gcd(int x, int y) { if(y == 0) return x; return Gcd(y, x % y); } void CoGcd(int m) { for(int i = 1; i <= m; i++) { for(int j = 1; j <= m; j++) { Gcd(i, j);
51nod 1964 陵陵曾玩的数论题
t56346的博客
08-21 130
题目描述： 腾腾是一名数论大师，经常在51nod上切高分数论题。凭借着高超的实力，拿到了51nod的#2。 陵陵是腾腾的好朋友。2017年1月，他还在苦恼一直都不会刷>=320分数论题。一天，他拿到了栋栋寄来的一道数论题。不幸的是，陵陵看错了这道数论题的题意，他感觉这太难了，于是只能拿去问腾腾：给定一个长度为n的序列a，对于一个区间[l,r]，如果al，al+1……ar的最小公倍数记作lcm，最大值记作max，这个区间[l,r]的价值定义为lcmmaxlcmmax（即lcm的max次幂）。...
51nod参考答案
Pierce110110的博客
12-28 4962
算术运算 2393 取余 #include <stdio.h> main() { int x,y; scanf("%d%d",&x,&y); printf("%d",x%y); } 2394 四则运算1 #include <stdio.h> int mai...
51nod3204分数运算
07-09
对于题目51nod3204，这是一个分数运算的问题。很抱歉，由于我无法直接访问外部网站或查找特定问题的详细答案，我不能为你提供该问题的具体解答。然而，我可以为您提供一般的分数运算知识和算法，以帮助您解决该问题...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值