高精度算法

高精度加法：

#include <iostream>

#include <string.h> 

using namespace std;

int main()

{

    char a[1001],b[1001];

    int la,lb,i;

    int x[1001]={0},y[1001]={0},he[1002]={0};

    int t=0;

    cin>>a>>b;

    la=strlen(a);

    lb=strlen(b);

    for (i=0;i<la;i++) x[la-i-1]=(int)a[i]-48;

    for (i=0;i<lb;i++) y[lb-i-1]=(int)b[i]-48;

    

    for (i=0;i<1001;i++)

    {

        t=x[i]+y[i]+he[i];

        if (t<10) he[i]=t;

        else

        {

            he[i+1]++;

            he[i]=t-10;

        }

    }

    for (i=1000;i>=0;i--) if (he[i]!=0) break;

    for (i=i;i>=0;i--) cout<<he[i];

    cout<<endl;

    system("pause");

    return 0;

}

高精度减法：

#include <iostream>

#include <math.h>

#include <string.h>

using namespace std;

char a[1001],b[1001],t[1001];

int main()

{

    

    int x[1001]={0},y[1001]={0},s[1002]={0};

    cin>>a>>b;

    int i,j,la,temp,lb;

    bool flag=false;

    la=strlen(a);lb=strlen(b);

    if (la<lb)

    {

        strcpy(t,a);

        strcpy(a,b);

        strcpy(b,t);

        flag=true;

    }

    else if (la==lb && strcmp(a,b)==-1)

         {

             strcpy(t,a);

             strcpy(a,b);

             strcpy(b,t);

             flag=true;

         }

    //如何逆序，并转换为数值 

    for(i=0;i<la;i++)x[la-i-1]=(int)a[i]-48; 

    for(i=0;i<lb;i++)y[lb-i-1]=(int)b[i]-48; 

    //数学加法原理 

    for(i=0;i<1001;i++)

    {

       temp=x[i]-y[i];

       if(temp>=0)

         s[i]=temp;

       else

         {

            s[i]=temp+10;

            x[i+1]--;

         }

    }

    //输出结果 

    for(i=1000;i>=0;i--)

    {

       if(s[i]!=0)

         break;

    }

    if (flag) cout<<"-";

    for(i=i;i>=0;i--)

       cout<<s[i];

    system("pause");

    return 0;

}

高精度乘法：

#include<cstdio>

#include<cstring>

#include<cmath>

#include<malloc.h>

int an,bn,fa=1,fb=1;

char b1[250],b2[250];/*纪录需要计算的两个高精度数据*/

void input(int a1[],int a2[])

{

int i,ai=1,bi=1;

scanf("%s%s",b1,b2);

an=strlen(b1);

bn=strlen(b2);

if(b1[0]==45){an--;fa=-1;ai=0;}

if(b2[0]==45){bn--;fb=-1;bi=0;}

for(i=0;i<an;i++,ai++){a1[i]=b1[an-ai]-'0';

}

for(i=0;i<bn;i++,bi++)a2[i]=b2[bn-bi]-'0';

return;

}

void subtraction(int a[250],int b[250],int q);

void multiplication(int a[],int b[])

{

  int i,j,k;

  int c[501];

  memset(c,0,sizeof c);

  k=an+bn-1;

  for(i=0;i<an;i++)

  for(j=0;j<bn;j++)

  {

   c[i+j]=a[i]*b[j]+c[i+j];

   c[i+j+1]=c[i+j]/10+c[i+j+1];

   c[i+j]=c[i+j]%10;

  }

  while(!c[k])k--;

  if(fa*fb<0) printf("-");

  for(i=k;i>=0;i--) printf("%d",c[i]);

}

  main()

{

  int a[250],b[250];

  memset(a,0,sizeof a);

  memset(b,0,sizeof b);

  input(a,b);

  

  multiplication(a,b);

  system("pause");

  //while(1);

}