高精度运算

1、 有关%运算

类似题目之前写过一次，但是这次写到又给忘记了(#￣～￣#)

接下来看题：



代码如下：

#include<stdio.h>
__int64 yanghui[1050][1050];
int main(){
int i,j;

//求杨辉三角形
for(i=0;i<1010;i++){
for(j=0;j<=i;j++){
if(j==0||j==i){
yanghui[i][j]=1;
}
else {
/*
由于yanghui[i][j]范围最终可能超过int
所以应当边%边+ 使yanghui[i][j]大小控制在一定范围内

同理如遇到乘法也应边%边*
*/
yanghui[i][j]=(yanghui[i-1][j-1]+yanghui[i-1][j])%707063423;
}
}
}

scanf("%d%d",&n,&m);
printf("%I64d",yanghui
[m]);

return 0;
}

2、 数组模拟乘法

运用数组模拟竖式乘法的过程来达到int范围不支持的高精度的运算

例题题目：



代码如下：

/**************************************************************
***竖式乘法运算***

1   2   3   4
×             5
--------------------
1   2
7   0  写下个位，十位以上往前进位
3×5+2=17
**************************************************************
*
数组模拟运算类似此法
每位数字先*i再进行进位 与 由低位到高位边乘i边进位效果相同
*
**************************************************************/

#include<stdio.h>
int a[3000]=1;//记录数字的每一位数。1000!大约有2500+位
int main(){
int i,j;
int n;
int length=1;//记录数字长度
scanf("%d",&n);

a[1]=1;
for(i=1;i<=n;i++){

//每位数*i
for(j=1;j<=length;j++){
a[j] *= i;
}

//进阶
for(j=1;j<=length;j++){
a[j+1] += a[j]/10;//进位
a[j] %= 10;//取个位
}

//此时j==length+1
/*拆分最高位
*
*若a[j]==0则说明原最高位*i之后没有进位；a[j]!=0时则有进位
*
*a[j]可能不是一位数字 故需对a[j]进行拆分
*并重新计算数字长度
*
*/
while(a[j] != 0){
length++;
a[j+1] += a[j]/10;
a[j] %= 10;
j++;
}
}

for(i=length;i>=1;i--){
printf("%d",a[i]);
}
printf("\n");

return 0;
}

以上代码中，a[i]仅存了一位数字。在内存和精度要求更高的情况下，可以使a[i]存三位四位甚至更多位的数字。只需要将以上代码中的10改成相应的10的倍数即可。