1. 程式人生 > >C/C++程式設計小練習 大數減法之浮點大數減法演算法

C/C++程式設計小練習 大數減法之浮點大數減法演算法

#include <iostream>
#include <cstring>
using namespace std;
void reverse_str(char *a,int size){
	for(int i=0;i<size/2;++i){
		char t=a[i];
		a[i]=a[size-i-1];
		a[size-i-1]=t;
	}
}

int f_to_i(char *a,char *c){//將浮點數字符串轉化為整形字串,並返回浮點數小數部分的長度
	int len=strlen(a);
	char *p1=a;
	char *p2=NULL;
	for(int i=0;i<len;++i){
		if(a[i]=='.'){
			a[i]='\0';
			p2=&a[i+1];
			break;
		}
	}
	strcat(p1,p2);
	strcpy(c,p1);
	return strlen(p2)+1;
}

//核心演算法和我的浮點大數加法一樣,就是按照小數點的位置做對齊處理,小數部分不夠的補0,然後將其當作大整數減法來處理
void large_float_sub(char a[],char b[],char c[]){
	char p1[20000]={0};
	char p2[20000]={0};
	int pos1=f_to_i(a,p1);
	int pos2=f_to_i(b,p2);
	int max_pos;
	if(pos1>pos2){
		int t=pos1-pos2;
		int len=strlen(p2);
		int i;
		for(i=0;i<t;++i){
			p2[i+len]='0';
		}
		p2[i+len]='\0';
		max_pos=pos1;
	}else{
		int t=pos2-pos1;
		int len=strlen(p1);
		int i;
		for(i=0;i<t;++i){
			p1[len+i]='0';
		}
		p1[i+len]='\0';
		max_pos=pos2;
	}

	int a_len=strlen(p1);
	int b_len=strlen(p2);
	reverse_str(p1,a_len);
	reverse_str(p2,b_len);
	int max_len=a_len>b_len?a_len:b_len;
	if(a_len>b_len){
		for(int i=0;i<a_len||i<b_len;++i){
			c[i]+=p1[i]-p2[i];
			if(i>=a_len||i>=b_len){
				c[i]-='0';
			}
			if(c[i]<0){
				c[i]+=10;
				c[i+1]-=1;
			}
		}
	}else if(a_len<b_len){
		cout<<"-";
		for(int i=0;i<a_len||i<b_len;++i){
			c[i]+=p2[i]-p1[i];
			if(i>=a_len||i>=b_len){
				c[i]-='0';
			}
			if(c[i]<0){
				c[i]+=10;
				c[i+1]-=1;
			}
		}
	}else{
		if(strcmp(a,b)>=0){
			for(int i=0;i<a_len||i<b_len;++i){
				c[i]+=p1[i]-p2[i];
				if(i>=a_len||i>=b_len){
					c[i]-='0';
				}
				if(c[i]<0){
					c[i]+=10;
					c[i+1]-=1;
				}
			}
		}else{
			cout<<"-";
			for(int i=0;i<a_len||i<b_len;++i){
				c[i]+=p2[i]-p1[i];
				if(i>=a_len||i>=b_len){
					c[i]-='0';
				}
				if(c[i]<0){
					c[i]+=10;
					c[i+1]-=1;
				}
			}
		}
	}

	int j;
	for(j=max_len;j>0;--j){//去掉前導0
		if(c[j]!=0){
			break;
		}
	}
	for(int i=j;i>=max_pos;--i){
		cout<<(int)c[i];
	}
	cout<<".";
	for(int i=max_pos-1;i>=0;--i){
		cout<<(int)c[i];
	}
	cout<<endl;
}

int main(){
	char a[20000]={0};
	char b[20000]={0};
	char c[40001]={0};
	cin>>a;
	cin>>b;
	large_float_sub(a,b,c);
	return 0;
}