高精度运算模板C语言

#include<iostream>
#include<string.h>
#include<algorithm>
#include<stdio.h>
#define maxx 100
using namespace std;
//高精度加法：
void add(char *s1,char *s2)
{
	int n,m,i,j,a[maxx+10]={0},b[maxx+10]={0};
	int l1=strlen(s1);
	int l2=strlen(s2);
	for(i=l1-1,j=maxx;i>=0;i--,j--)
		a[j]=s1[i]-'0';
	for(i=l2-1,j=maxx;i>=0;i--,j--)
		b[j]=s2[i]-'0';
	for(i=maxx;i>=0;i--)
	{
	a[i]+=b[i];
	if(a[i]>=10) a[i-1]+=a[i]/10,a[i]%=10;
	}
	int temp=0;
	for(i=0;i<=maxx;i++)
	{
	if(a[i]!=0 && temp==0) temp=1;
	if(temp==1)
	printf("%d",a[i]);
	}
	puts("");
}
//高精度减法
void sub(char *s1,char *s2)
{
	int n,m,i,j,a[maxx+10]={0},b[maxx+10]={0},c[maxx+10]={0};
	int l1=strlen(s1);
	int l2=strlen(s2);
	for(i=l1-1,j=maxx;i>=0;i--,j--)
	a[j]=s1[i]-'0';
	for(i=l2-1,j=maxx;i>=0;i--,j--)
	b[j]=s2[i]-'0';
	for(i=maxx;i>=0;i--)
	{
	c[i]+=a[i]-b[i];
	if(c[i]<0)
	{
	c[i]+=10;
	c[i-1]--;
	}
	}
	int temp=0;
	for(i=0;i<=maxx;i++)
	{
	if(c[i]!=0 && temp==0) temp=1;
	if(temp==1)
	printf("%d",c[i]);
	}
	puts("");
}
//高精度乘法
void mul(char *s1,char *s2)
{
	int n,m,i,j,i1,j1,a[maxx+10]={0},b[maxx+10]={0},c[2*maxx+10]={0};
	int l1=strlen(s1);
	int l2=strlen(s2); 
	for(i=l1-1,j=0;i>=0;i--)
		a[j++]=s1[i]-'0';
	for(i=l2-1,j=0;i>=0;i--)
		b[j++]=s2[i]-'0';
	for(i=0;i<l2;i++)
		for(j=0;j<l1;j++)
		c[i+j]+=b[i]*a[j];
	for(i=0;i<maxx*2;i++)
	if(c[i]>=10)
	{
	c[i+1]+=c[i]/10;
	c[i]%=10;
	}
	int t=0; 
	for(i=maxx*2;i>=0;i--)
	if(t) 
	printf("%d",c[i]);
	else if(c[i]){
	printf("%d",c[i]);
	t=1;
	}
	 if(!t)  printf("0");
	puts("");
}
//高精度除法
#include<iostream>
#include<string.h>
#include<algorithm>
#include<stdio.h>
#define maxx 100
using namespace std;
int an1[maxx+10],an2[maxx+10];
int ans[maxx+10];
	/*
	sun函数:长度为 l1的大整数p1减去长度为l2的大整数p2，减的结果放在p1里，返回值
	代表结果的长度 
	*/
int sub(int *p1,int *p2,int l1,int l2)
{
	int i;
	if(l1<l2)
	  return -1;
	//下面判断p1是否比p2大，如果不是，返回-1 
	bool temp=false;
	if(l1==l2){
	for(i=l1-1;i>=0;i--){
	if(p1[i]>p2[i]) 
	temp=true;
	else if(p1[i]<p2[i]){
		if(!temp)
	return -1;
	}
	}
	}
	for(i=0;i<l1;i++){ //做减法 
	p1[i]-=p2[i]; //要求调用本函数时给的参数能确保当i>=len2时，p2[i]=0 
	if(p1[i]<0){
	p1[i]+=10;
	p1[i+1]--;
	}
	}
	for(i=l1-1;i>=0;i--)
	if(p1[i]) 
	  return i+1; 
		return 0;
}	
int main()
{
	int i,j;
	char s1[maxx+10],s2[maxx+10];
	scanf("%s%s",s1,s2);
	int l1=strlen(s1),l2=strlen(s2);
	memset(an1,0,sizeof(an1));
	memset(an2,0,sizeof(an2));
	memset(ans,0,sizeof(ans));
	for(i=l1-1,j=0;i>=0;i--)
	an1[j++]=s1[i]-'0';
	for(i=l2-1,j=0;i>=0;i--)
	an2[j++]=s2[i]-'0';
	if(l1<l2){
	printf("0\n");
	return 0;
	}
	l1=sub(an1,an2,l1,l2);
	if(l1<0){ 
	printf("0\n");
	return 0;
	}
	else if(l1==0){
	printf("1\n");
	return 0;
	}
	ans[0]++;
	int ntimes=l1-l2;
	if(ntimes<0) 
	goto outputresult;
	else if(ntimes>0){
	for(i=l1-1;i>=0;i--){
	if(i>=ntimes)
	an2[i]=an2[i-ntimes];
	else
	an2[i]=0;
	}
	}
	l2=l1;
	for(j=0;j<=ntimes;j++){
	int ntmp;//一直减到不够为止，先减去若干个 ans*（10 的 ntimes 次方）
	//不够减了，再减去若干个ans*（10 的 ntimes 次方）
	while((ntmp=sub(an1,an2+j,l1,l2-j))>=0){
	l1=ntmp;
	ans[ntimes-j]++;
	}
	}
	outputresult:
	//下面的循环统一处理进位问题 
	for(i=0;i<maxx;i++){
	if(ans[i]>=10){
	ans[i+1]+=ans[i]/10;
	ans[i]%=10;
	}
	} 
	bool t=false;
	for(i=maxx;i>=0;i--)
	if(t) printf("%d",ans[i]);
	else if(ans[i]){
	printf("%d",ans[i]);
	t=true;
	}
	if(!t)
	printf("0\n");
	printf("\n");
	return 0;
}