标题:求一个十进制300位的质数
取消只看楼主
etsunny
Rank: 1
等 级:新手上路
帖 子:40
专家分:0
注 册:2007-10-11
 问题点数:0 回复次数:3 
求一个十进制300位的质数

我用的是数组,可是做除法的算法没有解决

搜索更多相关主题的帖子: 质数 十进制 
2007-11-08 22:52
etsunny
Rank: 1
等 级:新手上路
帖 子:40
专家分:0
注 册:2007-10-11
得分:0 

这是我写的源代码,再这里发下,太混乱了,不知道怎么说好,可能不是这样算的

#include<stdio.h>
#include<math.h>
#define N 10000
#define M N/2+1

main ()
{
int i=0,j=0,k=0,l=0,m=0,n=0;
int a[N]={0};
int b[M]={0};
int c[N]={0};
int jiayi(int b[M],int i);
int yushu(int a[N],int b[M],int n,int m);

/*输出到文件D:\iloveyou\1中*/
/*
for(l=0;l<N;l++)
{
a[l]=9;
fprintf(D:\iloveyou\1 * fp,"%d",a[N-1-l]);
} */

/*求N位的大质数,用数组表示和运算*/

for(a[N]=2;a[0]<=9;jiayi(a[N],N))
{
for(b[M]=1;b[0]<1;jiayi(b[M],M))
{
if (yushu(a[N],b[M],n,m)==0)
break;
}
if(b[0]>0)
for(l=0;l<N;l++)
{
c[l]=a[l];
}
}
for(i=0;i<N;i++)
{
printf("%d\n",c[i]);
}


}

/* i=m; 设b[0]为最高位,b[M]为个位,自增1;*/
jiayi(int b[M],int i)
{
b[i]++;
if (b[i]>9)
b[i]=0;
i--;
jiayi(b[M],i);
reaturn (b[M],i);
}

/*a[N]除以b[M],按各个位进行计算,即n位整数除以m位整数,取余*/
yushu(int a[N],int b[M],int n,int m)
{
/* int n=0,m=1; */
if (n<m)
{
if (a[n]>b[m])
{
a[n]=a[n]-b[m];
yushu(a[N],b[M],n,m);
}
else if (a[n]=b[m])
{
n++;
m++;
yushu(a[N],b[M],n,m);
}
else if (a[n]<b[m])
{
n++;
yushu(a[N],b[M],n,m);
}
}
else
{
if (a[n]>b[m])
{
a[n]=a[n]-b[m];
yushu(a[N],b[M],n,m);
}
else if (a[n]=b[m])
{
n++;
m++;
if (m<M)
yushu(a[N],b[M],n,m);
reaturn (0);
}
else if (a[n]<b[m])
{
reaturn (1);
}
}

}


圣火从希腊蔓延,带着普罗米修斯的祝愿
2007-11-10 17:36
etsunny
Rank: 1
等 级:新手上路
帖 子:40
专家分:0
注 册:2007-10-11
得分:0 
编译没有通过

圣火从希腊蔓延,带着普罗米修斯的祝愿
2007-11-10 17:36
etsunny
Rank: 1
等 级:新手上路
帖 子:40
专家分:0
注 册:2007-10-11
得分:0 

后来我在考虑32768进制,那样效率会不会高一点?


圣火从希腊蔓延,带着普罗米修斯的祝愿
2007-11-11 11:49



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-184661-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.361317 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved