C++复习（十五）——素数幻方

最新推荐文章于 2024-04-30 13:01:01 发布
最新推荐文章于 2024-04-30 13:01:01 发布
阅读量1k 收藏
点赞数 1
分类专栏： C++ 文章标签： C++ 素数幻方
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/QJM1995/article/details/87987957
版权

 

问题三十一：素数幻方 ，在一个4*4矩阵中，每一行，每一列，两条对角线上的四个数字组成的四位数均为可逆素数，求所有矩阵

/*
题目三十一：素数幻方 ，在一个4*4矩阵中，每一行，每一列，两条对角线上的四个数字组成的四位数均为可逆素数，求所有矩阵

分析：可逆素数可以判断，但是如果4*4方阵全部穷举，所需要的计算量极大
短时间内根本无法得出结果，需要化简问题 ，因为行数列数对角线均为可逆
素数，所以可以先求出所有的四位可逆素数进行穷举，计算量大大减少，因为是素数
所以每一行的末尾，每一列的末尾不可出现2的倍数和5倍数的数字，进一步缩减穷举
数据范围，在进行穷举则可得结果 
1、第一行，第一列不能有零
2、最后一行，最后一列不能有2,5的倍数
3、四位可逆素数
4、测试前两行 
5、测试前三行
6、输出结果 

*/
#include<iostream>
#include<cmath>
#define MAX 500
using namespace std;

int all_sum;
int last_all_sum;
int first_all_sum;

int ke_ni(int n){//返回逆数 
	int i,num = n,total = 0;
	int k,now = 1;
	for(i = 0;i<4;i++){
		for(k = 3;k>i;k--){
			now *= 10;
		}
		total += num%10*now;
		num /= 10;
		now = 1;
	}
	return total;
}

int su_shu(int n){//判断是否为素数 
	int i,now;
	if(n==1||n==0)return 0;
	if(n==2)return 1; 
	now = (int)sqrt((double)n) + 1;
	//now = n>>1;
	for(i = 2;i<=now;i++){
		if(n%i == 0)return 0;
	}
	return 1;
}

int get_n(int num,int i){//返回数n的第i位值 
	int last;
	for(;i > 1;i--){
		num /= 10;
	}
	last = num%10;
	return last;
}

int last_row_len(int n){
	int i;
	for(i=1;i<=4;i++){
		if(get_n(n,i)%2 == 0||get_n(n,i)%5 == 0){
			return 0;
		}
	}
	return 1;
}

int first_row_len(int n){
	int i;
	for(i=1;i<=4;i++){
		if(get_n(n,i) == 0){
			return 0;
		}
	}
	return 1;
}

int yes(int *all,int *first,int *last,int a,int b,int c,int d){
	int i;
	int flag = 0;
	for(i = 0;i<first_all_sum;i++){
		if(first[i] == get_n(a,4)*1000+get_n(b,4)*100+get_n(c,4)*10+get_n(d,4)){
			flag++;
			break;
		}
	}
	if(flag != 1)return 0;
	for(i = 0;i<all_sum;i++){
		if(all[i] == get_n(a,3)*1000+get_n(b,3)*100+get_n(c,3)*10+get_n(d,3)){
			flag++;
			break;
		}
	}
	if(flag != 2)return 0;
	for(i = 0;i<all_sum;i++){
		if(all[i] == get_n(a,2)*1000+get_n(b,2)*100+get_n(c,2)*10+get_n(d,2)){
			flag++;
			break;
		}
	}
	if(flag != 3)return 0;
	for(i = 0;i<last_all_sum;i++){
		if(last[i] == get_n(a,1)*1000+get_n(b,1)*100+get_n(c,1)*10+get_n(d,1)){
			flag++;
			break;
		}
	}
	if(flag != 4)return 0;
	for(i = 0;i<all_sum;i++){
		if(all[i] == get_n(a,4)*1000+get_n(b,3)*100+get_n(c,2)*10+get_n(d,1)){
			flag++;
			break;
		}
	}
	if(flag != 5)return 0;
	for(i = 0;i<all_sum;i++){
		if(all[i] == get_n(a,1)*1000+get_n(b,2)*100+get_n(c,3)*10+get_n(d,4)){
			flag++;
			break;
		}
	}
	if(flag != 6)return 0;
	return 1;
}

int yes_two(int *all,int *first,int *last,int a,int b){
	int i;
	int flag = 0;
	for(i = 0;i<first_all_sum;i++){
		if(get_n(first[i],4)*10+get_n(first[i],3) == get_n(a,4)*10+get_n(b,4)){
			flag++;
			break;
		}
	}
	if(flag != 1)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*10+get_n(all[i],3) == get_n(a,3)*10+get_n(b,3)){
			flag++;
			break;
		}
	}
	if(flag != 2)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*10+get_n(all[i],3) == get_n(a,2)*10+get_n(b,2)){
			flag++;
			break;
		}
	}
	if(flag != 3)return 0;
	for(i = 0;i<last_all_sum;i++){
		if(get_n(last[i],4)*10+get_n(last[i],3) == get_n(a,1)*10+get_n(b,1)){
			flag++;
			break;
		}
	}
	if(flag != 4)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*10+get_n(all[i],3) == get_n(a,4)*10+get_n(b,3)){
			flag++;
			break;
		}
	}
	if(flag != 5)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*10+get_n(all[i],3) == get_n(a,1)*10+get_n(b,2)){
			flag++;
			break;
		}
	}
	if(flag != 6)return 0;
	return 1;
} 

int yes_three(int *all,int *first,int *last,int a,int b,int c){
	int i;
	int flag = 0;
	for(i = 0;i<first_all_sum;i++){
		if(get_n(first[i],4)*100+get_n(first[i],3)*10+get_n(first[i],2) == get_n(a,4)*100+get_n(b,4)*10+get_n(c,4)){
			flag++;
			break;
		}
	}
	if(flag != 1)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*100+get_n(all[i],3)*10+get_n(all[i],2) == get_n(a,3)*100+get_n(b,3)*10+get_n(c,3)){
			flag++;
			break;
		}
	}
	if(flag != 2)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*100+get_n(all[i],3)*10+get_n(all[i],2) == get_n(a,2)*100+get_n(b,2)*10+get_n(c,2)){
			flag++;
			break;
		}
	}
	if(flag != 3)return 0;
	for(i = 0;i<last_all_sum;i++){
		if(get_n(last[i],4)*100+get_n(last[i],3)*10+get_n(last[i],2) == get_n(a,1)*100+get_n(b,1)*10+get_n(c,1)){
			flag++;
			break;
		}
	}
	if(flag != 4)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*100+get_n(all[i],3)*10+get_n(all[i],2) == get_n(a,4)*100+get_n(b,3)*10+get_n(c,2)){
			flag++;
			break;
		}
	}
	if(flag != 5)return 0;
	for(i = 0;i<all_sum;i++){
		if(get_n(all[i],4)*100+get_n(all[i],3)*10+get_n(all[i],2) == get_n(a,1)*100+get_n(b,2)*10+get_n(c,3)){
			flag++;
			break;
		}
	}
	if(flag != 6)return 0;
	return 1;
}

int main(){
	int i,count=0,j = 0;
	int all_su[MAX];
	for(i=1000;i<=9999;i++){
		if(su_shu(i)==1&&su_shu(ke_ni(i))==1){
			count++;
			all_su[j] = i;
			j++;
			//cout<<i<<" ";
		}
	}
	all_sum = j;
	int last[MAX],k=0;
	for(i=0;i<all_sum;i++){
		if(last_row_len(all_su[i]) == 1){
			last[k] = all_su[i];
			//cout<<last[k]<<" ";
			k++;
		}
	}
	last_all_sum = k;
	k = 0;
	int first[MAX];
	for(i=0;i<all_sum;i++){
		if(first_row_len(all_su[i]) == 1){
			first[k] = all_su[i];
			//cout<<first[k]<<" ";
			k++;
		}
	}
	first_all_sum = k;
	//方法一，速度有所提升，但仍然很慢 
	/*
	int all[MAX][MAX];
	int x;
	int a=0;
	for(i = 0;i<first_all_sum;i++){
		for(j = 0;j<all_sum;j++){
			for(k = 0;k<all_sum;k++){
				for(x = 0;x<last_all_sum;x++){
					if(yes(all_su,first,last,first[i],all_su[j],all_su[k],last[x])==1){
						cout<<first[i]<<'\n'<<all_su[j]<<'\n'<<all_su[k]<<'\n'<<last[x]<<'\n'<<endl;
						a++;
					}
				}
			}
		}
	}
	cout<<a<<endl;
	*/
	
	//方法二，速度极大提升 
	int x;
	int a=0;
	for(i = 0;i<first_all_sum;i++){
		for(j = 0;j<all_sum;j++){
			if(yes_two(all_su,first,last,first[i],all_su[j])==0)continue;
			for(k = 0;k<all_sum;k++){
				if(yes_three(all_su,first,last,first[i],all_su[j],all_su[k])==0)continue;
				for(x = 0;x<last_all_sum;x++){
					if(yes(all_su,first,last,first[i],all_su[j],all_su[k],last[x])==1){
						cout<<first[i]<<'\n'<<all_su[j]<<'\n'<<all_su[k]<<'\n'<<last[x]<<'\n'<<endl;
						a++;
					}
				}
			}
		}
	}
	cout<<a<<endl;
	return 0;
}

       这其中有不合适或者不正确的地方欢迎指正，我的QQ号码：2867221444（乔金明），谢谢，也可以相互交流下，备注信息随意，只要能看得出是开发者或者学习者即可。

《数据结构精讲与习题详解》C++实现版。章节——1.1
MarkBeY博客
11-05 247
《数据结构精讲与习题详解》C++实现版。章节——1.1
素数幻方
qq_42200028的博客
07-30 1476
求一个四阶的素数幻方，即在一个4*4的矩阵中，每一个数字，使每一行、每一列、和两个对角线上的4个数字所组成的四位数，均为可逆素数。 #include&lt;stdio.h&gt; #include&lt;math.h&gt; int number[210][5]; int select[110]; int array[4][5]; int count; int selecount;...
c++语言幻方
12-03
所谓幻方，就是一个n行n列的正方形，共有n2个格子，将1、2、3、……、n2这些数字放到这些格子里，使其每行的和、每列的和及两条对角线的和都是一个相同的数S，S称为幻和。当n为奇数时，称为奇数阶幻方，当n为偶数时，称为偶阶幻方。当n可被4整除时，称方为双偶幻方。当n不可被4整除时，称为单偶幻方。多少年来，许多数学家都在研究这个古老而有趣的问题，试图找出一般的解法，但一般都是针对当n是奇数和n是4的倍数的情况。当n是奇数时的算法：首先，将1放在第一行中间一个格子里。其次，依次将后一个数放到前一个数的右上格，如：将2放到1的右上格。将３放到２的右上格等等。可能出现下面的情况。①若右上格从上面超出，则将后一数放到与右上格同列的最后一行。②若右上格从右面超出，则将后一数放到与右上格同行的最后一列。③若右上格既从右面超出又从上面超出，则将后一数放到一数前一数的下面前一数的下面。④若右上格已被数字填充，则将后一数放到前一数的下面依以上法则，你可以很快的写出奇数阶幻方！当然，这种写法只是其中一个答案，而不是唯一答案。
通常的n阶幻方由1,2,3,4,...,n^2填入构成。而素数幻方全是由素数构成的各行、各列与两对角线之和均相等的方阵。 试使用c语言代码在一般区间[c,d]找出9个素数，构成一个三阶素数幻方，使得...
weixin_42587866的博客
02-18 353
答：可以使用C语言代码在一般区间[c,d]中搜索9个素数，构成一个三阶素数幻方，使得该方阵中3行、3列与两对角线上的3个数之和均相等。然后统计出素数幻方的个数，要注意的是转置后仍是素数幻方，但不再计入总数。 ...
穷举算法—素数幻方
午夜阳光
04-06 4047
素数幻方 1.  案例提出          通常的n阶幻方由1,2,...,n2构成的各行、各列与两对角线之和均相等n行n列方阵。素数幻方全是由素数构成的各行、各列与两对角线之和均相等方阵。 试寻求9个素数，构造一个3阶素数幻方，使得该素数方阵中3行、3列与两对角线上的3个素数之和均等于给定的整数s。 2. 设计要点      （1） 数学建模 设幻方正中间数为n，幻和(即每行,每列
SCAU软件开发基础C++复习
shuaiyoutiao的博客
12-08 674
C++语言程序设计（第4版） 文章目录C++语言程序设计（第4版）第一章 绪论1.1 面向对象的基本概念1.1.1 对象1.1.2 类（数据类型）1.1.3 封装1.1.4 继承1.1.5 多态性第二章 C++简单程序设计2.1 C++的特点2.2 基本数据类型2.3 变量的存储类型2.4 运算方式和表达式2.4.1 算数运算式与算数表达式2.4.2 赋值运算符与赋值表达式2.4.3 逗号运算和逗号表达式2.4.4 逻辑运算与逻辑表达式2.4.5 条件运算符与条件表达式2.4.6 sizeof运算符2.4.
复习笔记（九）——C++中的容器（STL容器）
别呀的博客
05-03 672
定义 容器：能够操作多种数据结构和算法的模板类和函数库。能够像现实中的容器一样存放各种类型的对象，一个容器中的所有对象必须是同一个类型的。 容器分为两大类：顺序容器、关联容器。 顺序容器 顺序容器有以下三种：vector（向量）、 deque（双队列）、list（列表）。 它们之所以被称为顺序容器，是因为元素在容器中的位置同元素的值无关，即容器不是排序的。将元素插入容器时，指定在什么位置（尾部、头部或中间某处）插入，元素就会位于什么位置。 1、vector 连续存储结构，每个元素在内存上是连续的； 支持高效
第十五届蓝桥杯软件赛C/C++大学A组个人省赛题解
最新发布
2301_78234743的博客
04-30 1762
** * 代码中的类名、方法名、参数名已经指定，请勿修改，直接返回方法规定的值即可 * * * @从三月初开始零碎的投递，到三月底开始海投，三月的时候面试的机会其实很少，到4月初面试机会才逐渐多了起。如题，具体地址是渝北区重庆涉外商务区B4幢附近，这附近有什么好的房源吗，以及有合租的小伙伴吗，预计五。#字节跳动内推#我的内推码：JNG3GB9我的内推码：JNG3GB9我的内推码：JNG3GB9我的内。
幻方问题C语言实现：设n为奇数，试在n n的表格内不重复地填上1, 2, …, n n数字，使得横向、纵向、两对角线上数之和都相等
06-14
幻方问题C语言实现：设n为奇数，试在n n的表格内不重复地填上1, 2, …, n n数字，使得横向、纵向、两对角线上数之和都相等
C语言幻方矩阵问题源代码
05-07
C语言幻方矩阵问题源代码，幻方矩阵是C语言的一个经典问题。
c三阶素数魔方阵
08-29
c语言编写的魔方阵，可以实现三阶素数魔方阵以及任意阶魔方阵，大一时候的c语言课程设计，当时编了好久才想出来。希望对大家有用。
HBUOJ--素数幻方
绀香零八的博客
11-30 1101
HBUOJ–素数幻方 求四阶的素数幻方。即在一个4*4的矩阵中，每一个格填入一个数字，使每一行、每一列和两条对角线上的四个数字所组成的四位数， 均为可逆素数。 【可逆素数：将某一素数的各位数字的顺序颠倒后得到的数仍是素数】 **输入格式要求："%d" 提示信息：“There are magic aquares with invertable primes as follow:\n” **输出格式要...
素数幻方 求四阶的素数幻方。即在一个4*4的矩阵中，每一个格填入一个数字，使每一行、每一列和两条对角线上的四个数字所组成的四位数，均为可逆素数（改数逆序后仍是素数）。 提示：这样的矩阵有很多。
listboy的博客
02-05 6074
素数幻方 求四阶的素数幻方。即在一个4*4的矩阵中，每一个格填入一个数字，使每一行、每一列和两条对角线上的四个数字所组成的四位数，均为可逆素数（改数逆序后仍是素数）。 提示：这样的矩阵有很多。 **输出提示信息：“There are magic aquares with invertable primes as follow:\n” **输出格式要求：“No.%d\n” “%6d” 程序运行示例如下： There are magic aquares with invertable primes as fol
c语言编程素数幻方,求素数幻方，这道题怎么做啊？(c语言)
weixin_28707365的博客
05-18 211
4阶的话，可以考虑直接递归搜索。我试了试有2992组解。1 2 11 124 9 8 57 10 3 146 13 16 151 2 11 124 9 8 513 10 3 146 7 16 151 2 11 124 15 8 57 16 3 146 13 10 91 2 ...
C++题目：验证幻方
SunnyLi1106的博客
01-03 569
幻方是一类有趣的数阵，它是由若干个整数组成的正方形。对于一个幻方，它每一行、每一列和两条对角线上的数相加之和均为同一个值。给定一个正方形的数阵，请编程判断它是否是一个幻方。
ZJGSU 2019 三阶素数幻方
2301_79821571的博客
12-11 468
试在一般区间[c,d]找出9个素数，构成一个三阶素数幻方，使得该方阵中3行、3列与两对角线上的3个数之和均相等。最后统计出素数幻方的个数。通常的n阶幻方由1,2,3,4,...,n^2填入构成。而素数幻方全是由素数构成的各行、各列与两对角线之和均相等的方阵。输入正整数c和d(c
C语言实现素数幻方
"C语言实现素数幻方的程序，通过找出1001到10000之间的素数，构建一个4x4的幻方，使得每行、每列、主对角线和副对角线上的数字之和都是素数。" 在编程领域，素数幻方是一种特殊的矩阵，其中每个单元格包含一个素数...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值