高精度数在c语言中实现的核心思想是
1、利用数组按位存储高精度数
2、利用循环等结构模拟基础的数学运算
1、高精度数加法
#include <iostream>
#include <string>
#include <vector>
#include <algorithm>
#include <stdexcept>
#include <string.h>
using namespace std;
#define MAXN 1002
int a[MAXN],b[MAXN];
int main(int argc, const char * argv[])
{
string str1,str2;//保存输入
long int len1,len2;
long int i,j,k;
int up;
/*输入流程*/
while(cin>>str1>>str2)
{
/*初始化各量*/
len1 = str1.length();
len2 = str2.length();
memset(a, 0, sizeof(a));
memset(b, 0, sizeof(b));
/*注意,必须倒着保存数据*/
for (i = len1 - 1, k = 0; i != -1; -- i)//倒序保存
{
a[k] = str1[i] - '0';
k++;
}
for (j = len2 - 1, k = 0; j != -1; -- j)
{
b[k] = str2[j] - '0';
k++;
}
for (i = 0, up = 0; i < MAXN; ++ i)//求和,结果保存在a中
{
a[i] = a[i] + b[i] + up;
up = a[i] / 10;
a[i] %= 10;
}
for (i = MAXN - 1; i != -1; -- i)//找到最高位
{
if (a[i])
{
break;
}
}
for (k = i; k != -1; --k)
{
cout<<a[k];
}
cout<<endl;
}
return 0;
}
2、大整数乘法(大整数*非大整数)
思路与加法类似,都是倒序保存,从低位不断相乘并进位(先),不再复述代码。
对于阶乘问题(大数乘法)先将一个数作为被乘数,同样按位存在一个整数数组,对于该大数所有位,用乘数乘以该位,并获取进位。
void mutex(int n)//ans初始化为0,下面的代码是计算n的阶乘的代码
{
int up=0;
int cnt=0;
while(n>1)
{
for(int i=0;i<=cnt;i++)//通过cnt记录当前的位数(大数*n,用cnt记录大数的位数)
{
ans[i] = ans[i]*n+up;
up = ans[i]/10;
ans[i]=ans[i]%10;
}
while(up)//当最高仍有进位,则cnt++
{
ans[++cnt]=up%10;
up/=10;
}//至此一次乘法结束
n--;
}
}
求阶乘的完整代码:
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
using namespace std;
int flag[101];//标记已经查询过的阶乘值
int ans[101][10000];//保存阶乘的结果
int mutex(int n)//ans初始化为0,下面的代码是计算n的阶乘的代码
{
int idx=n;
int up=0;
int cnt=0;
while(n>1)
{
for(int i=0;i<=cnt;i++)//通过cnt记录当前的位数(大数*n,用cnt记录大数的位数)
{
ans[idx][i] = ans[idx][i]*n+up;
up = ans[idx][i]/10;
ans[idx][i]=ans[idx][i]%10;
}
while(up)//当最高仍有进位,则cnt++
{
ans[idx][++cnt]=up%10;
up/=10;
}
n--;
}
return cnt;
}
int main()
{
int n;
memset(flag,-1,sizeof(flag));
for(int i=0;i<101;i++)
ans[i][0]=1;//初始化
while(cin>>n)
{
if(flag[n]<0){int cnt = mutex(n);flag[n]=cnt;}
for(int i=flag[n];i>=0;i--)
printf("%d",ans[n][i]);
cout<<endl;
}
return 0;
}
三、高精度除法:
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
int m;
char tmp[100];
int num[100];
memset(num,0,sizeof(num));
scanf("%s%d",tmp,&m);
int n = strlen(tmp);
for(int i=0;i<n;i++)//注意这里的保存顺序,其他运算都是倒序保存,由于除法是高位影响低位,所以正序保存方便一点
num[i]=tmp[i]-48;
int re=0;
for(int i=0;i<n;i++)/**除法核心**/
{
int t = num[i];
num[i]=(re*10+t)/m;//保留商
re=(re*10+t)%m;//保留余数
}
int i=0;while(!num[i]){i++;}//找到首位
for(i;i<n;i++)/**输出的核心**/
{
printf("%d",num[i]);
}
printf("\n余数为: %d",re);
return 0;
}
四、两高精度数相乘:通过start实现错位,同时先用sum保存结果,再按位取余相加。
int start;
for(int i=0;i<l1;i++)
{
start=i; //start实现错位
for(int j=0;j<l2;j++)
{
ans[start++]+=a[i]*b[j];
}
}
for(int i=0;i<start;i++)//这里比较特殊,要与普通乘法区分开,结果先完全保存在ans中,再对ans单独处理得到正确结果
{//start记录了两个大数的位数之和,结果数字的位数不会大于二者加1.
ans[i+1]+=ans[i]/10;
ans[i]=ans[i]%10;
}
五、大整数取模
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
char n[100];
int m;
scanf("%s%d",n,&m);//n为大数,m为取余号后的数(n%m)
int len = strlen(n);
int ans=0;
for(int i=0;i<len;i++)
{
ans=(int)(((long long)ans*10+n[i]-'0')%m);
}
printf("%d\n",ans);
return 0;
}
注意:除法与取余操作都是将大数正序保存在数组里(num[0]为最高位),其余(加 减 乘 除)需要倒叙保存