大数,即数字位数超乎想象(至少超乎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;
}