高精度计算时一般用一个数组来存一个数,数组的一个元素对应于数的一位,由于数组要仅为,可将数组的一个元素对应于数的一位,由于数计算时可能要进位,所以将数由低位开始存在数组,a[0]表示数的位数:
一、高精度的存储与加法:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[30],b[30];
int sum(int a[],int b[]){
int k=a[0]>b[0]?a[0]:b[0];
for(int i=1;i<=k;i++){
a[i+1]+=(a[i]+b[i])/10;//先进位
a[i]=(a[i]+b[i])%10;//计算当前位数字
}
if(a[k+1]>0)a[0]=k+1;
else a[0]=k;
return 0;
}
void print(int c[]){
for(int i=c[0];i>=1;i--){
cout<<c[i];
}
cout<<endl;
}
int main(){
string s1,s2;
cin>>s1>>s2;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
a[0]=s1.length();
for(int i=1;i<=a[0];i++){
a[i]=s1[a[0]-i]-'0';
}
b[0]=s2.length();
for(int i=1;i<=b[0];i++){
b[i]=s2[b[0]-i]-'0';
}
sum(a,b);
print(a);
}
二、减法
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[30],b[30];
void print(int c[]){
for(int i=c[0];i>=1;i--){
cout<<c[i];
}
cout<<endl;
}
int compare(int a[],int b[]){
if(a[0]>b[0])return 1;
if(a[0]<b[0])return -1;
for(int i=a[0];i>0;i--){
if(a[i]>b[i])return 1;
if(a[i]<b[i])return -1;
}
return 0;
}
int subtract(int a[],int b[]){
int flag=compare(a,b);
if(flag==0){memset(a,0,sizeof(a));return 0;}
if(flag==1){//a>b
for(int i=1;i<=a[0];i++){
if(a[i]<b[i]){a[i+1]--;a[i]+=10;}
a[i]=a[i]-b[i];
}
while(a[a[0]]==0)a[0]--;
return 0;
}
if(flag==-1){//b>a
for(int i=1;i<=b[0];i++){
if(b[i]<a[i]){b[i+1]--;b[i]+=10;}
a[i]=b[i]-a[i];
}
a[0]=b[0];
while(a[b[0]]==0)a[0]--;
return 0;
}
}
int main(){
string s1,s2;
cin>>s1>>s2;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
a[0]=s1.length();
for(int i=1;i<=a[0];i++){
a[i]=s1[a[0]-i]-'0';
}
b[0]=s2.length();
for(int i=1;i<=b[0];i++){
b[i]=s2[b[0]-i]-'0';
}
subtract(a,b);
print(a);
}
三、高精度乘单精度:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[30],b[30];
int multi1(int a[],int n){
for(int i=1;i<=a[0];i++){
a[i]*=n;
}
for(int i=1;i<=a[0];i++){
a[i+1]+=a[i]/10;
a[i]%=10;
}
int m=a[0]+1;
while(a[m]>0){a[m+1]=a[m]/10;a[m]%=10;m++;a[0]++;}
}
void print(int c[]){
for(int i=c[0];i>=1;i--){
cout<<c[i];
}
cout<<endl;
}
int main(){
string s1;
int d;
cin>>s1>>d;
memset(a,0,sizeof(a));
a[0]=s1.length();
for(int i=1;i<=a[0];i++){
a[i]=s1[a[0]-i]-'0';
}
multi1(a,d);
print(a);
}
四、高精度乘高精度
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[100]={0},b[100]={0},c[200]={0};
int main(){
string s1,s2;
cin>>s1>>s2;
int len1=s1.length(),len2=s2.length();
int k=len1+len2;
for(int i=0;i<len1;i++)a[i]=s1[len1-i-1]-'0';
for(int i=0;i<len2;i++)b[i]=s2[len2-i-1]-'0';
for(int i=0;i<len1;i++){
for(int j=0;j<len2;j++){
c[i+j]+=a[i]*b[j];
c[i+j+1]+=c[i+j]/10;
c[i+j]%=10;
}
}
while(!c[k])k--;
if(k<0)printf("0\n");
else{
for(int i=k;i>=0;i--)printf("%d",c[i]);
}
}
五、高精度除单精度
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[100]={0},c[100]={0};
int main(){
string s1;
int n,d=0;
cin>>s1>>n;
int k=s1.length();
for(int i=0;i<k;i++)a[i]=s1[k-i-1]-'0';
for(int i=k-1;i>=0;i--){
d=d*10+a[i];
c[i]=d/n;
d=d%n;
}
while(c[k-1]==0&&k>1)k--;
for(int i=k-1;i>=0;i--)printf("%d",c[i]);
printf("\n余数=%d\n",d);
}
六、大数取模(模数是底精度)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
#include<map>
#include<set>
#include<string>
#include<cstring>
using namespace std;
#define ll long long
const int mod=1e9+7;
int main(){
string s1;
cin>>s1;
ll ans=0;
for(int i=0;i<s1.length();i++){
ans=ans*10+s1[i]-'0';
ans%=mod;
}
printf("%lld\n",ans);
}
七、大数除法、取模。转: https://blog.csdn.net/lisp1995/article/details/52319348
#include<stdio.h>
#include<string.h>
char a[100],b[100];//用两个字符串用来输入两个大数
int x[100],y[100],z[100],m[100];//被除数 除数 商 余数
int digit;//大数的位数
void sub(int x[],int y[],int len1,int len2)//大数减法
{
int i;
for(i=0;i<len1;i++)
{
if(x[i]<y[i])
{
x[i]=x[i]+10-y[i];
x[i+1]--;
}
else
x[i]=x[i]-y[i];
}
for(i=len1-1;i>=0;i--)//判断减法结束之后,被除数的位数
{
if(x[i])
{
digit=i+1;
break;
}
}
}
int judge(int x[],int y[],int len1,int len2)
{
int i;
if(len1<len2)
return -1;
if(len1==len2)//若两个数位数相等
{
for(i=len1-1;i>=0;i--)
{
if(x[i]==y[i])//对应位的数相等
continue;
if(x[i]>y[i])//被除数 大于 除数,返回值为1
return 1;
if(x[i]<y[i])//被除数 小于 除数,返回值为-1
return -1;
}
return 0;//被除数 等于 除数,返回值为0
}
}
int main()
{
int i,j=0,k=0,temp;
int len1,len2,len;//len两个大数位数的差值
while(~scanf("%s %s",a,b))
{
len1=strlen(a);//被除数位数
len2=strlen(b);//除数位数
for(i=len1-1,j=0;i>=0;i--)//将字符串中各个元素倒序储存在数组中
x[j++]=a[i]-'0';
for(i=len2-1,k=0;i>=0;i--)
y[k++]=b[i]-'0';
if(len1<len2)//当被除数位数 小于 除数位数时
{
printf("商是:0\n");
printf("余数是:");
puts(a);
}
else //当被除数位数大于或者等于除数位数时
{
len=len1-len2;//两个大数位数的差值
for(i=len1-1;i>=0;i--)
//将除数后补零,使得两个大数位数相除。
//被除数:4541543329 除数:98745,加零后:9874500000
{
if(i>=len)
y[i]=y[i-len];
else
y[i]=0;
}
len2=len1;//将两个大数数位相同
digit=len1; //将原被除数位数赋值给digit
for(j=0;j<=len;j++)
{
z[len-j]=0;
while(((temp=judge(x,y,len1,len2))>=0)&&digit>=k)
//判断两个数之间的关系以及位数与除数原位数的关系
{
sub(x,y,len1,len2); //大数减法函数
z[len-j]++;//储存商的每一位
len1=digit;//重新修改被除数的长度
if(len1<len2&&y[len2-1]==0)
len2=len1;//将len1长度赋给len2;
}
if(temp<0)
//若被除数 小于 除数,除数减小一位。
//例如:被除数:4541543329 除数:(原)98745,(加零后)9874500000,后退一位后:0987450000
{
for(i=1;i<len2;i++)
y[i-1]=y[i];
y[i-1]=0;
if(len1<len2)
len2--;
}
}
printf("商是:");
for(i=len;i>0;i--)//去掉前缀0
{
if(z[i])
break;
}
for(;i>=0;i--)
printf("%d",z[i]);
printf("\n");
printf("余数是:");
for(i=len1;i>0;i--)
{
if(x[i])
break;
}
for(;i>=0;i--)
printf("%d",x[i]);
printf("\n");
}
}
return 0;
}
学习了这个博主的文章:点击打开链接