poj 1503 任意精度计算 求 R 的 n 次幂
//思路：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
//
将字符转化为数字
8
unsigned
int
change (
char
s[LENGTH], unsigned
int
s1[LENGTH
-
1
])
9
{
10
int ss[LENGTH - 1]; //ss 放未逆置的整数
11 memset (ss, 0, sizeof(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
32
void
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
50
int
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
2
#include
<
stdio.h
>
3
#include
<
stdlib.h
>
4
#include
<
string
.h
>
5
#define
LENGTH 6
//
小数的位数（含小数点）
6
7
//
将字符转化为数字
8
unsigned
int
change (
char
s[LENGTH], unsigned
int
s1[LENGTH
-
1
])9
{10
int ss[LENGTH - 1]; //ss 放未逆置的整数 11
memset (ss, 0, sizeof(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
32
void
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
50
int
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