#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define LENGTH 6 //小數的位數（含小數點）
//將字符轉化為數字
unsigned int change (char s[LENGTH], unsigned int s1[LENGTH - 1])
{
int ss[LENGTH - 1]; //ss 放未逆置的整數
memset (ss, 0, sizeof(ss));
int k = 0 ;
for (int i = 0; i < LENGTH && s[i]; ++i) //考慮特殊數據如：0.0001
{
if (s[i] != '.')
ss[k++] = s[i] - '0';
}
for (int j = 0;j < LENGTH - 1; j++)
{
s1[j] = ss[LENGTH - 2 - j];
}
int m = 0;
while ( (s[m] != '.') && s[m] )
++m;
return LENGTH - 1 - m; //小數點位數
}
//大數乘法運算
//函數返回 s2
void mu1 (unsigned int s1[LENGTH - 1],unsigned int s2[130])
{
int ss[130];
memset ( ss, 0, sizeof(ss) );
for ( int i = 0; i < LENGTH - 1; i++)
for (int j = 0;j < 130; j++) //難點:因為返回新的s2之后位數會增加 最多時 5* 25 = 125
ss [i + j] += s1[i] * s2[j];
//將 兩個大數相乘得的積ss中進行進位處理后放到s2 中
int c = 0;
for (int i = 0;i < 130;i++)
{
s2[i] = (c + ss[i]) % 10;
c = (c + ss[i]) / 10;
}
}
int main()
{
int n;
char s[LENGTH]; //要處理的冪 R
unsigned int s1[LENGTH - 1]; //將 R 轉化成數字
unsigned int s2[130];
while(scanf ("%s%d", s, &n) != EOF)
{
memset (s1, 0, sizeof (s1));
memset (s2, 0, sizeof (s2));
int j = change (s, s1); //得到小數點所在位置
change (s,s2); //得到s2 和 s1 進行冪運算
for ( int i = 1; i < n; i ++)
mu1 (s1,s2);
//在s2中前面的代表小數位，后面的代表整數位，
//所以關鍵是通過數值關系找到小數點的位置
//例：0.1010 * 0.1010 = 0.01020100
int m = 129;//去掉前導0
while ( (!s2[m]) && m)
--m;
int k = 0; //去掉尾0
while ( ( !s2[k] ) && (k < 130))
++k;
//輸出整數位
for (int i = m; i >= n * j; i--)
printf ("%d",s2[i]);
//輸出小數點
if ( j && n * j >= k + 1)
printf (".");
for (int i = n*j -1; i >= k; --i)
printf ("%d", s2[i]);
printf ("\n");
}
return 0;
// system ("pause");
}
posted on 2010-08-29 11:35
雪黛依夢 閱讀(1912)
評論(0) 編輯 收藏 引用 所屬分類:
大數