posts - 99,  comments - 8,  trackbacks - 0

//思路:1.调用函数change() 将小数 s 转化成去掉小数点的整数 s1 并且返回得到小数点的位置。
//      2.调用函数mu1()用大数乘法模拟整数运算,并将相乘的结果放在s2中,返回s2和是进行下面的 n -1次相乘
//      3.在主函数里面通过已知的小数点的位置,利用数值关系输出
//      4.注意小数点前没有前导0,小数点后面没有尾0。

//难点:for ( int i = 1; i < n; i ++)     mu1 (s1,s2);

 1             
 2#include <stdio.h>
 3#include <stdlib.h>
 4#include <string.h>
 5#define LENGTH 6  //小数的位数(含小数点) 
 6
 7//将字符转化为数字 
 8unsigned int change (char s[LENGTH], unsigned int s1[LENGTH - 1])
 9{
10     int ss[LENGTH - 1];  //ss 放未逆置的整数 
11     memset (ss, 0sizeof(ss));
12     
13     int k = 0 ;
14     for (int i = 0; i < LENGTH && s[i]; ++i)  //考虑特殊数据如:0.0001 
15     {
16         if (s[i] != '.')
17            ss[k++= s[i] - '0';
18     }
 
19     for (int j = 0;j < LENGTH - 1; j++)
20     {
21         s1[j] = ss[LENGTH - 2 - j];  
22     }

23     
24     int m = 0;
25     while ( (s[m] != '.'&& s[m] )
26           ++m; 
27           return LENGTH - 1 - m;   //小数点位数 
28}
 
29
30//大数乘法运算 
31//函数返回 s2 
32void mu1 (unsigned int s1[LENGTH - 1],unsigned int s2[130])
33{
34     int ss[130];
35     memset ( ss, 0, sizeof(ss) );
36     
37     for ( int i = 0; i < LENGTH - 1; i++)
38         for (int j = 0;j < 130; j++)  //难点:因为返回新的s2之后位数会增加 最多时 5* 25 = 125 
39         ss [i + j] += s1[i] * s2[j];
40     
41     //将 两个大数相乘得的积ss中进行进位处理后放到s2 中  
42     int c = 0;
43     for (int i = 0;i < 130;i++)
44     {
45         s2[i] = (c + ss[i]) % 10;
46         c = (c + ss[i]) / 10;
47     } 
48}
49
50int main()
51{
52    int n;
53    char s[LENGTH];  //要处理的幂 R 
54    unsigned int s1[LENGTH - 1];  //将 R 转化成数字 
55    unsigned int s2[130];
56     
57    while(scanf ("%s%d", s, &n) != EOF)
58    {
59        memset (s1, 0, sizeof (s1));
60        memset (s2, 0, sizeof (s2)); 
61        
62        int j = change (s, s1);      //得到小数点所在位置 
63        change (s,s2);              //得到s2 和 s1 进行幂运算 
64        for ( int i = 1; i < n; i ++)
65            mu1 (s1,s2); 
66        
67        //在s2中前面的代表小数位,后面的代表整数位,
68        //所以关键是通过数值关系找到小数点的位置
69      
70      
71         //例:0.1010  * 0.1010 = 0.01020100  
72        int m = 129;//去掉前导0 
73        while ( (!s2[m]) && m)
74        --m;
75        
76        int k = 0; //去掉尾0                       
77        while ( ( !s2[k] ) && (k < 130))                                                           
78        ++k;
79        
80        //输出整数位 
81        for (int i = m; i >= n * j; i--)
82            printf ("%d",s2[i]);
83            
84        //输出小数点
85        if ( j && n * j >= k + 1) 
86        printf (".");
87        
88        for (int i = n*j -1; i >= k; --i)
89        printf ("%d", s2[i]);
90        printf ("\n");
91    }
92    
93    return 0;
94   // system ("pause");   
95}

96
posted on 2010-08-09 13:21 雪黛依梦 阅读(600) 评论(0)  编辑 收藏 引用 所属分类: 大数

只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   知识库   博问   管理


<2011年3月>
272812345
6789101112
13141516171819
20212223242526
272829303112
3456789

常用链接

留言簿(4)

随笔分类

随笔档案

文章档案

搜索

  •  

最新评论

阅读排行榜

评论排行榜