高手们，帮忙啊，c语言程序设计

Posted 2023-04-20

素数幻方
求四阶的素数幻方。即在一个4X4 的矩阵中，每一个格填 入一个数字，使每一行、每一列和两条对角线上的4 个数字所组成的四位数，均为可逆素数。
要求（1） 请给出问题分析与算法设计，并画出流程图；
（2） 编写程序，并给出说明与注释；
简单点的最好····

*问题分析与算法设计
最简单的算法是：采用穷举法，设定4X4矩阵中每一个元素的值后，判断每一行、每一列和两条对角线上的4个数字组成的四位数是否都是可逆素数，若是则求出了满足题意的一个解。
这种算法在原理是对的，也一定可以求出满足题意的全部解。但是，按照这一思路编出的程序效率很低，在微机上几个小时也不会运行结束。这一算法致命的缺陷是：要穷举和判断的情况过多。
充分利用题目中的“每一个四位数都是可逆素数”这一条件，可以放弃对矩阵中每个元素进行的穷举的算法，先求出全部的四位可逆素数(204个)，以矩阵的行为单位，在四位可逆素数的范围内进行穷举，然后将穷举的四位整数分解为数字后，再进行列和对角线方向的条件判断，改进的算法与最初的算法相比，大大地减少了穷举的次数。
考虑矩阵的第一行和最后一行数字，它们分别是列方向四位数的第一个数字和最后一个数字，由于这些四位数也必须是可逆素数，所以矩阵的每一行和最后一行中的各个数字都不能为偶数或5。这样穷举矩阵的第一行和最后一行时，它们的取值范围是：所有位的数字均不是偶数或5的四位可逆数。由于符合这一条件的四位可逆素数很少，所以这一范围限制又一次减少了穷举的次数。
对算法的进一步研究会发现：当设定了第一和第二行的值后，就已经可以判断出当前的这种组合是否一定是错误的(尚不能肯定该组合一定是正确的)。若按列方向上的四个两位数与四位可逆数的前两位矛盾(不是其中的一种组合)，则第一、二行的取值一定是错误的。同理在设定了前三行数据后，可以立刻判断出当前的这种组合是否一定是错误的，若判断出矛盾情况，则可以立刻设置新的一组数据。这样就可以避免将四个数据全部设定好以后再进行判断所造成的低效。
根据以上分析，可以用伪语言描述以上改进的算法：
开始
找出全部四位的可逆素数；
确定全部出现在第一和最后一行的四位可逆素数；
在指定范围 内穷举第一行
在指定范围内穷举第二行
若第一、第二、三行已出现矛盾，则继续穷举下一个数；
在指定范围内穷举第四行
判断列和对角方向是否符合题意
若符合题意，则输出矩阵；
否则继续穷举下一个数；
结束
在实际编程中，采用了很多程序设计技巧，假如设置若干辅助数组，其目的就是要最大限度的提高程序的执行效率，缩短运行时间。下面的程序运行效率是比较高的。

*程序说明与注释
#include<stdio.h>
#include<math.h>
int number[210][5]; /*存放可逆素数及素数分解后的各位数字*/
int select[110]; /*可以放在矩阵第一行和最后一行的素数的下标*/
int array[4][5]; /*4X4的矩阵，每行0号元素存可逆素数对应的数组下标*/
int count; /*可逆素数的数目*/
int selecount; /*可以放在矩阵第一行和最后一行的可逆素数的数目*/
int larray[2][200]; /*存放素数前二、三位数的临时数组所对应的数量计数器*/
int lcount[2];
int num(int number);
int ok(int number);
void process(int i);
void copy_num(int i);
int comp_num(int n);
int find1(int i);
int find2(void);
int find0(int num);
void p_array(void);

int main()

int i,k,flag,cc=0,i1,i4;
printf("there are magic squares with invertable primes as follw:\n");
for(i=1001;i<9999;i+=2) /*求满足条件的可逆素数*/

k=i/1000;
if(k%2!=0&&k!=5&&num(i)) /*若可逆素数的第一位不是偶数或5*/

number[count][0]=i; /*存入数组*/
process(count++); /*分解素数的各位数字*/
if(number[count-1][2]%2!=0&& /*若可逆素数满足放在矩阵第一行*/
number[count-1][3]%2!=0&& /*和最后一行的条件，记录可逆素数的*/
number[count-1][2]!=5&& /*下标，计数器加1*/
number[count-1][3]!=5)
select[selecount++]=count-1;


larray[0][lcount[0]++]=number[0][0]/100; /*临时数组的第一行存前二位*/
larray[1][lcount[1]++]=number[0][0]/10; /*临时数组的第二行存前三位*/
for(i=1;i<count;i++) /*将素数不重复的前二、三位存入临时数组中*/

if(larray[0][lcount[0]-1]!=number[i][0]/100)
larray[0][lcount[0]++]=number[i][0]/100;
if(larray[1][lcount[1]-1]!=number[i][0]/10)
larray[1][lcount[1]++]=number[i][0]/10;

for(i1=0;i1<selecount;i1++) /*在第一行允许的汇聚围内穷举*/

array[0][0]=select[i1]; /*取对应的素数下标*/
copy_num(0); /*复制分解的素数*/
for(array[1][0]=0;array[1][0]<count;array[1][0]++) /*穷举第二行*/

copy_num(1); /*复制分解的数字*/
if(!comp_num(2))
continue; /*若每列的前两位的组成与素数相矛盾，则试探下一个数*/
for(array[2][0]=0;array[2][0]<count;array[2][0]++) /*穷举第三行*/

copy_num(2); /*复制分解的数字*/
if(!comp_num(3))
continue; /*若每列的前三位的组成与素数相矛盾，则试探下一个数*/
for(i4=0;i4<selecount;i4++) /*在最后一行允许的范围内穷举*/

array[3][0]=select[i4];
copy_num(3); /*复制分解的数字*/
for(flag=1,i=1;flag&&i<=4;i++) /*判断每列是否可逆素数*/
if(!find1(i))flag=0;
if(flag&&find2()) /*判断对角线是否为可逆素数*/
printf("No.%d\n",++cc); p_array(); /*输出幻方矩阵*/






int num(int number) /*判断是否可逆素数*/

int j;
if(!ok(number)) return 0;
for(j=0;number>0;number/=10) /*将素数变为反序数*/
j=j*10+number%10;
if(!ok(j)) return 0; /*判断反序数是否为素数*/
return 1;


int ok(int number) /*判断是否为素数*/

int i,j;
if(number%2==0) return 0;
j=sqrt((double)number)+1;
for(i=3;i<=j;i+=2)
if(number%i==0) return 0;
return 1;


void process(int i) /*将第i个整数分解为数字并存入数组*/

int j,num;
num=number[i][0];
for(j=4;j>=1;j--,num/=10)
number[i][j]=num%10;


void copy_num(int i) /*将array[i][0]指向的素数的各位数字复制到array[i]中*/

int j;
for(j=1;j<=4;j++)
array[i][j]=number[array[i][0]][j];


int comp_num(int n) /*判断array中每列的前n位是否与可逆素数允许的前n位矛盾*/

static int ii; /*用内部静态变量保存前一次查找到的元素下标*/
static int jj; /*ii:前一次查找前二位的下标，jj:前一次查找前三位的下标*/
int i,num,k,*p; /*p:指向对应的要使用的前一次下标ii或jj*/
int *pcount; /*pcount:指向要使用的临时数组数量的计数器*/
switch(n) /*根据n的值选择对应的一组控制变量*/
case 2:pcount=&lcount[0];p=ⅈbreak;
case 3:pcount=&lcount[1];p=&jj;break;
default:return 0;

for(i=1;i<=4;i++) /*对四列分别进行处理*/

for(num=0,k=0;k<n;k++) /*计算前n位数字代表的数值*/
num=num*10+array[k][i];
if(num<=larray[n-2][*p]) /*与前一次最后查找到的元素进行比较*/
for(;*p>=0&&num<larray[n-2][*p];(*p)--);/*若前次查找到的元素大，则向前找*/
else
for(;p<pcount&&num>larray[n-2][*p];(*p)++); /*否则向后找*/
if(*p<0||*p>=*pcount)

*p=0; return 0;

if(num!=larray[n-2][*p])
return 0; /*前n位不是可逆素数允许的值则返回0*/

return 1;


int find1(int i) /*判断列方向是否是可逆素数*/

int num,j;
for(num=0,j=0;j<4;j++)
num=num*10+array[j][i];
return find0(num);


int find2(void) /*判断对角线方向是否是可逆素数*/

int num1,num2,i,j;
for(num1=0,j=0;j<4;j++)
num1=num1*10+array[j][j+1];
for(num2=0,j=0,i=4;j<4;j++,i--)
num2=num2*10+array[j][i];
if(find0(num1)) return(find0(num2));
else return 0;


int find0(int num) /*查找是否为满足要求的可逆素数*/

static int j;
if(num<=number[j][0])for(;j>=0&&num<number[j][0];j--);
else for(;j<count&&num>number[j][0];j++);
if(j<0||j>=count) j=0;return 0;
if(num==number[j][0]) return 1;
else return 0;


void p_array(void) /*输出矩阵*/

int i,j;
for(i=0;i<4;i++)

for(j=1;j<=4;j++) printf("%d ",array[i][j]);
printf("\n");



*问题的进一步讨论
程序中大量技巧是用于尽早发现矛盾，减少循环次数，缩短运行时间。从实际效果看是相当不错的。但目前的程序仍然可以进一步优化。
当第四行设定了前三行后，尚未设定的行就没必要再使用穷举的方法，因为列方向设定好的三位数字已经限制了最后一个数字可能的取值，在可逆数中找出前三位数字与设定好的三位数字相同的素数。这些素数就是在这一列前面已设定好的三位数字的限制条件下可能的取值。此时每一列上只有不超过四个可能的取值。找出全部各列可能的取值(可能的四位可逆素数)，求出它们的交集。若交集为空，即没有共同的可能取值，则列间数据相互矛盾否满足则将交集中的数据填 入矩阵中就是题目的一个解。
算法可再进一步优化。先穷举一、二和四列的数据，然后用上面的算法来确定第三行的值，这样可进一步缩小穷举的范围，提高运行效率。
分析输出的结果。可以看出本题的基本解只有17种，每个解可通过旋转与反射获得同构的其它7个解，可以进一步改进程序，只输出17个基本解。 参考技术A 瞧好啦！别不懂！

#include<stdio.h>
#include<math.h>
int number[210][5]; /*存放可逆素数及素数分解后的各位数字*/
int select[110]; /*可以放在矩阵第一行和最后一行的素数的下标*/
int array[4][5]; /*4X4的矩阵，每行0号元素存可逆素数对应的数组下标*/
int count; /*可逆素数的数目*/
int selecount; /*可以放在矩阵第一行和最后一行的可逆素数的数目*/
int larray[2][200]; /*存放素数前二、三位数的临时数组所对应的数量计数器*/
int lcount[2];
int num(int number);
int ok(int number);
void process(int i);
void copy_num(int i);
int comp_num(int n);
int find1(int i);
int find2(void);
int find0(int num);
void p_array(void);

int main()

int i,k,flag,cc=0,i1,i4;
printf("there are magic squares with invertable primes as follw:\n");
for(i=1001;i<9999;i+=2) /*求满足条件的可逆素数*/

k=i/1000;
if(k%2!=0&&k!=5&&num(i)) /*若可逆素数的第一位不是偶数或5*/

number[count][0]=i; /*存入数组*/
process(count++); /*分解素数的各位数字*/
if(number[count-1][2]%2!=0&& /*若可逆素数满足放在矩阵第一行*/
number[count-1][3]%2!=0&& /*和最后一行的条件，记录可逆素数的*/
number[count-1][2]!=5&& /*下标，计数器加1*/
number[count-1][3]!=5)
select[selecount++]=count-1;


larray[0][lcount[0]++]=number[0][0]/100; /*临时数组的第一行存前二位*/
larray[1][lcount[1]++]=number[0][0]/10; /*临时数组的第二行存前三位*/
for(i=1;i<count;i++) /*将素数不重复的前二、三位存入临时数组中*/

if(larray[0][lcount[0]-1]!=number[i][0]/100)
larray[0][lcount[0]++]=number[i][0]/100;
if(larray[1][lcount[1]-1]!=number[i][0]/10)
larray[1][lcount[1]++]=number[i][0]/10;

for(i1=0;i1<selecount;i1++) /*在第一行允许的汇聚围内穷举*/

array[0][0]=select[i1]; /*取对应的素数下标*/
copy_num(0); /*复制分解的素数*/
for(array[1][0]=0;array[1][0]<count;array[1][0]++) /*穷举第二行*/

copy_num(1); /*复制分解的数字*/
if(!comp_num(2))
continue; /*若每列的前两位的组成与素数相矛盾，则试探下一个数*/
for(array[2][0]=0;array[2][0]<count;array[2][0]++) /*穷举第三行*/

copy_num(2); /*复制分解的数字*/
if(!comp_num(3))
continue; /*若每列的前三位的组成与素数相矛盾，则试探下一个数*/
for(i4=0;i4<selecount;i4++) /*在最后一行允许的范围内穷举*/

array[3][0]=select[i4];
copy_num(3); /*复制分解的数字*/
for(flag=1,i=1;flag&&i<=4;i++) /*判断每列是否可逆素数*/
if(!find1(i))flag=0;
if(flag&&find2()) /*判断对角线是否为可逆素数*/
printf("No.%d\n",++cc); p_array(); /*输出幻方矩阵*/






int num(int number) /*判断是否可逆素数*/

int j;
if(!ok(number)) return 0;
for(j=0;number>0;number/=10) /*将素数变为反序数*/
j=j*10+number%10;
if(!ok(j)) return 0; /*判断反序数是否为素数*/
return 1;


int ok(int number) /*判断是否为素数*/

int i,j;
if(number%2==0) return 0;
j=sqrt((double)number)+1;
for(i=3;i<=j;i+=2)
if(number%i==0) return 0;
return 1;


void process(int i) /*将第i个整数分解为数字并存入数组*/

int j,num;
num=number[i][0];
for(j=4;j>=1;j--,num/=10)
number[i][j]=num%10;


void copy_num(int i) /*将array[i][0]指向的素数的各位数字复制到array[i]中*/

int j;
for(j=1;j<=4;j++)
array[i][j]=number[array[i][0]][j];


int comp_num(int n) /*判断array中每列的前n位是否与可逆素数允许的前n位矛盾*/

static int ii; /*用内部静态变量保存前一次查找到的元素下标*/
static int jj; /*ii:前一次查找前二位的下标，jj:前一次查找前三位的下标*/
int i,num,k,*p; /*p:指向对应的要使用的前一次下标ii或jj*/
int *pcount; /*pcount:指向要使用的临时数组数量的计数器*/
switch(n) /*根据n的值选择对应的一组控制变量*/
case 2:pcount=&lcount[0];p=ⅈbreak;
case 3:pcount=&lcount[1];p=&jj;break;
default:return 0;

for(i=1;i<=4;i++) /*对四列分别进行处理*/

for(num=0,k=0;k<n;k++) /*计算前n位数字代表的数值*/
num=num*10+array[k][i];
if(num<=larray[n-2][*p]) /*与前一次最后查找到的元素进行比较*/
for(;*p>=0&&num<larray[n-2][*p];(*p)--);/*若前次查找到的元素大，则向前找*/
else
for(;p<pcount&&num>larray[n-2][*p];(*p)++); /*否则向后找*/
if(*p<0||*p>=*pcount)

*p=0; return 0;

if(num!=larray[n-2][*p])
return 0; /*前n位不是可逆素数允许的值则返回0*/

return 1;


int find1(int i) /*判断列方向是否是可逆素数*/

int num,j;
for(num=0,j=0;j<4;j++)
num=num*10+array[j][i];
return find0(num);


int find2(void) /*判断对角线方向是否是可逆素数*/

int num1,num2,i,j;
for(num1=0,j=0;j<4;j++)
num1=num1*10+array[j][j+1];
for(num2=0,j=0,i=4;j<4;j++,i--)
num2=num2*10+array[j][i];
if(find0(num1)) return(find0(num2));
else return 0;


int find0(int num) /*查找是否为满足要求的可逆素数*/

static int j;
if(num<=number[j][0])for(;j>=0&&num<number[j][0];j--);
else for(;j<count&&num>number[j][0];j++);
if(j<0||j>=count) j=0;return 0;
if(num==number[j][0]) return 1;
else return 0;


void p_array(void) /*输出矩阵*/

int i,j;
for(i=0;i<4;i++)

for(j=1;j<=4;j++) printf("%d ",array[i][j]);
printf("\n");

参考技术B 这种题太简单了，还在这问。 参考技术C 自己弄

C语言 四则运算 各位C语言高手帮忙做下这题！！谢谢谢谢谢谢

一、题目：四则运算

二、目的与要求

1． 目的：

通过编写四则运算程序，培养学生综合利用C语言进行程序设计的能力，熟悉字符串的操作方法，加强函数的运用，提高软件系统分析能力和程序文档建立、归纳总结的能力。

2． 基本要求：

1）要求用C语言编程，在Visual C++环境下调试完成；

2）要求按照程序功能分成几个功能模块来实现，各个功能模块分别使用函数来完成；

3）要求应用本课所讲授的程序设计语言知识来解决问题.

三、设计方法和基本原理

1． 课题功能描述

本程序的功能，就是实现数的加减乘除四则运算，如自动计算3+5*8的结果。

2． 问题详细描述

程序运行时，首先提示用户输入四则运算表达式(表达式中最多有3个运算符，数据范围为1-100的整数)；

输入结束后，程序自动进行计算并给出结果。

如：输入13-10+5/8时输出3.625。

3． 问题的解决方案

注意：问题的解决方案有很多，下面给出的仅供同学们参考。

该问题主要注意四则运算的优先级问题：乘除运算比加减运算优先级高，同级运算按从左到右的顺序运算。

本程序应采用模块化设计方法,设计几个功能模块。例如（仅供参考）：

字符串解析函数（将输入字符串分解成数和运算符）

将数据（字符数组中）转换成十进制数（long）
判读是否存在高优先级运算符，若存在首先计算其运算结果并保存。
同级运算按先后顺序进行。

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define maxsize 50
void trans(char str[],char exp[])/*将算术表达式str转换成后缀表达式exp*/

struct
char data[maxsize]; /*存放运算符*/
int top; /*栈指针*/
opr; /*定义运算符栈*/
char ch;
int i=0,t=0; /*t作为exp的下标，i作为str的下标*/
opr.top=-1; /*初始化设定top的值为负一*/
ch=str[i];i++; /*逐个读取字符串中的字符*/
while (ch!='\0') /*str表达式未扫描完时循环*/
switch(ch) /*判定*/

case '(':
opr.top++;opr.data[opr.top]=ch; /*判定为'('号，则将其入栈opr*/
break;
case ')':
while (opr.data[opr.top]!='(') /*判定为')'号*/
exp[t]=opr.data[opr.top]; /*将栈opr中'('以后的字符依次删除并存入数组exp中*/
opr.top--;
t++;

opr.top--; /*将左括号删除*/
break;
case '+': /*判定为加号或减号*/
case '-':
while (opr.top!=-1 &&opr.data[opr.top]!='(')
exp[t]=opr.data[opr.top]; /*将当前栈opr中(以前的所有字符依次删除并存入数组exp中*/
opr.top--;
t++;

opr.top++;opr.data[opr.top]=ch; /*将ch存入栈opr中*/
break;
case '*':
case '/':
while (opr.data[opr.top]=='*'||opr.data[opr.top]=='/'||opr.data[opr.top]=='^')
exp[t]=opr.data[opr.top]; /*将当前栈opr中连续的'*'或'/'或'^'依次删除并存入数组exp中*/
opr.top--;
t++;

opr.top++;opr.data[opr.top]=ch; /*将ch存入栈opr中*/
break;
case '^': /*判定为乘方号*/
while (opr.data[opr.top]=='^')
exp[t]=opr.data[opr.top]; /*将当前栈opr中连续的'^'依次删除并存入数组exp中*/
opr.top--;
t++;

opr.top++;opr.data[opr.top]=ch; /*将ch存入栈opr中*/
break;
case ' ': break; /*过滤掉空格*/
default:
while(ch>='0'&& ch<='9'||ch=='.') /*判定为数字*/
exp[t]=ch;t++; /*将后续数字依次存入数组中*/
ch=str[i];i++;

i--;
exp[t]='#';t++; /*用#标示一个数值串结束*/

ch=str[i];i++;

while (opr.top!=-1) /*此时str扫描完毕，栈不空时循环*/
exp[t]=opr.data[opr.top];
t++;opr.top--;

exp[t]='\0'; /*给exp表达式添加结束标示*/

float compvalue(char exp[]) /*计算后缀表达式的值*/

struct
float data[maxsize]; /*存放数值*/
int top; /*栈指针*/
st; /*定义数值栈*/
float d,d2;double po;
char ch;
int t=0,flag=1,i,count; /*t作为exp的下标*/
st.top=-1;
ch=exp[t];t++;
while (ch!='\0') /*exp字符串为扫描完时循环*/
switch(ch)

case '+':st.data[st.top-1]=st.data[st.top-1]+st.data[st.top]; /*执行两次退栈，并将计算结果入栈*/
st.top--;break;
case '-':st.data[st.top-1]=st.data[st.top-1]-st.data[st.top];
st.top--;break;
case '*':st.data[st.top-1]=st.data[st.top-1]*st.data[st.top];
st.top--;break;
case '/':
if(st.data[st.top]!=0)
st.data[st.top-1]=st.data[st.top-1]/st.data[st.top];
else
printf("\n除零错误！\n");
exit(0); /*除数为零，异常退出*/

st.top--;break;
case '^':
po=pow(st.data[st.top-1],st.data[st.top]); st.data[st.top-1]=(float)po;/*调用pow子函数进行乘方运算*/
st.top--;break;
default:
d=0; flag=1; d2=0; /*将数字字符转换成对应的数值存放到d中*/
while(ch>='0'&&ch<='9'&&flag) /*判定为数字字符*/
d=10*d+ch-'0';
ch=exp[t];t++;
if(ch=='.')
flag=0;


if(flag==0)
ch=exp[t];t++;count=0;
while(ch>='0'&&ch<='9') /*判定为数字字符*/
d2=10*d2+ch-'0';
ch=exp[t];t++;count++;

for(i=1;i<=count;i++)
d2=0.1*d2;

d+=d2;
st.top++;
st.data[st.top]=d;

ch=exp[t];t++;

return st.data[st.top];

int main()

char str[maxsize],exp[maxsize]; /*str存储原算术表达式，exp存储对应的后缀表达式*/
printf("the arithmetic expression is:\n");
gets(str);
trans(str,exp);
printf("the postfix expression is:%s\n",exp);
printf("the result is %g\n",compvalue(exp));
参考技术A 败笔技术论坛，中国最大的编程论坛，建议你上那去学习吧

参考资料：败笔网络安全小组

参考技术B 在网吧 1361019294 我的Q 加老回去给你说。机子上有。 参考技术C 我就喜欢害人：）
#include
#include
#include
#define
err
-1
#define
max
100
/*定义堆栈的大小*/
int
stack[max];
/*用一维数组定义堆栈*/
int
top=0;
/*定义堆栈指示*/
int
push(int
i)
/*存储运算数,入栈操作*/

if(top
评论
0
0
加载更多