其他分享
首页 > 其他分享> > 高精度整合

高精度整合

作者:互联网

9.7 函数整合

 

  1 #include<cstdio>
  2 #include<iostream>
  3 #include<cstring>
  4 #define ll long long
  5 #define MAXN 400000
  6 
  7 using namespace std;
  8 
  9 char sa[MAXN],sb[MAXN];
 10 int a[MAXN],b[MAXN],c[MAXN],d;
 11 
 12 template<typename type>
 13 inline void get(char sa[],type a[],bool mode)//读入一个数,按位存到a数组中,mode为0正序存,为1倒序存 
 14 {
 15     memset(a,0,sizeof(a));
 16     cin>>sa;
 17     int len=strlen(sa);
 18     if(mode) for(int i(1);i<=len;i++) a[i]=sa[len-i]^48;
 19     else for(int i(1);i<=len;i++) a[i]=sa[i-1]^48; 
 20 }
 21 
 22 template<typename type>
 23 inline void put(int len,type c[],bool mode)//按位输出一个长度为len的数组,mode为0正序输出,为1倒序输出
 24 {
 25     if(mode) for(int i(len);i>0;i--) cout<<c[i];
 26     else for(int i(1);i<=len;i++) cout<<c[i];
 27     puts("");
 28 }
 29 
 30 //高精度加法 
 31 template<typename type>
 32 inline void add(char sa[],char sb[],type a[],type b[],type c[])
 33 {
 34     memset(c,0,sizeof(c));
 35     get(sa,a,1);
 36     get(sb,b,1);
 37     int lena=strlen(sa),lenb=strlen(sb),lenc(0),carry(0);//carry存进位 
 38     while(lenc<=lena||lenc<=lenb)
 39     {
 40         c[lenc]=a[lenc]+b[lenc]+carry;
 41         carry=c[lenc]/10;
 42         c[lenc]%=10;
 43         lenc++;
 44     }
 45     c[lenc]=carry;
 46     while((!c[lenc])&&lenc>1) lenc--;
 47     put(lenc,c,1);
 48     return;
 49 }
 50 
 51 //高精度减法 
 52 template<typename type>
 53 inline void sub(char sa[],char sb[],type a[],type b[],type c[])
 54 {
 55     memset(a,0,sizeof(a));
 56     memset(b,0,sizeof(b));
 57     memset(c,0,sizeof(c));
 58     cin>>sa>>sb;
 59     int lena=strlen(sa),lenb=strlen(sb),lenc(0);
 60     char temp[MAXN];
 61     //"||"的优先级低于"&&"  //strcmp()为字符串比较函数,sa==sb返回0,s1>s2返回一个正整数,s1<s2返回一个负整数 
 62     if(lena<lenb || (lena==lenb)&&strcmp(sa,sb)<0)//处理被减数小于减数的情况 
 63     {
 64         strcpy(temp,sa);//将字符串sa复制到temp
 65         strcpy(sa,sb);
 66         strcpy(sb,temp);
 67         putchar('-');
 68         lena=strlen(sa),lenb=strlen(sb);
 69     }
 70     for(int i(1);i<=lena;i++) a[i]=sa[lena-i]^48;
 71     for(int i(1);i<=lenb;i++) b[i]=sb[lenb-i]^48;
 72     while(lenc<=lena||lenc<=lenb)
 73     {
 74         if(a[lenc]<b[lenc])
 75         {
 76             a[lenc]+=10;
 77             a[lenc+1]--;
 78         }
 79         c[lenc]=a[lenc]-b[lenc];
 80         lenc++;
 81     }
 82     while((!c[lenc])&&lenc>1) lenc--;
 83     put(lenc,c,1);
 84     return;
 85 }
 86 
 87 //高精度乘法 
 88 template<typename type>
 89 inline void mul(char sa[],char sb[],type a[],type b[],type c[])
 90 {
 91     memset(c,0,sizeof(c));
 92     get(sa,a,1);
 93     get(sb,b,1);
 94     int lena=strlen(sa),lenb=strlen(sb),lenc(0),carry(0);
 95     for(int i(1);i<=lena;i++)
 96     {
 97         carry=0;
 98         for(int j(1);j<=lenb;j++)
 99         {
100             c[i+j-1]=c[i+j-1]+a[i]*b[j]+carry;
101             carry=c[i+j-1]/10;
102             c[i+j-1]%=10;
103         }
104         c[i+lenb]=carry;
105     }
106     lenc=lena+lenb;
107     while((!c[lenc])&&lenc>1) lenc--;
108     put(lenc,c,1);
109     return;
110 }
111 
112 //高精度除法-高精除以低精(整除) 
113 template<typename type>
114 inline void div1(char sa[],type a[],type b,type c[])
115 {
116     memset(a,0,sizeof(a));
117     memset(c,0,sizeof(c));
118     cin>>sa>>b;
119     if(!b)
120     {
121         puts("0");
122         return;
123     }
124     int lena=strlen(sa),lenc(0);
125     ll carry(0);
126     for(int i(1);i<=lena;i++) a[i]=sa[i-1]^48;
127     for(int i(1);i<=lena;i++)//对被除数的每一位(包括余数)都除以除数 
128     {
129         c[i]=(carry*10+a[i])/b;
130         carry=(carry*10+a[i])%b;
131     }
132     lenc=1;
133     while((!c[lenc])&&lenc<lena) lenc++;
134     for(int i=lenc;i<=lena;i++) cout<<c[i];
135 }
136 
137 //高精度除法-高精除以高精(整除)
138 //用减法模拟除法,对被除数的每一位(包括余数)都减去除数,直到当前位置的数字(包括余数)小于除数 
139 template<typename type>
140 inline int comp(type a[],type b[])//比较a和b的大小关系,若a>b则为1,若a<b则为-1,若a=b则为0(类似于strcmp) 
141 {
142     if(a[0]>b[0]) return 1;
143     if(a[0]<b[0]) return -1;
144     for(int i=a[0];i>0;i--) //从最高位开始比较
145     {
146         if(a[i]>b[i]) return 1;
147         if(a[i]<b[i]) return -1;
148     }
149     return 0;
150 }
151 template<typename type>
152 inline void sub2(type a[],type b[])//计算a=a-b
153 {
154     int flag=comp(a,b);
155     if(!flag)
156     {
157         a[0]=0;
158         return;
159     }
160     if(flag=1)
161     {
162         for(int i=1;i<=a[0];i++)
163         {
164             if(a[i]<b[i])
165             {
166                 a[i]+=10;
167                 a[i+1]--;
168             }
169             a[i]-=b[i];
170         }
171         while((!a[a[0]])&&a[0]>0) a[0]--;
172         return;
173     }
174 }
175 template<typename type>
176 inline void div2(char sa[],char sb[],type a[],type b[],type c[])//在这里,x[0]存储的是x的位数 
177 {
178     memset(c,0,sizeof(c));
179     get(sa,a,1);//因为是逆序存,a[1]就是个位,a[lena]就是最高位 
180     get(sb,b,1);
181     int lena=strlen(sa),lenb=strlen(sb),lenc(0),carry(0),temp[MAXN];
182     if(lenb==1&&(!b))
183     {
184         puts("0");
185         return;
186     }
187     a[0]=lena,b[0]=lenb,c[0]=a[0]-b[0]+1;//商的位数不超过被除数的位数-除数的位数+1
188     for(int i=c[0];i>0;i--)
189     {
190         memset(temp,0,sizeof(temp));
191         for(int j=1;j<=b[0];j++) temp[j+i-1]=b[j];//复制b数组到temp数组从i开始的地方 
192         temp[0]=b[0]+i-1;
193         while(comp(a,temp)>=0)
194         {
195             c[i]++;
196             sub2(a,temp);
197         }
198     }
199     
200     while((!c[c[0]])&&c[0]>0) c[0]--;//处理前导0 
201     
202     if(!c[0]) puts("0");//输出商 
203     else
204     {
205         for(int i=c[0];i>0;i--) cout<<c[i];
206         puts("");
207     }
208     
209     if(!a[0]) puts("0");//输出余数 
210     else
211     {
212         for(int i=a[0];i>0;i--) cout<<a[i];
213         puts("");
214     }
215     return;
216 }
217 
218 signed main()
219 {
220     add(sa,sb,a,b,c);
221     sub(sa,sb,a,b,c);
222     mul(sa,sb,a,b,c);
223     div1(sa,a,d,c);
224     div2(sa,sb,a,b,c);
225     return 0;
226 }

 

标签:lenc,高精度,int,char,整合,sb,sa,type
来源: https://www.cnblogs.com/void-null/p/15239453.html