关于大数(C语言)

大数,即数字位数超乎想象(至少超乎long long想象)的数。因其位数过大,在c语言中不能对其进行直接运算;若想实现运算,则需要一些奇淫小技巧。

大数相加

两数相加的手算操作:
将两数右对齐,由右至左对位相加,结果上再加上一位的进位(如果有的话);保留该结果个位为这一位,再将多出的位数(也是如果有的话)进给下一位;反复进行上述操作,直至两数均消耗完。
懂的都懂,不懂的我也没有办法

模拟手算,即可实现大数相加

下面放代码

#include <stdio.h>

int main()
{
    char a[500], b[500];
    scanf("%s %s",a,b);//输入
    int da = 0, db = 0;//用于录每个数字的长度
    //将ascii码转换成数字
    int tmpa[500] = {0};
    for(int i = 0; a[i] != 0; i++)
    {
        tmpa[i+1] = a[i] - '0';
        da++;
    }
    for(int i = 1; i <= da; i++)
    {
        a[i] = tmpa[i];
    }
    int tmpb[500] = {0};
    for(int i = 0; b[i] != 0; i++)
    {
        tmpb[i+1] = b[i] - '0';
        db++;
    }
    for(int i = 1; i <= db; i++)
    {
        b[i] = tmpb[i];
    }
    //右对齐处理
    int d;//有效运算长度
    if(da > db)
    {
        d = da;
        int tmp[500] = {0};
        for(int i = d; i >= db - 1; i--)
        {
            tmp[i] = b[i - (da - db)];
        }
        for(int i = 1; i <= d; i++)
        {
            b[i] = tmp[i];
        }
    }
    if(db > da)
    {
        d = db;
        int tmp[500] = {0};
        for(int i = d; i >= da - 1; i--)
        {
            tmp[i] = a[i-(db-da)];
        }
        for(int i = 1; i <= d; i++)
        {
            a[i] = tmp[i];
        }
    }
    if(da == db) d = da;
    //模拟手算加法
    char c[510] = {0}, add = 0;
    for(int i = d; i > 0; i--)
    {
        int tmp = a[i] + b[i];
        c[i] = tmp % 10 + add;
        add = tmp / 10;
    }
    //输出
    if(add != 0) printf("%d",add);
    for(int i = 1; i <= d; i++)
    {
        printf("%d",c[i]);
    }
}

 
 

大数相乘

学废了大数相加,大数相乘不是so easy?即像手算乘法一样,将第一个数的每一位与第二个数相乘,得到的一堆数进行合理的错位相加,再处理一下进位即可。如图:

为了简化代码,我们也可以将两个数的各位相乘,根据结果的位权进行错位相加,再稍稍处理一下进位即可。如图:在这里插入图片描述
代码实现:

#include <stdio.h>
#include <string.h>

int main()
{
    char a1[10001], b1[10001];
    int a[10001], b[10001], x, len, c[10001];
    //输入
    gets(a1);
    gets(b1);
    int lena = strlen(a1);
    int lenb = strlen(b1);
    //将ascii码转换成数字
    for(int i = 1; i <= lena; i++)
    {
        a[i] = a1[lena-i] - '0';
    }
    for(int i = 1; i <= lenb; i++)
    {
        b[i] = b1[lenb-i]-'0';
    }
    //模拟手算乘法
	for(int i = 1; i <= lenb; i++)
    {
        for(int j = 1; j <= lena; j++)
        {
            c[i+j-1]+=a[j]*b[i];
        }
    }
    //结果处理
    for(int i = 1; i < lena + lenb; i++)
	{
        if(c[i] > 9)
	    {
		    c[i+1]+=c[i]/10;
		    c[i]%=10;
	    }
    }
    //判断结果位数
	len = lena + lenb;
    while(c[len]==0 && len>1) 
    {
        len--;
    }
    //输出
    for(int i = len; i >= 1; i--) 
    {
        printf("%d",c[i]);
    }
}

 
 

大数阶乘

按照大数乘法的方法模拟即可
这里为了偷懒 简化代码,将每6位数字存在一个int里,故只能计算1000以内阶乘

#include<stdio.h>
#include<math.h>

void factorial(int n);

int main()
{
    int n, d = 0;  //输入
    scanf("%d", &n);
    factorial(n);
}

void factorial(int n)
{
	int a[10000];
    int more = 0;
    int length = 1;
	int temp = 1; 
    int i, j;
    a[0] = 1;
    //阶乘计算
	for(i = n; i > 1; i--)
    {
		more = 0;
		for(j = 0; j < length; j++)
		{
			temp = a[j] * i + more;
			a[j] = temp % 1000000;
			more = temp / 1000000;
		}
		a[j] = more;
		if(more > 0) length++;
	} 
	//输出
	printf("%d", a[length-1]);
	for(i = length - 2; i >= 0; i--)
	    printf("%06d", a[i]);
}

 
 

大数相减

没什么好说的,直接放代码

#include <stdio.h>
#include <string.h>

int main()
{
    char a[10001], b[10001];
    scanf("%s%s",a,b);//输入
    int lena = strlen(a);
    int lenb = strlen(b);
    //将ascii码转换成数字
    int tmpa[10001] = {0};
    for(int i = 0; a[i] != 0; i++)
        tmpa[i+1] = a[i] - '0';
    for(int i = 1; i <= lena; i++)
        a[i] = tmpa[i];
    int tmpb[10001] = {0};
    for(int i = 0; b[i] != 0; i++)
        tmpb[i+1] = b[i] - '0';
    for(int i = 1; i <= lenb; i++)
        b[i] = tmpb[i];
    //右对齐处理
        int tmp[10001] = {0};
        for(int i = lena; i > lena - lenb; i--)
            tmp[i] = b[i - (lena - lenb)];
        for(int i = 1; i <= lena; i++)
            b[i] = tmp[i];

    //模拟手算减法
    char c[10001] = {0};
    for(int i = lena; i >= 1; i--)
    {
        if(a[i] >= b[i]) 
            c[i] = a[i] - b[i];
        else
        {
            c[i] = a[i] - b[i] + 10;
            a[i-1]--;
        }
    }
    //输出
    int i;
    for(i = 1; i <= lena; i++)
        if(c[i]!=0) 
        {
            printf("%d",c[i]);
            break;
        }
        i++;
    for(; i <= lena; i++)
            printf("%d",c[i]);
}

 
 

大数除法

这里根据大佬的思路整理了一下代码,具体思路可以前往大佬博客

#include<stdio.h>
#include<string.h>

char a[100], b[100];  //用两个字符串用来输入两个大数
int x[100], y[100], z[100];  //被除数, 除数, 商
int digit;  //大数的位数

void sub(int *x, int *y, int len1, int len2);  //大数减法
int judge(int *x, int *y, int len1, int len2);  //判断两数大小

int main()
{
    int i, j = 0, k = 0, temp;
    int len1, len2, d;  //d为两个大数位数的差值

    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("the result is:0\n");
        printf("the quotient is:");
        puts(a);
    }
    else //当被除数位数大于或者等于除数位数时
    {
        d = len1 - len2;  //两个大数位数的差值
        for(i = len1 - 1; i >= 0; i--)  //将除数后补零,使得两个大数位数相同
        {
            if(i >= d) y[i] = y[i-d];
            else y[i] = 0;
        }
        len2 = len1;  //将两个大数数位相同

        digit = len1;  //将原被除数位数赋值给digit

        for(j = 0; j <= d; j++)
        {
            z[d-j] = 0;

            while( ( ( temp = judge(x, y, len1, len2) ) >= 0) && digit >= k)  //判断两个数的大小以及被除数位数与除数原位数的关系
            {
                sub(x, y, len1, len2);  //大数减法函数

                z[d-j]++;  //储存商的每一位

                len1 = digit;  //重新修改被除数的长度

                if(len1 < len2 && y[len2-1] == 0)
                    len2 = len1;  //将len1长度赋给len2
            }

            if(temp < 0)  //若被除数小于除数,除数减小一位
            {
                for(i = 1; i < len2; i++)
                    y[i-1] = y[i];
                y[i-1] = 0;
                if(len1 < len2) len2--;
            }
        }

        printf("the result is:");
        for(i = d; i > 0; i--)  //去掉前缀0
            if(z[i]) break;

        for(; i >= 0; i--)
            printf("%d", z[i]);
        printf("\n");

        printf("the quotient is:");
        for(i = len1; i > 0; i--)  //去掉前缀0
            if(x[i]) break;

        for(; i >= 0; i--)
            printf("%d", x[i]);
        printf("\n");
    }
}

void sub(int *x, int *y, int len1, int len2)//大数减法
{
    for(int i = 0; i < len1; i++)
    {
        if(x[i] < y[i])
        {
            x[i] = x[i] - y[i] + 10;
            x[i+1]--;
        }
        else
            x[i] = x[i] - y[i];
    }
}

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
    }
}

 
 

快速幂取余

快速幂,即幂运算的优化。在数学中有 (m)2n=(m2)n ,于是我们便可以通过这种方法对幂运算做一些优化。
eg:

216 = (22)8 = (42)4 = (162)2 = 2562 = 65536

95 = 9*(92)2 = 9*6561 = 59049

当指数较大时,快速幂的时间复杂度能在循环乘法上大大优化。
关于快速幂取余,我们先来介绍一个引理:

(a * b) mod n = ( ( a mod n ) * ( b mod n ) ) mod n

当对幂运算的结果进行取余的时候,便可以用这种方法优化。
eg:

(410) mod 3 = ( (4 mod 3)10 )mod 3 = (110) mod 3 = 1

下面是代码:

#include <stdio.h>

int fastPower(int base, int power, int a);

int main()
{
	int x, y, a;  //基数,幂次,取余值
	scanf("%d %d %d", &x, &y, &a);
	printf("%d", fastPower(x, y, a));
}

int fastPower(int base, int power, int a) 
{
    int result = 1;
    while (power > 0) 
	{
        if (power & 1)
		{
            result = result * base % a;
        }
        power >>= 1;
        base = (base * base) % a;
    }
    return result;
}
  • 5
    点赞
  • 35
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值