C/C++程式設計小練習 大數減法之浮點大數減法演算法
阿新 • • 發佈:2019-02-08
#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; }