部分高精度算法c语言

1、部分高精度算法 高精度除以低精度； 算法：按照从高位到低位的顺序，逐位相除。在除到第j位时，该位在接受了来自第j+1位的余数后与除数相除，如果最高位为零，则商的长度减一。源程序如下： #include #define N 500 main() int aN = 0, cN = 0; int i, k, d, b; char a1N; printf(Input 除数:); scanf(%d, &b); printf(Input 被除数:); scanf(%s, a1); k = strlen(a1); for(i = 0; i = 0 ; i-) d = d * 10 + ai; ci = d

2、/ b; d = d % b; while(ck - 1 = 0 & k 1) k-; printf(商=); for(i = k - 1; i = 0; i-) printf(%d, ci); printf(n余数=%d, d); 高精度乘以高精度（要求用尽可能少的存储单元）； 算法：用数组保存两个高精度数，然后逐位相乘，注意考虑进位和总位数。源程序如下： #include main() int a240 = 0, b240 = 0, c480 = 0; int i, j, ka, kb, k; char a1240, b1240; gets(a1); ka = strlen(a1); ge

3、ts(b1); kb = strlen(b1); k = ka + kb; for(i = 0; i ka; i+) ai = a1ka-i-1 - 0; for(i = 0; i kb; i+) bi = b1kb-i-1 - 0; for(i = 0; i ka; i+) for(j = 0; j = 0; i-) printf(%d, ci); 高精度除以高精度（要求用尽可能少的存储单元）； 算法：用计算机模拟手算除法，把除法试商转化为连减。 #include #define N 500 int bj(int a, int b, int k1, int k2) /*比较大小函数*/ in

4、t i, t, flag; /*flag作标志位*/ if(k1 k2) flag = 1; /*被除数大于除数返回1*/ else /*被除数和除数位数相等则逐位进行比较*/ i = k1; t = 0; while(t = 0 & i 0) if(ai bi) t = 1; flag = 1; else if(ai = bi) i-; else t = 1; flag = 0; if(i = 0 & t = 0) flag = 2; /*被除数等于除数返回2*/ return flag; int jf(int a, int b, int k1, int k2) /*减法运算*/ int i

5、, k, dN; for(i = 0; i k2; i+) di = bi; /*把除数赋给数组d*/ for(i = k2; i N; i+) di = 0; /*d数组无数据的高位置0*/ k = k1 - k2 - 1; /*计算减法起始位置*/ if(k 0) for(i = k2 - 1; i = 0; i-) di + k = di; /*移动减数位数与被减数对齐*/ for(i = 0; i k; i+) di = 0; /*移动后的其余位置0*/ for(i = 0; i = di) ai -= di; else ai + 1 = ai + 1 - 1; ai = 10 + a

6、i - di; return k; main() int aN = 0, bN = 0, cN = 0, dN = 0; int i, ka, kb, m, t, t1, t2, k, x, kd, kk; char a1N, b1N; printf(Input 被除数:); scanf(%s, a1); ka = strlen(a1); for(i = 0; i ka; i+) ai = a1ka - i -1 - 0; printf(Input 除数:); scanf(%s, b1); kb = strlen(b1); for(i = 0; i = 1) k = jf(a, b, ka,

7、kb); ck+; if(k m) m = k; t1 = 0; for(i = k; i 0) m+; cm = t1; while(t = 1); if(t2 = 0) printf(商=0); printf(n余数=); for(i = kd - 1; i = 0; i-) printf(%d, ai); exit(1); if(t2 = 2) printf(商 = 1); printf(n余数 = 0); exit(1); kk = kd; while(!ckd - 1) kd-; printf(商 = ); for(i = kd - 1; i = 0; i-) printf(%d,

8、ci); while(!akk) kk-; printf(n余数 = ); if(kk = 0; i-) printf(%d, ai); ！，要求精确到位（。 算法：结果用数组a保存，开始时a0=1，依次乘以数组中各位，注意进位和数组长度的变化。源程序如下： #include #define M 1000 main() int aM, i, n, j, flag = 1; printf(n=); scanf(%d,&n); printf(n!=); a0 = 1; for(i = 1; i M; i+) ai = 0; for(j = 2; j = n; j+) for(i = 0; i fl

9、ag; i+) ai *= j; for(i = 0; i = 10) ai+1 += ai/10; ai = ai % 10; if(i = flag-1) flag+; for(j = flag - 1; j = 0; j-) printf(%d, aj); 问题5. 以字符串方式由键盘输入两个高精度的进制正整数，串长小于，以第一个数为被除数，第二个数为除数，进行高精度除法运算，并显示商和余数。 算法： 问题6. 麦森数 【问题描述】形如2P-1的素数称为麦森数，这时P一定也是个素数。但反过来不一定，即如果P是个素数，2P-1不一定也是素数。到1998年底，人们已找到了37个麦森数。最大的

10、一个是P=3021377，它有909526位。麦森数有许多重要应用，它与完全数密切相关。 任务：从文件中输入P（1000P3100000），计算2P-1的位数和最后500位数字（用十进制高精度数表示） 【输入格式】 文件中只包含一个整数P（1000P3100000） 【输出格式】 第一行：十进制高精度数2P-1的位数。 第2-11行：十进制高精度数2P-1的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0） 不必验证2P-1与P是否为素数。 【输入样例】 1279 【输出样例】 386 000000000000000000000000000000000000000000

11、00000000 00000000000000000000000000000000000000000000000000 00000000000000104079321946643990819252403273640855 38615262247266704805319112350403608059673360298012 23944173232418484242161395428100779138356624832346 49081399066056773207629241295093892203457731833496 615835504729594205476898112116936771

12、47548478866962 50138443826029173234888531116082853841658502825560 46662248318909188018470682222031405210266984354887 32958028878050869736186900714720710555703168729087 算法：2的幂可以转化成左移运算，为了提高运算速度，可每次左移10位，即每次乘210。对于个位单独考虑，每次左移一位。源程序如下： #include #include #define MAX 100000 main() int p; int i, j; scanf(

13、%d, &p); printf(%dn, (int)(p * log10(2.0) + 1); long store110 = 0; store0 = 1; int left = p % 10; p /= 10; for(i = 1; i = p; i+) for(j = 0; j = 100; j+) storej = 10; for(j = 0; j = MAX) storej + 1 += storej / MAX; storej %= MAX; for(i = 1; i = left; i+) for(j = 0; j = 100; j+) storej = 1; for(j = 0;

14、 j = MAX) storej + 1 += storej / MAX; storej %= MAX; store0 -= 1; for(i = 1; i 100; i+) if(storei - 1 = 0; i-) printf(%05d, storei); if(100 - i) % 10 = 0) printf(n); 问题7. 精确计算：1n+2n+3n+mN的后K位数，其中1=n=9999，1=m=9999，1=K=99。N，m，K从键盘输入。 例：输入：n=3 m=3 K=3 输出：36 问题8. 有一个正整数N（N可能达到120位），它是由若干个不大于65535的正整数相乘而

15、得到的。请把这个数分解成素数因子（质因子）的乘积。 输入：输入文件只有一行为N的值。 输出：（1）素数因子由小到大分行输出； （2）每一行输出一个素数因子和该素数因子的个数，用一个空格分开； （3）如果正整数N的分解中有一个以上的大于65535的素数，请按照（1）、（2）的要求输出分解中的小于65535的素数后，在下一行输出“DATA ERROR！”。 算法：先将2到65535之间的所有素数保存在数组中，用这个数去除数组中的每一个数，得到一个质因数就打印出来。源程序如下： #include #include int length, temp120; int sushu(int a) int i

16、, j, k = 0, m; for(i = 2; i = 65537; i+) m = sqrt(i); for(j = 2; j m) ak = i; k+; return k; int divide(int a, int k) int i, d = 0; for(i = length - 1; i = 0; i-) d = d * 10 + ai; tempi = d / k; d = d % k; if(!d) while(templength - 1 = 0 & length 1) length-; for(i = 0; i length; i+) ai = tempi; tempi

17、 = 0; for(i = length; i 120; i+) ai = 0; else for(i = 0; i length; i+) tempi = 0; return d; main() int i, k, s, d; /*s计数器; d余数*/ int a6600, b120 = 0, c120 = 0; char b1120; gets(b1); length = strlen(b1); for(i = 0; i length; i+) bi = b1length - i - 1 - 0; k = sushu(a); for(i = 0; i k; i+) s = 0; d =

18、divide(b, ai); while(!d) s+; d = divide(b, ai); if(i = k - 1) printf(Data Error!); break; 附录资料:不需要的可以自行删除Abstract: Based on the comprehensive analysis on the plastic parts structure service requirement, mounding quality and mould menu factoring cost. A corresponding injection mould of internal side

19、core pulling was designed. By adopting the multi-direction and multi-combination core-pulling. A corresponding injection mould of internal side core pulling was designed, the working process of the mould was introducedC语言详解 - 枚举类型注：以下全部代码的执行环境为VC+ 6.0在程序中，可能需要为某些整数定义一个别名，我们可以利用预处理指令#define来完成这项工作，您的

20、代码可能是：#define MON 1#define TUE 2#define WED 3#define THU 4#define FRI 5#define SAT 6#define SUN 7在此，我们定义一种新的数据类型，希望它能完成同样的工作。这种新的数据类型叫枚举型。1. 定义一种新的数据类型 - 枚举型 以下代码定义了这种新的数据类型 - 枚举型enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN;(1) 枚举型是一个集合，集合中的元素(枚举成员)是一些命名的整型常量，元素之间用逗号,隔开。(2) DAY是一个标识符，可以看成这个集合的名字，是一个

21、可选项，即是可有可无的项。(3) 第一个枚举成员的默认值为整型的0，后续枚举成员的值在前一个成员上加1。(4) 可以人为设定枚举成员的值，从而自定义某个范围内的整数。(5) 枚举型是预处理指令#define的替代。(6) 类型定义以分号;结束。2. 使用枚举类型对变量进行声明新的数据类型定义完成后，它就可以使用了。我们已经见过最基本的数据类型，如：整型int, 单精度浮点型float, 双精度浮点型double, 字符型char, 短整型short等等。用这些基本数据类型声明变量通常是这样：char a; /变量a的类型均为字符型charchar letter;int x, y, z; /变量

22、x,y和z的类型均为整型intint number;double m, n;double result; /变量result的类型为双精度浮点型double既然枚举也是一种数据类型，那么它和基本数据类型一样也可以对变量进行声明。方法一：枚举类型的定义和变量的声明分开enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN;enum DAY yesterday;enum DAY today;enum DAY tomorrow; /变量tomorrow的类型为枚举型enum DAYenum DAY good_day, bad_day; /变量good_day和bad

23、_day的类型均为枚举型enum DAY方法二：类型定义与变量声明同时进行：enum /跟第一个定义不同的是，此处的标号DAY省略，这是允许的。 saturday, sunday = 0, monday, tuesday, wednesday, thursday, friday workday; /变量workday的类型为枚举型enum DAYenum week Mon=1, Tue, Wed, Thu, Fri Sat, Sun days; /变量days的类型为枚举型enum weekenum BOOLEAN false, true end_flag, match_flag; /定义枚举

24、类型并声明了两个枚举型变量方法三：用typedef关键字将枚举类型定义成别名，并利用该别名进行变量声明：typedef enum workday saturday, sunday = 0, monday, tuesday, wednesday, thursday, friday workday; /此处的workday为枚举型enum workday的别名workday today, tomorrow; /变量today和tomorrow的类型为枚举型workday，也即enum workdayenum workday中的workday可以省略：typedef enum saturday, s

25、unday = 0, monday, tuesday, wednesday, thursday, friday workday; /此处的workday为枚举型enum workday的别名workday today, tomorrow; /变量today和tomorrow的类型为枚举型workday，也即enum workday也可以用这种方式：typedef enum workday saturday, sunday = 0, monday, tuesday, wednesday, thursday, friday;workday today, tomorrow; /变量today和tom

26、orrow的类型为枚举型workday，也即enum workday注意：同一个程序中不能定义同名的枚举类型，不同的枚举类型中也不能存在同名的命名常量。错误示例如下所示：错误声明一：存在同名的枚举类型typedef enum wednesday, thursday, friday workday;typedef enum WEEK saturday, sunday = 0, monday, workday; 错误声明二：存在同名的枚举成员typedef enum wednesday, thursday, friday workday_1;typedef enum WEEK wednesday,

27、 sunday = 0, monday, workday_2;3. 使用枚举类型的变量3.1 对枚举型的变量赋值。实例将枚举类型的赋值与基本数据类型的赋值进行了对比：方法一：先声明变量，再对变量赋值#include/* 定义枚举类型 */enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN ;void main() /* 使用基本数据类型声明变量，然后对变量赋值 */ int x, y, z; x = 10; y = 20; z = 30; /* 使用枚举类型声明变量，再对枚举型变量赋值 */ enum DAY yesterday, today, tomor

28、row; yesterday = MON; today = TUE; tomorrow = WED; printf(%d %d %d n, yesterday, today, tomorrow);方法二：声明变量的同时赋初值#include /* 定义枚举类型 */enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN ;void main() /* 使用基本数据类型声明变量同时对变量赋初值 */ int x=10, y=20, z=30; /* 使用枚举类型声明变量同时对枚举型变量赋初值 */ enum DAY yesterday = MON, today

29、= TUE, tomorrow = WED; printf(%d %d %d n, yesterday, today, tomorrow);方法三：定义类型的同时声明变量，然后对变量赋值。#include /* 定义枚举类型，同时声明该类型的三个变量，它们都为全局变量 */enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN yesterday, today, tomorrow;/* 定义三个具有基本数据类型的变量，它们都为全局变量 */int x, y, z;void main() /* 对基本数据类型的变量赋值 */ x = 10; y = 20; z

30、= 30; /* 对枚举型的变量赋值 */ yesterday = MON; today = TUE; tomorrow = WED; printf(%d %d %d n, x, y, z); /输出：10 20 30 printf(%d %d %d n, yesterday, today, tomorrow); /输出：1 2 3方法四：类型定义，变量声明，赋初值同时进行。#include /* 定义枚举类型，同时声明该类型的三个变量，并赋初值。它们都为全局变量 */enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN yesterday = MON, t

31、oday = TUE, tomorrow = WED;/* 定义三个具有基本数据类型的变量，并赋初值。它们都为全局变量 */int x = 10, y = 20, z = 30;void main() printf(%d %d %d n, x, y, z); /输出：10 20 30 printf(%d %d %d n, yesterday, today, tomorrow); /输出：1 2 33.2 对枚举型的变量赋整数值时，需要进行类型转换。#include enum DAY MON=1, TUE, WED, THU, FRI, SAT, SUN ;void main() enum DA

32、Y yesterday, today, tomorrow; yesterday = TUE; today = (enum DAY) (yesterday + 1); /类型转换 tomorrow = (enum DAY) 30; /类型转换 /tomorrow = 3; /错误 printf(%d %d %d n, yesterday, today, tomorrow); /输出：2 3 303.3 使用枚举型变量#includeenum BELL = a, BACKSPACE = b, HTAB = t, RETURN = r, NEWLINE = n, VTAB = v, SPACE =

33、;enum BOOLEAN FALSE = 0, TRUE match_flag;void main() int index = 0; int count_of_letter = 0; int count_of_space = 0; char str = Im Ely efod; match_flag = FALSE; for(; strindex != 0; index+) if( SPACE != strindex ) count_of_letter+; else match_flag = (enum BOOLEAN) 1; count_of_space+; printf(%s %d ti

34、mes %c, match_flag ? match : not match, count_of_space, NEWLINE); printf(count of letters: %d %c%c, count_of_letter, NEWLINE, RETURN);输出：match 2 timescount of letters: 10Press any key to continue4. 枚举类型与sizeof运算符#include enum escapes BELL = a, BACKSPACE = b, HTAB = t, RETURN = r, NEWLINE = n, VTAB = v, SPACE = ;enum BOOLEAN FALSE = 0, TRUE match_flag;void main() printf(%d bytes n, sizeof(enum escapes); /4 b