大数计算 总结

万能开头

贴一下代码

#include<bits/stdc++.h> 
using namespace std;

大数主要思路

1、对数据进行判断处理 swap strlen strcmp
2、进行反向存储
3、进行计算，进位 借位（一般可以不与进制数进行比较，直接进位，因 为前置零可以去除），注意最后一位的处理
进位：
阶乘：a[n]=a[n]*j;
a[n+1]+=a[n]/n;
a[n]%=n;
乘法：a[i+j]+=a[i]*a[j];
a[i+j+1]+=a[i+j]/n;
a[i+j]%=n;
借位：
a[i]+=10; a[i+1]–; c[i]=a[i]-b[i];
c[i]=(x*10+a[i])/b;
x=(x*10+a[i])%b;
4、前置零清除：while(s[s[0]]==0)–s[0];

大数加法

主要思路就是数字用字符数组反向读入，然后最好判断一下长度，把长的放在前面然后做读入。做完之后从后面依次进位,
进位的时候要注意最后一位的处理！！！
以后尽量不用字符串第一个a[0].

//交换字符串
void swap(char *p,char *q)
{
char a[20];
strcpy(a,p);
strcpy(p,q);
strcpy(q,a);
}

//用a[0]读字符数组的长度，反向读入字符数组将他变成整型方便计算
//尽量把这个写成函数
a[0] = strlen(s);
for(i=1;i<=a[0];++i)
a[i]=s[s[0]-i]-48;

//进位
while(i<=l1||i<=l2)
{
    c[i]=a[i]+b[i]+x;
    x=c[i]/10;
    c[i]%=10;
    i++;

}
c[i]=x;//很重要！！！！
while(c[i]==0)i--;//去除前置零
for(j=i;j>=1;++j)cout<<c[j];

大数减法

总得来说和大数加法差不多 这里是借位 然后最重要的是首先要判断那个哪个大。同样要注意最后一位

//判断字符串
if(strlen(s1<strlen(s2)||
strlen(s1)==strlen(s2)&&strcmp(s1,s2)<0){
    swap(s1,s2);
    cout<<"-";
}

//借位
t2=min(strlen(s1),strlen(s2));
while(i<=t2){
        x=aa[i]-bb[i];
        if(x<0){aa[i]=aa[i]-bb[i]+10;

        --aa[++i];
        }
        else {aa[i]=aa[i]-bb[i];
        ++i; 
        }

    }
    while(aa[t1]==0) --t1;  
    for(i=t1;i>=1;--i)cout<<aa[i];

高精度乘法

最重要的就是进位思想是首先反向存入，然后就是进位处理

for(i=0;i<la;++i)
    for(j=0;j<lb;++j){
        cc[i+j]+=aa[i]*bb[j];
        cc[i+j+1]+=cc[i+j]/10;
        cc[i+j]%=10;

    }
    k=la+lb;
    while(cc[k]==0)--k;

高精度除以低精度

而且其实这里不需要反向存入，直接从前面开始除就好了
但是需要反向输出
就是一个公式

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


    for(i=la-1;i>=0;--i){
        k=k*10+aa[i];//每位依次
        cc[i]=k/b;
        k%=b;   
        //cc[i]=(x*10+aa[i])/b;
        //x=(x*10+aa[i])%b;
    }

    while(cc[la-1]==0&&la>1)--la;
    for(i=la-1;i>=0;--i) cout<<cc[i]; 
    cout<<endl;
    cout<<k;//x

大整数的因子

是高精度除以低精度除法的一个变式，借助算因子的那个公式

for(j=2;j<=9;++j){
        t=0;
        for(i=0;i<l;++i){
            if(a[i]-48+t>=j)
            t=((a[i]-48+t)%j)*10;
            else
            t=(a[i]-48)*10;
        }
        if(t==0){cout<<j<<" ";flag=1;
        }
    }

t=t*10+c[i]-48;
t%=n;

幂运算 阶乘 阶乘和

幂运算和阶乘运算一致，先进行进位，最后处理最高位。乘一次变一次

//引入第三个数组
for(y=1;y<=n;++y)
    {
    memset(b,0,sizeof(b));
      b[0]=a[0];
        for(j=1;j<=b[0];++j)
        {
            b[j]+=a[j]*y;
            b[j+1]=b[j]/10;
            b[j]%=10; 
        }
        while(b[b[0]+1]>0)
        {
            b[b[0]+2]=b[b[0]+1]/10;
            b[b[0]+1]%=10;
            b[0]++;
        }
        for(i=1;i<=b[0];++i)a[i]=b[i];
        a[0]=b[0];


//引入余数
for(j=1;j<=n;++j){

        k=0;//记得清零
        for(i=0;i<count;++i){
            t=c[i]*j+k;//变这个j就行了
            c[i]=t%10;
            k=t/10;
        }
        while(k)//记录进位 
        {
            c[count++]=k%10;
            k/=10;
        }
    }

#include<bits/stdc++.h>
using namespace std;
int a[100000010];
char c [60];
int main() {
//从1乘到n
    int i,len,n;
    scanf("%d",&n);
    a[1]=1;
    len=1;
    for(int k=2; k<=n; k++) {
        for(i=1; i<=len; i++) {
            a[i]=a[i]*k;
        }
        for(i=1; i<=len; i++) {
            if(a[i]>=10) {
                a[i+1]+=a[i]/10;
                a[i]%=10;
            }
        }
        while(a[len+1]>0) {
            len++;
            a[len+1]=a[len]/10;
            a[len]=a[len]%10;
        }
    }
    for(i=len; i>=1; i--)
        printf("%d",a[i]);

}

c[j]=c[j]*i+k;
k=c[j]/10；
c[j]%=10;
if(k>0&&j>=len)++len;

阶乘和就是运用两个阶乘相加

#include<bits/stdc++.h> 
using namespace std;
int a[50000]={1,1};
int s[55000]={1,0};
int b[50000]={1,1};
int r[55000]={1,0};
int main(){

    int i,j,n,k,y;
    cin>>n;
    for(y=1;y<=n;++y)
    {
    memset(b,0,sizeof(b));
      b[0]=a[0];
        for(j=1;j<=b[0];++j)
        {
            b[j]+=a[j]*y;
            b[j+1]=b[j]/10;
            b[j]%=10; 
        }
        while(b[b[0]+1]>0)
        {
            b[b[0]+2]=b[b[0]+1]/10;
            b[b[0]+1]%=10;
            b[0]++;
        }
        for(i=1;i<=b[0];++i)a[i]=b[i];
        a[0]=b[0];
    memset(r,0,sizeof(r));
    r[0]= max(a[0],s[0]);   
        for(k=1;k<=r[0];++k)
        {
            r[k]+=s[k]+a[k];
            r[k+1]=r[k]/10;
            r[k]%=10;
        }
        while(r[r[0]+1]>0)
        {
            r[r[0]+2]=r[r[0]+1]/10;
            r[r[0]+1]%=10;
            r[0]++;
        }

        s[0]=r[0];

    for(i=1;i<=r[0];++i)s[i]=r[i];


        }
while(s[s[0]]==0)--s[0];
for(i=s[0];i>=1;--i)cout<<s[i]; 
        return 0;   
}